维纳滤波及图像去噪

1、-维纳滤波与图像去噪摘要首先选取对图像降噪比较有代表性的维纳滤波，在加有高斯噪声、椒盐噪声和乘性噪声的图像上进展处理，再将维纳滤波与中值滤波和均值滤波抑制噪声的效果进展比较，通过实验仿真及其处理效果，详细分析维纳滤波在图像去噪中的特点及各自作用的利弊。关键词维纳滤波；中值滤波；均值滤波；图像去噪Wiener filtering and image denoisingLIMeng,ZHAOQi*ian University of Posts and Telemunications, School of munication and information engineering，*ian7100

2、00, ChinaAbstract：Select the first is a representative of wiener filtering for image noise reduction with gauss noise and salt and pepper noise and multiplicative noise of image processing, then wiener filtering and median filtering and mean filtering noise effect parison, through the e*perimental s

3、imulation and the treatment effect, detailed analysis of wiener filtering in image denoising, the characteristics and the pros and cons of each role.Keywords：Wiener filtering,Median filtering,Mean filtering,Image denoising. z-0 引言图像在成像、传输、转换或存储的过程中会受到各种随机干扰信号即噪声的影响，从而会使画面变得粗糙、质量下降、特征淹没。为了减弱噪声、复原真实的画

4、面，就需要用到降噪滤波器对图像数据进展处理1。通过选取对图像降噪比较有代表性的维纳滤波对加有高斯噪声、椒盐噪声和乘性噪声的图像进展了处理，并将维纳滤波与中值滤波和均值滤波抑制噪声的效果进展比较，结合实验仿真及其处理效果，详细分析维纳滤波在图像去噪中的特点及各自作用的利弊。1图像去噪图像去噪是数字图像处理中的重要环节和步骤。去噪效果的好坏直接影响到后续的图像处理工作如图像分割、边缘检测等。图像信号在产生、传输过程中都可能会受到噪声的污染，一般数字图像系统中的常见噪声主要有：高斯噪声主要由阻性元器件部产生、椒盐噪声主要是图像切割引起的黑图像上的白点噪声或光电转换过程中产生的泊松噪声等。目前比较经典

5、的图像去噪算法主要有以下三种2：均值滤波：也称线性滤波，主要思想为邻域平均法，即用几个像素灰度的平均值来代替每个像素的灰度。有效抑制加性噪声，但容易引起图像模糊，可以对其进展改进，主要避开对景物边缘的平滑处理3。中值滤波：基于排序统计理论的一种能有效抑制噪声的非线性平滑滤波信号处理技术。中值滤波的特点即是首先确定一个以*个像素为中心点的邻域，一般为方形邻域，也可以为圆形、十字形等等，然后将邻域中各像素的灰度值排序，取其中间值作为中心像素灰度的新值，这里领域被称为窗口，当窗口移动时，利用中值滤波可以对图像进展平滑处理。其算法简单，时间复杂度低，但其对点、线和尖顶多的图像不宜采用中值滤波，很容易自

6、适应化4。Wiener维纳滤波：最典型的一种滤波，20世纪40年代，维纳奠定了最正确滤波器研究的根底，即假定输入是有用信号和噪声信号的合成，并且它们都是广义平稳过程，而且他们的二阶统计特性都。维纳根据最小均方准则即滤波器的输出信号与需要信号的均方值最小，求得了最正确线性滤波器的参数，这种滤波器被称为维纳滤波器。这是一种自适应滤波器，根据局部方差来调整滤波器效果，对于去除高斯噪声效果明显。由于基于维纳滤波器的图像复原效果比较好，具有一定的抑制噪声能力，近年来被广泛的应用到图象复原领域，维纳滤波算法得到不断的改进和开展。现在，许多有效的图像复原算法都是在此根底形成的5。2维纳滤波原理维纳滤波wie

7、ner filtering) 一种基于最小均方误差准则、对平稳过程的最优估计器。这种滤波器的输出与期望输出之间的均方误差为最小，因此，它是一个最正确滤波系统。它可用于提取被平稳噪声所污染的信号6。维纳滤波综合了退化函数和噪声统计特性两个方面进展复原处理。维纳滤波建立在最小化统计准则的根底上，它所得的结果只是平均意义上的最优。从退化图像复原出原图像的估计值，噪声为。用向量表示，为对的线性算子最小二乘方问题可看成是使形式为的函数服从约束条件的最小化问题。也就是说，在约束条件下求的最小化而得到f的最正确估计。这种有条件的极值问题可以用拉格朗日乘数法来处理。用拉格朗日法建立目标函数：其中为一常数，是拉

