孫偉峰　戴永壽（中國(guó)石油大學(xué)（華東）信息與控制工程學(xué)院青島266580）

?

采用多級(jí)殘差濾波的非局部均值圖像去噪方法

孫偉峰*戴永壽
（中國(guó)石油大學(xué)（華東）信息與控制工程學(xué)院青島266580）

為充分利用殘差中的圖像信息以提升非局部均值算法的去噪性能，該文提出一種多級(jí)殘差圖像濾波新方法。首先對(duì)含噪圖像進(jìn)行非局部均值濾波得到初始的去噪圖像和權(quán)值分布矩陣，然后對(duì)殘差圖像進(jìn)行固定權(quán)值非局部均值濾波來(lái)提取圖像結(jié)構(gòu)信息，將提取的信息經(jīng)高斯平滑抑噪后作為補(bǔ)償圖像，與去噪圖像相加得到增強(qiáng)的恢復(fù)圖像。針對(duì)上述方法提出一種多級(jí)濾波的實(shí)現(xiàn)方案，從理論上推導(dǎo)證明了該方法的原理及可行性，并提出一種無(wú)需參考圖像的迭代停止準(zhǔn)則來(lái)自適應(yīng)地優(yōu)選濾波級(jí)數(shù)。實(shí)驗(yàn)結(jié)果表明，提出的迭代停止準(zhǔn)則能夠達(dá)到與峰值信噪比一致的優(yōu)選結(jié)果；與經(jīng)典的非局部均值算法相比，在計(jì)算效率相當(dāng)?shù)那闆r下，所提方法能夠顯著地提升其去噪性能，峰值信噪比平均可以提高1.2 dB，且具有更好的細(xì)節(jié)保持能力。

圖像去噪；非局部均值；殘差濾波；停止準(zhǔn)則

1　引言

數(shù)字圖像在獲取、傳輸?shù)冗^(guò)程中不可避免地會(huì)受到噪聲的影響，噪聲的引入會(huì)降低圖像質(zhì)量，進(jìn)而影響圖像分割、配準(zhǔn)、識(shí)別等圖像處理任務(wù)的性能。因此，作為一種重要的預(yù)處理手段，圖像去噪一直是圖像處理領(lǐng)域中一個(gè)基礎(chǔ)且重要的研究課題，其難點(diǎn)在于如何在噪聲去除與圖像特征保留之間達(dá)到平衡。

非局部均值（Non-Local Means，NLM）算法［1］是空域?yàn)V波方法的新發(fā)展，被廣泛應(yīng)用于圖像與視頻降噪［24］-。以此方法為基礎(chǔ)，發(fā)展了一系列基于圖像分塊（patch-based）表示的圖像去噪方法［510］-，取得了先進(jìn)的去噪結(jié)果。非局部均值算法利用圖像分塊間的相似性質(zhì)去除噪聲，具有良好的邊緣保持能力；然而，由于算法參數(shù)選擇不合適、相似性權(quán)值計(jì)算不準(zhǔn)確等不足，在濾除噪聲的同時(shí)，也會(huì)丟失一些圖像細(xì)節(jié)成分，這部分信息將存在于殘差圖像中，目前，大多數(shù)圖像去噪方法都致力于設(shè)計(jì)算法使得去噪圖像更逼近原始圖像，而很少考慮如何令殘差圖像逼近純?cè)肼暎瑳](méi)有很好地利用殘差圖像提供的信息。

