基于高阶奇异值分解的磁共振图像去噪方法技术

一种基于高阶奇异值分解的磁共振图像去噪方法，包括：（1）对原始图像进行方差稳定变换；（2）对变换后的图像进行第一次高阶奇异值分解去噪，具体是：（a）对每一个标的像素对应的参考块，逐个寻找该参考块的相似块组成高维数组；（b）对高维数组进行高阶奇异值分解变换得到系数和自适应基；（c）通过第一阈值操作；（d）进行高阶奇异值分解逆变换；（3）加权平均进行像素合并；（4）加权平均得到第一次加权平均后的图像；（5）对第一次加权平均后的图像进行第二次高阶奇异值分解去噪获得第二次去噪后的图像；（6）对第二次去噪后的图像进行方差稳定逆变换得到滤波图像。本发明专利技术的方法能够有效地抑制磁共振图像噪声，提高图像质量。

【技术实现步骤摘要】
基于高阶奇异值分解的磁共振图像去噪方法
本专利技术属于医学图像处理的
，具体地说属于一种基于高阶奇异值分解的磁共振图像去噪方法。
技术介绍
磁共振成像是当前临床医学影像学的重要检查手段之一。然而，由于成像机制的影响，成像过程中不可避免的会引入噪声。图像中的噪声会大大降低图像的质量，使得图像边缘变得模糊，细微结构难以辨认，从而影响临床诊断，降低分析任务的可信度，例如图像配准，图像分割和一些相关组织的参数测量(如灌注图像和弛豫时间的相关参数)。因此，对于临床诊断和图像分析来说，降低噪声是十分必要的。降低噪声可以通过两种方法来实现，一种是增加采集次数，平均多次采集的数据来提高图像的信噪比。这种方法是以增加采集时间来提高信噪比的。另一种方法是通过图像后处理滤波技术来提高信噪比。图像滤波技术不需要增加采集时间，因而在临床研究中得到广泛应用。近期，有研究者采用高阶奇异值分解方法(Rajwade，A.，Rangarajan，A.，Banerjee，A..Imagedenoisingusingthehigherordersingularvaluedecomposition.IEEETransPatternAnalMach2013，Intell35，849-862)对加有高斯噪声的自然图像进行去噪。该方法是一种简单的，基于块匹配，机器学习的稀疏去噪方法。该方法是将相似块组成高维数组，利用相似块之间的稀疏性进行去噪。该方法针对不同的高维数组，通过训练学习得到不同的自适应基，能够更好地表达图像内容。不同于奇异值分解，高阶奇异值分解不需要将高维数组展开成矩阵来进行分解，不会破坏块内部的相关信息。然而，这种通过图像本身学习得到的基，易受图像中噪声的影响尤其是当噪声严重的时候。因此，针对现有技术不足，提供一种对加有莱斯噪声的磁共振图像进行去噪处理的方法以克服现有技术不足甚为必要。
技术实现思路
本专利技术的目的在于针对现有技术不足，提供一种基于高阶奇异值分解的磁共振图像去噪方法，具有去噪性好，且能够提高去噪处理后的图像质量。本专利技术的上述目的通过如下技术手段实现。一种基于高阶奇异值分解的磁共振图像去噪方法，依次包括如下步骤：(1)对加有莱斯噪声的原始磁共振图像进行方差稳定变换，将依赖信号分布的莱斯噪声变成独立于信号分布的噪声，得到变换后的噪声图像；(2)将变换后的噪声图像中的每一个像素作为标的像素或者按照一定距离间隔取像素作为标的像素，取标的像素及其周围的像素作为参考块，执行以下操作进行第一次高阶奇异值分解去噪，具体是：(a)通过k最近邻方法依次逐个寻找参考块的相似块，组成高维数组；(b)对高维数组进行高阶奇异值分解变换，得到对应的系数和自适应基；(c)通过第一阈值操作，将绝对值小于第一阈值的系数置零得到修正后的系数；(d)根据自适应基与修正后的系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的估计值；(3)对步骤(2)得到的估计值加权平均进行像素合并，得到像素合并后的估计值，对每一个像素进行像素合并构成第一次去噪后的图像；(4)对方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均，得到第一次加权平均后的图像；(5)对第一次加权平均后的图像进行第二次高阶奇异值分解去噪，具体是：(5.1)将第一次加权平均后的图像中的每一个像素作为二次标的像素或者按照一定距离取像素作为二次标的像素，取二次标的像素及其周围的像素作为二次参考块，执行以下操作：(a)以第一次去噪后的图像进行相似块判断，通过k最近邻方法依次逐个寻找二次参考块的二次相似块，组成二次高维数组；(b)对二次高维数组进行高阶奇异值分解变换，得到对应的二次系数和二次自适应基；(c)通过二次阈值操作，将绝对值小于二次阈值的二次系数置零得到修正后的二次系数；(d)根据二次自适应基与修正后的二次系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的二次估计值；(5.2)对步骤(5.1)得到的二次估计值加权平均进行像素合并，得到像素合并后的二次估计值，对每一个像素进行像素合并构成第二次去噪后的图像；(6)对第二次去噪后的图像进行方差稳定逆变换，得到最终的滤波图像。优选的，上述步骤(2)中(c)分步骤的第一阈值大小其中p为块大小，d为图像维度，K为相似块的个数，σ为方差稳定变换后图像的标准差。优选的，上述步骤(4)具体是采用迭代正则化方法将方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均。优选的，上述步骤(4)对方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均具体为：其中表示第一次去噪后的图像，Y表示方差稳定变化后的噪声图像，ξ∈[0，1]是松弛参数，表示噪声图像Y对图像Yw的贡献。优选的，上述步骤(5)中(c)分步骤的二次阈值大小其中γ∈[0，1]为尺度参数，控制图像Yw的标准差σw大小。优选的，上述步骤(2)的(a)中，通过k最近邻方法依次逐个寻找参考块的相似块，相似块的块大小和相似块个数的取值根据噪声水平而定；当噪声水平不大于2％时，相似块的块大小为3，相似块个数为30-70；当噪声水平大于2％而小于等于7％时，相似块的块大小为4，相似块个数为70-80；当噪声水平大于7％而小于等于15％时，相似块的块大小为5，相似块个数为100-130。优选的，上述步骤(5)的(a)中，通过k最近邻方法依次逐个寻找二次参考块的二次相似块，二次相似块的块大小和二次相似块个数的取值根据噪声水平而定；当噪声水平不大于2％时，二次相似块的块大小为3，二次相似块个数为30-70；当噪声水平大于2％而小于等于7％时，二次相似块的块大小为4，二次相似块个数为70-80；当噪声水平大于7％而小于等于15％时，二次相似块的块大小为5，二次相似块个数为100-130。本专利技术基于高阶奇异值分解的磁共振图像去噪方法，该方法通过两步高阶奇异值分解去噪。第一步去噪后，对方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均，对加权平均后的图像进行第二步高阶奇异值分解去噪。与原始噪声图像相比，加权平均后的图像含有的噪声变小，第二步去噪可以减少噪声对高阶奇异值分解所得基的影响，使得高阶奇异值分解后的基能够更好地表达图像内容。故本专利技术方法在去噪的同时可以很好的保留图像细节，提高图像质量。附图说明利用附图对本专利技术作进一步的说明，但附图中的内容不构成对本专利技术的任何限制。图1为本专利技术基于高阶奇异值分解的磁共振图像去噪方法的流程示意图；图2为加入噪声水平为15％莱斯噪声的T1加权、T2加权和质子密度(PD)加权图像使用本专利技术方法去噪后的结果示意图；图3为对噪声水平为3％的真实数据，采用本专利技术方法去噪后的结果示意图；图4是对应图3的局部放大图；图5为对噪声水平为4.5％的真实数据，采用本专利技术方法去噪后的结果示意图；图6是对应图5的局部放大图。具体实施方式结合以下实施例对本专利技术作进一步描述。实施例1。如图1所示，一种基于高阶奇异值分解的磁共振图像去噪方法，依次包括如下步骤。(1)对加有莱斯噪声的原始磁共振图像进行方差稳定变换，将依赖信号分布的莱斯噪声变成独立于信号分布的噪声，得到变换后的噪声图像。(2)将变换后的噪声图像中的每一个像素作为标的像素或本文档来自技高网...