8、格朗日乘数。加上约束条件后，就可以按一般求极小值的方法进展求解。将上式两边对微分，并另其结果为零,得:求解，有： (1)式中可以调节，以满足约束条件。1式即为维纳滤波复原方法的根底。设和分别为 f 和 n 的相关矩阵，即：的第 ij 个元素是，代表 f 的第 i 个和第 j个元素的相关。因为 f 和 n 中的元素都是实数，所以和都是实对称矩阵，对于大多数图像而言，相邻像素之间相关性很强，在 20 30 个像素之外趋于零。在此条件下, 和可以近似为分块循环矩阵，并进展对角化处理，有:式中 A和 B 为对角阵，W 为酉阵，A 和 B 中的元素对应和中的相关元素的傅里叶变换。这些相关元素的傅叶变换称

9、为图像和噪声的功率谱。假设用来代替，则(1)式变为 (2)由循环矩阵对角化的知识可知，分块循环矩阵：其中D为对角矩阵，其元素正是H的本征值，D* 是D的复共轭。因而，2式变为：上式两边同乘以,得：写成频率域形式为:(3)上式称为维纳滤波，括号中的项组成的滤波器通常称为维纳滤波器或最小均方误差滤波器。其中，是退化图像的傅里叶变换，是退化函数，是噪声功率谱, 是原始图像功率谱。由上面原理推导可知，维纳滤波器的传递函数为：如果噪声是零，则噪声的功率谱消失，并且维纳滤波退化为逆滤波，所以说逆滤波是维纳滤波的特例。3维纳滤波原理个人推导其中，是原来的信号高清图像，是模糊退化矩阵，是参加的噪声，是模糊下采

10、样后的图像，是滤波系统，是估计信号恢复的图像。根据正交原则误差信号与进入估计的信号正交，可得：其中，是误差信号，是进入估计的信号。其中，g的自相关函数：所以，对上式，两边取傅里叶变换，得：因为信号f、噪声n互不相关，所以进一步得到：所以可以得到：因为在频域可写为令所以所以为了到达自适应的目的，参加一个调整系数s，得到：思路：将退化加噪后的信号作用于一个“白化滤波器，得到一个输出，再将此输出用滤波器进展滤波，得到最正确估计,即。实际系统都是因果的，可以根据信号模糊图像的功率谱采用谱分解定理求出信号模型的传输函数H，根据谱分解的结果可以确定。这样，就可以根据得到信号的最正确估计。4实验仿真与结果分

11、析程序见附录在仿真实验中，主要利用MATLAB实验平台，在MATLAB中可以按照维纳滤波的原理和公式来编写语句进展滤波，但由于此种方法较为复杂，同时MATLAB也有自带的维纳滤波器的函数，因此本课题中使用MATLAB自带的函数进展维纳滤波。在MATLAB中与维纳滤波有关的函数有wiener2()和deconvwnr()，这两个函数都能够完成维纳滤波的功能，deconvwnr()强调图象复原方面，wiener2()强调图像空间域锐化的作用7。其中wiener2()函数只支持二维滤波，由于此处选的是一Lena的灰度图片，使用wiener2()函数。以下四个实验中，均采用默认的3 3窗口维纳滤波器。

12、一实验一、维纳滤波对不同噪声的滤波效果结论：从图中可以看到维纳滤波对高斯噪声、乘性噪声都有明显的抑制作用。二实验二：给图像参加高斯噪声，作均值、中值和维纳滤波结论：对于高斯噪声，维纳滤波对噪声的抑制效果更好,容易失去边缘信息。三实验三：给图像参加椒盐噪声，作均值、中值和维纳滤波结论：中值滤波对于去除椒盐噪声效果最好，而维纳滤波去除效果则最差。四实验四：给图像加乘性噪声，作均值、中值和维纳滤波结论：对于乘性噪声，维纳滤波对噪声的抑制效果更好，但也容易失去边缘信息。5结论维纳滤波对高斯噪声、乘性噪声都有明显的抑制作用，正常情况下，图像恢复效果比均值滤波和中值滤波效果好，但维纳滤波容易失去边缘信息。

13、而对于椒盐噪声几乎没有抑制效果。维纳滤波适应面广，去噪效果好，但算法复杂度高，处理效率不高。维纳滤波复原法存在着几个实质性的局限。第一，最有标准是基于最小均方误差的且对所有误差等权处理，这个标准在数学上可以承受，但却是个不适合人眼的方式，原因在于人类对复原错误的感知在具有一致灰度和亮度的区域中更为严重，而对于出现在暗的和高梯度区域的误差敏感性差得多。第二，空间可变的退化不能用维纳滤波复原法复原，而这样的退化是常见的。第三，维纳滤波不能处理非平稳信号和噪声。通过MATLAB仿真实验，使我们更加深刻地了解到维纳滤波的原理及其在图像处理方面的应用。在图像恢复处理中使用的方法还有很多，应该根据具体情况

