专利名称：基于图像域叠加的propeller磁共振数据重建方法
技术领域：
本发明涉及磁共振数据重建技术领域，具体来说涉及一种基于图像域叠加的 PROPELLER磁共振数据重建新方法。
背景技术：
磁共振成像(Magnetic Resonance Imaging, MRI)是当前临床医学影像学的重 要检查手段之一，由于数据采集时间比较长，成像易受运动干扰，如何有效地克服运动对成 像的影响，一直以来是磁共振成像研究的热点与技术难题之一。PROPELLER (PeriodicalIy Rotated Overlapping ParallEL Lines with Enhanced Reconstruction,螺方宠奖)(图 1) 方法能够利用K空间中心重叠采样区域的数据来估计采集过程中受检查者的运动信息，从 而实现K空间条的运动补偿，对于运动伪影有很好的消除效果。完成运动补偿后，PROPELLER数据的重建是个典型的非笛卡尔数据的重建问题，非 笛卡尔数据的重建一直是磁共振重建研究的热点问题之一，直接求和(DirectSummation) 傅立叶变换，一般被认为可以较高精度的实现图像重建，但是由于计算复杂度很高，很难 推广到实际应用中。目前应用于PROPELLER采样数据重建的主要算法为卷积插值网格化 (Gridding)算法，但由于这种算法的点扩展函数含有旁瓣，导致重建结果中含有混叠现象。根据PROPELLER数据采样的独特性和傅里叶变换的线性性质，并参考FBP重建算 法，提出本发明所涉及的基于图像域叠加的PROPELLER重建算法，这种算法与Gridding算 法的本质区别在于，Gridding算法在K空间域进行插值，而本文算法在图像域进行插值，所 以该算法简称为iPROPELLER算法，i代表图像空间(image)。

发明内容
本发明的目的在于提供一种基于图像域叠加的PROPELLER磁共振数据重建新方 法，该磁共振数据重建新方法可以有效的避免混迭伪影的干扰，获得更高的图像重建精度。本发明的目的可通过以下的技术措施来实现一种基于图像域叠加的PROPELLER磁共振数据重建新方法，其特征在于包括以下 步骤(1)基于PROPELLER采集方式采集磁共振数据；(2)对每个K空间条进行密度补偿后进行二维逆快速傅里叶变换重建出临时图 像；(3)在图像域对临时图像进行运动补偿和旋转校正；(4)将校正后的临时图像叠加得到最后的重建结果。所述步骤⑴中的K空间采用个数和每个K空间条采样点数设定为每10度采集 一个K空间条，共采集18个，每个K空间条采集24行相位编码线，每条相位编码线采样256 个数据，信号叠加平均次数为1。本基于图像域叠加的PROPELLER磁共振数据重建新方法与现有技术相比，能够明显减少混迭现象，使重建图像更为清晰。


