[发明专利]基于迭代投影的MRI图像重构方法

权利要求书

1.一种基于迭代投影的核磁共振MRI图像重构方法，包括如下步骤：

(1)对MRI图像x进行傅里叶变换，得到变换系数：x_f＝DFT(x)，用变密度欠采样观测矩阵Φ对图像x的傅里叶变换系数x_f进行观测采样，获得傅里叶域的观测值y＝Φx_f，对观测值y进行傅里叶反变换得到MRI图像x的初始值：x⁰＝DFT^-1(y)；

(2)结合MRI图像局部光滑特性和小波域稀疏的先验知识，将MRI图像x的初始值x⁰用全变差方法进行滤波，然后用小波域双变量阈值方法对该滤波结果x^TV进行优化，得到MRI图像x的优化结果x′；

(3)对优化结果x′进行投影，得到投影结果x^*＝x′+Φ^T(y-Φx′)，并判断x^*是否满足迭代终止条件：|D⁽ⁱ⁺¹⁾-D⁽ⁱ⁾|＜10^-5，其中N是图像的像素个数，i表示第i次迭代，若满足条件，将投影结果x^*作为MRI图像x的重构结果并输出，否则返回步骤(2)继续执行。

2.根据权利要求1所述的基于迭代投影的MRI图像重构方法，其中步骤(2)所述的用小波域双变量阈值方法对全变差滤波结果x^TV进行优化，按如下步骤进行：

2.1)选取双变量阈值为：Threshold(f,λ)=(f2+fP2-λ3σ(i)σξ)f2+fP2·f,]]>λ是控制算法收敛性的常数，取λ＝20，f＝Ψx^TV是x^TV的小波变换系数，Ψ表示小波变换，i表示低i次迭代，f_p是f的父系数，如果则取阈值为0，是系数ξ3×3邻域的边缘方差估计，σ⁽ⁱ⁾是中值估计函数，

2.2)根据确定的双变量阈值Threshold(f，λ)选择比阈值大的小波变换系数f_s保留，其余的系数置0，再对保留系数f_s进行小波反变换就得到对x^TV的优化结果x′＝Ψ^-1f_s。