【技术保护点】
一种基于高阶奇异值分解的磁共振图像去噪方法，依次包括如下步骤：（1）对加有莱斯噪声的原始磁共振图像进行方差稳定变换，将依赖信号分布的莱斯噪声变成独立于信号分布的噪声，得到变换后的噪声图像；（2）将变换后的噪声图像中的每一个像素作为标的像素或者按照一定距离间隔取像素作为标的像素，取标的像素及其周围的像素作为参考块，执行以下操作进行第一次高阶奇异值分解去噪，具体是：（a）通过k最近邻方法依次逐个寻找参考块的相似块，组成高维数组；（b）对高维数组进行高阶奇异值分解变换，得到对应的系数和自适应基；（c）通过第一阈值操作，将绝对值小于第一阈值的系数置零得到修正后的系数；（d）根据自适应基与修正后的系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的估计值；（3）对步骤（2）得到的估计值加权平均进行像素合并，得到像素合并后的估计值，对每一个像素进行像素合并得到第一次去噪后的图像；（4）对方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均，得到第一次加权平均后的图像；（5）对第一次加权平均后的图像进行第二次高阶奇异值分解去噪，具体是：（5.1）将第一次加权平均后的图像中的每一个像素作为二次标的像素或者按照一定距离取像素作为二次标的像素，取二次标的像素及其周围的像素作为二次参考块，执行以下操作：（a）以第一次去噪后的图像进行相似块判断，通过k最近邻方法依次逐个寻找二次参考块的二次相似块，组成二次高维数组；（b）对二次高维数组进行高阶奇异值分解变换，得到对应的二次系数和二次自适应基；（c）通过二次阈值操作，将绝对值小于二次阈值的二次系数置零得到修正后的二次系数；（d）根据二次自适应基与修正后的二次系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的二次估计值；（5.2）对步骤（5.1）得到的二次估计值加权平均进行像素合并，得到像素合并后的二次估计值，对每一个像素进行像素合并得到第二次去噪后的图像；（6）对第二次去噪后的图像进行方差稳定逆变换，得到最终的滤波图像。...

