汪美玲 周先春石蘭芳

(1.南京信息工程大學(xué)電子與信息工程學(xué)院，南京，210044；2.南京信息工程大學(xué)，江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心,南京,210044；3.南京信息工程大學(xué),江蘇省氣象探測與信息處理重點實驗室,南京,210044；4.南京信息工程大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院，南京，210044)

雙向增強擴散濾波的圖像去噪模型*

汪美玲1,2,3周先春1,2,3石蘭芳4

(1.南京信息工程大學(xué)電子與信息工程學(xué)院，南京，210044；2.南京信息工程大學(xué)，江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心,南京,210044；3.南京信息工程大學(xué),江蘇省氣象探測與信息處理重點實驗室,南京,210044；4.南京信息工程大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院，南京，210044)

提出一種雙向增強擴散濾波的圖像去噪模型。簡化擴散方程建立雙向擴散系數(shù)，使模型在擴散過程中能夠?qū)崿F(xiàn)平滑與銳化的雙向過程，為加強平滑和銳化強度，用小波變換增強圖像，使整體圖像輪廓得到增強和局部圖像紋理特征得到弱化。然后，對閾值進行了自適應(yīng)設(shè)計和改進，使其根據(jù)圖像的最大灰度值和迭代次數(shù)自動控制閾值，進一步保留圖像邊緣和細節(jié)特征。實驗仿真和可行性的驗證結(jié)果表明，新模型去噪效果較理想，不但能抑制噪聲，而且能保護細節(jié)信息，峰值信噪比得到了有效的提高，性能更優(yōu)越。

雙向擴散系數(shù)；自適應(yīng)閾值；圖像去噪

引 言

當(dāng)今，許多科學(xué)領(lǐng)域都包含著數(shù)字圖像處理技術(shù)，而在偏微分方程(Partial differential equation,PDE)圖像去噪算法中，各向異性擴散已成為當(dāng)前研究熱點[1-4]。Perona和Malik[5]在1990年首次令傳導(dǎo)系數(shù)依賴圖像梯度的二階PDE，建立PM模型。Rudin等[6]在1992年從能量泛函角度建立了能更好地保持圖像邊緣紋理等細節(jié)特征的全變分(Total variation,TV)模型，是PDE在圖像去噪中的重要體現(xiàn)。Guy Gilboa 等[7]提出了復(fù)模型，是從實數(shù)域延伸到復(fù)數(shù)域的PM模型的改進，該模型能夠同時實現(xiàn)前后向擴散，擴大了去除噪聲的范圍。自從提出PM模型后，近些年來，PDE的圖像去噪方法取得了很大的進展[8-10]。最近，出現(xiàn)了許多較新的圖像去噪算法，如Buades等提出了非局部平均(Non-local means,NLM)濾波器，該濾波器能夠保護圖像的紋理特征，但計算復(fù)雜、運行效率慢[11]。文獻 [12] 提出了塊匹配三維協(xié)同(Block matching 3D cooperative,BM3D)方法，該方法的峰值信噪比較高，視覺效果很好，但有復(fù)雜的時間處理過程。文獻[13]提出了基于概率修復(fù)(Probablistic patch-based, PPB )濾波器，該濾波器對合成孔徑雷達(Synthetic aperture radar, SAR)圖像濾波效果較好，但處理效率較低，有過長的濾波過程。針對該問題，又提出了非迭代的PPB濾波方法，該方法雖然解決了處理效率的問題，但該濾波器對圖像紋理特征的保持仍有欠缺。本文在經(jīng)典的擴散模型基礎(chǔ)上，建立局部坐標系，對擴散模型進行簡化與分析，建立擴散濾波過程，使圖像在擴散過程中實現(xiàn)平滑與銳化的雙向過程。為了更加突出圖像的整體結(jié)構(gòu)，增強邊緣的銳化作用，作一些特殊處理，用小波對圖像進行分解，在頻域里增強圖像，對圖像進行系數(shù)的分解，從而達到圖像整體結(jié)構(gòu)的突出和圖像細節(jié)特征的弱化，然后，對閾值進行了自適應(yīng)設(shè)計，使其根據(jù)圖像的最大灰度值和迭代次數(shù)自動控制閾值。

1 PM和TV模型

1.1 PM模型

Perona和Malik首次提出了PM模型，其能量泛函為

(1)

(2)

(3)

(4)

式中：k為圖像的梯度閾值。

1.2TV模型

TV模型是另一種經(jīng)典的去噪模型，是將其構(gòu)造的能量泛函轉(zhuǎn)化為偏微分方程來求解，其能量泛函公式為