近年來(lái)，越來(lái)越多的學(xué)者認(rèn)識(shí)到殘差圖像中蘊(yùn)含信息的重要性，提出了許多殘差圖像處理方法來(lái)提高去噪算法的性能。例如，文獻(xiàn)［11］提出采用自適應(yīng)維納濾波從殘差圖像中提取圖像的結(jié)構(gòu)信息。文獻(xiàn)［12］設(shè)計(jì)了一種多邊濾波器對(duì)殘差圖像進(jìn)行濾波；文獻(xiàn)［13］將該方法取代BM 3D算法中的維納濾波環(huán)節(jié)，更好地保持了圖像的結(jié)構(gòu)特征。文獻(xiàn)［14］利用殘差圖像的信息修正NLM算法中的權(quán)值表達(dá)，提高了權(quán)值計(jì)算的準(zhǔn)確性。文獻(xiàn)［15］采用小波閾值的方法對(duì)殘差圖像進(jìn)行處理，提升了NLM算法的去噪性能。文獻(xiàn)［16］提出一種殘差圖像后處理方式來(lái)增強(qiáng)K-SVD算法的去噪能力，并發(fā)展了一種通用的增強(qiáng)圖像去噪算法性能的遞歸算法［17］?？梢?jiàn)，殘差中圖像信息的利用為提升去噪算法的性能提供了一種新的途徑。

為進(jìn)一步提升NLM算法的去噪性能，本文提出一種多級(jí)殘差圖像濾波新方法。該方法以固定權(quán)值非局部均值濾波與高斯濾波作為基本的處理模塊，通過(guò)對(duì)殘差中圖像信息的逐級(jí)提取與補(bǔ)償，提高去噪圖像質(zhì)量；基于殘差圖像提出了一種無(wú)需參考圖像的圖像質(zhì)量評(píng)價(jià)準(zhǔn)則進(jìn)行最優(yōu)濾波級(jí)數(shù)選擇，兼顧了噪聲去除與圖像細(xì)節(jié)保留。從理論推導(dǎo)分析與實(shí)驗(yàn)結(jié)果驗(yàn)證兩個(gè)方面證明了提出方法的有效性。

2　非局部均值去噪算法與殘差圖像

2.1非局部均值去噪算法

給定一幅含噪圖像y，對(duì)于其中的任何一個(gè)像素i，NLM算法利用搜索鄰域內(nèi)所有像素值的加權(quán)平均作為該點(diǎn)的估計(jì)值（）x i?，即：

其中，Ω表示以像素i為中心的搜索鄰域；（，）i jω為像素i與j之間的相似性權(quán)值，定義為

式（2）中，iN表示以像素i為中心的固定大小的圖像塊，稱為相似鄰域，為歸一化常數(shù)，h為濾波參數(shù)。文中將由所有（，）i jω構(gòu)成的矩陣W稱為NLM算法的權(quán)值分布矩陣，將采用已知權(quán)值分布矩陣的NLM算法稱為固定權(quán)值非局部均值濾波方法。

2.2殘差圖像

殘差圖像也被稱為“方法噪聲”［1］，定義為含噪圖像與去噪圖像之間的差值。理想情況下，去噪過(guò)程中應(yīng)該僅有噪聲被濾除，因而對(duì)于一個(gè)好的去噪方法，其殘差圖像應(yīng)該類似白噪聲。而實(shí)際上，即使對(duì)于目前最先進(jìn)的去噪方法，其殘差圖像中仍然包含圖像細(xì)節(jié)信息［15］。

在Lena圖像中加入標(biāo)準(zhǔn)差為25的高斯白噪聲，利用非局部均值算法進(jìn)行去噪，得到的去噪圖像及相應(yīng)的殘差圖像如圖1所示。可見(jiàn)，殘差圖像中包含明顯的圖像結(jié)構(gòu)信息，若能夠有效提取這些信息并將其反饋至去噪圖像中，必然能夠增強(qiáng)去噪圖像中的細(xì)節(jié)特征?；谶@個(gè)思路，為了充分提取殘差中的圖像結(jié)構(gòu)信息，本文提出了一種多級(jí)殘差圖像濾波方法，用以增強(qiáng)NLM算法的去噪性能。

3　采用多級(jí)殘差濾波的非局部均值圖像去噪方法

