何 坤，鄭秀清，琚生根，張永來

（1. 四川大學(xué)計算機學(xué)院 成都 610065；2. 四川師范大學(xué)信息技術(shù)學(xué)院 成都 618300）

改進全變分的圖像去噪

何 坤1，鄭秀清2，琚生根1，張永來2

（1. 四川大學(xué)計算機學(xué)院 成都 610065；2. 四川師范大學(xué)信息技術(shù)學(xué)院 成都 618300）

為了彌補各向同性擴散去噪的非保邊性、異性擴散的耗時性，分析了各向同性和異性擴散的機理，根據(jù)噪聲對像素變化的影響，設(shè)計了新的擴散函數(shù)，理論上分析了該函數(shù)的擴散性能：對平滑區(qū)各向同性擴散，邊緣區(qū)實現(xiàn)各向異性擴散。在傳統(tǒng)全變分去噪的基礎(chǔ)上，提出了改進全變分的圖像去噪模型，運用固定點代算法設(shè)計了相應(yīng)的離散迭代函數(shù)。實驗結(jié)果表明，該算法在圖像平滑區(qū)進行各向同性擴散，繼承了各向同性的優(yōu)點，降低了傳統(tǒng)全變分的運行時間；在邊緣區(qū)實現(xiàn)了各向異性擴散保護了圖像邊緣，提高了圖像的峰值信噪比和視覺效果。

各向異性； 擴散函數(shù)； 圖像去噪； 各向同性； 全變分

人眼對圖像簡化、提取其結(jié)構(gòu)信息，實現(xiàn)對圖像的認知和理解。但圖像在獲取、保存及傳輸過程中不可避免地受到噪聲污染。噪聲和結(jié)構(gòu)信息均屬于高頻信息［1］，降低了視覺對圖像簡化的質(zhì)量，影響了其分析理解［2］。為了提高機器視覺對圖像分析理解的準確性，需對圖像進行保結(jié)構(gòu)信息的去噪預(yù)處理［3］。

最簡單的圖像去噪處理是根據(jù)噪聲和圖像的主要內(nèi)容屬于不同頻段，對圖像進行濾波（如高斯平滑）處理。該類算法實質(zhì)上是對圖像中任意像素進行各向（切向和法向）同性擴散［4］，切、法向等速擴散使得系統(tǒng)運行速度較快，對圖像平滑區(qū)去噪效果較好［5］，但法向擴散導(dǎo)致圖像結(jié)構(gòu)信息模糊［6］。

為了保護結(jié)構(gòu)信息，傳統(tǒng)全變分（TV）算法［7］從像素角度出發(fā)，分析了噪聲對像素的影響——加速像素變化。根據(jù)此性質(zhì)以像素變化的1e-范數(shù)建立了圖像去噪模型［8］。該算法是對任意像素僅沿切線方向擴散［9］，保護了邊緣信息，但平滑區(qū)的單向擴散使得系統(tǒng)運行效率較低［10］。

為了實現(xiàn)快速地圖像保邊去噪，本文分析了各向同性和異性擴散［11］的機理，將兩種擴散方式有機結(jié)合，繼承各自優(yōu)點彌補其不足，從像素變化角度出發(fā)，以梯度為變量設(shè)計了新的擴散函數(shù)。理論上分析了該擴散函數(shù)的局部擴散性能：圖像平滑區(qū)沿切、法線方向以不同的速度同時擴散，提高了該區(qū)域的平滑速度，彌補了傳統(tǒng)TV單向擴散［12］的不足；圖像邊緣和紋理區(qū)域沿切線方向擴散保護了邊緣。在此基礎(chǔ)上，本文結(jié)合傳統(tǒng)全變分圖像去噪模型，提出了改進的全變分圖像去噪：首先根據(jù)圖像邊緣形成的光學(xué)原理分析圖像梯度；其次根據(jù)梯度幅度將圖像分為平滑區(qū)和邊緣紋理區(qū)，建立相應(yīng)的偏微分方程；最后運用固定點迭代算法分析偏微分方程，設(shè)計了離散迭代函數(shù)。實驗結(jié)果表明，本文擴散函數(shù)在圖像邊緣紋理區(qū)各向異性擴散，繼承了傳統(tǒng)全變分的優(yōu)點，保護了邊緣；在圖像平滑區(qū)進行各向同性擴散，降低了系統(tǒng)運行時間，提高了圖像峰值信噪比（PSNR）和視覺效果。

1 擴散函數(shù)

從信號的頻率來看，圖像結(jié)構(gòu)信息和高斯噪聲表現(xiàn)為相同的特性——高頻信號，但結(jié)構(gòu)信息是分析理解圖像的基層信息。為了抑制噪聲提高分析理解的準確性，在圖像去噪過程中應(yīng)保護其結(jié)構(gòu)信息。假設(shè)原始圖像u（x， y）被零均值方差為σ2的高斯噪聲n（x， y）污染的圖像為：

