[发明专利]一种消除磁共振图像中尖峰噪声的方法

权利要求书

1.一种消除磁共振图像中尖峰噪声的方法，其特征在于，包括：

步骤一：定位采集到的磁共振图像的k空间数据中的所有尖峰噪声所在的位置，设置所述尖峰噪声所在位置的信号强度值作为初始值，并设定算法迭代参数，包括能量函数的阈值、能量函数变化量的阈值以及最大迭代次数；

步骤二：将所述磁共振图像的K空间数据转换为图像域数据；

步骤三：对所述图像域数据进行小波变换得到小波域数据，并对所述图像域数据进行有限差分运算获得全变差数据；

步骤四：根据所述小波域数据与所述全变差数据计算能量函数值；

步骤五：判断所述能量函数值，若所述能量函数值大于所述能量函数值的阈值，且所述能量函数值的变化量大于所述能量函数值变化量的阈值，且迭代次数没有达到所述的最大迭代次数阈值，则利用优化算法改变所述尖峰噪声所在位置的信号强度的数值，并重复执行所述步骤二至步骤五进行迭代；否则就将当前的所述尖峰噪声所在位置的信号强度值作为修复值，完成尖峰噪声去除过程。

2.如权利要求1所述的消除磁共振图像中尖峰噪声的方法，其特征在于，利用图像的稀疏性来进行尖峰噪声的消除，所述能量函数值的计算方法如下：

E(I)＝||Ψ(I)||₁+αTV(I)

式中，I代表图像域数据，Ψ(I)表示所述小波域数据，TV(I)表示所述全变差数据，α表示权重分配。

3.如权利要求1所述的消除磁共振图像中尖峰噪声的方法，其特征在于，所述步骤五中，所述优化算法包括共轭梯度法。

4.如权利要求1所述的消除磁共振图像中尖峰噪声的方法，其特征在于，所述步骤五中进一步包括，当所述能量函数小于所述能量函数的阈值，或者当所述能量函数的变化量小于所述能量函数变化量的阈值，或者当迭代的次数超过所述最大迭代次数时，则停止迭代计算过程，将当前所述尖峰噪声所在位置的信号强度值作为修复所述的尖峰噪声所在位置信号强度的修复值。