图像滤波去噪处理

摘要图像是信息社会人们获取信息的重要来源之一。在通过图像传感器将现实世界中的有用图像信号进行采集、量化、编码、传输、恢复的过程中，存在大量影响图像质量的因素。因此图像在进行使用之前，一般都要经过严格的预处理如去噪、量化、压缩编码等。噪声的污染直接影响着对图像边缘检测、特征提取、图像分割、模式识别等处理，使人们不得不从各种角度进行探索以提高图像的质量。所以采用适当的方法尽量消除噪声是图像处理中一个非常重要的预处理步骤。图像处理技术在20世纪首先应用于图像的远距离传送，而改善图像质量的应用开始于1964年美国喷气动力实验室用计算机对“徘徊者七号”太空船发回的月球照片进行处理，并获得巨大成功。现在图像处理技术已深入到科学研究、军事技术、工农业生产、医学、气象及天文学等领域。科学家利用人造卫星可以获得地球资源照片、气象情况；医生可以通过X射线或CT对人体各部位的断层图像进行分析。但在许多情况下图像信息会受到各种各样噪声的影响，严重时会影响图像中的有用信息，所以对图像的噪声处理就显得十分重要。因此我选择图像去噪方面进行了解及研究，现将自己已了解的知识进行汇总。目录摘要2一、图像滤波的应用4二、均值滤波51、 均值滤波的思想1、 均值滤波的算法1、 均值滤波的实验结果三、中值滤波71中值滤波的思想1中值滤波的算法1中值滤波的实验结果四、维纳滤波88维纳滤波的思想8维纳滤波的算法8维纳滤波的实验结果五、小波变换9小波变换滤波的思想小波变换滤波的算法小波变换滤波的实验结果六、Contourlet变换的图像去噪11Contourlet变换的基本思想Contourlet变换的算法七、全变差正则化的Shearlet收缩去噪12Shearlet收缩去噪原理简介Shearlet收缩去噪算法八、结果分析及自己的收获12结果分析自己的收获参考文献13图像滤波的应用图像是信息社会人们获取信息的重要来源之一。在现实世界中，存在大量影响图像质量的因素。因此图像在进行使用之前，一般都要经过严格的预处理如去噪、量化、压缩编码等。噪声的污染直接影响着对图像边缘检测、特征提取、图像分割、模式识别等处理，因此我们不得不从各种角度进行探索以提高图像的质量。采用适当的方法尽量消除噪声是图像处理中一个非常重要的预处理步骤。现在图像处理技术已深入到科学研究、军事技术、工农业生产、医学、气象及天文学等领域。科学家利用人造卫星可以获得地球资源照片、气象情况；医生可以通过X射线或CT对人体各部位的断层图像进行分析。但在许多情况下图像信息会受到各种各样噪声的影响，严重时会影响图像中的有用信息，所以对图像的噪声处理就显得十分重要。为了改善图像质量，从图像中提取有效信息，必须对图像进行去噪预处理。根据噪声的频谱分布的规律和统计特征以及图像的特点，出现了多种多样的去噪方法。经典去噪方法有：空间域合成法、频域合成法和最优线性合成法等。与之相适应的出现了许多应用方法：如均值滤波器、中值滤波器、低通滤波器、维纳滤波器、最小失真法等。这些方法广泛应用，促进了数字信号处理的极大发展。显著提高了图像质量。近年来，小波变换去除噪声的方法得到广泛的应用。与传统的去噪方法相比。它利用的是非线性域值，在时间域和频率域同时具有良好的局部化性质，而且时窗和频窗的宽度可以调节。对高频成分采用逐渐精细的时域或空域取样步长，从而可以聚焦到对象的任意细节。因此可以提高散斑高散射特性的噪声对比度，很好的消除散斑噪声9。小波变换去除噪声的方法在不断地发展，去噪方法很多，如非线性小波变换阈值法去噪、小波变换模极大值去噪及基于小波变换域的尺度相关性去噪法均值滤波均值滤波的基本思想是将某像素邻域内的各点的灰度平均值来代替该像素原来的灰度级。通常邻域都取成N*N的方形窗口,其降噪平滑后的图像为取力s(ij)其中，s是点(x,y)邻域内的点集，M是点集S中的总点数。将算术均值滤波器作一改进，将某像素邻域内的各点的灰度加权平均值来代替该像素原来的灰度值，得到加权均值滤波。从权值上看，灰度越接近中心像素其权值越大。加权平均的算法可表示为gfey)=O-jlES其中，w(i,j)是权值，表示其所起作用的大小。均值滤波的算法用MATLA段计均值滤波算法为:(1)高斯噪声sample=imread(lenna.jpg);sgray=rgb2gray(sample);I=imnoise(sample,gaussian,0,0.005);Igray=rgb2gray(I);subplot(3,2,1);imshow(sgray);title(subplot(3,2,2);imshow(Igray);title(K1=filter2(fspecial(averageK2=filter2(fspecial(averageK3=filter2(fspecial(averageK4=filter2(fspecial(averagesubplot(3,2,3);imshow(K1);title(subplot(3,2,4);imshow(K2);title(subplot(3,2,5);imshow(K3);title(subplot(3,2,6);imshow(K4);title(2)椒盐噪声sample=imread(lenna.jpg);sgray=rgb2gray(sample);原始图像);加入高斯噪声后的图像);,3),Igray)/255;,5),Igray)/255;,7),Igray)/255;,9),Igray)/255;模板尺寸为3*3的滤波后图像);模板尺寸为5*5的滤波后图像);模板尺寸为7*7的滤波后图像);模板尺寸为9*9的滤波后图像);I=imnoise(sample,salt&pepper,0.02);Igray=rgb2gray(I);subplot(3,2,1);imshow(sgray);title(subplot(3,2,2);imshow(Igray);title(原始图像);加入椒盐噪声后的图像);K1=filter2(fspecial(average,3),Igray)/255;K2=filter2(fspecial(average,5),Igray)/255;K3=filter2(fspecial(average,7),Igray)/255;K4=filter2(fspecial(average,9),Igray)/255;subplot(3,2,3);imshow(K1);title(subplot(3,2,4);imshow(K2);title(subplot(3,2,5);imshow(K3);title(subplot(3,2,6);imshow(K4);title(2.3均值滤波的实验结果模板尺寸为3*3的滤波后图像);模板尺寸为5*5的滤波后图像);模板尺寸为7*7的滤波后图像);模板尺寸为9*9的滤波后图像);(1)高斯噪声的滤波结果j;u 收 Jul*1 茹本诺更互叫 wdULnimjwTextGo=iRndlselEi(raj=rib2(Jubplot(3j2aubplot(3f2yl=filter2(f2=filter2(f3=filter2(f4=filter2(fubplot(Sj2,ubplot(3,工ubplot。2Tubplot%2j:erencepagedocfilter2结果分析：通过通过图像的识别率可以看出：模板越大去噪效果越好，但图像模糊度也随之增加。(2)椒盐噪声去噪结果,Editor-G:2C3iaoyanzjaosheng.m由图可得：模板越大，去噪效果越好，因此使用均值滤波去噪时选用的模板尺寸（邻域半径）越大效果越好。三、中值滤波中值滤波的基本思想此方法是用该像素的相邻像素的灰度中值来代替像素值。是一种典型的排序滤波器。1=median其中，是邻域中各点的灰度值。所以，中值滤波后的图像g(x,y)是以(x,y)为中心的N*N窗口的各像素的灰度中间值。中值滤波的MATLAB法用MATLA破置中值滤波算法为：sample=imread(lenna.jpg);sgray=rgb2gray(sample);I1=imnoise(sample,salt&pepper,0.02);Igray1=rgb2gray(I1);I2=imnoise(sample,gaussian,0,0.005);Igray2=rgb2gray(I2);K1=medfilt2(Igray1);K2=medfilt2(Igray2);subplot(2,3,1);imshow(sgray);title(原始图形);subplot(2,3,2);imshow(Igray1);title(加椒盐噪声后的图像);fileEditText:唱唱-juzhiLm*:OVt?1.0sanj)le=LK2=filteK4=filtEquzIgray=r1subplot(sgray=iiI=i3mDL5subplot(2,3,3);imshow(Igray2);title(加高斯噪声后的图像,);subplot(2,3,4);imshow(K1);title(中值滤波1(椒盐);subplot(2,3,5);imshow(K2);title(中值滤波2(高斯),);3.3中值滤波的实验结果通过结果可以看出，中值滤波对于椒盐噪声的去噪能力比对高斯噪声的去噪能力好。四、维纳滤波维纳滤波的基本思想假定线性滤波器的输入为有用信号和噪声之和，两者均为广义平稳过程且知它们的二阶统计特性，维纳根据最小均方误差准则(滤波器的输出信号与需要信号之差的均方值最小)，求得了最佳线性滤波器的参数，这种滤波器被称为维纳滤波器。维纳滤波器的优点是适应面较广，无论平稳随机过程是连续的还是离散的，是标量的还是向量的，都可应用。对某些问题，还可求出滤波器传递函数的显式解，并进而采用由简单的物理元件组成的网络构成维纳滤波器。维纳滤波器的缺点是，要求得到半无限时间区间内的全部观察数据的条件很难满足，同时它也不能用于噪声为非平稳的随机过程的情况，对于向量情况应用也不方便。因此，维纳滤波在实际问题中应用不多。实现维纳滤波的要求是：(1)输入过程是广义平稳的；(2)输入过程的统计特性是已知的。维纳滤波的算法维纳滤波的MATLABT法为：sample=imread(lenna.jpg);sgray=rgb2gray(sample);I1=imnoise(sample,salt&pepper,0.02);Igray1=rgb2gray(I1);I2=imnoise(sample,gaussian,0,0.005);Igray2=rgb2gray(I2);小波变换小波变换滤波的基本思想近年来，小波理论得到了非常迅速的发展，由于其良好的时频特性，实际应用非常广泛将小波用于去噪也得到很多学者的重视。小波去噪利用了不同中心频率的带通滤波器对信号滤波，把主要反映噪声频率的那些尺度的系数去掉，再把剩余各尺度的系数综合起来做反变换，从而使噪声得到较好的抑制。小波变换具有时频局域化特性，能够检测到局部突变的边缘特性，而且可将图像结构和纹理分别表现在不同分辨率层次上。其中，小波域值算法是利用信号和噪声小波系数幅值上的差异，通过选择合适的域值对小波系数进行处理，以达到去除噪声又保留有用信号的目的。此方法中，小波系数域的处理方法及阈值的估计是两个关键技术。在此方法基础上又开发出了，基于Bayes方法的小波去噪和基于假设检验的小波去噪。小波变换滤波的算法小波变换的MATLA眸法为:X=imread(lenna.jpg);I=rgb2gray(X);K1=wiener2(Igray1,33);K2=wiener2(Igray2,33);subplot(2,3,1);imshow(sgray);title(subplot(2,3,2);imshow(Igray1);title(subplot(2,3,3);imshow(Igray2);title(subplot(2,3,4);imshow(K1);title(subplot(2,3,5);imshow(K2);title(4.3维纳滤波的实验结果原始图像）;加椒盐噪声后的图像）;加高斯噪声后的图像）;维纳滤波1（椒盐）；维纳滤波2（高斯）；通过实验结果可以看出：维纳滤波具有较大的局限性，但其对高斯噪声的去噪效果相对较好。J=imnoise(I,gaussian,0,0.005);c,l=wavedec2(J,2,sym4J1=wrcoef2(a,c,l,sym4J2=wrcoef2(a,c,l,sym4thr,sorh,keepapp=ddencmp(J3=wdencmp(gbl,J,sym4J4=wdencmp(gbl,J,sym4J5=medfilt2(J3);J6=medfilt2(J4);figure;subplot(2,3,1);imshow(I);title(subplot(2,3,2);imshow(J);title(subplot(2,3,3);image(J1);title(subplot(2,3,4);image(J2);title(subplot(2,3,5);image(J3);title(subplot(2,3,6);image(J4);title(figure;);,1);,2);den,wv,J);,1,thr,sorh,keepapp);,2,thr,sorh,keepapp);原始图像）；加入高斯噪声后的图像）第一次wrcoef2去噪）；第二次wrcoef2去噪）；第一次阈值去噪图像）；第二次阈值去噪图像）；原始图像）；加入高斯噪声后的图像）第二次滤波1）;第二次滤波2）;+国 r*9国 k 目3-9-11-12-13-库-17-1S-19-20-21-subplot(2,3,1);imshow(I);title(subplot(2,3,2);imshow(J);title(subplot(2,3,4);image(J5);title(subplot(2,3,5);image(J6);title(5.3通过小波变换滤波的实验结果通过小波分解的图像以及全局阈值滤波后的图像如下所示:FilflEditTextFilflEditText1.01.0subplot(2subplot1 12 2subplot(2(2subplot(2figure;subplot(2subplot2subplot22subplot(2(2J5=n?dfilJ6=Piedfilfigure：subplot(2Contourlet变换的图像去噪Contourlet变换的基本思想Contourlet变换可以满足曲线的各向异性尺度关系，能够很好地抓住图像的几何结构。Contourlet变换是一种新的图像二维表示方法，具有多分辨率、局部定位、多方向性、近邻界采样和各向异性等性质，其基函数分布于多尺度、多方向上，少量系数即可有效地捕捉图像中的边缘轮廓，而边缘轮廓正是自然图像中的主要特征。Contourlet变换的基本思想是首先用一个类似小波的多尺度分解捕捉边缘奇异点，再根据方向信息将位置相近的奇异点汇集成轮廓段。选用Burr和Adelson于1983年提出的拉普拉斯塔式滤波器结构(LP)对图像多分辨率分解来捕捉奇异点。Contourlet变换的算法由于Contourlet变换缺乏平移不变性因此在应用它进行去噪时会带来人为的视觉效果。为克服这些人为视觉效果这里引入了基于循环平移的平移不变Contourlet去噪算法。假设含噪图像f(i,j),0=i,j=N。定义F为平移n位的平移算子Fnf(i,j)=f(in,jn)mod(N)1由于F是一一X应的所以其反变换因子可以设为F_n=Fn，一幅图像中包含许多奇异点对某个奇异点来说是效果最佳的平移量”而对另一个奇异点来说效果可能很差。为了解决这个问题在一定范围内进行循环平移运算再将平移运算的所有结果求平均值。将经过小波变换并经过全局阈值处理的图像再用中值滤波进行第二次滤波，以加强滤波效果。结果如下所示：以Figure2 2I1=1FileEditViewInsertToolsDesktopWindowHelp丁/I相)“十|日国|口Editor-Editor-6 6 旌口转三内废像房03|笛勺里J图唱-1QI+|j-J3=wdencmpCgbl-J4=wdencmp(jbl-J5=medfilt2(J33-J6=medfilt2(J4J-figure;-subplot(2j3j1)-subplot(2j3)；15-subplot(,2,(,2,3j3)；ia-subplot03,4);17-subplot(2j3j5):13-subplot3j6);1 电figur#;-subplot(231)-subplot32)：e esubplot3,4)：-subplot3j5)：加入高斯噪声后的图像quzaochuli,m*zhong.原始图像具体算法描述如下：1)对含噪图像f(i,j)进行循环平移；2)对每次平移后的图像做离散Contourlet变换，可以得到各尺度各方向上的Contourlet变换系数c。3)Contourlet变换系数c进行阈值处理，得到去噪后的系数cAo4)利用cA进行离散Contourlet反变换，得到去噪后图像fA。全变差正则化的Shearlet收缩去噪Shearlet收缩去噪原理简介收缩法是变换域中最为广泛的图像去噪方法。阈值收缩主要基于以下事实，比较大的系数一般都是以实际信号为主，而比较小的系数则很大程度上是噪声。即在系数中，低频分量中含有大量的信息，应该给予保留；同时在高频分量中，一些绝对值大的重要的系数并不是噪声，而是边缘信息，也应保留。硬阈值收缩法的收缩函数为：式中r为收缩的阈值。硬阈值是将绝对值小于阈值的系数置为零，而将绝对值大于阈值的系数不加任何处理给予保留。硬阈值的优点是运算速度快，能取得较好的去噪效果。但由于其收缩函数是不连续的，所以在图像中会产生许多“人为”噪声，图像会出现振铃、伪吉布斯效应等视觉失真。图3给出了硬阈值函数示意图。图中的小正方形是代表绝对值小于阈值的系数，因为在硬阈值处理中将小于阈值的系数直接置零，这样就丢失了图像的一部分信息。伪吉布斯现象可以视为一种整体振荡，而全变差正则化对整体振荡有很好的抑制作用。为了避免硬阈值收缩去噪丢失的图像信息，以及更好地消除伪吉布斯现象，本文引入全变差正则化。Shearlet收缩去噪算法1)Shearlet硬阈值估计AA对含噪图像进行Shearlet变换，得到shearlet系数djl(f)和尺度系数Cjl(f);J,lJ,l(2)Shearlet系数进行硬阈值去噪，用愀表示保留系数的指标集。用Shearlet反变换得到图像的初始估计f。2)迭代修正初始化lx=1,N代表迭代次数，%=Lmax,dj=LmaX/N。计算Jp(f)的次梯度2Jp(f)；(2)将17fJp(f)投影到V空间彳导到R(fJp(f)；利用式(13)计算f(nW=fnKnFV(VfJp(f)；(4)n=n+l,九0书=4-d九，当儿n之。时转至步骤2)中第(1)步循环执行，否则结束迭代退出循环。结果分析及自己的收获结果分析通过所有滤波设计及结果对比，我们可以得到:第一次wrco日12去瞟50005000在处理服从高斯分布的一类噪声时，维纳滤波与中值滤波去除效果较好一些；中值滤波对于去除椒盐噪声效果好，而维纳滤波去除效果差，中值滤波对于去除椒盐噪声效果明显，是因为椒盐噪声只在画面上的部分点随机出现，而中值滤波根据数据排序，将未被污染的点代替噪声点的值的概率较大，所以抑制效果好。对点、线和尖顶较多的图像不宜采用中值滤波，因为一些细节点可能被当成噪声点。而且均值滤波，中值滤波都在一定程度上模糊了图像边缘，提高滤波效果是以模糊边缘为代价的。小波分解可以把图像分层次按照小波基展开，并且可以根据图像的性质及给定的处理标准确定展开到哪一级为止，还可以把细节分量和近似分量分开。由图可得：wrcoef2和wpdencm薛函数可以有效地进行去噪处理。而且在小波变换的基础上再利用中值滤波可以进一步增加去噪效果。因此对边缘具有很好的保留，不会模糊边缘。因此滤波效果最好。在分析过程中可以看出，针对不同类型的噪声需要相应的滤波去噪算法才能取得较好的均值滤,波图像100200效果，才能使后续的图像处理工作得以更加优质的进行。自己的收获通过这次对图像滤波去噪处理的了解，以及对所有方法基本思想的了解，进一步激发了自己对该课题的兴趣，也拓宽了自己的知识面，也巩固了自己所学的数字信号处理方面相关的知识。在整个课题中，通过对不同的去噪算法的分析，以及结果的比较，更好地了解了在不同情况下应利用不同的方法来解决问题，了解每一个滤波器的优缺点，在遇到不同问题时，可以根据各个滤波方法的优缺点选择合适的滤波方法，从而使我们的效果达到最好，在了解滤波器的原理的过程中，我了解到许多数字信号处理的知识。对于数字图像信号，噪声表现为或大或小的极值，这些极值通过加减作用于图像象素的真实灰度值上，在图像造成亮、暗点干扰，极大降低了图像质量，影响图像复原、分割、特征提取、图识别等后继工作的进行。要构造一种有效抑制噪声的滤波机必须考虑两个基本问题能有效地去除目标和背景中的噪声;同时，能很好地护图像目标的形状、大小及特定的几何和拓扑结构特征。图像频率域去噪方法是对图像进行某种变换，将图像从空间域转换到频率域，对频率域中的变换系数进行处理，再进行反变换将图像从频率域转换到空间域来达到去除图像噪声的目的，例如小波变换。通过这次综述，对图像的处理有了进一步的理解，并学会了在网上查找自己所需的资料，可以很好的有利于自己理解自己所做的课题以及各种算法的利弊。参考文献：1胡海智，孙辉，邓承志，陈习，柳枝华.全变差正则化的she盯let收缩去噪J.中国图象图形学报，2011,16(2):l醴-1732郭蕾，田松，许悦雷，等.一种小波自适应比例萎缩去噪改进算法J.电视技术，2012,36(11)3YiS;LabateD;Eas1eyGRAshearletapproachtoedgeanalysisanddetection2009(05)4胡海智.孙辉.邓承志.陈习.柳枝华.占惠星.HUHai-zhi.SUNHui.DENGCheng-zhi.CHENXi.LIUZhi-hua.ZHAN,Hui-xing基于Shearlet变换的图像去噪算法期刊论文-计算机应用2010,30(6)5夏良正.数字图像处理M南京：东南大学出版社，19996JeongB.LeeYH.DesignofweightedorderstatisticfllterusingtheperceptionaIgorilhmJ.IEEETransactionsonSignalProcessing,19947谢杰成，张大力，徐文立.小波去噪综述J.中国图象图形学报，2002,7(3)8费佩燕，郭宝龙.基于多小波的去噪技术研究J.中国图象图形学报，2005,10(1)9LeventSendur,I.vanW.Selesnick.Bivariateshrinkagefunctionsforwavelet2baseddenoisingexploitinginterscaledependencyJ.IEEETrans.SignalProcessing,2002,50(11):2744275610梁栋；沈敏；高清维一种基于Contourlet递归CycleSpinning的图像去噪方法期刊论文-电子学报2005(11)