非局部均值算法的權(quán)值分布矩陣反映了像素間的結(jié)構(gòu)相似性，觀察圖1（b）可見(jiàn)，殘差圖像與含噪圖像中包含相似的圖像結(jié)構(gòu)，可以采用相同的權(quán)值分布矩陣對(duì)殘差圖像進(jìn)行非局部均值濾波來(lái)提取圖像結(jié)構(gòu)信息。由于濾波過(guò)程在強(qiáng)噪聲背景下進(jìn)行，提取的圖像信息中必然包含噪聲成分，為了削弱噪聲的影響，在反饋至去噪圖像之前，本文采用高斯濾波進(jìn)行降噪處理。這樣，就構(gòu)成了本文提出的殘差圖像濾波方法的兩個(gè)核心模塊：固定權(quán)值非局部均值濾波與高斯平滑。

圖1　去噪圖像與殘差圖像示例

為了更加充分地提取殘差中的圖像信息，對(duì)上述方法提出一種多級(jí)濾波的實(shí)現(xiàn)方案，使圖像信息提取與補(bǔ)償環(huán)節(jié)交織進(jìn)行；此外，利用去噪圖像、殘差圖像與含噪圖像之間的相關(guān)性提出了一種迭代停止準(zhǔn)則，能夠?qū)V波級(jí)數(shù)進(jìn)行自適應(yīng)優(yōu)選。由此，發(fā)展了一種新的多級(jí)殘差圖像濾波方法。

3.1算法原理

W為N N×的正定矩陣，對(duì)其進(jìn)行特征值分解得=W V d。其中，其對(duì)角線上包含由小到大排列的特征值為正交矩陣，其每一列為對(duì)應(yīng)特征值λ1，λ2，…，λN的特征向量。利用矩陣V的列向量構(gòu)成的空間對(duì)圖像x進(jìn)行分解得=x V b。

基于上述表示，對(duì)含噪圖像進(jìn)行非局部均值濾波得到初始去噪圖像對(duì)于本文提出的基于固定權(quán)值非局部均值濾波的多級(jí)殘差濾波方法，第k次迭代得到的去噪圖像kx?可以表示為

將式（4）展開(kāi)表示成含噪圖像y的函數(shù)的形式為

式（5）中，I表示N N×的單位矩陣。為了分析提出方法的有效性，對(duì)迭代過(guò)程中所得降噪圖像的均方誤差值隨迭代次數(shù)的變化規(guī)律進(jìn)行討論。第k次迭代時(shí)，所得估計(jì)圖像的偏差為

其中，E（）?表示求數(shù)學(xué)期望。偏差模值的平方表示為

其中，tr（）?與Cov（）?分別表示求矩陣的跡與協(xié)方差。由此，第k次迭代后所得去噪圖像的均方誤差可以表示為

由式（7）～式（9）可知，隨著迭代次數(shù)k的增加，所得估計(jì)圖像的偏差將逐漸減小，而方差將逐漸增大，當(dāng)k→∞時(shí)，所得去噪圖像的方差將趨于，即將收斂到含噪圖像y。

令式（10）等于零，得到對(duì)第p個(gè)像素點(diǎn)估計(jì)時(shí)的最優(yōu)迭代次數(shù)

上述分析證明，在一定的迭代次數(shù)范圍內(nèi)，對(duì)殘差圖像進(jìn)行多級(jí)固定權(quán)值非局部均值濾波，利用提取的信息對(duì)去噪圖像進(jìn)行補(bǔ)償，能夠在均方誤差意義上提高去噪圖像的質(zhì)量，從而證明了本文方法在原理上的可行性。然而，在實(shí)際的去噪過(guò)程中，很難對(duì)每一個(gè)像素點(diǎn)處的信號(hào)能量進(jìn)行準(zhǔn)確估計(jì)，因而無(wú)法通過(guò)直接計(jì)算得到。注意到，盡管對(duì)于不同的像素點(diǎn)，也不盡相同，很難找到一個(gè)統(tǒng)一的迭代次數(shù)，使得對(duì)每個(gè)像素的估計(jì)都能達(dá)到均方誤差意義上的最優(yōu)，但是，必然存在一個(gè)迭代次數(shù)，能夠從整體上優(yōu)化去噪效果，本文在3.3節(jié)提出一種停止準(zhǔn)則來(lái)確定該“最優(yōu)”的迭代次數(shù)。

