[发明专利]用于执行并行磁共振成像的方法

权利要求书

1.一种身体的并行磁共振成像的方法，包括：

-从具有已知的或估算的灵敏度分布和噪声协方差矩阵的各个接收天线中获取所述身体的一组基础磁共振图像，所述基础图像在k-空间中欠采样；以及

-执行所述身体的磁共振图像的正则化重建；

其中，执行磁共振图像的正则化重建的所述步骤通过最小化成本函数而在离散的框架空间中进行，所述成本函数包括：

-误差项，在给定所获取的基础图像的情况下，所述误差项表示重建图像的似然；和

-框架惩罚项，所述框架惩罚项表示所述重建图像的框架系数的实际统计分布与所述系数的先验分布之间的偏差；

基于所述身体的辅助的磁共振图像估算所述重建图像的框架系数的先验分布。

2.根据权利要求1所述的方法，其中，在给定所述所获取的基础图像的情况下，所述误差项表示所述重建图像的负对数似然。

3.根据权利要求1或2所述的方法，其中，在给定所述所获取的基础磁共振图像和所述框架系数的先验分布的情况下，通过在所述框架空间中最大化限定所述身体的图像的一组框架系数的完全后验分布，而进行所述执行所述身体的磁共振图像的正则化重建的步骤。

4.根据前述任一项权利要求所述的方法，其中，根据所述所获取的基础磁共振图像重建所述身体的所述辅助的磁共振图像。

5.根据权利要求4所述的方法，其中，使用灵敏度编码-SENSE算法，根据所述所获取的基础磁共振图像重建所述身体的所述辅助的磁共振图像。

6.根据权利要求5所述的方法，其中，使用选自以下算法中的一种算法，根据所述所获取的基础磁共振图像重建所述身体的所述辅助的磁共振图像：

-未正则化的SENSE算法；

-在图像空间中正则化的SENSE算法；和

-在k-空间中正则化的SENSE算法。

7.根据前述任一项权利要求所述的方法，其中，假设所述框架系数的广义Gauss-Laplace先验统计分布，使用最大似然估算量或后验均值估算量，基于所述身体的所述辅助的磁共振图像估算所述分布的参数。

8.根据前述任一项权利要求所述的方法，其中，所述误差项为二次平均误差项。

9.根据前述任一项权利要求所述的方法，其中，所述所获取的基础图像为三维图像，且其中，在离散的三维框架空间中进行所述执行磁共振图像的正则化重建的步骤。

10.根据权利要求9所述的方法，其中，通过叠加待成像的对象的切片的二维基础图像来获得所述所获取的三维基础图像。

11.根据前述任一项权利要求所述的方法，其中，所述执行磁共振图像的正则化重建的步骤基于冗余小波框架表示。

12.根据权利要求1至10中任一项所述的方法，其中，所述执行磁共振图像的正则化重建的步骤基于非冗余小波表示。

13.根据前述任一项权利要求所述的方法，其中，所述成本函数还包括至少一个选自以下项的空间域惩罚项：

-所述重建图像的全变差范数；和

-凸约束。

14.一种执行身体的动态的并行磁共振成像的方法，包括：

-从具有已知的或估算的灵敏度分布和噪声协方差矩阵的各个接收天线中获取所述身体的一组时间序列的基础磁共振图像，所述基础图像在k-空间中欠采样；以及

-执行所述身体的时间序列的磁共振图像的正则化重建；

其中，通过最小化成本函数进行所述执行时间序列的基础磁共振图像的正则化重建的步骤，所述成本函数包括：

-误差项，在给定相应的所获取的基础图像的情况下，所述误差项表示各重建图像的似然；和

-时间惩罚项，所述时间惩罚项表示在所述序列的连续图像之间的像素间或体素间的差异。

15.根据权利要求14所述的方法，其中，所述时间惩罚项基于边缘保持函数。

16.根据权利要求15所述的方法，其中，所述时间惩罚项基于凸边缘保持函数。

17.根据权利要求16所述的方法，其中，所述时间惩罚项基于L_p范数，其中，p≥1且优选地1≤p<1.5。