基于奇异值权重函数的非局部TV模型图像去噪方法技术

一种基于奇异值权重函数的非局部TV模型图像去噪方法。步骤如下：(1)首先输入噪声图像。(2)设置算法的相关参数，包括非局部搜索窗口大小N

Image denoising method of non local TV model based on singular value weighting function

An image denoising method based on singular value weighting function for non local TV model. The steps are as follows: (1) input noise image first. (2) setting the parameters of the algorithm, including the size of the nonlocal search window N

【技术实现步骤摘要】
基于奇异值权重函数的非局部TV模型图像去噪方法一、
本专利技术属于图像处理
，具体涉及去除加性噪声的图像去噪领域和改进的非局部TV模型的图像去噪方法。二、
技术介绍
图像去噪旨在通过对被噪声污染的图像进行某种处理，以降低噪声对原始有用信息的影响，尽可能地还原出被噪声污染前的图像。GUYGILBOA和STANLEYOSHER提出的非局部TV模型(NonlocalTotalVariation,NLTV)是将GUYGILBOA和STANLEYOSHER提出的非局部算子(见文献：NONLOCALOPERATORSWITHAPPLICATIONSTOIMAGEPROCESSING.SIAMMultiscaleModelingandSimulation.Vol.7,No.3,pp.1005–1028)引入到Rudin-Osher-Fatemi提出的总变分(TotalVariation，TV)模型中而得到的，NLTV模型具有在去除图像噪声的同时能较好地保留图像的纹理细节的特点，该模型包含保真项和正则项。保真项表示观测图像(即待去噪的图像)与去噪后图像的接近程度；正则项表达去噪后图像光滑程度的先验知识。在正则项中，像素相似度权重函数具有很重要的作用，其表达式与含噪图像的图像块中各像素灰度值直接相关，但是由于噪声的存在，噪声会对像素灰度值产生干扰，进而影响像素相似度权重。奇异值分解是一种基于特征向量的矩阵变换方法，在信号处理、模式识别、数字水印技术等方面都得到了应用。由于图像具有矩阵结构，因此本专利技术将其用于提取图像块的主要特征，将提取的图像主要特征用于构建新的像素相似度权重函数，以降低噪声对像素相似度权重的影响，再将构建的新的像素相似度权重函数应用到NLTV模型，得到新的NLTV模型，用于图像去噪。本文中所提到的NLTV模型是指GUYGILBOA和STANLEYOSHER提出的非局部TV模型。三、
技术实现思路
本专利技术的目的是克服目前现有基于非局部理论的图像去噪技术中，在图像受到噪声干扰后，图像的像素值被噪声污染而改变，像素相似度权重赋值会受到噪声干扰的缺陷。本专利技术通过构建一种新型的像素相似度权重函数来减少噪声的干扰，提高像素相似度权重赋值的准确性。本专利技术所述的基于奇异值权重函数的非局部TV模型图像去噪方法包括如下步骤：(1)首先输入N0×N0大小的待去噪的图像f；(2)设置本专利技术方法的相关参数，包括非局部搜索窗口大小N1×N1、邻域窗口大小N2×N2、像素相似度权重函数的参数h和j，高斯核的标准差σ，分裂的Bregman迭代辅助变量bk的初始值b0、平滑参数θ，以及保真参数λ；(3)通过奇异值分解方法获得图像f中各领域图像块(图像块大小为N2×N2)的奇异值矩阵和最大奇异值；(4)构建基于步骤(3)得到的最大奇异值的图像f的像素相似度权重函数；(5)应用步骤(4)构建的权重函数，建立基于该权重函数的用于对图像f进行去噪的新的非局部TV模型；(6)对步骤(5)建立的新的非局部TV模型，采用分裂的Bregman算法，通过逐次迭代逼近来求解；并设迭代计数变量初值k＝0；(7)进行分裂的Bregman算法数值迭代运算，获得本次迭代后的输出图像uk+1；(8)每次完成步骤(7)的迭代后，计算出迭代结果uk+1的评价值，即uk+1的峰值信噪比(PSNR)；如果本次迭代后的输出图像uk+1的峰值信噪比大于上一次迭代输出图像uk的峰值信噪比，则不满足迭代停止条件，令k＝k+1，并返回步骤(7)，继续迭代运算；如果本次迭代结果uk+1的峰值信噪比小于等于上一次迭代输出图像uk的峰值信噪比，则满足迭代停止条件，将上一次迭代输出图像uk作为最优值ufinal输出，即令ufinal＝uk，并转到步骤(9)；(9)将结果ufinal作为最后去噪结果图像输出。本专利技术的优点是建立新的像素相似度权重函数。将图像块进行奇异值分解并选取最大奇异值，将选取的最大奇异值用于像素相似度权重函数的构建，减少了噪声对权重函数的干扰，能更准确的对像素相似度权重进行赋值。将由此构建的基于图像块的像素相似度权重函数应用到NLTV模型，用于图像去噪，提高图像去噪的效果。本专利技术方法适合用于所含噪声模型为加性噪声的含噪图像的去噪。四、附图说明图1为本专利技术方法的流程图。图2是原始图像未加噪的图像，其中，图2(11)是Peppers图，图2(12)是Cameraman图，图2(13)是Lena图。图3为本专利技术对Lena图去噪仿真结果，其中，图3(11)是待去噪的图像(是在原始图像上叠加了均值为零标准差为10的高斯噪声后的结果)，图3(12)非局部TV方法去噪，图3(13)本专利技术方法去噪。五、具体实施方式本专利技术所述的基于奇异值权重函数的非局部TV模型图像去噪方法，具体步骤如下：(1)首先输入N0×N0大小的待去噪图像f；(2)设置本专利技术方法的相关参数，包括非局部搜索窗口大小N1×N1、邻域窗口大小N2×N2、像素相似度权重函数的参数h和j，高斯核的标准差σ，分裂的Bregman迭代辅助变量bk的初始值b0、平滑参数θ和保真参数λ；(3)设Mx是以步骤(1)输入的含噪图像f中像素点x∈Ω为中心的大小为N2×N2的图像块像素灰度值矩阵,Ω为f的图像空间，不同的像素点对应不同的图像块。将各图像块Mx进行奇异值分解：Mx＝UxΛxVxT。式中Ux、Vx分别是Mx的左奇异矩阵和右奇异矩阵，大小都是N2×N2。Λx是Mx的奇异值矩阵，它的大小也是N2×N2，它的对角元素不为零，其他元素全为零，Λx对角元素一共有N2个：按从大到小排列为：Λx的对角元素就是Mx的奇异值，它包含了Mx的全部特征。奇异值的定义：对于m×n阶矩阵A，ATA的n个特征值的非负平方根叫作A的奇异值；(4)利用步骤(3)得到的图像块的奇异值构建新的像素相似度权重函数。由于奇异值包含了图像块的主要特征，所以相似的图像块之间的奇异值是相近的；在图像块的奇异值中，大的奇异值包含了图像块的主要特征，小的奇异值包含了图像块的次要特征；并且在含有噪声的图像中，噪声不是图像的主要特征，所以在构建新的像素相似度权重函数时，只选用图像块的最大奇异值，通过图像块的主要特征来判断图像块之间的相似性，这样就减少了噪声的干扰；构建图像f中两个像素点x和y的相似度权重函数：其中x为当前像素点，y是以x为中心的搜索窗口内的一点，ax是以x为中心、大小为N2×N2的邻域内的像素灰度值矩阵Mx的最大奇异值，ay是以y为中心、大小为N2×N2的邻域的像素灰度值矩阵My的最大奇异值；表示在图像f中以x为中心的大小为N2×N2的图像块与以y为中心的大小为N2×N2的图像块之间的高斯加权距离，表示求和范围是以x或y为中心的大小为N2×N2的邻域内的每一像素点(不包括x或y自身)，共N2×N2-1项，Gσ(c)是标准差为σ的高斯核函数；h和j是常数，通过干预指数函数的衰减速度来控制权重函数ω的大小,h和j的取值越大，权重函数ω的值越接近1，算法收敛速度快，但是难以达到最优值，h和j的取值越小，权重函数ω的值越接近0，算法经过多次迭代可以收敛到最优值，但是耗费时间多，h和j的取值原则上要综合以上两点因素，取大小合适的值；(5)建立非局部TV模型：其中J(u本文档来自技高网...

【技术保护点】
基于奇异值权重函数的非局部TV模型图像去噪方法，包括如下步骤：(1)首先输入N

【技术特征摘要】
1.基于奇异值权重函数的非局部TV模型图像去噪方法，包括如下步骤：(1)首先输入N0×N0大小的待去噪图像f；(2)设置本发明方法的相关参数，包括非局部搜索窗口大小N1×N1、邻域窗口大小N2×N2、像素相似度权重函数的参数h和j，高斯核的标准差σ，分裂的Bregman迭代辅助变量bk的初始值b0、平滑参数θ和保真参数λ；(3)设Mx是以步骤(1)输入的含噪图像f中像素点x∈Ω为中心的大小为N2×N2的图像块像素灰度值矩阵,Ω为f的图像空间，不同的像素点对应不同的图像块。将各图像块Mx进行奇异值分解：Mx＝UxΛxVxT。式中Ux、Vx分别是Mx的左奇异矩阵和右奇异矩阵，大小都是N2×N2。Λx是Mx的奇异值矩阵，它的大小也是N2×N2，它的对角元素不为零，其他元素全为零，Λx对角元素一共有N2个：按从大到小排列为：Λx的对角元素就是Mx的奇异值，它包含了Mx的全部特征。奇异值的定义：对于m×n阶矩阵A，ATA的n个特征值的非负平方根叫作A的奇异值；(4)利用步骤(3)得到的图像块的奇异值构建新的像素相似度权重函数。由于奇异值包含了图像块的主要特征，所以相似的图像块之间的奇异值是相近的；在图像块的奇异值中，大的奇异值包含了图像块的主要特征，小的奇异值包含了图像块的次要特征；并且在含有噪声的图像中，噪声不是图像的主要特征，所以在构建新的像素相似度权重函数时，只选用图像块的最大奇异值，通过图像块的主要特征来判断图像块之间的相似性，这样就减少了噪声的干扰；构建图像f中两个像素点x和y的相似度权重函数：其中x为当前像素点，y是以x为中心的搜索窗口内的一点，ax是以x为中心、大小为N2×N2的邻域内的像素灰度值矩阵Mx的最大奇异值，ay是以y为中心、大小为N2×N2的邻域的像素灰度值矩阵My的最大奇异值；表示在图像f中以x为中心的大小为N2×N2的图像块与以y为中心的大小为N2×N2的图像块之间的高斯加权距离，表示求和范围是以x或y为中心的大小为N2×N2的邻域内的每一像素点(不包括x或y自身)，共N2×N2-1项，Gσ(c)是标准差为σ的高斯核函数；h和j是常数，通过干预指数函数的衰减速度来控制权重函数ω的大小,h和j的取值越大，权重函数ω的值越接近1，算法收敛速度快，但是难以达到最优值，h和j的取值越小，权重函数ω的值越接近0，算法经过多次迭代可以收敛到最优值，但是耗费时间多，h和j的取值原则上要综合以上两点因素，取大小合适的值；(5)建立非局部TV模型：其中J(u)是模型的目标函数，λ是保真参数，f是输入的含噪图像，u是去噪后得到的图像，Ω是f、u的图像空间，是GUYGILBOA和STANLEYOSHER提出的非局部梯度算子，u(y)和u(x)分别是图像u上的点x和y的像素灰度值，(6)对步骤(5)建立的NLTV模型，采用分裂的Bregman算法迭代逼近来求解，迭代过程分成三步，降低了迭代的...

【专利技术属性】
技术研发人员：金燕，蒋文宇，万宇，赵羿，杜伟龙，王雪丽，
申请(专利权)人：浙江工业大学，
类型：发明
国别省市：浙江,33