图1是本发明的基于图像域叠加的PROPELLER磁共振数据重建新方法的图像处理 过程的示意图；图2是PR0PELLERMRI的K空间采样轨迹图；图3 (a)是本发明方法进行图像处理的仿真的Sh印p-Logan体模图像；图3(b)是采用传统DFT算法重建的图像显示结果；图3(c)是采用传统Gridding算法重建的图像显示结果；图3(d)是采用本发明方法的重建图像显示结果；图3(e)是图3(b)与图3(a)的差值图像；图3(f)是图3(c)与图3(a)的差值图像；图3(g)是图3(d)与图3(a)的差值图像。
具体实施例方式本发明的基于图像域叠加的PROPELLER磁共振数据重建新方法的具体实施过程 如图1所示，包括四个步骤如下步骤1数据采集。按照图2的设计采样方式，本发明的实验中K空间采样个数和 每个K空间条采样点数设定为每10度采集一个K空间条，共采集18个，每个K空间条采 集24行相位编码线，每条相位编码线采样256个数据，信号叠加平均次数为1。将采样数据 存入一 3维数组。步骤2用采集的K空间条数据分别重建临时图像。将采集到的18个水平K空间 条分别进行采样密度补偿，然后补零到笛卡尔坐标空间中，即将24*256的矩阵放在K空间 的中间位置，周围补零，成为256*256的K空间矩阵，然后通过逆FFT变换重建出临时图像。步骤3图像域运动补偿和旋转校正。通过频域相关或图像域优化搜索等方法估计 采集每个K空间条时，成像物体的层内旋转运动参数Δ θ与水平与竖直方向上平移运动参 数(ΔΧ，Δ y)，在图像域直接进行相应补偿。对于旋转补偿，将重建出的临时图像沿相反方 向旋转Δ θ即可。对于平移补偿，只需根据(ΔΧ，Δ y)，将旋转补偿后的子图像做相反的 方向平行移动，最终使该子图像与参考图像保证空间位置上一一对准关系即可。层内运动 补偿后，还可能存在一些由于层间运动，或弹性变形造成的不能补偿的运动，可通过相关加 权方法予以补偿，相关系数的确定可以通过计算子图像之间的相关系数获得。另外采集到的各个K空间条都要放在第一个K空间条的位置上然后补零重建，相 当于频域旋转了一定的角度，重建的临时图像需要逆旋转同样的角度进行补偿。以采集到 的第二个K空间条为例，设它与作为参考的第一个K空间条的角度为θ，如果把第二个K空 间条的数据放在第一个K空间条的位置上，相当于把它沿顺时针方向旋转了 θ，进行傅立 叶变换后，根据傅立叶变换旋转不变定理(频域的旋转等同与图像域同样角度的旋转)，要 将重建出的临时图像沿逆时针方向旋转θ，才能进入下一步骤的处理中。步骤4图像域叠加。将步骤3得到的校正后的临时图像直接叠加起来，就得到了 最终的重建结果。
为了将传统图像重建方法与本方法做比较，因此对图3(a)所示的体模图像分别 采用DFT算法、Gridding算法和本方法的进行三种算法的重建处理，最后分析比较各种重 建方法的效果。图3(b)至图3(d)示出了 PROPELLER采样轨迹参数为(18，24，256)时，分别应用 DFT、Gridding与本发明算法的最终重建图像。图(a)为用于仿真的Sh印p-Logan体模图 像；图(b)为DFT算法重建结果，显然图像周围存在明显的混迭；图(c)为Gridding算法 重建结果，图像内部也可见混迭现象；图(d)本发明所涉及iPROPELLER算法重建结果，很明 显，基本上没有任何混迭现象。图3(e)、图3(f)、图3(g)分别给出了三种重建结果与真实 参考图像图3(a)之间的差值图像，对于重建中图像内部混迭现象的显示更为清楚。本发明的实施方式不限于此，在本发明上述基本技术思想前提下，按照本领域的 普通技术知识和惯用手段对本发明内容所做出其它多种形式的修改、替换或变更，均落在 本发明权利保护范围之内。
权利要求
一种基于图像域叠加的PROPELLER磁共振数据重建方法，其特征在于包括以下步骤(1)基于PROPELLER采集方式采集磁共振数据；(2)对每个K空间条进行密度补偿后进行二维逆快速傅里叶变换重建出临时图像；(3)在图像域对临时图像进行运动补偿和旋转校正；(4)将校正后的临时图像叠加得到最后的重建结果。
2.根据权利要求1所述的基于图像域叠加的PROPELLER磁共振数据重建方法，其特征 在于所述步骤(1)中的K空间采样个数和每个K空间条采样点数设定为每10度采集一 个K空间条，共采集18个，每个K空间条采集24行相位编码线，每条相位编码线采样256 个数据。
全文摘要
本发明公开了一种基于图像域叠加的PROPELLER磁共振数据重建方法，包括以下步骤(1)基于PROPELLER采集方式采集磁共振数据；(2)对每个K空间条进行密度补偿后进行二维逆快速傅里叶变换重建出临时图像；(3)在图像域对临时图像进行运动补偿和旋转校正；(4)将校正后的临时图像叠加得到最后的重建结果。的特点。该磁共振数据重建新方法可以有效的避免混迭伪影的干扰，获得更高的图像重建精度。
文档编号G01R33/56GK101995561SQ201010517528
公开日2011年3月30日 申请日期2010年10月22日 优先权日2010年10月22日
发明者冯衍秋, 陈武凡 申请人:南方医科大学