[发明专利]一种磁共振去除运动伪影的方法及装置

权利要求书

1.一种磁共振去除运动伪影的方法，其特征在于，包括以下步骤：

1)获得若干通道全采集的原始K空间；

2)利用所述全采集的原始K空间，计算获得合并系数，利用所述合并系数，计算获得K₁空间；

3)将K₁空间的数据减去原始K空间的数据得到K₂空间，将K₂空间变换至图像域得到S₂图像，提取在计算K₁空间时由于算法不理想而引入的误差，得到S_2’图像；

4)将所述S_2’图像变换至K空间域，得到K₃空间数据，将K₁空间的数据减去K₃空间的数据，得到K₄空间数据；

5)判断是否进行下一次迭代计算，若是，则将K₄空间数据作为新的原始K空间数据，回到步骤2)；若否，则执行步骤6)；

6)将K₄空间变换至图像域，得到最终的图像。

2.如权利要求1所述的磁共振去除运动伪影的方法，其特征在于，所述步骤3)中提取在计算K₁空间时由于算法不理想而引入的误差的方法具体为，将所述S₂图像中所有大于第一阈值的像素点的值替换为0，得到S_2’图像；或将所述S₂图像中所有大于第一阈值的像素点的值替换为0，再进行均值滤波，得到S_2’图像。

3.如权利要求1所述的磁共振去除运动伪影的方法，其特征在于，所述步骤2)具体包括以下子步骤：

21)利用原始K空间数据求得并行采集重建方法的合并系数；

22)对每一通道的原始K空间按照相同的规则进行分组，获得若干分组K空间，所述分组K空间为与所述原始K空间相同大小的矩阵，所述分组K空间包括原始K空间中的部分已采集数据，K空间中的其他数据点作为待填补数据；

23)利用所述合并系数，以及分组K空间中的已采集数据，使用并行采集重建方法对分组K空间中的待填补数据进行填补；

24)将填补完成后的各分组K空间数据进行合并，获得K₁空间数据。

4.如权利要求3所述的磁共振去除运动伪影的方法，其特征在于，所述步骤22)具体为：对每一通道的原始K空间进行预处理分组，分为两个分组K空间，所述分组K空间为与所述原始K空间相同大小的矩阵，其中一个分组K空间的第奇数条数据线为原始K空间中第奇数条数据线，第偶数条数据线为待填补数据；另一个分组K空间的第偶数条数据线为原始K空间中第偶数条数据线，第奇数条数据线为待填补数据。

5.如权利要求3所述的磁共振去除运动伪影的方法，其特征在于，所述步骤24)具体为,计算填补完成后的各分组K空间数据的平均值，作为K₁空间数据。

6.如权利要求2所述的磁共振去除运动伪影的方法，其特征在于，所述第一阈值为，该次迭代计算获得的S₂图像中像素值的平均值。

7.如权利要求1所述的磁共振去除运动伪影的方法，其特征在于，所述步骤5)中，通过以下方法判断是否进行迭代计算：预先设定迭代次数，若没有达到预先设定的迭代次数，则继续进行迭代计算，将K₄空间数据作为新的原始K空间数据，回到步骤2)；若达到预先设定的迭代次数，则停止进行迭代计算，执行所述步骤6)。

8.如权利要求1所述的磁共振去除运动伪影的方法，其特征在于，所述步骤5)中，通过以下方法判断是否进行下一次迭代计算：设定第二阈值，若该次迭代获得的S₂图像的像素值的标准差大于等于所述第二阈值时，则继续进行迭代计算，将K₄空间数据作为新的原始K空间数据，回到步骤2)；若该次迭代获得的S₂图像的像素值的标准差小于所述第二阈值时，则停止进行迭代计算，执行所述步骤6)。

9.如权利要求8所述的磁共振去除运动伪影的方法，其特征在于，所述第二阈值为，第一次迭代计算时获得的S₂图像的像素值的标准差的30％。

10.一种磁共振去除运动伪影的装置，其特征在于，包括：

采集单元，用于获得若干通道的全采集的原始K空间；

K₁空间获得单元，与所述采集单元相连，用于利用所述全采集的原始K空间计算获得合并系数，利用所述合并系数，计算获得K₁空间；S_2’图像获得单元，与所述K₁空间获得单元相连，用于将K₁空间的数据减去原始K空间的数据，得到K₂空间，将K₂空间变换至图像域，得到S₂图像，提取在计算K₁空间时由于算法不理想而引入的误差，得到S_2’图像；

K₄空间获得单元，与所述S_2’图像获得单元相连，用于将所述S_2’图像变换至K空间域，得到K₃空间数据，将K₁空间的数据减去K₃空间的数据，得到K₄空间数据；

迭代计算判断单元，与所述K₁空间获得单元和K₄空间获得单元分别相连，用于判断是否进行下一次迭代计算；

图像获得单元，与所述迭代计算判断单元相连，用于将K4空间数据变换至图像域，得到最终的图像。