K空间重建方法及磁共振成像方法发明专利.docxVIP

K空间重建方法及磁共振成像方法 技术领域 本发明涉及磁共振成像领域，尤其涉及一种K空间重建方法及磁共振成像方法。 背景技术 在磁共振成像技术中，成像的速度是衡量成像方法的一个很重要标准。限制成像 速度的很重要因素是数据采集以及k空间数据填充。一般的数据采集方式要采满k空间后才 能进行重建得到图像，磁共振并行采集重建技术利用线圈重组合并的方式，采集部分K空间 数据，然后对欠采样的K空间进行数据填补，最后利用填补完整的k空间数据进行图像重建。 采用并行采集技术，无需对整个K空间的数据进行采集，耗时少，提高了磁共振成像的速度。 图1为4个通道(线圈1,2,3,4)拟合的GRAPPA方法的采样图。传统的GRAPPA算法如 图1，黑色点代表实际采集的K空间数据，也可称为应采集数据，白色点代表欠采样需要填补 的数据，灰色点代表为了计算线圈合并系数，而适量全采的K空间数据，也可称为校准数据。 图中任意一个白色点可以表示为周围黑色点的线性叠加，相当于对多个线圈的数据进行了 合并，而合并系数n?ij(表示第i个线圈，第j个位置)可以通过黑色点拟合灰色点来确定，图1 中从某个灰色点辐射出去到其他每个黑色点之间的线表示拟合灰色点的关系。线圈合并系 数确定后，白色点可根据求得的线圈合并系数将各线圈的数据合并后计算得到。 虽然这种填补方式可以将成像速度提高，但是由于合并系数的精度会受到多种因 素的影响(包括噪音)，导致拟合的结果有一定的误差，使得填补的数据与实际值产生偏差， 导致重建后图像产生卷折等伪影，影响成像质量(在加速因子较大时更为明显)。所以，如果 能够克服噪声等因素的影响，减小线圈合并系数的误差，会对成像的质量带来很大的改进。 发明内容 本发明要解决的问题是提供一种K空间重建方法及装置，以以减小线圈合并系数 计算过程中，噪音等不理想因素的干扰。 为解决上述问题，本发明提供的K空间重建方法，包括以下步骤： a，并行采集K空间数据，所采集数据包括应采集数据和校准数据； b，提取所述校准数据形成第一数据组，对所述第一数据组执行像素合并操作获得 第二数据组，所述第一数据组与第二数据组的扫描视野相同； c，利用所述第二数据组计算线圈合并系数，基于所述线圈合并系数对欠采样的K 空间进行填补，以获得满采样的K空间。 进一步地，所述步骤b中，对所述第一数据组执行像素合并操作获得第二数据组, 所述第一数据组与第二数据组的扫描视野相同的步骤，进一步包括： b11,在所述第一数据组中，选取至少一个维度； b12,将所述第一数据组变换至所选各维度的时域，在各维度的时域中，对所述第 一数据组执行卷积操作，得到第二数据组； b13,将所述第二数据组由时域变换至K空间域。 进一步地，所述步骤b12和步骤b13中，所述变换为傅里叶变换。 进一步地，所述步骤b12之前，还包括：在各维度的时域中，对所述第一数据组的相 位执行去线性操作。 进一步地，所述在各维度的时域中，对所述第一数据组的相位执行去线性操作的 步骤之前，还包括：在各维度的时域中，对所述第一数据组取模值。 进一步地，所述步骤b中，对所述第一数据组执行像素合并操作获得第二数据组, 所述第一数据组与第二数据组的扫描视野相同的步骤，进一步包括： b21,在所述第一数据组中，选取至少一个维度； b22，在所选维度的K空间域，对所述第一数据组执行点乘操作以得到第二数据组。 基于上述K空间重建方法，本发明还提供了一种磁共振成像方法，包括以下步骤： 以上述的K空间重建方法对K空间进行重建； 将重建之后的K空间变换到图像域以获得图像。 本发明提供的K空间重建方法，在利用校准数据计算线圈合并系数之前，对校准数 据执行不改变扫描视野的像素合并操作，提高了数据的信噪比，优化了线圈合并系数的计 算精度，减小了由拟合误差引起的图像伪影。 附图说明： 图1是4个通道(线圈1,2,3,4)拟合的GRAPPA方法的采样示意图； 图2是本发明提供的K空间重建方法的流程示意图； 图3是本发明实施例一的K空间数据采样示意图； 图4是利用本发明提供的K空间重建方法进行K空间重建时的数据变化示意图； 图5是本发明实施例一中一种C?N操作的实现过程的示意图； 图6a是利用传统GRAPPA的方法，在加速因子为4时，填补缺失k空间数据后重建的 图像； 图6b是利用本发明实施例一的方法提供的K空间重建方法，在加速因子为4时，填 补缺失k空间数据后重建的图像。 具体实施方式： 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和实施例 对本发明的具体实施方式做详细的说明。 图2是本发明提供的K空间重建方法的流程示意图。请参考图2，所述K空间重建方 法，包括以下步骤： a,并行采集K空间数据，所采集数据包括应采集数据和校准数据； b,提取所述校准数据形成第