[发明专利]在多梯度回波MRI中对场波动的回顾性校正

权利要求书

1.一种存储指令的非瞬态存储介质，所述指令能由计算机(22)读取和运行以执行磁共振(MR)图像重建方法，所述磁共振图像重建方法包括：

执行对多梯度回波MR成像数据的迭代重建以生成重建的MR图像(36)，所述迭代重建使用将所述MR成像数据与所述重建的MR图像进行关联的模型(30)，其中，所述模型包括参数化磁场波动分量(32)；并且

在执行所述迭代重建期间，更新所述模型的所述参数化磁场波动分量的参数以优化成本函数(40)，所述成本函数取决于所述重建的MR图像相对于所述模型的所述参数化磁场波动分量的所述参数的偏导数。

2.根据权利要求1所述的非瞬态存储介质，其中，所述参数化磁场波动分量(32)包括其中，T_E是在k空间位置k(t)处的所述MR成像数据的回波时间，并且ω(r，k(t))包括所述参数化磁场波动分量的参数。

3.根据权利要求2所述的非瞬态存储介质，其中，将所述MR成像数据与重建的MR图像进行关联的所述模型(30)包括：

其中，f(r，T_E)是所述重建的MR图像，s_c(r)及其复共轭对MR接收线圈敏感性进行建模，FT表示傅里叶变换，FT^-1表示傅里叶逆变换，c是MR成像系统的通道的数目，并且是由磁场波动损坏的所述重建的MR图像。

4.根据权利要求1-3中的任一项所述的非瞬态存储介质，其中，所述参数化磁场波动分量(32)将磁场波动的空间变化建模为沿着三个正交空间轴中的每个轴的线性变化。

5.根据权利要求1所述的非瞬态存储介质，其中，所述参数化磁场波动分量(32)将所述磁场波动建模为空间均匀的时间磁场波动，并且所述参数化磁场波动分量包括其中，T_E是在k空间位置k处的所述MR成像数据的回波时间，并且ω(k)包括所述参数化磁场波动分量的参数。

6.根据权利要求5所述的非瞬态存储介质，其中，将所述MR成像数据与重建的MR图像进行关联的所述模型(30)包括：

其中，是所述MR成像数据，对MR接收线圈敏感性进行建模，FT^-1表示傅里叶逆变换，c是MR成像系统的通道的数目，并且f(r，T_E，ω(k))是由磁场波动损坏的所述重建的MR图像。

7.根据权利要求1-3中的任一项所述的非瞬态存储介质，其中，所述参数化磁场波动分量(32)将所述磁场波动建模为在生成所述多梯度回波MR成像数据的MR成像数据采集的单个重复时间(TR)内在时间上恒定。

8.根据权利要求1-3中的任一项所述的非瞬态存储介质，其中，所述成本函数(40)包括总变差(TV)成本函数。

9.根据权利要求1-3中的任一项所述的非瞬态存储介质，其中，所述更新包括应用Broyden-Fletcher-Goldfarb-Shanno(BFGS)算法来更新所述模型(30)的所述参数化磁场波动分量(32)的所述参数，从而优化所述成本函数(40)。

10.根据权利要求1-3中的任一项所述的非瞬态存储介质，其中，所存储的指令还能由所述计算机(22)读取和运行以执行以下操作中的一项或多项：

使用所述重建的MR图像(36)来执行R₂*映射以生成R₂*映射图(50)，

使用所述重建的MR图像来执行磁敏感加权成像(SWI)以生成SWI图像(52)，以及

使用所述重建的MR图像来执行定量磁敏感映射(QSM)以生成QSM映射图(54)。