【技术特征摘要】
1.一种基于高阶奇异值分解的磁共振图像去噪方法，依次包括如下步骤：(1)对加有莱斯噪声的原始磁共振图像进行方差稳定变换，将依赖信号分布的莱斯噪声变成独立于信号分布的噪声，得到变换后的噪声图像；(2)将变换后的噪声图像中的每一个像素作为标的像素或者按照一定距离间隔取像素作为标的像素，取标的像素及其周围的像素作为参考块，执行以下操作进行第一次高阶奇异值分解去噪，具体是：(a)通过k最近邻方法依次逐个寻找参考块的相似块，组成高维数组；(b)对高维数组进行高阶奇异值分解变换，得到对应的系数和自适应基；(c)通过第一阈值操作，将绝对值小于第一阈值的系数置零得到修正后的系数；(d)根据自适应基与修正后的系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的估计值；(3)对步骤(2)得到的估计值加权平均进行像素合并，得到像素合并后的估计值，对每一个像素进行像素合并得到第一次去噪后的图像；(4)对方差稳定变换后的噪声图像和第一次去噪后的图像进行加权平均，得到第一次加权平均后的图像；(5)对第一次加权平均后的图像进行第二次高阶奇异值分解去噪，具体是：(5.1)将第一次加权平均后的图像中的每一个像素作为二次标的像素或者按照一定距离取像素作为二次标的像素，取二次标的像素及其周围的像素作为二次参考块，执行以下操作：(a)以第一次去噪后的图像进行相似块判断，通过k最近邻方法依次逐个寻找二次参考块的二次相似块，组成二次高维数组；(b)对二次高维数组进行高阶奇异值分解变换，得到对应的二次系数和二次自适应基；(c)通过二次阈值操作，将绝对值小于二次阈值的二次系数置零得到修正后的二次系数；(d)根据二次自适应基与修正后的二次系数，进行高阶奇异值分解逆变换，将逆变换后的结果作为该高维数组中所有图像块去除噪声后的二次估计值；(5.2)对步骤(5.1)得到的二次估计值加权平均进行像素合并，得到像素合并后的二次估计值，对每一个像素进行像素合并得到第二次去噪后的图像；(6)对第二次去噪后的图像进行方差稳定逆变换，得到最终的滤波图像。2.根据权利要求1所述的基于高阶奇异值分解的磁共振图像去噪方法，其特征在于：...

【专利技术属性】
技术研发人员：冯衍秋，张鑫媛，徐中标，陈武凡，
申请(专利权)人：南方医科大学，
类型：发明
国别省市：广东;44