(5)

式中:I0為原始圖像，I為變化中的灰度圖像，I為圖像像素的梯度，λ為拉格朗日算子，其Euler-lagrange方程為

(6)

用梯度下降流解式(6)，則可得TV模型

(7)

2 改進的模型

2.1 經(jīng)典擴散模型的簡化分析

圖1 圖像中某點局部坐標 Fig．1 Local coordinate of a point of the image

(8)

(9)

(10)

則式(2)最終表達為

(11)

(12)

2.2 新模型建立

根據(jù)以上分析對式(11)做進一步簡化

(13)

設(shè)

(14)

(15)

(16)

由式(16)可知，沿梯度方向逆擴散，銳化圖像邊緣，改善了邊緣模糊的現(xiàn)象，能很好地保護圖像的邊緣。

(17)

圖2 g(·)函數(shù)和 Iηη的擴散系數(shù)Fig．2 Function g(·) and diffusion coefficient curve of Iηη

由式(17)可知，在圖像的平滑區(qū)域，該模型轉(zhuǎn)化為了各向同性去噪方程，即在圖像的邊緣和梯度方向能夠達到同時平滑的效果，去噪能力較強。觀察梯度方向上Iηη前擴散系數(shù)的曲線變化圖2(b)所示。由圖2(b)可以看出，隨著梯度增大，當(dāng)梯度達到某一值時，擴散系數(shù)出現(xiàn)了負值，符合以上的分析過程。本文的整個擴散過程，實質(zhì)上是通過圖像梯度與閾值的大小進行比較實現(xiàn)的，通常在擴散過程中，會預(yù)先設(shè)定一個梯度閾值k，k一般是個常量，即不論像素點所在區(qū)域與相鄰區(qū)域的灰度值存在怎樣的關(guān)系，都采用這個預(yù)先設(shè)定好的常量梯度閾值，從而不能保持圖像某些區(qū)域的細節(jié)特征。隨著擴散時間的不斷增加，該部分像素點的梯度閾值逐漸減小，為了使邊緣和噪聲得以區(qū)分，當(dāng)擴散時間增加時，梯度閾值k則相應(yīng)減小，從而就能在小范圍區(qū)域內(nèi)判斷邊緣位置，能更好地保持圖像特征。

綜上分析，將梯度閾值k設(shè)計為隨擴散時間和擴散次數(shù)變化的函數(shù)，即有

(18)

在濾波過程中為了不讓圖像失真，給圖像添加一項保真項[14]，表示為

(19)

設(shè)計的雙向增強擴散模型為

(20)

圖3 實驗原始圖像與增強圖像 Fig．3 Original and enhanced images of experiment

為了解雙向增強擴散模型，采用中心差分法[15]，則雙向增強擴散模型的離散化形式為

(21)

式中:n=0,1,2,……,表示時間步長。其中

(22)

(23)

(24)

3 實驗結(jié)果與分析

(25)

(26)

式中：圖像的分辨率為W×H，I和I0分別為原始圖像和去噪后的圖像，MSE越小越好，PSNR越大越好。

首先對Nuist圖像用各算法和本文新模型進行平滑。其中，Δt=0.2，n=30，PM模型中的擴散系數(shù)取文章中的式(3)，10為取定的閾值，TV模型中，λ=0.1，新模型中，λ*=5。平滑結(jié)果、評價指標和局部放大圖分別如圖4、表1和圖5所示。為了更好地顯示濾波前后圖像的細節(jié)信息情況，采用Canny算子[16]檢測各模型濾波結(jié)果的圖像，實驗結(jié)果如圖6所示。

分別從圖像整體與局部放大的可視效果來看，如圖4(b),(c)和圖5(b),(c)，PM模型處理后的圖像中仍含有大量噪聲沒有被去除，且產(chǎn)生了“塊狀”效應(yīng)；TV模型去噪后的圖像中出現(xiàn)了一些孔洞，這是由

圖4 各模型去噪后圖像Fig．4 Denoised image of each model

類別加噪圖像PM模型[4]TV模型[5]NLM算法[10]BM3D算法[11]Non-iterativePPB算法[12]IterativePPB算法[12]新模型MSE385.2109222.970377.342571.35527.6887×10-466.138765.31110.5474PSNR22.273824.648329.246629.596531.141529.932529.985950.7474

圖5 各模型去噪后局部放大圖像Fig．5 Local amplification of denoised images of all models

圖6 各模型去噪后邊緣提取圖像Fig．6 Edge detection of denoised images of all models

