基于混合變分模型的圖像去噪

周 楊，溫興平，張麗娟，王 軍

昆明理工大學(xué) 國土資源工程學(xué)院，昆明 650093

基于混合變分模型的圖像去噪

周 楊，溫興平，張麗娟，王 軍

昆明理工大學(xué) 國土資源工程學(xué)院，昆明 650093

1 引言

圖像在獲取、傳輸、轉(zhuǎn)換等過程中不可避免地會受到不同程度的噪聲影響，導(dǎo)致圖像分割、圖像識別、圖像分析、圖像理解等處理結(jié)果出現(xiàn)較大誤差，所以圖像去噪成為圖像處理研究領(lǐng)域中的熱點問題[1]。

為了在去除噪聲的同時較好地保持圖像的清晰度、對比度、紋理信息等，大量的去噪方法相繼出現(xiàn)[2]。傳統(tǒng)圖像去噪分為空間域處理和頻率域處理：空間域去噪的主要思想是對噪聲圖像使用平滑濾波進行平滑處理，主要方法包含中值濾波、均值濾波和非局部平均濾波等[3-5]。頻率域去噪主要通過Fourier變換、小波變換等將空間域信號轉(zhuǎn)化為頻率域信號，然后對頻率域信號進行低通濾波，最后再轉(zhuǎn)化為空間域信號[6-7]，但基于空間域和頻率域的圖像去噪不能保持較好的幾何機構(gòu)。近年來，基于偏微分方程（Partial Differential Equation，PDE）的圖像處理方法發(fā)展迅速，在圖像去噪和去模糊等方面有著良好的應(yīng)用效果[8]。其中基于全變分（Total Variation，TV），其是通過建立數(shù)學(xué)模型得到能量泛函，選取適當?shù)膮?shù)，對圖像進行反復(fù)迭代，通過解方程得到最優(yōu)解從而得到很好的去噪效果。但其得到的圖像無法較好地保持邊緣，會產(chǎn)生“階梯”效應(yīng)。因此，為了改進算法缺陷，針對這兩類經(jīng)典模型存在的不足，人們提出了很多改進方法，文獻[9]提出一種廣義TV去噪模型，該模型能克服假邊緣的產(chǎn)生，在去噪的同時保持了邊緣細節(jié)，但該模型參數(shù)為常數(shù)，其取值對實驗結(jié)果比較敏感。為克服該不足，文獻[10]提出將梯度作為圖像邊緣的檢測算子計算參數(shù)，文獻[11]采用差曲率代替梯度作為邊緣判別算子計算參數(shù)，文獻[12]將非局部均值濾波引入?yún)?shù)中，這些改進方法都可根據(jù)圖像的局部特征自適應(yīng)選取參數(shù)的值，在不同程度上改善了圖像去噪效果。然而，這些算法也存在缺陷，雖然能夠很好地保持圖像的邊緣信息，但圖像中細小的邊緣會被過度平滑掉。

為了有效抑制噪聲，獲得更好的視覺效果，提出一種基于混合變分模型的圖像去噪方法。首先將調(diào)和模型和全變分模型進行融合，增強模型的去噪性能，然后自適應(yīng)選取加權(quán)參數(shù)，在有效減少噪聲的同時，保護圖像的邊緣細節(jié)，最后采用多幅含噪圖像進行仿真對比實驗。結(jié)果表明，相對于其他去噪模型，相同條件下，本文模型不僅獲得了較理想的視覺效果，提高了圖像質(zhì)量，而且獲得了更高的峰值信噪比，有效降低了均方根誤差，驗證了算法的有效性。

2 混合變分模型

2.1 全變分模型

全變分模型去噪的基本思想是對含噪圖像建立數(shù)學(xué)模型，得到相應(yīng)的能量泛函，通過反復(fù)迭代使得能量泛函趨于極小值，這時可以認為噪聲含量達到極小值，得到此時的解，即為去噪后的圖像。設(shè) f(x)∶Ω→R為圖像空間Ω上的灰度圖像，u(x)∶Ω→R為去除噪聲后的圖像，相似度函數(shù)定義如下：

基于上述定義，全變分模型的能量泛函定義為：

采用交替優(yōu)化方法得到關(guān)于μ(x)的Euler-Lagrange方程和ωk+1的廣義軟閾值公式，μ(x)即為所得的去噪后圖像。

2.2 調(diào)和模型

式中，?u表示圖像的梯度，ω為圖像的定義域，λ為Lagrange乘子。

假定μ滿足Neumann邊界條件Ё?u/?ν=0，v表示?Ω的外法向量，則調(diào)和模型存在唯一解且解滿足最優(yōu)性條件為：

2.3 混合變分模型

式中，ω為權(quán)函數(shù)，且0≤ω＜1。