3.2算法描述

采用多級(jí)殘差濾波的非局部均值圖像去噪算法的流程圖如圖2所示。

圖2　算法流程圖

具體步驟如下：

（1）令k=0，對(duì)含噪圖像y利用NLM算法進(jìn)行去噪，得到初始去噪圖像及權(quán)值分布矩陣W；

（3）令k=k+1，采用權(quán)值分布矩陣W對(duì)殘差圖像進(jìn)行固定權(quán)值NLM濾波，再經(jīng)過(guò)高斯濾波抑制噪聲后，得到補(bǔ)償圖像Ck；

（5）重復(fù)步驟（2）～步驟（4），直至達(dá)到最優(yōu)濾波級(jí)數(shù)k*；

（6）將得到的?xk*作為最終的去噪結(jié)果。

3.3停止準(zhǔn)則

多級(jí)殘差濾波過(guò)程中，每一級(jí)濾波得到的去噪圖像與殘差圖像都在變化，算法將最終收斂到含噪圖像。因此，停止準(zhǔn)則的選擇尤為重要，決定著算法的性能與計(jì)算復(fù)雜度。本文基于去噪圖像、殘差圖像與含噪圖像之間的相關(guān)性［18］，提出了一種無(wú)參考的圖像質(zhì)量評(píng)價(jià)方法，用來(lái)選擇最優(yōu)的濾波級(jí)數(shù)。根據(jù)皮爾遜相關(guān)系數(shù)，按式（12）計(jì)算兩幅大小為

基于上述考慮，本文提出采用內(nèi)積的方式實(shí)現(xiàn)兩者之間的線性相關(guān)，定義為

其中，me表示e的第m個(gè)元素。

4　實(shí)驗(yàn)結(jié)果與分析

利用大小為512×512的Lena，Barbara，Boat以及256×256的House，Peppers5幅標(biāo)準(zhǔn)測(cè)試圖像及一幅實(shí)際的含噪圖像開(kāi)展去噪實(shí)驗(yàn)，對(duì)本文方法的有效性進(jìn)行驗(yàn)證。

4.1停止準(zhǔn)則實(shí)驗(yàn)驗(yàn)證

以峰值信噪比（Peak Signal to Noise Ratio，PSNR）作為評(píng)價(jià)標(biāo)準(zhǔn)，利用提出的去噪方法進(jìn)行固定迭代次數(shù)去噪實(shí)驗(yàn)，迭代次數(shù)取10，對(duì)3.3節(jié)中提出的停止準(zhǔn)則進(jìn)行驗(yàn)證、評(píng)價(jià)。

首先，在Lena圖像中加入標(biāo)準(zhǔn)差為10～50的高斯白噪聲，采用PSNR準(zhǔn)則與本文提出的e值準(zhǔn)則判定最優(yōu)的迭代次數(shù)，結(jié)果如圖3所示。

圖3　最優(yōu)迭代次數(shù)比較

由圖3中的結(jié)果可見(jiàn)，由PSNR準(zhǔn)則與e值準(zhǔn)則確定的最優(yōu)迭代次數(shù)具有高度的一致性，一般經(jīng)過(guò)2～3次迭代即可取得最優(yōu)的結(jié)果。噪聲標(biāo)準(zhǔn)差為40時(shí)，由PSNR準(zhǔn)則與e值準(zhǔn)則確定的最優(yōu)迭代次數(shù)分別為2和3，此時(shí)，利用PSNR準(zhǔn)則得到的峰值信噪比僅比e值準(zhǔn)則得到的高0.05 dB。