于TV模型不滿足形態(tài)學(xué)的原則，從而使圖像遭到破壞，但該模型有一定的去噪效果，圖像的細節(jié)特征保持較好，如圖6(c)所示。由4(d)和表1可以看出，盡管NLM算法能較好地去除噪聲，但從圖5(d)、圖6(d)可以看出該算法不能很好地保護原圖像的結(jié)構(gòu)信息。圖4(e)、圖5(e) 和表1顯示了BM3D方法有較好的視覺效果和較高的PSNR值，但從圖6(e)可知，此方法不能有效地保護圖像細節(jié)特征。在表1中，Iterative PPB方法去噪性能較好，但在運行過程中花費的時間很長，從而不能有效地應(yīng)用到實際中，而 Non-iterative PPB方法解決了這個問題，但是該方法不能很好地保護原圖像的結(jié)構(gòu)信息，如圖6(f)所示。由圖4(h)、圖5(h)可以看到雙向增強擴散模型的可視性最好，主要是由于該模型建立平滑與銳化的雙向擴散過程，比較好地考慮圖像的局部特征，兼顧了去噪與紋理特征的保護。由表1可知，雙向增強擴散模型各項性能最好，與仿真圖的效果相同。

用不同強度噪聲反映各模型對噪聲的去除效果和圖像特征保護狀況。現(xiàn)對圖像添加不同方差的噪聲，實驗結(jié)果如圖7所示，由圖7可以看出新模型在所有模型中有最高的峰值信噪比和最低的均方差。

圖7 不同噪聲方差時MSE和PSNR的變化Fig．7 PSNR and MSE with different variance noise for each model

4 結(jié)束語

本文在擴散方程的基礎(chǔ)上，構(gòu)建雙向擴散系數(shù)，建立雙向增強擴散濾波模型。該模型在擴散過程中能夠?qū)崿F(xiàn)平滑與銳化的雙向過程，為了更加突出圖像的整體結(jié)構(gòu)，增強對邊緣銳化的作用，做一些特殊處理。首先用小波分解圖像，在頻域里增強圖像，處理分解系數(shù)，突出圖像的輪廓和整體結(jié)構(gòu)，弱化細節(jié)，然后對閾值進行了自適應(yīng)設(shè)計和改進，使其根據(jù)圖像的最大灰度值和迭代次數(shù)自動控制閾值，進一步保留圖像邊緣和細節(jié)特征。仿真實驗結(jié)果表明，本文新模型能夠兼顧圖像噪聲的去除和細節(jié)特征的保護，大幅提高了PSNR，較經(jīng)典模型去噪性能更優(yōu)越。

[1] Niang O,Thioune A, Gueirea M C, et al. Partial differential equation-based approach for empirical mode decomposition:Application on image analysis[J]. IEEE Transactions on Image Processing,2012,21(9):3991-3999.

[2] 陳婷婷, 沈民奮, 楊金耀. 一種復(fù)合各向異性擴散的圖像去噪算法[J].數(shù)據(jù)采集與處理, 2014, 29 (5): 757-763.

Chen Tingting, Shen Minfen, Yang Jinyao. Algorithm based on compound anisotropic diffusion for image denoising[J]. Journal of Data Acquisition and Processing, 2014, 29 (5): 757-763.

