K空间重建方法及装置发明专利.docxVIP

K空间重建方法及装置 技术领域 本发明涉及磁共振成像领域，尤其涉及一种K空间重建方法及装置。 背景技术 平面回波成像(echo planar imaging，EPI)是基于方向相反的频率读出梯度交替 采集MR信号的奇、偶回波。由于成像物体的磁敏感性、化学位移及磁场不均匀性等因 素引起的偏振效应，导致K空间数据在相位编码上产生相位调制。信号接收回路以及滤 波器的相位偏移会导致采集的奇、偶回波之间出现相位移动。此外，高速梯度的开关 场导致的涡流也会影响回波时序的准确性。所有这些因素导致经过傅立叶变换重建的 EPI图像在相位编码方向存在明显的伪影，严重破坏了图像的质量。参见图1，预想采 集数据如图1(a)所示，实际采集到的数据如图1(b)所示。在图1(b)中，粗线和细线分 成了两组，相互之间有一个移位。伪影可以通过优化扫描参数、改善滤波器的性能、 提高静磁场的均匀性、减小梯度涡流等措施得以控制。但是由于系统硬件性能的提高 是有限的，低强度的伪影几乎总是出现在EPI图像中。 为了克服这一缺点，人们提出了各种各样的基于后处理的伪影消除方法。 Buonocore在J.Magn Resorn Med2001，45(1)：96～108中提出了基于奇、偶重建的伪 影消除方法。但是该方法要求伪影与目标图像不完全重叠，而实际图像中目标和伪影 往往是重叠的，因此该方法在实际应用中效果并不十分理想。此外，该方法还假设奇、 偶相位差是按线性变化的，这进一步限制了它的应用。K.J Lee在J.Magn Reson Med 2002，47(4)：812～817中提出了基于凸集投影的方法来校正相位，但他假设目标图像 和含有伪影的图像相位差是按线性变化的，同时假设每一行的相位差是按同一参数变 化的。此外还有人提出用二阶矩的方法，但该方法要求图像是对称的，同时要求奇偶 相位差按线性变化，且所有行的相位差按同一参数变化。 目前比较常用的方法是采用相位恢复的方法消除伪影。在采集K空间数据时，共采 集两部分数据：应采集数据和校准数据，通过对采集到的校准数据进行分析，对应采 集数据进行校正，以消除伪影。但是，采集校准数据增加了额外的时间损耗，对于高 倍数加速的采集，采集校准数据所耗费的时间会大大限制成像速度。 发明内容 本发明要解决的问题是提供一种K空间重建方法及装置，以避免额外采集校准数 据，利用校准数据对采集的K空间数据进行校正。 为解决上述问题，本发明技术方案提供的K空间重建的方法，包括以下步骤： a，多次激发平面回波成像序列，并行采集K空间数据，所采集的数据包含至少两 个数据组，各数据组在频率编码方向上互成一定角度，且各数据组之间有重叠部分； b，将所述每个数据组内的相位编码线按照采集方向不同进一步划分为第一小组和 第二小组，所述第一小组内的相位编码线与所述第二小组内的相位编码线采集方向相 反； c，将各小组加入K空间固定，每小组加入时寻找与已加入K空间的其余小组相关 性最强的位置固定； d，利用上述步骤c的K空间数据进行欠采样数据填充，实现K空间的重建。 进一步地，所述步骤a中，所述并行采集K空间数据，所采集的数据包含至少两 个数据组的步骤中，所述各数据组对应的加速因子的值相同或不同。 进一步地，所述步骤a中，所述各数据组在频率编码方向上互成一定角度的步骤 中，所述角度为大于0度且小于180度。 进一步地，所述步骤a中，所述各数据组在频率编码方向上互成一定角度的步骤 中，所述角度为90度。 进一步地，所述步骤c中，所述每小组加入时寻找与已加入K空间的其余小组相 关性最强的位置固定，进一步包括以下步骤： c1，将一小组放入K空间中，横向和／或纵向移动该小组，使该小组与已加入K 空间的其余小组之间有重叠部分； c2，分别取出该小组与已加入K空间的其余小组之间的重叠部分，并将其统一插 值到同一公共网格空间，计算所述重叠部分的相似性值； c3，重复步骤c1和c2，以获得该小组在K空间中多个位置处所对应的与已加入K 空间的其余小组之间的重叠部分的相似性值； c4，将获得的多个相似性值进行对比，确定该小组与已加入K空间的其余小组相 关性最强的位置，将该小组加入K空间固定。 进一步地，所述步骤c2中，所述计算所述重叠部分的相似性值的步骤具体为：利 用方差估计方法或相关性估计方法计算所述重叠部分的相似性值。 进一步地，所述步骤c2中，计算所述重叠部分的相似性值之前，还包括以下步骤： 对所述取出的该小组与已加入K空间的其余小组之间的重叠部分分别取模值。 为了解决上述问题，本发明还提供了一种K空间成像装置，包括： 采集单元，适于通过多次激发平面回波成像序列，并行采集K空间数据，所采集 的数据包含至少两个数据组，各数据组在频率编码方向上互成一定角度，且各数据组 之间有重叠部分； 分组