然后，對(duì)PSNR準(zhǔn)則與e值準(zhǔn)則在不同迭代次數(shù)下的變化規(guī)律進(jìn)行比較。在Lena圖像中加入標(biāo)準(zhǔn)差為20的高斯白噪聲，不同迭代次數(shù)下，去噪圖像與含噪圖像、殘差圖像與含噪圖像之間的相關(guān)系數(shù)曲線如圖4（a）所示，PSNR與e值曲線如圖4（b）所示，其中，橫坐標(biāo)表示迭代次數(shù)，0值對(duì)應(yīng)初始NLM算法的去噪結(jié)果?？梢?jiàn)，隨著迭代次數(shù)的增加，去噪圖像與含噪圖像之間的相關(guān)系數(shù)首先逐漸增大，當(dāng)?shù)螖?shù)為7時(shí)增至最大；而殘差圖像中由于不斷有圖像信息被抽取出來(lái)，其與含噪圖像之間的相關(guān)性逐漸減弱。e值曲線與PSNR曲線呈現(xiàn)出一致的變化規(guī)律，當(dāng)?shù)螖?shù)為3時(shí)，兩者都達(dá)到最大值。

迭代次數(shù)取0～3時(shí)所得的去噪圖像及對(duì)應(yīng)殘差圖像的局部圖如圖5所示，可以看出，在達(dá)到最優(yōu)的迭代次數(shù)之前，隨著迭代的進(jìn)行，所得去噪圖像中的細(xì)節(jié)信息逐漸豐富，如帽子上的紋理部分，同時(shí)，對(duì)應(yīng)殘差圖像中包含的圖像結(jié)構(gòu)信息逐漸減弱。

4.2峰值信噪比與結(jié)構(gòu)相似性比較實(shí)驗(yàn)

在5幅測(cè)試圖像中分別加入標(biāo)準(zhǔn)差為10～50的高斯白噪聲，利用原始NLM［1］，NLFMT［15］，PPB 25it［6］和本文方法開(kāi)展去噪實(shí)驗(yàn)。所有方法中，相似性鄰域與搜索鄰域的大小分別取7×7與21×21［1］，NLM算法的濾波參數(shù)h =1.0×7×7×σ2；高斯濾波器大小為3×3，標(biāo)準(zhǔn)差為1.0。采用峰值信噪比與結(jié)構(gòu)相似指數(shù)準(zhǔn)則（Structural Sim ilarity Index Metric，SSIM）［19］對(duì)不同方法的去噪結(jié)果進(jìn)行定量比較。所得PSNR與SSIM的結(jié)果分別列于表1～表2中，每個(gè)單元格內(nèi)，左上、右上、左下、右下4個(gè)數(shù)據(jù)分別代表原始NLM，NLFMT，PPB 25it以及本文方法的去噪結(jié)果，最好的結(jié)果用加粗字體表示。

圖4　不同迭代次數(shù)下相關(guān)系數(shù)、e值與PSNR曲線

圖5　不同迭代次數(shù)下的去噪圖像與殘差圖像

表1　峰值信噪比比較結(jié)果

由表1中的數(shù)據(jù)可以看出，對(duì)于5幅測(cè)試圖像，在噪聲標(biāo)準(zhǔn)差小于30時(shí)，本文方法能夠得到最高的PSNR值；在噪聲水平較高時(shí)，由于原始NLM算法的權(quán)值分布矩陣計(jì)算不準(zhǔn)確，影響了本文方法的去噪效果，但仍然能夠取得與PPB 25it接近的峰值信噪比；與原始NLM算法相比，本文方法在峰值信噪比方面平均可以獲得1.2 dB的增益，最多可以提高1.54 dB。

由表2中的結(jié)果可以看出，本文方法在較低的噪聲水平下能夠獲得最高的結(jié)構(gòu)相似性，說(shuō)明本文方法具有良好的結(jié)構(gòu)保持能力。在噪聲水平較高時(shí)，本文方法會(huì)在某些去噪圖像中引入偽結(jié)構(gòu)信息，使得結(jié)構(gòu)相似性有所降低。