噪聲分量 n（x，y）加劇圖像像素變化，使 I（x，y）的整體梯度比原始圖像大。因此，為了消除噪聲保護結(jié)構(gòu)信息，利用梯度幅度函數(shù)［13］建立的目標函數(shù)為：

運用變分法得到式（2）的歐拉-拉格朗日方程：

為了在去噪過程中保護結(jié)構(gòu)信息提高運行效率，擴散函數(shù)應(yīng)滿足兩個條件：1） 為了減少系統(tǒng)運行時間，平滑區(qū)應(yīng)各向同性擴散，擴散函數(shù)的變化速度接近于2） 為了保護結(jié)構(gòu)信息，邊緣和紋理區(qū)域應(yīng)各向異性處理。本文繼承傳統(tǒng)各向同性和異性擴散的優(yōu)點，以梯度幅度為變量設(shè)計了新的擴散函數(shù)：

圖1 不同擴散函數(shù)

該函數(shù)在切、法向擴散速度分別為：

切、法向擴散速度之比值k 為：

為了繼承各向同性和異性處理的優(yōu)點，本文利用梯度幅度將圖像分為平滑區(qū)和邊緣紋理區(qū)，如果梯度小于閾值T時，該像素位于平滑區(qū)應(yīng)各向同性處理，此時當像素梯度大于閾值T時，該像素應(yīng)各向異性處理，此時不同區(qū)域相應(yīng)的歐拉方程為：

2 離散迭代

當像素 u（i， j）位于平滑區(qū)域，本文假設(shè)其四鄰域也位于平滑區(qū)域，其鄰域的擴散權(quán)重為：

當像素 u（i， j）位于邊緣紋理區(qū)域，其四鄰域也位于該區(qū)域，對應(yīng)像素的擴散權(quán)重為：

不同區(qū)域的擴散系數(shù)如圖2所示。

為了在圖像去噪的過程中保護邊緣紋理信息，式（8）運用了噪聲圖像信息。當位于平滑區(qū)，鄰域像素梯度較小，對應(yīng)的擴散系數(shù)ω（p）遠大于權(quán)重λ，去噪后 u（i，j）主要由其鄰域像素加權(quán)決定；若像素（i，j）位于邊緣紋理處，對應(yīng)的梯度較大，相應(yīng)的擴散系數(shù)ω（p）趨于0，去噪后 u（i，j）主要由圖像 I（i，j）決定。

圖2 不同區(qū)域的擴散系數(shù)

為了得到式（2）的最優(yōu)解，該文將式（8）進行迭代運算，第n步迭代表達式為：

噪聲圖像 I（i，j）中包含了原始圖像邊緣紋理信息，每步迭代都運用了噪聲圖像的相關(guān)信息，使邊緣紋理信息不會隨著迭代次數(shù)增加而模糊。當滿足下式停止迭代，有：

3 實驗結(jié)果分析

為了快速地實現(xiàn)圖像保邊去噪，本文將各向同性和異性擴散有機結(jié)合，設(shè)計了新的圖像保邊去噪擴散函數(shù)。該擴散函數(shù)對圖像平滑區(qū)沿切法向同時擴散，減少了該區(qū)域的運行時間；對邊緣紋理區(qū)域法向擴散速度幾乎為0，保護了圖像的邊緣紋理。根據(jù)圖像邊緣紋理形成的光學(xué)機理，邊緣紋理區(qū)像素變化幅度比平滑區(qū)大，本文運用梯度固定閾值T法，將圖像分為平滑和邊緣紋理兩個區(qū)域，如果閾值T取值較大，圖像中像素變化較慢形成的弱邊緣紋理誤作為平滑區(qū)，去噪后圖像的弱邊緣紋理進行各向同性擴散，該信息保護較差；反之，將平滑區(qū)的噪聲誤判為邊緣紋理，對該區(qū)域的噪聲進行保護，殘余噪聲較大。

對含有加性高斯噪聲的peppers圖像運用不同閾值的部分去噪結(jié)果如圖3所示。當閾值0.2T=時，去噪后圖像的整體質(zhì)量PSNR為15.13 dB，由于閾值較小，圖像邊緣紋理保護較好，但是圖像平滑區(qū)的強噪聲（幅度較大）被作為邊緣紋理而得以保護，所以去噪后殘余噪聲較大；增大閾值使T=0.5，平滑區(qū)殘余噪聲較小，邊緣紋理得以較好的保護，去噪后圖像PSNR為21.3 dB；若繼續(xù)增大閾值，導(dǎo)致將紋理作為平滑區(qū)進行各向同性擴散，導(dǎo)致紋理損失而過度平滑，甚至出現(xiàn)虛假邊緣，去噪后PSNR下降，但保護圖像強邊緣。當T=1.1時，雖然圖像中各個目標輪廓清晰（強邊緣），但由于紋理損失較大，圖像PSNR為12.91 dB。為了保護圖像邊緣紋理信息，同時殘余噪聲較小，本文閾值T取為0.5。

