圖像濾波去噪處理

1、摘要圖像是信息社會人們獲取信息的重要來源之一。在通過圖像傳感器將現(xiàn)實世界中的有用圖像信號進行采集、量化、編碼、傳輸、恢復(fù)的過程中，存在大量影響圖像質(zhì)量的因素。因此圖像在進行使用之前，一般都要經(jīng)過嚴格的預(yù)處理如去噪、量化、壓縮編碼等。噪聲的污染直接影響著對圖像邊緣檢測、特征提取、圖像分割、模式識別等處理，使人們不得不從各種角度進行探索以提高圖像的質(zhì)量。所以采用適當(dāng)?shù)姆椒ūM量消除噪聲是圖像處理中一個非常重要的預(yù)處理步驟。圖像處理技術(shù)在20世紀首先應(yīng)用于圖像的遠距離傳送，而改善圖像質(zhì)量的應(yīng)用開始于1964年美國噴氣動力實驗室用計算機對“徘徊者七號”太空船發(fā)回的月球照片進行處理，并獲得巨大成功?，F(xiàn)在圖

2、像處理技術(shù)已深入到科學(xué)研究、軍事技術(shù)、工農(nóng)業(yè)生產(chǎn)、醫(yī)學(xué)、氣象及天文學(xué)等領(lǐng)域?？茖W(xué)家利用人造衛(wèi)星可以獲得地球資源照片、氣象情況；醫(yī)生可以通過X射線或CT對人體各部位的斷層圖像進行分析。但在許多情況下圖像信息會受到各種各樣噪聲的影響，嚴重時會影響圖像中的有用信息，所以對圖像的噪聲處理就顯得十分重要。因此我選擇圖像去噪方面進行了解及研究，現(xiàn)將自己已了解的知識進行匯總。目錄摘要2一、圖像濾波的應(yīng)用4二、均值濾波51、 均值濾波的思想1、 均值濾波的算法1、 均值濾波的實驗結(jié)果三、中值濾波71中值濾波的思想1中值濾波的算法1中值濾波的實驗結(jié)果四、維納濾波88維納濾波的思想8維納濾波的算法8維納濾波的實驗

3、結(jié)果五、小波變換9小波變換濾波的思想小波變換濾波的算法小波變換濾波的實驗結(jié)果六、Contourlet變換的圖像去噪11Contourlet變換的基本思想Contourlet變換的算法七、全變差正則化的Shearlet收縮去噪12Shearlet收縮去噪原理簡介Shearlet收縮去噪算法八、結(jié)果分析及自己的收獲12結(jié)果分析自己的收獲參考文獻13圖像濾波的應(yīng)用圖像是信息社會人們獲取信息的重要來源之一。在現(xiàn)實世界中，存在大量影響圖像質(zhì)量的因素。因此圖像在進行使用之前，一般都要經(jīng)過嚴格的預(yù)處理如去噪、量化、壓縮編碼等。噪聲的污染直接影響著對圖像邊緣檢測、特征提取、圖像分割、模式識別等處理，因此我們不

4、得不從各種角度進行探索以提高圖像的質(zhì)量。采用適當(dāng)?shù)姆椒ūM量消除噪聲是圖像處理中一個非常重要的預(yù)處理步驟。現(xiàn)在圖像處理技術(shù)已深入到科學(xué)研究、軍事技術(shù)、工農(nóng)業(yè)生產(chǎn)、醫(yī)學(xué)、氣象及天文學(xué)等領(lǐng)域?？茖W(xué)家利用人造衛(wèi)星可以獲得地球資源照片、氣象情況；醫(yī)生可以通過X射線或CT對人體各部位的斷層圖像進行分析。但在許多情況下圖像信息會受到各種各樣噪聲的影響，嚴重時會影響圖像中的有用信息，所以對圖像的噪聲處理就顯得十分重要。為了改善圖像質(zhì)量，從圖像中提取有效信息，必須對圖像進行去噪預(yù)處理。根據(jù)噪聲的頻譜分布的規(guī)律和統(tǒng)計特征以及圖像的特點，出現(xiàn)了多種多樣的去噪方法。經(jīng)典去噪方法有：空間域合成法、頻域合成法和最優(yōu)線性合

5、成法等。與之相適應(yīng)的出現(xiàn)了許多應(yīng)用方法：如均值濾波器、中值濾波器、低通濾波器、維納濾波器、最小失真法等。這些方法廣泛應(yīng)用，促進了數(shù)字信號處理的極大發(fā)展。顯著提高了圖像質(zhì)量。近年來，小波變換去除噪聲的方法得到廣泛的應(yīng)用。與傳統(tǒng)的去噪方法相比。它利用的是非線性域值，在時間域和頻率域同時具有良好的局部化性質(zhì)，而且時窗和頻窗的寬度可以調(diào)節(jié)。對高頻成分采用逐漸精細的時域或空域取樣步長，從而可以聚焦到對象的任意細節(jié)。因此可以提高散斑高散射特性的噪聲對比度，很好的消除散斑噪聲9。小波變換去除噪聲的方法在不斷地發(fā)展，去噪方法很多，如非線性小波變換閾值法去噪、小波變換模極大值去噪及基于小波變換域的尺度相關(guān)性去噪