實(shí)驗(yàn)發(fā)現(xiàn)，本文方法一般需要迭代2～3次即可取得最優(yōu)的去噪結(jié)果。與取得先進(jìn)去噪效果的PPB 25it方法相比，本文方法在噪聲水平較低時(shí)，能夠得到優(yōu)于PPB 25it方法的峰值信噪比與結(jié)構(gòu)相似指數(shù)；在較高的噪聲水平下，本文方法可以取得與PPB 25it相近的去噪效果，但本文方法的計(jì)算復(fù)雜度遠(yuǎn)遠(yuǎn)低于該方法。4種方法的平均運(yùn)行時(shí)間對(duì)比見(jiàn)表3所示，所有方法都是在matlab R2013a平臺(tái)上運(yùn)行，計(jì)算機(jī)配置為3.20 GHz Intel Core i5 CPU，4.0 GB內(nèi)存。由于非局部均值算法的計(jì)算量主要集中在相似性權(quán)值的計(jì)算上，而本文方法采用固定權(quán)值非局部均值濾波，且只需進(jìn)行2～3次迭代，計(jì)算效率仍然較高。

表2　結(jié)構(gòu)相似指數(shù)比較結(jié)果

表3　 4種算法平均運(yùn)行時(shí)間對(duì)比（s）

4.3方法噪聲比較實(shí)驗(yàn)

以Lena圖像為例，對(duì)其加入標(biāo)準(zhǔn)差為20的高斯白噪聲，采用上述4種方法得到的去噪圖像與方法噪聲如圖6所示。由圖6中的結(jié)果可以看出，原始NLM算法所得去噪圖像被過(guò)度平滑，方法噪聲中包含明顯的圖像結(jié)構(gòu)信息；NLFMT算法所得方法噪聲中的圖像結(jié)構(gòu)信息明顯減弱，但仍清晰可見(jiàn)，所得去噪圖像中含有小波基函數(shù)帶來(lái)的偽影，視覺(jué)效果較差；PPB 25it與本文方法所得方法噪聲中已看不出明顯的圖像結(jié)構(gòu)特征，所得去噪圖像的邊緣及紋理結(jié)構(gòu)清晰，如帽沿及帽穗部分，本文方法所得去噪圖像中的細(xì)節(jié)特征更加自然。

4.4實(shí)際圖像去噪實(shí)驗(yàn)

本文方法不依賴于原始圖像，且噪聲的方差可以由估計(jì)方法得到［20］，因此，該方法可以應(yīng)用于對(duì)實(shí)際含噪圖像的降噪處理。圖7（a）所示是一幅實(shí)際的含噪圖像，分別利用NLM及本文方法對(duì)其進(jìn)行去噪，所得降噪圖像如圖7（b）和圖7（c）所示。

由圖7中的去噪結(jié)果可以看出，兩種方法都能夠有效地濾除噪聲；與原始的NLM算法相比，本文方法能夠更好地保留圖像的細(xì)節(jié)信息，如樹(shù)木的枝葉、瓦片的紋理、窗戶的邊緣等部分，如圖7（c）中虛線框所標(biāo)示，圖7（d）和圖7（e）給出了樹(shù)木枝葉的局部放大圖像。

圖6　去噪圖像與方法噪聲比較

圖7　實(shí)際圖像去噪結(jié)果

5　結(jié)束語(yǔ)

本文提出了一種適用于NLM算法的殘差圖像多級(jí)濾波方法，該方法將殘差圖像濾波與NLM算法有機(jī)結(jié)合，能夠充分地挖掘殘差中的圖像結(jié)構(gòu)信息；基于去噪圖像、殘差圖像與含噪圖像之間的相關(guān)性，提出了一種無(wú)需參考圖像的迭代停止準(zhǔn)則，能夠自適應(yīng)地選擇最佳濾波級(jí)數(shù)，得到與PSNR準(zhǔn)則類似的評(píng)價(jià)效果；去噪過(guò)程無(wú)需參考圖像，可應(yīng)用于對(duì)實(shí)際含噪圖像的降噪處理。實(shí)驗(yàn)結(jié)果表明，本文方法能夠顯著地提升非局部均值算法的去噪性能，且具有更好的細(xì)節(jié)保持能力。

本文從圖像信息提取與去噪效果評(píng)價(jià)兩個(gè)方面有效地利用了殘差圖像提供的信息。由此可見(jiàn)，殘差圖像在提升去噪算法的性能方面具有良好的潛力，探索更加有效且通用的殘差圖像濾波方法，是下一步需要開(kāi)展的工作。