從式（18）可知，當ω=0時，相當全變分模型，其可以較好地保持邊緣等細節(jié)信息；當ω=1時，相當于平滑較好的調(diào)和模型；當0＜ω＜1時，采用全變分模型處理圖像的邊緣細節(jié)，并采用調(diào)和模型處理平滑區(qū)域。

3 仿真實驗

3.1 圖像來源及仿真環(huán)境

為了驗證本文模型的有效性，選擇三種類型的圖像作為仿真對象，具體如圖1所示。仿真實驗均在Intel?CoreTM2.85 GHz CPU，2 GB RAM，Windows 7操作系統(tǒng)，采用Matlab 2012編程實現(xiàn)仿真實驗。

圖1 含噪圖像

3.2 結(jié)果與分析

分別采用調(diào)和模型、全變分模型和本文模型對圖1中圖像進行去噪處理，得到結(jié)果如圖2～4所示。對圖2～4進行分析可以得到如下結(jié)論：

（1）傳統(tǒng)調(diào)和模型在去除圖像中的噪聲時，會去除圖像的部分細節(jié)信息。

（2）相對于調(diào)和模型，全變分模型可以得到更優(yōu)的去噪效果，但在平坦區(qū)域仍存在不同程度的“階梯”效應(yīng)。

（3）相對于對比模型，本文模型不僅減少了平滑區(qū)域的殘余噪聲，降低了邊緣的模糊程度，同時可以較好保護圖像的邊緣細節(jié)，處理后圖像整體質(zhì)量較高，獲得更加理想的視覺效果。

圖2 lena圖像的去噪結(jié)果

圖3 Peppers圖像的去噪結(jié)果

圖4 Barbara圖像的去噪結(jié)果

3.3 客觀評價

在評價圖像去噪模型性能時，有許多客觀評價標準，本文采用峰值信噪比（PSNR）、均方根誤差（MSE）和去噪時間作為模型去噪性能的客觀評價指標[13]，設(shè)M×N大小的帶噪圖像為 f(m，n)，處理后的圖像為 f′(m，n)，其中m=1，2，…，M，n=1，2，…，N，則其計算方法如下：

其中，Q是量化級數(shù)。

調(diào)和模型、全變分模型和本文模型的PSNR和MSE如表1所示。從表1可知，本文模型具有較高的峰值信噪比，有效降低了均方根誤差，因此在保護細節(jié)信息的同時。一定程度上緩解了圖像的梯度效果，降噪性能比對比模型有了一定程度的改善。

表1 各種模型的客觀評價結(jié)果對比

調(diào)和模型、全變分模型和本文模型的去噪時間如圖5所示，從圖5可知，相對于全變分模型，調(diào)和模型的去噪時間相對較少，而本文模型的去噪時間要少于全變分模型，略高于調(diào)和模型，但是可以較高質(zhì)量地恢復(fù)圖像，因此綜合性能最優(yōu)。

圖5 不同模型的去噪速度對比

3.4 與經(jīng)典去噪模型的性能對比

為了進一步測試本文模型的性能，采用當前經(jīng)典圖像去噪模型進行對比實驗，它們分別為：Contourlet變換、小波閾值去噪、BM3D、BiShrink以及SURE-LET[13-16]，所有模型的去噪效果評價指標如表2所示。從表2可知，相對于當前經(jīng)典去噪模型，本文模型的PSNR值得到明顯提高，MSE的值有所降低，對比結(jié)果再次證明了本文算法的有效性和優(yōu)越性。

表2 與經(jīng)典模型的性能對比

4 結(jié)束語

為了提高圖像質(zhì)量，有效去除噪聲，針對全變分模型存在的不足，提出了基于混合變分的圖像去噪模型，并通過仿真對模型的性能進行測試。仿真實驗結(jié)果表明，本文模型具有良好的去噪性能，提高了圖像質(zhì)量，獲得了更理想的視覺效果。

致謝 本文得到昆明理工大學(xué)地質(zhì)過程與礦產(chǎn)資源省創(chuàng)新團隊和昆明理工大學(xué)遙感地球化學(xué)學(xué)科方向團隊聯(lián)合資助，在此表示感謝。

[1]Al-amri S S，Kalyankar N V，Khamitkar S D A.Comparative study of removal noise from remote sensing image[J]. IJCSI International Journal of Computer Science Issues，2010，7（1）：32-36.

[2]Mohammad R H，Ahmad M O，Wang Chunyan.An edgeadaptive laplacian kernel for nonlinear diffusion filters[J]. IEEE Transactions on Image Processing，2012，21（4）：1561-1572.

[3]牛和明，杜茜，張建勛.一種自適應(yīng)全變分圖像去噪算法[J].模式識別與人工智能，2011，24（6）：798-803.

[4]侯榆青，張歡，史晶，等.一種改進的全變分自適應(yīng)圖像去噪模型[J].西北大學(xué)學(xué)報，2008，38（3）：371-373.