14、做具体分析，维纳滤波是假设图像信号可以近似看成平稳随机过程的前提下，使输入图像和恢复图像之间的均方误差到达最小的准则函数来实现图像恢复的方法。如果图像的统计特性，则用维纳滤波来进展图像处理的效果还是可观的8。我们在上述实验过程中假定的噪声是高斯的且是加性的，噪声和信号相互独立，然而在实际测量中很多图像的噪声往往是非加性的，因此我们需要进一步将维纳滤波器推广到更复杂的情况中，由于时间和水平有限，在此就不加表达。参考文献1朝晖、弘编.数字图像处理及应用M.:机械工业，200406.2卓亚.图像去噪技术综述J.科技信息,2021,(15):317.3吴建华,迟生,周卫星.中值滤波与均值滤波的去噪性能

15、比较J.大学学报(工科版),1998,(01):33-36.4鸿林,忠民,羿宗琪.中值滤波技术在图像处理中的应用J.信息技术,2004,(07):26-27.5晓杰.维纳滤波的应用研究D.长安大学,2006.6丁玉美、阔永红、高新波编.数字信号处理M.:机械工业，200406.7曾敬枫.基于MATLAB维纳滤波算法在图像复原实验的应用J.现代计算机(专业版),2021,(34):3-5.8友凎.基于MATLAB的维纳滤波器仿真研究.中国科技论文网9田小平,娜,吴成茂.多尺度动静态联合彩色图像滤波处理.西邮学报10田小平;乔森;吴成茂.一种改进的图像模糊比照度增强算法.邮电大学学报. z-附录：

16、给图像加噪声：I=imread('C:UsersAdministratorDesktoplena.jpg');J=imnoise(I,'gaussian',0,0.005);subplot(2,3,1);imshow(I);title('原始图像');subplot(2,3,2); imshow(J);title('参加高斯噪声之后的图像');K1=filter2(fspecial('average',3),J)/255; K2=filter2(fspecial('average',5),J)/25

17、5;K3=filter2(fspecial('average',7),J)/255;K4= filter2(fspecial('average',9),J)/255; subplot(2,3,3);imshow(K1);title('改进后的图像1');subplot(2,3,4); imshow(K2);title('改进后的图像2');subplot(2,3,5);imshow(K3);title('改进后的图像3');subplot(2,3,6);imshow(K4);title('改进后的图像

18、9;);实验一：维纳滤波对不同噪声的滤波效果I=imread('C:UsersAdministratorDesktoplena.jpg');subplot(3,3,1);imshow(I);title('原始图像');J1=imnoise(I,'gaussian',0,0.02);subplot(3,3,2);imshow(J1);title('加高斯噪声图');K1=wiener2(J1,3 3); subplot(3,3,3);imshow(K1);title('维纳滤波后的图像');I=imread('

19、;C:UsersAdministratorDesktoplena.jpg');subplot(3,3,4);imshow(I);title('原始图像');J2=imnoise(I,'salt & pepper',0.02);subplot(3,3,5);imshow(J2);title('加椒盐噪声图');K1=wiener2(J2,3 3);subplot(3,3,6);imshow(K1);title('维纳滤波后的图像');I=imread('C:UsersAdministratorDesktopl

20、ena.jpg');subplot(3,3,7);imshow(I);title('原始图像');J3=imnoise(I,'speckle',0.02);subplot(3,3,8);imshow(J3);title('加乘性噪声图');K1=wiener2(J3,3 3); subplot(3,3,9);imshow(K1);title('维纳滤波后的图');实验二：参加高斯噪声，作均值、中值和维纳滤波I=imread('C:UsersAdministratorDesktoplena.jpg');J1=

21、imnoise(I,'gaussian',0,0.02);K2=filter2(fspecial('average',3),J)/255; K= medfilt2(J);¨K1=wiener2(J,3 3);¨subplot(2,3,1);imshow(I);title('原始图像');subplot(2,3,2);imshow(J);title('加噪图像');subplot(2,3,3);imshow(K2);title('均值滤波后的图像');subplot(2,3,4);imshow(K

22、);title('中值滤波后的图像');subplot(2,3,5);imshow(K1);title('维纳滤波后的图像');实验三：参加椒盐噪声，作均值、中值和维纳滤波I=imread('C:UsersAdministratorDesktoplena.jpg');J2=imnoise(I,'salt & pepper',0.02);K2=filter2(fspecial('average',3),J)/255; K= medfilt2(J);K1=wiener2(J,3 3); subplot(2,3,1);imshow(I);title('原始图像');subpl