［1］BUADESA，COLL B，and MOREL JM.A review of image denoising algorithm s，w ith a new one［J］.Multiscale Modeling and Simulation（SIAM Interdisciplanary Journal），2005，4（2）: 490-530.doi:10.1137/040616024.

［2］鐘瑩，楊學(xué)志，唐益明，等.采用結(jié)構(gòu)自適應(yīng)塊匹配的非局部均值去噪算法［J］.電子與信息學(xué)報(bào)，2013，35（12）:2908-2915. doi:10.3724/SP.J.1146.2013.00099.

ZHONG Ying，YANG Xuezhi，TANG Yim ing，etal.Non-local means denoising derived from structure-adap ted block matching［J］.JournalofElectronics&Information Technology，2013，35（12）:2908-2915.doi:10.3724/SP.J.1146.2013.00099.

［3］SUTOUR C，DELEDALLE C A，and AUJOL J F.Adaptive regularization of the n l-means:app lication to image and video denoising［J］.IEEE Transactions on Image Processing，2014，23（8）:3506-3521.doi:10.1109/TIP.2014.2329448.

［4］LU Lu，JIN W eiqi，and WANG Xia.Non-localm eans im age denoising w ith a soft threshold［J］.IEEE Signal Processing Letters，2015，22（7）:833-837.doi:10.1109/LSP.2014.2371332.

［5］DABOV K，F(xiàn)OI A，KATKOVNIK V，et al.Im age denoising by sparse 3D transform-domain collaborative filtering［J］. IEEE Transactions on Im age Processing，2007，16（8）: 2080-2095.doi:10.1109/TIP.2007.901238.

［6］DELEDALLE C A，DENISL，and TUPIN F.Iterativeweighted maximum likelihood denoising w ith probabilistic patch-based weights［J］.IEEE Transactions on Image Processing，2009，18（12）: 2661-2672.doi:10.1109/TIP.2009.2029593.

［7］羅亮，馮象初，張選德，等.基于非局部雙邊隨機(jī)投影低秩逼近圖像去噪算法［J］.電子與信息學(xué)報(bào)，2013，35（1）:99-105.doi: 10.3724/SP.J.1146.2012.00819.

LUO Liang，F(xiàn)ENG Xiangchu，ZHANG Xuande，et al.An image denoisingmethod based on non-local two-side random projection and low rank approximation［J］.Journal of Electronics&Information Technology，2013，35（1）:99-105. doi:10.3724/SP.J.1146.2012.00819.

［8］TALEBI H and M ILANFAR P.Global image denoising［J］. IEEE Transactions on Im age Processing，2014，23（2）:755-768. doi:10.1109/TIP.2013.2293425.

［9］LIU Ganchao，ZHONG Hua，and JIAO Licheng.Com paring noisy patches for image denoising:a double noise sim ilarity model［J］.IEEE Transactions on Image Processing，2015，24（3）: 862-872.doi:10.1109/TIP.2014.2387390.

［10］FENG Jianzhou，SONG Li，HUO Xiaom ing，et al.An optim ized pixel-w ise weighting approach for patch-based image denoising［J］.IEEE Signal Processing Letters，2015，22（1）:115-119.doi:10.1109/LSP.2014.2350032.

［11］BRUNET D，VRSCAY E R，and WANG Z.The use of residuals in image denoising［C］.Proceedings of the International Conference on Image Analysis and Recognition，Halifax，2009:1-12.doi:10.1007/978-3-642-02611-9_1.

［12］CHEN J，TANG C K，and WANG J.Noise brush:interactive high quality image noise separation［J］.ACM Transactions on Graphics，2009，28（5）:146:1-10.doi:10.1145/1618452.1618492.

［13］PYO Y，PARK RH，and CHANG S.Noise reduction in high-iso images using 3-d collaborative filtering and structure extraction from residual blocks［J］.IEEE Transactions on Consum er Electron ics，2011，57（2）:687-695.doi:10.1109/ TCE.2011.5955209.

