磁共振动态成像方法技术

一种磁共振动态成像方法，在当前帧图像，根据Ｌ级发射尺度激励波形和小波激励波形并计算尺度能量值，且根据小波激励波形计算Ｌ级小波能量值；产生队列，根据当前帧图像的Ｌ级小波能量值所对应的Ｌ级位置值预测Ｌ－１级位置值，把所述位置值放置所述队列；根据所述队列的位置值所对应的下一帧图像的位置值，获取预测小波能量值以及根据Ｌ级发射尺度激励波形并获得Ｌ级预测尺度能量值；小波能量值与预测小波能量值之差，判断是否大于预设判断值，是，则把所述位置值加入队列，否，则删除所述位置值；根据尺度能量值和小波能量值进行图像重建，并生成图像。采用上述的方法，能够通过预测的方法减少数据的采集量，减少成像的时间，达到快速动态成像的目的。

【技术实现步骤摘要】

本专利技术涉及成像领域，特别是涉及一种。
技术介绍
磁共振成像(MRI)是一种被广泛采用的成像诊断技术，可同时获得检查部位的 形态信息和功能信息，具有其它技术所无可比拟的优势，成为当今医学影像检查的重要 手段。然而磁共振成像技术对物体运动的敏感性较高，成像过程中被检查者身体的心 脏跳动、呼吸以及移动等生理性运动都会使影像模糊和对比度失真，以致无法得到具有 临床价值的诊断图像。可见传统的的成像速度不够快，无法快速的对运动的物体 进行成像。
技术实现思路
基于此，有必要提供一种快速成像的。一种，包括以下步骤S10:在当前帧图像，根据所需图 像分辨率确定的L级发射尺度激励波形和小波激励波形；S20:根据所述的尺度激励波形 采集L级尺度系数并计算当前帧图像的L级尺度能量值，且根据所述小波激励波形采集 L级小波系数并计算当前帧图像的L级小波能量值；S30:产生队列，根据所述当前帧图 像的L级小波能量值所对应的L级位置值预测L-I级位置值，把所述位置值放置所述队 列；S40:根据所述队列的位置值所对应的下一帧图像的位置值，发射小波激励波形并 获取预测小波能量值，以及在下一帧图像的根据所需图像分辨率确定的L级发射尺度激 励波形，并获得L级预测尺度能量值；S50:当前帧图像的小波能量值与下一帧图像的预 测小波能量值之差，判断是否大于预设判断值，是，则把所述位置值加入队列，否，则 删除所述位置值；S60:根据当前帧图像的L级尺度能量值和各级小波能量值，以及下一 帧图像的L级预测尺度能量值和各级预测小波能量值，进行图像重建，并生成当前帧图 像和下一帧图像。优选地，在所述当前帧图像和下一帧图像的L级上发射个尺度激励波形和 小波激励波形，所述M为相位编码时的采集数。优选地，所述步骤S30包括幻10:产生队列，判断所述L级小波能量值是否 大于预设L级门限值，所述L级门限值为氏，是，则L级小波能量值的位置值进入队列， 所述位置值存储形式为(m，Lth)，m表示位置，Lth表示所在级数；S330 :根据所述队 列的位置值(m，Lth)按照{2m，(L-I) th}和{2m+l，(L-I) th}的方式预测在L_1级上 两个小波能量值的位置值；S350 :在所述{2m，(L-l)th}和{2m+l，(L_l)th}的位置值 处发射小波激励波形并获取L-I级小波能量值，判断所述L-I小波能量值是否大于L-I级 门限值Eu，是，则所述L-I小波能量值的位置值进入所述队列，否，则从所述队列删除位置值。优选地，所述步骤S50之后还包括判断队列是否为空，是，则整理获取的当 前帧图像的L级尺度能量值和各级小波能量值，以及下一帧图像的L级预测尺度能量值和 各级预测小波能量值，否，则返回步骤S50。优选地，步骤S30包括S320:产生队列，判断所述L级小波能量值是否大于 预设L级门限值，所述L级门限值为是，则L级小波能量值的位置值进入队列，所 述位置值存储形式为(i，L，k，m)，i为1、2、3和4，k和m表示所在位置值；S340 根据所述队列的位置值(i，L，k，m)按照(i，L-1, 2k，2m)、(i，L-1, 2k，2m+l)、 (i，L-1, 2k+l, 2m)和(i，L-1, 2k+l, 2m+l)的方式预测四个在L_1级小波能量值 的采集位置值；S360 在所述(i，L-1, 2k，2m)、(i，L-1, 2k，2m+l)、(i，L-I, 2k+l，2m)和(i，L-1, 2k+l, 2m+l)的位置值处发射小波激励波形并获取L_1级小波 能量值，判断所述L-I小波能量值是否大于L-I级门限值τ κ，是，则进入所述队列， 否，则从所述队列删除位置值。采用上述的方法，能够通过预测的方法减少数据的采集量，减少成像的时间， 达到快速动态成像的目的。附图说明图1为一实施例的磁共振动态成像的方法流程图2为一实施例的磁共振动态成像中的二维动态成像方法流程图3为另一实施例的磁共振动态成像中的三维动态成像的方法流程图。具体实施方式为了提高磁共振动态成像速度，提出了一种对图像基于小波编码的方法来实现 快速的磁共振成像。请参阅附图1 2，为磁共振动态成像的方法流程图。SlO 在当前帧，根据所需图像分辨率确定的L级发射尺度激励波形和小波 激励波形。对图像编码进行分级，L级是根据所需要成像的分辨率大小确定，例如 320X240, 640X480、800 X 600或IOMX 768等，然后在L级上发射尺度激励波形和小波 激励波形。具体地，把L级分为低频区域和高频区域，在低频区域发射尺度激励波形，且为“(力^ ；在高频区域发射小波激励波形，且为以(Mt1 '其中小(y)为尺度函数，ψ (y)为与尺度函数相对应的小波函数，因为不同级的尺度函数和小波函数都是 在相对应的母函数基础上经过伸缩和平移得到。然后在确定的低频区域及高频区域对应 的发射尺度激励波形和小波激励波形，能够精确的获得特定区域的尺度系数并获得尺度 能量值和小波系数并获得小波能量值，根据划分的区域精确的预测需要采集的位置，提 高采集的速度和准确度。在二维成像的实施例中，在L级上发射个尺度激励波形和 小波激励波形，M为传统的相位编码时的采集数。在其它实施例中，在三维成像的实施 例中，小波激励波形为二维小波激励波形，该二维小波激励波形的个数为 个，N表示传统的相位编码时的采集数。S20根据所述尺度激励波形采集L级尺度系数并计算L级尺度能量值，根 据小波激励波形采集L级小波系数并计算L级小波能量值。根据发射尺度激励波形{ΦΜ-LAy)}^ 所获得的接收信号为 Sk(t) = / / ρ (χ, y) ΦM-L, k(y)elxGxtdxdy 并计算其能量值；发射小波激励波形·( —“(>01所获得的接收信号为Sk(t) =Sf ρ (χ, y) k(y)elxextdxdy并计算其能量值，其中，M为传统的相位编码时的采集数，L为所 在级数，k位置值。在其它实施例中，二维小波激励波形<_(>^)所获得的接收信号为 sUKm (0 = \\\p{x,y,z)^]Xm (y,z)e^dxdydz ,并计算其能量值。其中，M为传统的相位编码时的采集数，其中，j取1、2、3、4，M为传统的相位编码时的采集数，L为所在级 数，k和m位置值。S30 产生队列，根据当前帧图像的L级小波能量值所对应的L级位置值预测 L-I级位置值，把位置值放置该队列。具体地，S310:产生队列。首先设置一个队列，用于放置在各级所采集的小波 能量值的位置值数据。判断L级小波能量值是否大于预设L级门限值，该L级门限值为 El,是，则L级小波能量值的位置值进入队列，该位置值存储形式为(m，Lth)，m表示 位置，Lth表示所在级数；否，则小于门限值的^的位置值不进入队列。具体地，根据 在L级所获取的小波能量值，然后把该值与预设的L级门限值进行比较，同时把L级的 门限值用&表示。若所采集的L级小波能量值大于该门限值氏，则把L级小波能量值 的位置值按照规定的模式存储在步骤幻10所设置好的队列中。S330 根据所述队列的位置值(m，Lth)按照{2m，(L_l)th}和{2m+l，(L-I) th}的方式预测在(L-本文档来自技高网...

【技术保护点】
一种磁共振动态成像方法，其特征在于，包括以下步骤：Ｓ１０：在当前帧图像，根据所需图像分辨率确定的Ｌ级发射尺度激励波形和小波激励波形；Ｓ２０：根据所述的尺度激励波形采集Ｌ级尺度系数并计算当前帧图像的Ｌ级尺度能量值，且根据所述小波激励波形采集Ｌ级小波系数并计算当前帧图像的Ｌ级小波能量值；Ｓ３０：产生队列，根据所述当前帧图像的Ｌ级小波能量值所对应的Ｌ级位置值预测Ｌ－１级位置值，把所述位置值放置所述队列；Ｓ４０：根据所述队列的位置值所对应的下一帧图像的位置值，发射小波激励波形并获取预测小波能量值，以及在下一帧图像的根据所需图像分辨率确定的Ｌ级发射尺度激励波形，并获得Ｌ级预测尺度能量值；Ｓ５０：当前帧图像的小波能量值与下一帧图像的预测小波能量值之差，判断是否大于预设判断值，是，则把所述位置值加入队列，否，则删除所述位置值；Ｓ６０：根据当前帧图像的Ｌ级尺度能量值和各级小波能量值，以及下一帧图像的Ｌ级预测尺度能量值和各级预测小波能量值，进行图像重建，并生成当前帧图像和下一帧图像。

【技术特征摘要】

【专利技术属性】
技术研发人员：寇波，刘新，郑海荣，谢国喜，邱本胜，吴垠，戴睿彬，
申请(专利权)人：中国科学院深圳先进技术研究院，
类型：发明
国别省市：94[中国|深圳]