[3] Liu Peng, Huang Fang, Li Guoqing. Remote-sensing image de-noising using partial differential equations and auxiliary images as priors[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9 (3): 358-361.

[4] 周先春, 汪美玲, 石蘭芳. 閾值尋優(yōu)的高保真各向異性濾波模型[J]. 計算機輔助設(shè)計與圖形學(xué)學(xué)報, 2016, 28 (9): 1550-1559.

Zhou Xianchun, Wang Meiling, Shi Lanfang. Anisotropic filtering model of high-fidelity based on threshold optimization[J]. Journal of Computer-Aided Design & Computer Graphics, 2016, 28 (9): 1550-1559.

[5] Perona P, Malik J. Scale-space and edge detection using anisotropic diffusion[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1990, 12 (7): 629-639.

[6] Rudin L, Oshers, Fatemi E. Nonlinear total variation based noise removal algorithms [J]. Physic D, 1992, 60 (1-4): 259-268.

[7] Gilboa G. Nonlinear scale space with spatially varying stopping time[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2008, 30 (12): 2175-2186.

[8] Chumchob N. Vectorial total variation-based regularization for variational image registration[J]. IEEE Transactions on Image Processing, 2013, 22 (11): 4551-4558.

[9] 陳永彬，宋海英. 基于圖像亮度局部約束的總變分圖像去噪法[J]. 數(shù)據(jù)采集與處理, 2011, 26 (5): 564-568.

Chen Yongbin, Song Haiying. Total variation image denoising based on local intensity bound constraints[J]. Journal of Data Acquisition and Processing, 2011, 26 (5): 564-568.

[10] Zhong Junmei, Sun Huifang. Wavelet-based multiscale anisotropic diffusion with adaptive statistical analysis for image restoration[J]. IEEE Transactions on Circuits and System—I: Regular Papers, 2008, 55 (9):2716-2724.

[11] Buades A. A non-local algorithm for image denoising[C]∥ IEEE Computer Society Conference on Computer Vision and Pattern Recognition. San Diego:[s.n.], 2005, 60-65.

[12] Dabov K, Foi A, Katkovnik V, et al. Image denoising by sparse 3-D transform-domain collaborative filtering[J]. IEEE Transactions on Image Processing, 2007, 16 (8): 2080-2095.

[13] Deledalle C A, Denis L, Tupin F. Iterative weighted maximum likelihood denoising with probabilistic patch-based weights[J]. IEEE Transactions on Image Processing, 2009, 18 (12): 2661-2672.

[14] 朱立新，王平安，夏德深. 非線性擴散圖像去噪中的耦合自適應(yīng)保真項研究[J]. 計算機輔助設(shè)計與圖形學(xué)學(xué)報，2006，18 (10)：1519-1524.

Zhu Lixin, Wang Pingan, Xia Deshen. A study on the coupling adaptive fidelity term in anisotropic diffusion based image denoising[J]. Journal of Computered-Aided Design & Computer Graphics, 2006, 18 (10): 1519-1524.

[15] Ghita O, llea D E, Whelan P F. Adaptive noise removal approach for restoration of digital images corrupted by multimodal noise[J]. IET Image Process let, 2012, 6(8): 1148-1160.

[16] Xu Qian, Varadarajan S, Chakrabarti C, et al. A distributed canny edge detector: Algorithm and FPGA implementation[J]. IEEE Transactions on Image Processing, 2014, 23 (7): 2944-2959.

Image De-noising Model based on Bidirectional Enhanced Diffusion Filter

Wang Meiling1,2,3, Zhou Xianchun1,2,3, Shi Lanfang4

(1. School of Electronic and Information Engineering, Nanjing University of Information Science & Technology, Nanjing, 210044, China; 2. Jiangsu Collaborative Innovation Center of Atmospheric Environment and Equipment Technology, Nanjing University of Information Science and Technology, Nanjing, 210044, China; 3. Jiangsu Key Laboratory of Meteorological Observation and Information Processing, Nanjing University of Information Science and Technology, Nanjing, 210044, China; 4. School of Mathematics and Statistics, Nanjing University of Information Science & Technology, Nanjing, 210044, China)

A bidirectional enhanced diffusion filter image de-noising model is presented. The diffusion equation is firstly simplified and analyzed to establish bidirectional diffusion coefficient. Hence, the two-way process of smoothing and sharpening can be achieved by the model in the diffusion process, To further enhance the strength of the smoothing and sharpening, image enhancement is used to enhance the overall outline of the image using wavelet transform, thus weakening texture detail of the image. Then, the threshold will be designed and improved, and it will be automatically controlled by maximum image gray value and iterative times, which can retain the image edge and detail features. The proposed model is be simulated. The experimental result shows that the new model is ideal, and it can improve the performance of de-noising and the protection of edge. The texture detail information is satisfactory. The peak signal to noise ratio is promoted drastically. Therefore,the performance is better than classical algorithms.

bidirectional diffusion coefficient; adaptive threshold; image de-noising

國家自然科學(xué)基金(11202106, 61201444)資助項目；教育部高等學(xué)校博士學(xué)科點專項科研基金(20123228120005)資助項目；江蘇省“信息與通信工程”優(yōu)勢學(xué)科建設(shè)項目;江蘇省自然科學(xué)基金(BK20131005)資助項目；江蘇省青藍工程和江蘇省高校自然科學(xué)研究項目(13KJB170016)資助項目。

2014-12-30；

2015-10-12

TP391

A

汪美玲(1989-)，女，碩士研究生，研究方向:圖像處理、模式識別，E-mail:Mely nuaa@163.com。

周先春(1974-)，男，博士，副教授，研究方向：信號與信息處理。

石蘭芳(1976-)，女，博士，副教授，研究方向：非線性分析。