［14］ZHONG H，YANG C，and ZHANG X H.A new weight for nonlocalm eans denoising using m ethod noise［J］.IEEE Signal Processing Letters，2012，19（8）:535-538.doi:10.1109/ LSP.2012.2205566.

［15］KUMAR B K S.Im age denoising based on non-local m eans filter and itsm ethod noise th resholding［J］.Signal，Im age and Video Processing，2012，7（6）:1211-1227.doi:10.1007/ s11760-012-0372-7.

［16］ROMANO Y and ELAD M.Imp roving K-SVD denoising by post-processing its method noise［C］.IEEE International Conference on Image Processing，Melbourne，2013:435-439. doi:10.1109/ICIP.2013.6738090.

［17］ROMANO Y and ELAD M.Boosting of image denoising algorithm s［J］.SIAM Journal on Imaging Sciences，2015，8（2）: 1187-1219.doi:10.1137/140990978.

［18］KONG X F，LI K，YANG Q X，et al.A new image quality metric for im age au to-denoising［C］.14th IEEE International Conference on Com puter Vision，Sydney，2013:2888-2895.doi: 10.1109/ICCV.2013.359.

［19］WANG Z，BOVIK A C，SHEIKH H R，et al.Image quality assessment:From error visibility to structural sim ilarity［J］. IEEE Transactions on Image Processing，2004，13（4）:600-612. doi:10.1109/TIP.2003.819861.

［20］CHANG SG，YU B，and VETTERLIB.Adaptive wavelet thresholding for image denoising and com pression［J］.IEEE Transactions on Image Processing，2000，9（9）:1532-1546.doi: 10.1109/83.862633.

孫偉峰：男，1982年生，副教授，研究方向?yàn)閳D像處理、目標(biāo)跟蹤.

戴永壽：男，1963年生，教授，研究方向?yàn)楣收显\斷與智能信息處理、計(jì)算機(jī)測(cè)控技術(shù).

Non-local Means Image Denoising with Multi-stage Residual Filtering

SUN Weifeng DAIYongshou
（Co llege of Inform ation and Con trol Engineering，China University of Petroleum，Qingdao 266580，China）

In order to sufficiently exploit the image in formation residing in the residual image for boosting the denoising performance of the Non-local Means（NLM）algorithm，a novelmulti-stage residual filtering method is proposed.Firstly，the Non-LocalMeans algorithm is applied to a noisy image to p roduce an initial denoised image and a weight distributingmatrix.Then the fixed-weight NLM algorithm is app lied to the residual image followed by a Gaussian filtering p rocess，which can extract the im age conten t ou t from the residualas a com pensation im age. The com pensation image is then added back to the denoised im age to generate an enhanced restored im age.An iterative scheme，whose p rincip le and feasibility are derived and p roved theoretically，is developed for the above filtering p rocedure；meanwhile a novelstopping criterion with no reference image required isp roposed to determ ine the op timal number of iterations adaptively.Experimental results demonstrate that the proposed stopping criterion behaves sim ilarly as the PSNR ru le，and compared w ith theoriginalNLM approach，the p roposed method can boost the denoising performance significantly w ith 1.2 dB PSNR gains achieved on average and m ore detail inform ation p reserved，w hile the com putational com p lexity is not apparently increased.

Image denoising；Non-LocalMeans（NLM）；Residual filtering；Stopping criterion

TN911.73

A

1009-5896（2016）08-1999-08

10.11999/JEIT 151227

2015-11-03；改回日期：2016-03-04；網(wǎng)絡(luò)出版：2016-04-14

孫偉峰sw f0217@163.com

國(guó)家自然科學(xué)基金（61501520），山東省自然科學(xué)基金（ZR2013FL035），中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金（14CX 02083A）Foundation Item s:The National Natu ral Science Foundation of China（61501520），Shandong Provincial Natural Science Foundation（ZR2013FL035），The Fundamental Research Funds for the Central Universities（14CX 02083A）