6、法均值濾波均值濾波的基本思想是將某像素鄰域內(nèi)的各點的灰度平均值來代替該像素原來的灰度級。通常鄰域都取成N*N的方形窗口,其降噪平滑后的圖像為取力s(ij)其中，s是點(x,y)鄰域內(nèi)的點集，M是點集S中的總點數(shù)。將算術(shù)均值濾波器作一改進，將某像素鄰域內(nèi)的各點的灰度加權(quán)平均值來代替該像素原來的灰度值，得到加權(quán)均值濾波。從權(quán)值上看，灰度越接近中心像素其權(quán)值越大。加權(quán)平均的算法可表示為gfey)=O-jlES其中，w(i,j)是權(quán)值，表示其所起作用的大小。均值濾波的算法用MATLA段計均值濾波算法為:(1)高斯噪聲sample=imread(lenna.jpg);sgray=rgb2gray(sam

7、ple);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(

8、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);I

9、gray=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);t

10、itle(subplot(3,2,4);imshow(K2);title(subplot(3,2,5);imshow(K3);title(subplot(3,2,6);imshow(K4);title(2.3均值濾波的實驗結(jié)果模板尺寸為3*3的濾波后圖像);模板尺寸為5*5的濾波后圖像);模板尺寸為7*7的濾波后圖像);模板尺寸為9*9的濾波后圖像);(1)高斯噪聲的濾波結(jié)果j;u 收 Jul*1 茹本諾更互叫 wdULnimjwTextGo=iRndlselEi(raj=rib2(Jubplot(3j2aubplot(3f2yl=filter2(f2=filter2(f3=filter2(f

11、4=filter2(fubplot(Sj2,ubplot(3,工ubplot。2Tubplot%2j:erencepagedocfilter2結(jié)果分析：通過通過圖像的識別率可以看出：模板越大去噪效果越好，但圖像模糊度也隨之增加。(2)椒鹽噪聲去噪結(jié)果,Editor-G:2C3iaoyanzjaosheng.m由圖可得：模板越大，去噪效果越好，因此使用均值濾波去噪時選用的模板尺寸（鄰域半徑）越大效果越好。三、中值濾波中值濾波的基本思想此方法是用該像素的相鄰像素的灰度中值來代替像素值。是一種典型的排序濾波器。1=median其中，是鄰域中各點的灰度值。所以，中值濾波后的圖像g(x,y)是以(x,y

12、)為中心的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(原始圖形);subplo

13、t(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中值濾波的實驗結(jié)果通過結(jié)果可以看出，中值濾波對于椒鹽噪聲的去

14、噪能力比對高斯噪聲的去噪能力好。四、維納濾波維納濾波的基本思想假定線性濾波器的輸入為有用信號和噪聲之和，兩者均為廣義平穩(wěn)過程且知它們的二階統(tǒng)計特性，維納根據(jù)最小均方誤差準則(濾波器的輸出信號與需要信號之差的均方值最小)，求得了最佳線性濾波器的參數(shù)，這種濾波器被稱為維納濾波器。維納濾波器的優(yōu)點是適應(yīng)面較廣，無論平穩(wěn)隨機過程是連續(xù)的還是離散的，是標(biāo)量的還是向量的，都可應(yīng)用。對某些問題，還可求出濾波器傳遞函數(shù)的顯式解，并進而采用由簡單的物理元件組成的網(wǎng)絡(luò)構(gòu)成維納濾波器。維納濾波器的缺點是，要求得到半無限時間區(qū)間內(nèi)的全部觀察數(shù)據(jù)的條件很難滿足，同時它也不能用于噪聲為非平穩(wěn)的隨機過程的情況，對于向量情況

15、應(yīng)用也不方便。因此，維納濾波在實際問題中應(yīng)用不多。實現(xiàn)維納濾波的要求是：(1)輸入過程是廣義平穩(wěn)的；(2)輸入過程的統(tǒng)計特性是已知的。維納濾波的算法維納濾波的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);小波變換小波變換濾波的基本思想近年來，小波理論得到了非常迅速的發(fā)展，由于其良好的時頻特性，實際

16、應(yīng)用非常廣泛將小波用于去噪也得到很多學(xué)者的重視。小波去噪利用了不同中心頻率的帶通濾波器對信號濾波，把主要反映噪聲頻率的那些尺度的系數(shù)去掉，再把剩余各尺度的系數(shù)綜合起來做反變換，從而使噪聲得到較好的抑制。小波變換具有時頻局域化特性，能夠檢測到局部突變的邊緣特性，而且可將圖像結(jié)構(gòu)和紋理分別表現(xiàn)在不同分辨率層次上。其中，小波域值算法是利用信號和噪聲小波系數(shù)幅值上的差異，通過選擇合適的域值對小波系數(shù)進行處理，以達到去除噪聲又保留有用信號的目的。此方法中，小波系數(shù)域的處理方法及閾值的估計是兩個關(guān)鍵技術(shù)。在此方法基礎(chǔ)上又開發(fā)出了，基于Bayes方法的小波去噪和基于假設(shè)檢驗的小波去噪。小波變換濾波的算法小波

17、變換的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維納濾波的實驗結(jié)果原始圖像）;加椒鹽噪聲后的圖像）;加高斯噪聲后的圖像）

18、;維納濾波1（椒鹽）；維納濾波2（高斯）；通過實驗結(jié)果可以看出：維納濾波具有較大的局限性，但其對高斯噪聲的去噪效果相對較好。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(

19、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去噪）；第一次閾值去噪圖像）；第二次閾值去噪圖像）；原始圖像）；加入高斯噪聲后的圖像）

20、第二次濾波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通過小波變換濾波的實驗結(jié)果通過小波分解的圖像以及全局閾值濾波后的圖像如下所示:FilflEditTextFilflEditText1.01.0subplot(2subplot1 12 2subplot(2(2subplot(2fig

21、ure;subplot(2subplot2subplot22subplot(2(2J5=n?dfilJ6=Piedfilfigure：subplot(2Contourlet變換的圖像去噪Contourlet變換的基本思想Contourlet變換可以滿足曲線的各向異性尺度關(guān)系，能夠很好地抓住圖像的幾何結(jié)構(gòu)。Contourlet變換是一種新的圖像二維表示方法，具有多分辨率、局部定位、多方向性、近鄰界采樣和各向異性等性質(zhì)，其基函數(shù)分布于多尺度、多方向上，少量系數(shù)即可有效地捕捉圖像中的邊緣輪廓，而邊緣輪廓正是自然圖像中的主要特征。Contourlet變換的基本思想是首先用一個類似小波的多尺度分解捕捉邊

22、緣奇異點，再根據(jù)方向信息將位置相近的奇異點匯集成輪廓段。選用Burr和Adelson于1983年提出的拉普拉斯塔式濾波器結(jié)構(gòu)(LP)對圖像多分辨率分解來捕捉奇異點。Contourlet變換的算法由于Contourlet變換缺乏平移不變性因此在應(yīng)用它進行去噪時會帶來人為的視覺效果。為克服這些人為視覺效果這里引入了基于循環(huán)平移的平移不變Contourlet去噪算法。假設(shè)含噪圖像f(i,j),0=i,j=N。定義F為平移n位的平移算子Fnf(i,j)=f(in,jn)mod(N)1由于F是一一X應(yīng)的所以其反變換因子可以設(shè)為F_n=Fn，一幅圖像中包含許多奇異點對某個奇異點來說是效果最佳的平移量”而對

23、另一個奇異點來說效果可能很差。為了解決這個問題在一定范圍內(nèi)進行循環(huán)平移運算再將平移運算的所有結(jié)果求平均值。將經(jīng)過小波變換并經(jīng)過全局閾值處理的圖像再用中值濾波進行第二次濾波，以加強濾波效果。結(jié)果如下所示：以Figure2 2I1=1FileEditViewInsertToolsDesktopWindowHelp丁/I相)“十|日國|口Editor-Editor-6 6 旌口轉(zhuǎn)三內(nèi)廢像房03|笛勺里J圖唱-1QI+|j-J3=wdencmpCgbl-J4=wdencmp(jbl-J5=medfilt2(J33-J6=medfilt2(J4J-figure;-subplot(2j3j1)-subpl

24、ot(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)進行循環(huán)平移；2)對每次平移后的圖像做離散Contourlet變換，可以得到各尺度各方向上的Contourlet變換系數(shù)c。3)Contourlet變換系數(shù)c進行閾值處理，得到去噪后的系數(shù)cAo4)利用c

25、A進行離散Contourlet反變換，得到去噪后圖像fA。全變差正則化的Shearlet收縮去噪Shearlet收縮去噪原理簡介收縮法是變換域中最為廣泛的圖像去噪方法。閾值收縮主要基于以下事實，比較大的系數(shù)一般都是以實際信號為主，而比較小的系數(shù)則很大程度上是噪聲。即在系數(shù)中，低頻分量中含有大量的信息，應(yīng)該給予保留；同時在高頻分量中，一些絕對值大的重要的系數(shù)并不是噪聲，而是邊緣信息，也應(yīng)保留。硬閾值收縮法的收縮函數(shù)為：式中r為收縮的閾值。硬閾值是將絕對值小于閾值的系數(shù)置為零，而將絕對值大于閾值的系數(shù)不加任何處理給予保留。硬閾值的優(yōu)點是運算速度快，能取得較好的去噪效果。但由于其收縮函數(shù)是不連續(xù)的，

26、所以在圖像中會產(chǎn)生許多“人為”噪聲，圖像會出現(xiàn)振鈴、偽吉布斯效應(yīng)等視覺失真。圖3給出了硬閾值函數(shù)示意圖。圖中的小正方形是代表絕對值小于閾值的系數(shù)，因為在硬閾值處理中將小于閾值的系數(shù)直接置零，這樣就丟失了圖像的一部分信息。偽吉布斯現(xiàn)象可以視為一種整體振蕩，而全變差正則化對整體振蕩有很好的抑制作用。為了避免硬閾值收縮去噪丟失的圖像信息，以及更好地消除偽吉布斯現(xiàn)象，本文引入全變差正則化。Shearlet收縮去噪算法1)Shearlet硬閾值估計AA對含噪圖像進行Shearlet變換，得到shearlet系數(shù)djl(f)和尺度系數(shù)Cjl(f);J,lJ,l(2)Shearlet系數(shù)進行硬閾值去噪，用愀

27、表示保留系數(shù)的指標(biāo)集。用Shearlet反變換得到圖像的初始估計f。2)迭代修正初始化lx=1,N代表迭代次數(shù)，%=Lmax,dj=LmaX/N。計算Jp(f)的次梯度2Jp(f)；(2)將17fJp(f)投影到V空間彳導(dǎo)到R(fJp(f)；利用式(13)計算f(nW=fnKnFV(VfJp(f)；(4)n=n+l,九0書=4-d九，當(dāng)兒n之。時轉(zhuǎn)至步驟2)中第(1)步循環(huán)執(zhí)行，否則結(jié)束迭代退出循環(huán)。結(jié)果分析及自己的收獲結(jié)果分析通過所有濾波設(shè)計及結(jié)果對比，我們可以得到:第一次wrco日12去瞟50005000在處理服從高斯分布的一類噪聲時，維納濾波與中值濾波去除效果較好一些；中值濾波對于去除椒

28、鹽噪聲效果好，而維納濾波去除效果差，中值濾波對于去除椒鹽噪聲效果明顯，是因為椒鹽噪聲只在畫面上的部分點隨機出現(xiàn)，而中值濾波根據(jù)數(shù)據(jù)排序，將未被污染的點代替噪聲點的值的概率較大，所以抑制效果好。對點、線和尖頂較多的圖像不宜采用中值濾波，因為一些細節(jié)點可能被當(dāng)成噪聲點。而且均值濾波，中值濾波都在一定程度上模糊了圖像邊緣，提高濾波效果是以模糊邊緣為代價的。小波分解可以把圖像分層次按照小波基展開，并且可以根據(jù)圖像的性質(zhì)及給定的處理標(biāo)準確定展開到哪一級為止，還可以把細節(jié)分量和近似分量分開。由圖可得：wrcoef2和wpdencm薛函數(shù)可以有效地進行去噪處理。而且在小波變換的基礎(chǔ)上再利用中值濾波可以進一步

29、增加去噪效果。因此對邊緣具有很好的保留，不會模糊邊緣。因此濾波效果最好。在分析過程中可以看出，針對不同類型的噪聲需要相應(yīng)的濾波去噪算法才能取得較好的均值濾,波圖像100200效果，才能使后續(xù)的圖像處理工作得以更加優(yōu)質(zhì)的進行。自己的收獲通過這次對圖像濾波去噪處理的了解，以及對所有方法基本思想的了解，進一步激發(fā)了自己對該課題的興趣，也拓寬了自己的知識面，也鞏固了自己所學(xué)的數(shù)字信號處理方面相關(guān)的知識。在整個課題中，通過對不同的去噪算法的分析，以及結(jié)果的比較，更好地了解了在不同情況下應(yīng)利用不同的方法來解決問題，了解每一個濾波器的優(yōu)缺點，在遇到不同問題時，可以根據(jù)各個濾波方法的優(yōu)缺點選擇合適的濾波方法，從而使我們的效果達到最好，在了解濾波器的原理的過程中，我了解到許多數(shù)字信號處理的知識。對于數(shù)字圖像信號，噪聲表現(xiàn)為或大或小的極值，這些極值通過加減作用于圖像象素的真實灰度值上，在圖像造成亮、暗點干擾，極大降低了圖像質(zhì)量，影響圖像復(fù)原、分割、特征提取、圖識別等后繼工作的進行。要構(gòu)造一種有效抑制噪聲的