圖3 不同閾值T的去噪結(jié)果

噪聲圖像包含了邊緣紋理信息，為了保護圖像邊緣紋理，在擴散處理中依權(quán)重λ考慮噪聲圖像的相關(guān)信息。權(quán)重λ大小在一定程度上決定了邊緣紋理的保護能力，當λ相對于鄰域權(quán)重較小時，去噪后殘余噪聲較小，但邊緣紋理模糊嚴重；反之，保邊性較好，殘余噪聲較大。對含有高斯噪聲的Lena圖像運用不同λ處理的結(jié)果如圖4所示，當1.0λ=時，PSNR為19.14 dB，圖像的強邊緣保護較好，但紋理區(qū)存在一定程度的模糊，平滑區(qū)存在過度平滑，偽邊緣嚴重；增大λ，平滑區(qū)的偽邊緣減少，圖像的整體質(zhì)量提高；當λ=2.0時，PSNR為27.67 dB；如繼續(xù)增大λ，圖像邊緣保護較好，但平滑區(qū)殘余噪聲增大使得PSNR降低，當λ=3.0時，PSNR為19.16 dB。為了保護圖像邊緣紋理信息，本文權(quán)重λ取為2.0。

圖4 不同λ的去噪結(jié)果

在處理器為Intel-Core i5CPU @3.40 GHz，RAM為4 GB，VC6.0的環(huán)境下，本文算法對圖3a的噪聲圖像迭代40次收斂，CPU的運行時間為1.52 s。其迭代次數(shù)與圖像PSNR之間的關(guān)系如圖5a所示，前20次迭代圖像的PSNR提高較快，將PSNR從12.03 dB提高到20.15 dB；隨迭代次數(shù)增加，圖像質(zhì)量改善速度下降，后20次僅僅提高了1.15 dB。傳統(tǒng)TV［14］（λ=3.0）需迭代445次收斂，CPU的運行時間為12.8 s，其迭代次數(shù)與PSNR之間的關(guān)系如圖5b所示。

為了驗證算法的有效性，對疊加不同高斯噪聲的圖像分別運用本文算法、高斯平滑（為了提高去噪效果，不同含噪圖像高斯核參數(shù)σ不同：Lena圖像σ=0.8；peppers圖像σ=1.4；baboon圖像σ=2.0）、和傳統(tǒng)TV［14］（λ=3.0）進行去噪處理，其部分結(jié)果如圖6所示，去噪后圖像的PSNR和SSIM［15］如表1所示。從PSNR的角度看，該算法比高斯平滑去噪效果好，比傳統(tǒng)全變分去噪質(zhì)量略高。從視覺效果上看，該算法的保邊性與傳統(tǒng)全變分相當，兩者SSIM差異不大。不同算法CPU的運行時間如表2所示。從表2可見，高斯平滑的系統(tǒng)運行時間最少，該算法對平滑區(qū)域進行各向同性擴散，其運行時間比傳統(tǒng)TV短。同時隨著噪聲圖像的質(zhì)量下降（PSNR變?。?，該算法的去噪時間變長。

圖5 迭代次數(shù)與PSNR之間的關(guān)系

圖6 不同算法比較

表1 不同去噪算法的質(zhì)量評價

表2 不同去噪算法的CPU時間 單位：s

4 結(jié) 束 語

為了實現(xiàn)快速地保邊去噪，本文將各向同性和異性擴散方式有機結(jié)合，從噪聲對像素變化的影響出發(fā)設(shè)計了新的擴散函數(shù)，結(jié)合傳統(tǒng)全變分圖像去噪模型，提出了改進的全變分圖像去噪。以梯度幅度將圖像分為平滑區(qū)和邊緣紋理區(qū)，對不同區(qū)域建立相應(yīng)的偏微分方程，利用固定點迭代算法設(shè)計了相應(yīng)的離散迭代函數(shù)。實驗結(jié)果表明，該擴散函數(shù)在邊緣紋理區(qū)各向同異擴散，保護了圖像邊緣；在平滑區(qū)各向同性擴散，降低了系統(tǒng)運行時間，彌補傳統(tǒng)全變分去噪算法單向擴散的不足。

［1］ GAO Lei-fu， FU xin. An adaptive TV model for image denoising［C］//Proceedings-International Conference on Natural Computation. Shenyang， China： ICNC， 2013.

［2］ ALAMRI 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.

［3］ LIU Chan-juan， QIAN Xu， LI Cai-xia. A texture image denoising model using the combination of tensor voting and total variation minimization［J］. Journal of Theoretical and Applied Information Technology， 2012， 46（1）： 294-302.

