數(shù)字圖像去噪典型算法仿真與分析

數(shù)字圖像去噪典型算法仿真與分析數(shù)字圖像去噪典型算法仿真與分析數(shù)字圖像去噪典型算法仿真與分析數(shù)字圖像去噪典型算法仿真與分析編制僅供參考審核批準(zhǔn)生效日期地址：電話：傳真:郵編:數(shù)字圖像去噪典型算法仿真與分析個人信息*********摘要：圖像去噪是數(shù)字圖像處理中的重要環(huán)節(jié)和步驟。本文首先介紹了常見的圖像噪聲；然后，在介紹圖像去噪的基本方法和原理的基礎(chǔ)上，討論了均值濾波、中值濾波和維納濾波三種典型的圖像去噪方法；最后，對包含有高斯噪聲和椒鹽等噪聲的圖像進(jìn)行去噪，并對其去噪效果進(jìn)行了仿真和分析比較，得出了三種方法各自的適用性特點。關(guān)鍵詞：圖像去噪；均值濾波；中值濾波；維納濾波SimulationandAnalysisofImageDe-noisingMethodsinDigitalImageName:***(個人信息****)Abstract:Imagedenoisingisoneofthemostimportantpartsandstepsofimageprocessing.Firstly,thepaperintroducesthecommonimagenoise.Then,basedontheprincipleandmethodsofeliminatingimagenoise,itdiscussesmeanfiltering,medianfiltering,andWienerfilteringwhicharetypicalimagedonoising.Finally,itusesthesemethodstoeliminateimagenoisewhichcontainsGaussiannoiseandsalt&peppernoise.Andthroughcomparingandanalyzingtheeffectofthesemethods,itconcludestheapplicabilityofeachmethodindifferentapplication.Keywords:imagedenoising;meanfiltering;medianfiltering;Wienerfiltering0引言數(shù)字圖像是現(xiàn)代人們獲取信息的主要來源。由于成像系統(tǒng)、傳輸介質(zhì)和記錄設(shè)備等的不完善，數(shù)字圖像在其形成、傳輸記錄過程中往往會收到多種噪聲的污染。一般來說，現(xiàn)實中的圖像都是帶噪圖像。噪聲使圖像變得模糊，甚至淹沒圖像特征，給后面圖像區(qū)域分割、分析判斷等工作帶來了困難。因此，在邊緣檢測、圖像分割、特征提取、模式識別等高層次處理之前，選用適當(dāng)?shù)姆椒ūM量地去除噪聲干擾是一個非常重要的環(huán)節(jié)和步驟，也一直是圖像處理研究領(lǐng)域進(jìn)行的主要課題之一。在圖像去噪的研究方面，國內(nèi)外的一些學(xué)者提出了大量的算法，如一些典型的圖像去噪算法，均值濾波、中值濾波、維納濾波、小波變換等[1-5]。近年來，一些改進(jìn)的、新型的算法也被許多學(xué)者提出。文[6]提出了一種改進(jìn)的均值濾波算法，該算法針對均值濾波在抑制噪聲的過程中會損失圖像的邊緣等細(xì)節(jié)信息問題，在計算局部窗口內(nèi)中心像素灰度均值時，既考慮了窗口內(nèi)各像素與中心像素間的灰度差異，又顧及了窗口內(nèi)各像素與中心像素間的距離；但是容易導(dǎo)致圖像細(xì)節(jié)因過度平滑而變得模糊。對于復(fù)雜的噪聲圖像，文[7]提出了一種基于同性質(zhì)點個數(shù)的噪聲點檢測算法，這是一種改進(jìn)的中值濾波算法，該算法對于脈沖噪聲來說不僅在濾除噪聲方便有較好的效果，保持圖像細(xì)節(jié)信息方面也取得了一定的成就；但對未知噪聲類型的圖像進(jìn)行濾除時其效果就不明顯。針對經(jīng)典維納濾波器存在的不足，文[8]提出了一種新的自適應(yīng)維納濾波器，該濾波器能夠根據(jù)不同的圖像特性在給定的多個模板之間自適應(yīng)的選擇模板，使得濾波效果更加理想；但是對于不太復(fù)雜的圖像，新的自適應(yīng)維納濾波和普通維納濾波相比，改善空間不是很大，效果不是很明顯。綜上，現(xiàn)有的圖像去噪方法大致可以分為兩類：一類是空間域方法，另一類是頻率域方法。本文主要圍繞著空間域的去噪方法，分析幾種典型的空間域去噪方法，研究總結(jié)各算法的優(yōu)缺點、適用性及處理效率等，并通過MATLAB進(jìn)行仿真和分析比較。1噪聲噪聲可以理解為妨礙人的視覺器官或系統(tǒng)傳感器對所接收圖像源信息進(jìn)行理解或分析的各種因素。噪聲在理論上可以定義為“不可預(yù)測，只能用概率統(tǒng)計方法來認(rèn)識的隨機誤差”[9]。圖像系統(tǒng)中的常見噪聲依據(jù)噪聲產(chǎn)生的原因，將經(jīng)常影響圖像質(zhì)量的噪聲源分為三類：阻性元器件內(nèi)部產(chǎn)生的高斯噪聲；光電轉(zhuǎn)換過程中的泊松噪聲（椒鹽噪聲）；感光過程中產(chǎn)生的顆粒噪聲。根據(jù)噪聲和信號的關(guān)系也可以將其分為兩種形式：一類是加性噪聲，另一類是乘性噪聲。圖像去噪的意義噪聲在圖像處理中是一個非常重要的問題，它對圖像的輸入、采集、處理的各個環(huán)節(jié)以及最終的輸出結(jié)果都會產(chǎn)生一定的影響。特別是在圖像的輸入、采集過程中，噪聲是個十分關(guān)鍵的問題，若輸入伴有較大噪，必然影響之后的處理以及最終的處理效果。因此，任何一個良好的圖像處理系統(tǒng)，無論是模擬處理還是計算機處理無不把減少最前一級的噪聲作為主攻目標(biāo)。去噪處理已成為圖像處理中極其重要的環(huán)節(jié)和步驟。2圖像去噪典型算法現(xiàn)有的圖像去噪方法大致可以劃分為兩類：一類是空間域方法，主要采用各種圖像平滑模板對圖像進(jìn)行卷積處理，以達(dá)到壓抑或消除噪聲的目的；另一類是頻率域方法，主要通過對圖像進(jìn)行變換后，選用適當(dāng)?shù)念l率帶通濾波器進(jìn)行濾波處理，再經(jīng)過反變換得到去噪后圖像。本文的工作主要圍繞著空間域方法，對三種典型的空間域去噪方法，均值濾波、中值濾波和維納濾波，進(jìn)行討論、仿真與分析。均值濾波均值濾波也稱為線性濾波，其采用的主要方法為鄰域平均法[1-4]。均值濾波的基本思想是用領(lǐng)域幾個像素值的均值來代替原圖像中的各個像素值，其領(lǐng)域的選取通常為像素的4領(lǐng)域和8領(lǐng)域。假定有一幅由L*H個像素組成的圖像f(x,y)，令Sxy表示中心在(x,y)點，尺寸為M*N的滑動模板窗口。均值濾波過程就是計算模板中所有像素的均值g(x,y)，由公式1給出；然后將g(x,y)代替該點(x,y)原來的像素值f(x,y)，得到去噪后圖像F(x,y)，如公式2所示。(1)(2)這個操作可以用其系數(shù)1/M*N的卷積模板來實現(xiàn)。由上式可知，經(jīng)過均值處理后，噪聲的均值不變，方差變小，說明噪聲的強度減弱了，即噪聲得到了抑制。當(dāng)然，這種方法在平滑噪聲的同時，也會模糊信號的細(xì)節(jié)和邊緣，即在清除噪聲的同時也會對圖像的高頻細(xì)節(jié)成分造成破壞和損失，使圖像模糊。中值濾波中值濾波是一種非線性信號處理方法，它的基本原理是把數(shù)字圖像或數(shù)字序列中的一點的值用該點的一個鄰域中的各點值的中值代替。通俗地講中值濾波就是用一個活動窗口沿圖象移動，窗口中心位置的象素灰度用窗口內(nèi)所有象素灰度的中值來代替。其算法步驟為：首先確定一個以某個像素(x,y)為中心點的鄰域，常用的有方形、十字形和圓形等鄰域；然后將鄰域中的各個像素的灰度值進(jìn)行排序，取其中間值作為中心點像素灰度的新值，這里的鄰域通常被稱為窗口，這個窗口可以是3*3、5*5等；當(dāng)窗口在圖像中上下左右進(jìn)行移動后，利用中值濾波算法可以很好地對圖像進(jìn)行平滑處理[1-4]。它是一種鄰域運算，類似于卷積，但計算的不是加權(quán)求和，而是把鄰域中的像素按灰度級進(jìn)行排序然后選擇該組的中間值作為輸出像素值。即用該像素的相鄰像素的灰度中值來代替該像素的值，由公式3給出。(3)領(lǐng)域的大小決定在多少個數(shù)值中求中值，窗口的形狀決定在什么樣的幾何空間中去元素計算中值。窗口的大小和形狀有時對濾波效果影響很大。維納濾波維納(Wiener)濾波器是一種自適應(yīng)濾波器，它根據(jù)圖像的局部方差調(diào)整濾波器的輸出，它的最終目標(biāo)是使恢復(fù)圖像g(x,y)與原始圖像f(x,y)的均方差最小。所謂自適應(yīng)濾波，就是利用潛意識可以獲得的濾波器參數(shù)等結(jié)果，自動的調(diào)節(jié)現(xiàn)時刻的濾波器參數(shù)，以適應(yīng)信號的噪聲未知的或隨時間變換的統(tǒng)計特性。Wiener算法[1][5]首先估計出像素的局部矩陣均值和方差：(4)然后，對于每一個像素利用Wiener濾波器估計出其像素值：(5)這里，v2是圖像中噪聲的方差。維納濾波去噪方法根據(jù)圖像的局部方差來調(diào)整濾波器的輸出，當(dāng)局部方差大時，濾波器的平滑效果較弱；當(dāng)局部方差小時，濾波器平滑效果較強。它比線性濾波器具有更好的選擇性，可以更好地保存圖像的邊緣和高頻細(xì)節(jié)信息。3仿真結(jié)果及分析利用MATLAB圖像處理工具箱的函數(shù)可以非常方便地實現(xiàn)以上三種濾波方法。但是，為了加深對這三種濾波方法的理解和運用，本文采用了自編的MATLAB程序來進(jìn)行仿真，其中Wiener濾波器采用MATLAB自帶函數(shù)。為了比較以上幾種不同圖像去噪的方法的優(yōu)劣，本文對同一幅圖像利用MATLAB工具箱人為添加不同的圖像噪聲；然后分別采用均值濾波、中值濾波、維納濾波等平滑處理辦法對圖像進(jìn)行去噪復(fù)原，得到了豐富的實驗結(jié)果。限于篇幅，本文只給出了對添加有均值為零的高斯噪聲和椒鹽噪聲的Leda圖像進(jìn)行三種濾波算法處理的結(jié)果，如圖1、圖2、圖3所示。圖1Leda原始圖像和添加噪聲圖像圖2Leda添加高斯噪聲的仿真結(jié)果圖3Leda添加椒鹽噪聲的仿真結(jié)果通過對實驗結(jié)果進(jìn)行分析比較，從主觀上可以得出以下結(jié)論：1)對于均值濾波，由圖2可以看出，圖像中噪聲含量很少，說明均值濾波對均值為零的高斯噪聲有比較好的抑制作用。由圖3可以看出均值濾波對椒鹽噪聲的去除效果不明顯，噪聲雖然得到了一定的抑制，但是圖像邊緣變得比較模糊，而且隨著濾波器尺寸的增加，圖像的細(xì)節(jié)銳化程度相應(yīng)降低。因為均值濾波只是將某點出現(xiàn)的噪聲強度，讓周圍的數(shù)據(jù)平均分擔(dān)了，所以得到的結(jié)果是噪聲幅值減小，但是噪聲點的顆粒面積同時變大，而椒鹽噪聲均值不為零。2)對于中值濾波，由圖2可以看出，中值濾波對高斯噪聲的抑制效果不明顯，因為高斯噪聲使用隨機大小的幅值污染所有的點，因此無論怎樣進(jìn)行數(shù)據(jù)選擇得到的始終還是被污染的值。由圖3可以看出中值濾波對去除椒鹽噪聲效果明顯，因為椒鹽噪聲只在畫面中的部分點上隨機出現(xiàn)，根據(jù)中值濾波原理可知，通過數(shù)據(jù)排序的方法將圖像中未被噪聲污染的點代替噪聲點的值的概率比較大，因此噪聲的抑制效果很好。中值濾波與均值濾波相比，對畫面清晰度的保持效果明顯，缺點是因為涉及大量排序運算，運算速度較慢，對圖像的實時處理有影響，對一些細(xì)節(jié)多，特別是點、線、尖頂細(xì)節(jié)較多的圖像不宜采用中值濾波的方法，因為細(xì)節(jié)點有可能被當(dāng)成了噪聲點。3)對于維納濾波，由圖(1-3)的對比可以看出，維納濾波器去除高斯白噪聲的效果較好，去除椒鹽噪聲的效果較差，從圖中可以清晰的看到一些明顯的噪聲點。維納濾波能較好保存圖像的邊緣和高頻細(xì)節(jié)信息，但去噪后的圖像仍略顯模糊。為了更客觀地說明各算法的濾波性能，采用峰值信噪比PSNR(PeakSignalNoiseRatio)作為客觀評價的尺度，PSNR值越大，表示恢復(fù)圖像與原始圖像越接近，PSNR定義為：(6)其中，M,N表示圖像的尺寸，f(x,y)表示原始圖像的灰度值，g(x,y)表示濾波后的圖像灰度值。各算法對噪聲圖像處理后的PSNR如表1所示：表1Leda峰值信噪比測試結(jié)果單位：dBPSNR噪聲類型均值濾波中值濾波維納濾波高斯噪聲椒鹽噪聲由表中PSNR數(shù)據(jù)也可以得出，均值濾波對高斯噪聲抑制效果較好，中值濾波對椒鹽噪聲有明顯的抑制作用，維納濾波對高斯噪聲有較好的抑制作用，而對椒鹽噪聲的抑制作用比較差。4結(jié)束語本文只研究了空間域的幾種典型的去噪方法。本文在MATLAB環(huán)境下，結(jié)合不同的噪聲類型，分析討論了均值濾波、中值濾波和維納濾波三種典型的圖像去噪方法原理和應(yīng)用范圍，以及它們在去除幾種典型噪聲方面的性能優(yōu)劣和視覺效果對比，得出了結(jié)論。抑制噪聲和保持圖像中的細(xì)節(jié)往往是一對矛盾，也是圖像處理中尚未很好解決的問題。在實際應(yīng)用中，處理噪聲圖像前應(yīng)明確以下兩點：1)圖像受到何種類型的噪聲干擾；2)受噪聲干擾的程度。然后選擇適當(dāng)?shù)娜コ肼暭夹g(shù)，減少圖像去噪過程的盲目性。如何將各種算法進(jìn)行改進(jìn)，得到一種新的兼有細(xì)節(jié)保護(hù)和噪聲抑制的優(yōu)良特性的算法，使其去噪效果更為理想，應(yīng)用更為廣泛，這是圖像處理工作的目標(biāo)，也是今后的研究方向。參考文獻(xiàn)：賽地瓦爾地·買買提.基于Matlab的幾種圖像去噪方法研究[J].河南科學(xué),2013,31(9):1387-1390王文慶,晏婷.基于空間域的圖像去噪方法比較研究[J].西安郵電學(xué)院學(xué)報,2012,17(2):75-79胡蕾,張偉,覃慶炎.幾種去噪算法的應(yīng)用分析[J].信息技術(shù),2007(7):81-83.盛仲飆.基于Matlab的圖像去噪算法研究