[发明专利]一种MRI图像重构方法

权利要求书

1.一种磁共振成像图像重构方法，包括：

S1、采用非均匀傅立叶变换对磁共振成像信号进行插值处理；

S2、采用l₁范数和全变分混合正则化项建立压缩感知模型；

其中，‖m‖_TV表示全变分正则化项，m表示磁共振成像图像，S表示磁共振成像信号，表示经欠采样所得的测量数据，m表示由恢复得到的图像；α和β为正则化参数，表示平衡两个正则化项在目标函数中所占的比重；Φ为感知矩阵；

S3、基于交替方向乘子法重构磁共振成像图像；具体为：由步骤S2得到的优化问题为：

其中，‖·‖₁表示求l₁范数，表示求l₂范数，z＝Dm表示图像的空间域的有限差分，D为全变分算子，ρ₁、ρ₂是添加的惩罚项的惩罚参数，u₁、u₂是拉格朗日乘子的函数；

通过交替最小化方法求解，分别优化变量m、z、u₁以及u₂，得到重构磁共振成像图像。

2.根据权利要求1所述的磁共振成像图像重构方法，其特征在于，所述步骤S3所述分别优化m，z，u₁和u₂变量，通过交替最小化方法求解得到重构磁共振成像图像，具体为：

A1、通过固定变量z，u₁和u₂，得到求解变量m第i+1次的值的子优化问题：

通过固定变量m，u₁和u₂，得到求解变量z第i+1次的值的子优化问题：

A2、根据下式利用得到的mⁱ⁺¹和zⁱ⁺¹更新变量u₁和u₂：

最后，通过循环执行步骤A1、A2，直到满足迭代终止条件，得到重建图像信号m。

3.根据权利要求2所述的磁共振成像图像重构方法，其特征在于，

所述迭代终止条件为：

|SSIM(mⁱ)-SSIM(m^i-1)|＜ε₃；

其中，SSIM(·)表示求结构相似度值运算，ε₃表示迭代终止阈值。