一种磁共振去除运动伪影的方法及装置的制造方法
【专利说明】一种磁共振去除运动伪影的方法及装置 【技术领域】
[0001] 本发明涉及磁共振成像（MRl，Magnetic Resonance Imaging)技术领域，尤其涉及 一种磁共振去除运动伪影的方法及装置。 【【背景技术】】
[0002] 在磁共振成像技术中，磁共振信号空间称为K空间，即为傅里叶变换空间，将磁共 振扫描采集获得的磁共振信号填入到K空间中，将K空间数据经傅里叶逆变换即可得到磁 共振图像。
[0003] 运动干扰一直以来是制约磁共振成像质量的一个棘手问题，特别是在对某些部位 的扫描中，例如腹部扫描、心脏扫描，人体的呼吸或心脏的跳动会造成运动伪影，从而影响 磁共振成像质量。
[0004] 当前有一些磁共振去除运动伪影的方法，如PROPELLER方法，可以在一定程度上 缓解运动对成像造成的干扰，去除运动伪影，但是这些方法都是针对平面内的刚体运动，对 于有些部位的扫描，影响成像质量的主要因素是非刚体运动或平面外的运动，这些方法无 能为力。
[0005] 还有一些方法，如流动补偿技术，可以补偿血液流动造成的伪影，但是这些方法局 限于某一种特定的情况下才能应用（即血液流动造成的伪影）。
[0006] 总之，现有的磁共振成像抑制运动伪影的方法均有较大的局限性。
[0007] 因此，需要提出一种新的磁共振去除运动伪影的方法及装置，可以有效抑制各种 情况造成的运动伪影，提高磁共振成像质量。 【
【发明内容】
】
[0008] 本发明解决的问题是现有的磁共振去除运动伪影的方法无法有效抑制各种情况 造成的运动伪影的问题。
[0009] 为了解决上述问题，本发明提出一种磁共振去除运动伪影的方法，包括以下步 骤：
[0010] 1)获得若干通道全采集的原始K空间；
[0011] 2)利用所述全采集的原始K空间，计算获得合并系数，利用所述合并系数，计算获 得心空间；
[0012] 3)将Ki空间的数据减去原始K空间的数据得到K2空间，将K2空间变换至图像域 得到S 2图像，提取部分在计算&空间时由于算法不理想而引入的误差，得到S2，图像；
[0013] 4)将所述S2，图像变换至K空间域，得到K3空间数据，将Ki空间的数据减去K 3空 间的数据，得到K4空间数据；
[0014] 5)判断是否进行下一次迭代计算？若是，则将K4空间数据作为新的原始K空间数 据，回到步骤2);若否，则执行步骤6);
[0015] 6)将1(4空间变换至图像域，得到最终的图像。
[0016] 可选地，所述步骤3)中提取部分在计算1空间时由于算法不理想而引入的误差 的方法具体为，将所述s 2图像中所有大于第一阈值的像素点的值替换为0,得到s2，图像；或 将所述s2图像中所有大于第一阈值的像素点的值替换为0,再进行均值滤波，得到s 2，图像。
[0017] 可选地，所述步骤1)具体为，采集获得若干通道的全采集的原始K空间；或采集获 得若干通道的欠采集K空间，并利用相关方法获得全采集的原始K空间；
[0018] 可选地，所述相关方法为并行采集重建方法或半傅里叶重建方法。
[0019] 可选地，所述步骤2)具体包括以下子步骤：
[0020] 21)利用原始K空间数据求得合并系数；
[0021] 22)对每一通道的原始K空间按照相同的规则进行分组，获得若干分组K空间，所 述分组K空间为与所述原始K空间相同大小的矩阵，所述每一分组K空间都包括原始K空 间中的部分已采集数据，K空间中的其他数据点作为待填补数据；
[0022] 23)利用所述合并系数，以及分组K空间中的已采集数据，对每一分组K空间中的 待填补数据进行填补；
[0023] 24)将填补完成后的各分组K空间数据进行合并，获得&空间数据。
[0024] 可选地，所述步骤21)中，所述合并系数为并行采集重建方法的合并系数。
[0025] 可选地，所述步骤23)中，使用并行采集重建方法对每一分组K空间中的待填补数 据进行填补。
[0026] 可选地，所述并行采集重建方法为GRAPPA方法或SPIRiT方法。
[0027] 可选地，所述步骤22)具体为：对每一通道的原始K空间进行预处理分组，分为两 个分组K空间，所述分组K空间为与所述原始K空间相同大小的矩阵，其中一个分组K空间 的第奇数条数据线为原始K空间中第奇数条数据线，第偶数条数据线为待填补数据；另一 个分组K空间的第偶数条数据线为原始K空间中第偶数条数据线，第奇数条数据线为待填 补数据。
[0028] 可选地，所述步骤24)具体为：计算填补完成后的各分组K空间数据的平均值，作 为I空间数据。
[0029] 可选地，所述第一阈值为，该次迭代计算获得的S2图像中像素值的平均值。
[0030] 可选地，所述步骤5)中，通过以下方法判断是否进行迭代计算：预先设定迭代次 数，若没有达到预先设定的迭代次数，则继续进行迭代计算，将K 4空间数据作为新的原始K 空间数据，回到步骤2);若达到预先设定的迭代次数，则停止进行迭代计算，执行所述步骤 6)。
[0031] 可选地，所述步骤5)中，通过以下方法判断是否进行下一次迭代计算：设定第二 阈值，若该次迭代获得的S 2图像的像素值的标准差大于等于所述第二阈值时，则继续进行 迭代计算，将K4空间数据作为新的原始K空间数据，回到步骤2);若该次迭代获得的S2图 像的像素值的标准差小于所述第二阈值时，则停止进行迭代计算，执行所述步骤6)。
[0032] 可选地，所述第二阈值为，第一次迭代计算时获得的S2图像的像素值的标准差的 30%。
[0033] 可选地，通过傅里叶逆变换将K空间变换至图像域，通过傅里叶变换将图像变换 至K空间域。
[0034] 本发明还提出一种磁共振去除运动伪影的装置，包括：
[0035] 采集单元，用于获得若干通道的全采集的原始K空间；
[0036] &空间获得单元，与所述采集单元相连，用于利用所述全采集的原始K空间计算获 得合并系数，利用所述合并系数，计算获得I空间；
[0037] S2，图像获得单元，与所述&空间获得单元相连，用于将&空间的数据减去原始K 空间的数据，得到K2空间，将K2空间变换至图像域，得到S2图像，提取部分在计算Ki空间时 由于算法不理想而引入的误差，得到S2，图像；
[0038] K4空间获得单元，与所述&，图像获得单元相连，用于将所述S2，图像变换至K空间 域，得到K3空间数据，将&空间的数据减去K3空间的数据，得到K4空间数据；
[0039] 迭代计算判断单元，与所述&空间获得单元和K4空间获得单元分别相连，
[0040] 用于判断是否进行下一次迭代计算；
[0041] 图像获得单元，与所述迭代计算判断单元相连，用于将K4空间数据变换至图像 域，得到最终的图像。
[0042] 本发明对比现有技术有如下的有益效果：本发明利用迭代逼近的思想，将K空间 中数据线之间由于运动造成的差异进行抑制甚至去除，留下数据线之间共有的特性，去除 数据线中由于运动造成的差异，使得重建获得的图像中的运动伪影得到抑制，并且对各种 情况造成的运动伪影均有较好的抑制效果，从而提高磁共振成像质量。 【【附图说明】】
[0043]图1是本发明的磁共振去除运动伪影的方法的流程图；
[0044] 图2是GRAPPA方法的示意图；
[0045] 图3是本发明的磁共振去除运动伪影的方法中步骤S02的具体过程的示意图；
[0046] 图4是未使用本发明的方法重建图像和使用本发明的方法重建图像的对比；
[0047] 图5是本发明的磁共振去除运动伪影的装置的示意图。 【【具体实施方式】】
[0048] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明 的【具体实施方式】做详细的说明。
[0049]如图1所示，本发明的磁共振去除运动伪影的方法具体包括以下步骤：
[0050] 执行步骤S01，获得若干报道的全采集的原始