[5]Cho D，Bui T D.Multivariate statistical modeling for image denoising using wavelet transforms[J].IEEE Transactions on Signal Processing Image Communication，2005，20（1）：77-89.

[6]戴維，于盛林，孫栓.基于Contourlet變換自適應(yīng)閾值的圖像去噪算法[J].電子學(xué)報，2007，35（10）：1939-1943.

[7]劉文，吳傳生，許田.自適應(yīng)全變分圖像去噪模型及其快速求解[J].計算機應(yīng)用研究，2011，28（12）：4797-4800.

[8]Beck A，Teboulle M.Fast gradient-based algorithms for constrained total variation image denoising and deblurring problems[J].IEEE Transactions on Image Processing，2009，18（11）：2419-2434.

[9]Yin Ming，Liu Wei，Shui Jun，et al.Quaternion wavelet analysis and application in image denoising[J].Mathematical Problems in Engineering，2012，2012：1-21.

[10]蔡良師，鄭南山.基于小波系數(shù)相關(guān)性的圖像去噪[J].測繪科學(xué)，2012，37（1）：94-95.

[11]仵冀穎，阮秋琦.偏微分方程在圖像去噪中的應(yīng)用[J].計算機工程與應(yīng)用，2006，42（22）：69-71.

[12]高國榮，許錄平，馮冬竹.利用非抽樣Shearlet域GSM模型進行圖像去噪[J].武漢大學(xué)學(xué)報，2013，38（7）：778-782.

[13]賈建，焦李成，魏玲.基于概率模型的非下采樣Contourlet變換圖像去噪[J].西北大學(xué)學(xué)報：自然科學(xué)版，2009，39（1）：3-18.

[14]郭強，郁松年.基于三變量模型的剪切波去噪方法[J].自動化學(xué)報，2010，36（8）：1062-1072.

[15]Luisier F，Blu T，Unser M.A new SURE approach to image denoising：interscale orthonormal wavelet thresholding[J].IEEE Trans on Image Processing，2007，16（3）：593-606.

[16]王凱，劉甲甲，袁建英，等.基于小波子帶系數(shù)能量的優(yōu)化權(quán)值降噪[J].計算機應(yīng)用，2013，33（8）：2341-2345.

ZHOU Yang,WEN Xingping,ZHANG Lijuan,WANG Jun

Engineering Institute of Land and Resources,Kunming University of Science and Technology,Kunming 650093,China

In order to restrain the noise and obtain better visual effect,a novel image denoising method is proposed based on hybrid variation model.The harmonic model and the total variation model are integrated to enhance the filtering performance,and then the parameters are adaptively selected.When combinatorial coefficient is large,hybrid model is close to total variation model,and when combinatorial coefficient is small,hybrid model is close to harmonic model to effectively remove the staircase effect to protect the edge detail.The simulation experiment is carried out to test performance of model. The results show that,compared with other denoising models,the proposed model has better denoising effect and can better protect edge information under the same conditions.

image denoising;harmonic model;total variational model;adaptive denoising

為了有效抑制噪聲，獲得更好的視覺效果，提出了一種基于混合變分模型的圖像去噪方法。將調(diào)和模型和全變分模型進行融合，增強模型的去噪性能，根據(jù)自適應(yīng)選取組合系數(shù)，組合系數(shù)較大時偏向于全變分模型，較小時偏向于調(diào)和模型，這樣不僅可以有效去除階梯效應(yīng)，同時保護邊緣細節(jié)，采用仿真對比實驗以測試模型性能。結(jié)果表明，相對其他去噪模型，相同條件下，該模型取得更優(yōu)的去噪聲效果，提高了圖像的質(zhì)量。

圖像去噪；調(diào)和模型；全變分模型；自適應(yīng)去噪

A

TP317

10.3778/j.issn.1002-8331.1402-0417

ZHOU Yang,WEN Xingping,ZHANG Lijuan,et al.Image denoising based on hybrid variational filter model.Computer Engineering and Applications,2014,50（24）：183-186.

國家自然科學(xué)基金（No.41101343）；云南省高技術(shù)產(chǎn)業(yè)發(fā)展項目計劃。

周楊（1989—），女，碩士生，主要從事遙感地質(zhì)工作；溫興平（1970—），通訊作者，男，博士，教授，主要從事地球化學(xué)與遙感地質(zhì)研究工作；張麗娟（1988—），女，碩士生，主要從事遙感地質(zhì)工作；王軍（1987—），男，碩士生，主要從事遙感地質(zhì)工作。E-mail：wfxyp2008@gmail.com

2014-03-03

2014-05-13

1002-8331（2014）24-0183-04

CNKI網(wǎng)絡(luò)優(yōu)先出版：2014-07-16，http∶//www.cnki.net/kcms/doi/10.3778/j.issn.1002-8331.1402-0417.html