［4］ RUSSO F. New method for performance evaluation of grayscale image denoising filters［J］. IEEE Signal Processing Letters， 2010， 17（5）： 417-420.

［5］ LI Hong-juang， WU Ming-ni. Image noise reduction using wiener filtering with pseudo-inverse［J］. Measurement：Journal of the International Measurement Confederation，2010， 43（10）： 1649-1655.

［6］ HAJIABOL M R. A self-governing fourth-order nonlinear diffusion filter for image noise removal［J］. IPSJ Transactions on Computer Vision and Applications， 2011：92（2）： 177-191.

［7］ PAN Zhen-kuan， GAO Lei， WEI Wei-bo， et al. The dual Bregman algorithms of generalized TV models for image denoising［C］//Proceedings of the 2012 International Conference on Image Processing， Computer Vision， and Pattern Recognition. Las Vegas， USA： ［s.n.］， 2012.

［8］ CHO S I， KIM Y H. Edge connectivity-based image denoising for digital TV systems［C］//IEEE International Conference on Consumer Electronics. Las Vegas， USA：［s.n.］， 2013.

［9］ LI Dong-ming， ZHANG li-juan. Study on image denoising method based on an adaptive total variation model［C］//2011 International Conference on Transportation， Mechanical，and Electrical Engineering. Changchun， China： ICNC， 2011.

［10］ 洪丹楓， 魏偉波， 潘振寬， 等. 基于curvelet變換和非局部TV模型的圖像去噪［J］. 青島大學(xué)學(xué)報（自然科學(xué)版），2014， 27（1）： 59-65. HONG Dan-feng， WEI Wei-bo， PAN Zhen-kuan， et al. Image denosing based on curvelet transform and non-local TV model［J］. Journal of Qingdao University （Natural Science Edition）， 2014， 27（1）： 59-65.

［11］ 張春曉， 趙剡， 連遠鋒. 各向異性擴散圖像去噪現(xiàn)狀研究［J］. 電光與控制， 2012， 19（5）： 51-55. ZHANG Chun-xiao， ZHAO Yan， LIAN Yuan-feng. Research on anisotropic diffusion denoising development［J］. Electronics optics&Control， 2012， 19（5）： 51-55.

［12］ 王旭東， 馮象初， 霍雷剛. 去除乘性噪聲的重加權(quán)各向異性全變差模型［J］. 自動化學(xué)報， 2012， 38（3）： 444-451. WANG Xu-dong， FENG Xiang-chu， HUO Lei-gang. Iteratively reweighted anisotropic-TV based multiplicative noise removal model［J］. Acta Automatica Sinica， 2012，38（3）： 444-451.

［13］ 端金鳴. 潘振寬， 臺雪成. 彩色紋理圖像恢復(fù)的非局部TV模型［J］. 中國圖象圖形學(xué)報， 2013， 18（7）： 753-760. DUAN Jin-ming， PAN Zhen-kuan， TAI Xue-cheng. Non-local TV models for restoration of color texture images［J］. Journal of Image and Graphics， 2013， 18（7）：753-760.

［14］ AUBERT G， KORNPROBST P. Mathematical problem in image processing： Partial differential equations and the calculus of variations［M］. New York， USA： Springer Science Business Media， 2009.

［15］ ZHOU W， BOVIK A C. Modern image quality assessment［M］. New York， USA： Morgan & Claypool Publishers，2006.

編 輯 黃 莘

Image Denoising on Improved Total Variation

HE Kun1， ZHENG Xiu-qing2， JU Sheng-gen1， and ZHANG Yong-lai2
（1. School of Computer Science， Sichuan University Chengdu 610065；2. College of Information and Technology， Sichuan Normal University Chengdu 618300）

By analyzing the mechanism of isotropic and anisotropic diffusion in image denoising， a new diffusion function based on the effect of noise on the image pixel variation is designed in this paper. The diffusion performance of this function is isotropic diffusion in the smoothing sub-region and anisotropic diffusion on the edge. Then， an improved total variation denoising model is proposed based on the traditional total variation （TV），and the corresponding discrete iterative function is designed by using the fixed point iteration algorithm. The algorithm uses isotropic diffusion on the smooth region， inheriting the advantages of the isotropic diffusion，reducing the running time of the traditional TV， and protecting the edge by using anisotropic diffusion. Experimental results show that this algorithm can improve image peak signal to noise ratio （PSNR） and visual effects.

anisotropic； diffusion function； image denoising； isotropic； total variation

TP394.1

A

10.3969/j.issn.1001-0548.2016.02.026

2014 - 07 - 16；

2016 - 01 - 19

四川省科技支撐計劃（2013SZ0157）

何坤（1972 - ），男，博士，副教授，主要從事圖像處理、模式識別等方面的研究.