姚培娟 張志利 張亞娟

（黃河科技學(xué)院工學(xué)部 鄭州 450061）

1 引言

現(xiàn)有圖像去噪算法可以分為基于圖像稀疏性的字典學(xué)習(xí)去噪方法，變換域去噪算法以及空間域去噪算法等［1］。基于圖像稀疏性的字典方法有全局字典算法，K-SVD（K-Singular Value Decomposition）算法等［2］。變換域算法主要包括傅里葉變換去噪［3］、余弦變換去噪［4］、K-L 變換去噪［5］、小波變換去噪算法［6］以及三維塊匹配去噪算法［7］（Block matching and 3D filtering algorithm，BM3D）等?？臻g域算法主要有鄰域加權(quán)濾波算法［8］和雙邊濾波算法［9］（Bilateral Filters，BF），以及非局部均值算法（Non-local Means，NLM）等。

NLM 算法是應(yīng)用較廣的一種去噪方法，但是由于方法原理的限制，其去噪性能依然存在較大的提升空間。針對上述情況，國內(nèi)外學(xué)者先后提出基于圖像塊均值以及梯度方向的NLM 算法［10］，基于圖像塊均值和方差的NLM 算法等［11］。上述算法效率有所提升，但是去噪效果依然不夠理想。隨后Tasdizen 等提出主鄰域字典的NLM 算法［12］，以及采用極大似然估計(jì)的NLM 算法，基于相似性驗(yàn)證的非局部均值算法等（Similarity Validation Based Non-local Means，NLM-SVB）［13］，上述算法中引入?yún)?shù)較多，具有較大的不確定性，導(dǎo)致圖像細(xì)節(jié)容易損傷，從而影響圖像的去噪效果。

為解決上述NLM 算法的不足，本文提出一種改進(jìn)自適應(yīng)NLM 算法。該算法采用改進(jìn)的邊緣檢測算子提取圖像邊緣特征，將其應(yīng)用于鄰域相似度權(quán)重計(jì)算；然后聯(lián)合基于角度的鄰域相似性和基于距離的邊緣相似性來度量圖像塊之間的相似性，根據(jù)目標(biāo)圖像塊的局部結(jié)構(gòu)信息，對目標(biāo)圖像塊內(nèi)部相鄰像素加權(quán)進(jìn)行中心像素估計(jì)，從而有效保留局部細(xì)節(jié)信息。經(jīng)仿真模擬驗(yàn)證，改進(jìn)的NLM 算法具有較好的去噪效果和邊緣紋理細(xì)節(jié)保護(hù)能力，BM3D算法具有均等圖像質(zhì)量，計(jì)算效率高。

2 常規(guī)NLM算法

NLM 算法通過圖像鄰域間的歐式距離計(jì)算相似性權(quán)值，對相似像素值進(jìn)行加權(quán)平均實(shí)現(xiàn)去噪目的［14］，其公式為

其中，y(i)表示原始含噪圖像像素，x?(j)表示去噪后的像素值，Sj表示以像素j 為中心的搜索窗口，wNL(i,j)表示以像素點(diǎn)i和像素點(diǎn)j為中心的兩鄰域相似性權(quán)重，其表達(dá)式為

其中，y(Nj)表示以j 為中心的鄰域圖像塊像素，y(Ni)表示以i 為中心的鄰域圖像塊像素，h 表示衰減系數(shù)，用于控制濾波強(qiáng)度。

然而，上述權(quán)重計(jì)算公式在高比例噪聲情況下很難準(zhǔn)確得到不同鄰域圖像塊之間的相似性，導(dǎo)致該方法的去噪能力有待提升。

3 邊緣約束的改進(jìn)NLM算法

3.1 邊緣檢測算子

目前，按照邊緣檢測微分算子的使用階數(shù)，傳統(tǒng)的邊緣檢測算子可以分為一階算子和二階算子，其中一階算子包括Roberts算子，Prewitt算子，Sobel算子，Kirsch 算子等；二階算子包括Marr-Hildreth算子，Canny 算子，Laplacian 算子等［15］。然而上述常規(guī)算子在含噪圖像中，很難準(zhǔn)確定位和識(shí)別邊緣特征，針對上述情況，本文采用改進(jìn)的多尺度小波變換方法提取圖像的邊緣特征。

多尺度小波變換通過在每個(gè)尺度下沿梯度方向?qū)ふ夷５臉O大值，利用閾值法得到圖像邊緣特征信息，從而有效均衡噪聲去除和邊緣定位之間的協(xié)同［16］。將二維高斯函數(shù)作為尺度函數(shù)，其表達(dá)式為

則X，Y兩個(gè)方向上的一階導(dǎo)數(shù)分別為

其中，x，y 表示坐標(biāo)位置，σ為高斯函數(shù)的標(biāo)準(zhǔn)差（尺度系數(shù)）。對圖像進(jìn)行邊緣檢測時(shí)，尺度選取尤為重要，尺度較大容易模糊邊緣信息，尺度較小容易遭受噪聲干擾提取到偽邊緣信息。因此，本文提出一種自適應(yīng)尺度參數(shù)計(jì)算方式提取邊緣特征信息，其示意圖如圖1所示。

圖1 三種尺寸鄰域示意圖

圖1為一個(gè)M×N圖像的三種鄰域尺寸示意圖，其自適應(yīng)求取最優(yōu)尺寸步驟如下，首先對圖1（a）鄰域中的頂端像素點(diǎn)（x-1，y-1）與其他像素值作差，將所有結(jié)果加權(quán)求取平均值；然后對圖1（b）鄰域中的邊界像素點(diǎn)（x-1，y）與其他像素值作差，并加權(quán)求取平均值；最后對圖1（c）鄰域中的中間像素點(diǎn)（x，y）和其他的像素點(diǎn)作差，加權(quán)求取均值。遍歷圖像的所有像素后，得到一個(gè)M×N 尺寸的矩陣，對矩陣求取均值就得到當(dāng)前圖像的最優(yōu)變換尺度。

計(jì)算最優(yōu)尺度參數(shù)σ后，再次選取σ-1 和σ+1分別作為另外兩個(gè)尺度。然后在三個(gè)尺度下，將式（5）和原始圖像進(jìn)行卷積運(yùn)算實(shí)現(xiàn)多尺度分解，計(jì)算每個(gè)尺度中像素點(diǎn)的梯度幅值A(chǔ)grd和相角α，公式如下：

由式（7）可知，相角范圍為[-π/2,π/2]，將其分為四個(gè)區(qū)間，不同區(qū)間分別對應(yīng)不同的梯度方向（0°（水平）方向，45°方向，90°（垂直）方向，135°方向），其具體示意圖如圖2所示。

圖2 梯度方向劃分示意圖

在不同尺度下沿不同梯度方向?qū)ふ耶?dāng)前位置的模極大值，通過直方圖統(tǒng)計(jì)邊緣數(shù)目，得到高低閾值，采用雙閾值連接模值得到邊緣特征［24］。

3.2 邊緣約束的鄰域相似性度量

本文應(yīng)用主成分分析權(quán)重函數(shù)對噪聲表現(xiàn)的魯棒性特點(diǎn)，定義基于邊緣特征約束的相似性權(quán)重度量公式為

其中，fd表示主成分分析的子空間，SI( )Ni表示以像素i 為中心的邊緣像素集，則表示分別以像素i 和j 為中心的圖像塊的邊緣相似度。h和h′分別表示濾波參數(shù)，用以調(diào)節(jié)邊緣特征和平滑特征對相似度權(quán)重的影響。

為提升算法對邊緣/紋理細(xì)節(jié)的保護(hù)，本文將相似性權(quán)重度量公式進(jìn)行改進(jìn)，則公式為

3.3 濾波參數(shù)的自適應(yīng)選取

鄰域中心像素疊加高斯噪聲之后，主成分分析的子空間的相似性產(chǎn)生差異，因此，根據(jù)圖像結(jié)構(gòu)特征自適應(yīng)設(shè)置濾波參數(shù)［25］。根據(jù)高斯分布特征，應(yīng)該有一定的像素分布在不超過噪聲標(biāo)準(zhǔn)差β倍的范圍內(nèi)，即：

則濾波參數(shù)h 和噪聲方差的近似關(guān)系可以表示為

不同的圖像細(xì)節(jié)結(jié)構(gòu)，式（12）中的β選取不同，因此采用最大類間方差法確定自適應(yīng)閾值th，將邊緣特征像素值和自適應(yīng)閾值映射至區(qū)間［0，6］，得到紋理結(jié)構(gòu)圖IV和閾值re，則β可以表示為

其中，re是映射到區(qū)間［0，6］之間的閾值，由式（13）可知，當(dāng)前圖像塊平滑，選取較大的β1，在紋理細(xì)節(jié)區(qū)域時(shí)，選取較小的β2。

3.4 目標(biāo)像素去噪

為提高圖像邊緣/紋理細(xì)節(jié)的保護(hù)能力，本文對中心像素以及周圍的鄰域像素進(jìn)行加權(quán)平均生成目標(biāo)圖像塊，目標(biāo)圖像塊的建立方式如圖3所示。

圖3 目標(biāo)圖像塊搜索建立執(zhí)行過程示意圖

圖3中，中心的3×3圖像塊為目標(biāo)圖像塊，左下角為搜索鄰域圖像塊，虛線箭頭表示搜索路徑。將圖像塊作為輸入?yún)?shù)，式（9）變?yōu)?/p>

其中，Ni,Nj分別表示以i，j 為中心的圖像塊，結(jié)合式（12）和式（13），則包含目標(biāo)像素i 的目標(biāo)圖像塊更新為

其中，W(Ni,Nj)表示式（14）定義的數(shù)據(jù)塊塊相似性權(quán)重公式，?(Ni),y(Nj)分別表示目標(biāo)圖像塊的像素集和原始含噪鄰域圖像塊像素集。

確定目標(biāo)圖像塊之后，本文根據(jù)多窗口劃分思想實(shí)現(xiàn)中心像素值去噪處理，窗口劃分如圖4 所示。

圖4 多窗口劃分示意圖

由圖4 可知，在5×5 窗口中，包含中心像素（黑色實(shí)心圓）的窗口一共有4個(gè)六邊形，4個(gè)五邊形，1個(gè)正方形。因此根據(jù)3×3 的中心目標(biāo)圖像塊與周圍8 個(gè)多邊形鄰域圖像像素塊的相似性進(jìn)行加權(quán)，則最終的去噪結(jié)果表達(dá)式為

其中，W(i,jk)表示式（9）定義的相似性權(quán)重公式，?(i),y(jk)分別表示去噪后像素和原始含噪像素。

4 方法測試及結(jié)果分析

4.1 標(biāo)準(zhǔn)灰度圖像測試

標(biāo)準(zhǔn)灰度圖像Baboon 圖像添加均值為0、方差為20 的高斯白噪，則含噪圖像和去噪結(jié)果局部放大圖如圖5 所示，在本次仿真測試中，搜索窗口大小為13×13。

由圖5 可知，常規(guī)NLM 方法去噪結(jié)果圖像模糊，存在部分“虛假紋理”；文獻(xiàn)［10］方法相對較好，但是圖像質(zhì)量差，細(xì)節(jié)信息損傷嚴(yán)重；文獻(xiàn)［11］方法細(xì)節(jié)保持較好，但是去噪能力有限；文獻(xiàn)［12］方法去噪能力相對較好，但是邊緣紋理等細(xì)節(jié)模糊；文獻(xiàn)［13］方法相對以上算法效果較好，但是其去噪結(jié)果仍舊有待提高；BM3D 去噪結(jié)果細(xì)節(jié)最為清晰，圖像的“鼻梁”和“胡須”輪廓清楚。本文方法的去噪結(jié)果和BM3D結(jié)果近似，優(yōu)于另外五種去噪結(jié)果。

圖5 方差為20的含噪測試圖像及其去噪結(jié)果

為進(jìn)一步定量分析不同方法的去噪效果，采用峰值信噪比（PSNR）和平均結(jié)構(gòu)相似性（Mean Structural Similarity，MSSIM）作為評(píng)價(jià)標(biāo)準(zhǔn)，將不同測試結(jié)果的PSNR 和MSSIM 進(jìn)行統(tǒng)計(jì)，結(jié)果如表1所示。

由表1 可知，常規(guī)NLM 算法PSNR 平均提升43.7%，平均結(jié)構(gòu)相似性平均值為0.75；文獻(xiàn)［10］算法PSNR 值平均提升46.8%，平均結(jié)構(gòu)相似性平均值為0.78；文獻(xiàn)［11］算法平均提升48.6%，平均結(jié)構(gòu)相似性平均值為0.79；文獻(xiàn)［12］算法的PSNR 平均提升50.3%，平均結(jié)構(gòu)相似性平均值為0.80；文獻(xiàn)［13］算法的PSNR 平均提升52.5%，平均結(jié)構(gòu)相似性平均值為0.805；BM3D 算法的PSNR 平均提升55.3%，平均結(jié)構(gòu)相似性平均值為0.83；而本文算法的PSNR 值平均提升54.9%，平均結(jié)構(gòu)相似性平均值為0.82。由此可知，本文算法的去噪效果和BM3D 算法的去噪效果幾乎一致，證實(shí)本文算法的有效性和合理性。

表1 不同方法的Baboon測試結(jié)果PSNR（dB）-MSSIM

4.2 實(shí)際B超圖像驗(yàn)證

常規(guī)去噪算法在進(jìn)行B 超圖像噪聲去除時(shí)，圖像的邊緣/細(xì)節(jié)紋理和去噪效果之間很難同時(shí)兼顧，為進(jìn)一步驗(yàn)證本文算法的測試效果，將本文算法和其他六種算法的B 超圖像去噪效果進(jìn)行對比，其效果如圖6、圖7所示。

圖6 不同方法的B超圖像去噪結(jié)果

圖7 不同算法的去噪殘差

由圖6 和圖7 可知，常規(guī)NLM 算法結(jié)果紋理細(xì)節(jié)較為模糊，其去噪殘差中存在較多的連續(xù)結(jié)構(gòu)信息，表明對圖像的細(xì)節(jié)損傷較為嚴(yán)重；文獻(xiàn)［10］和文獻(xiàn)［11］結(jié)果相對較好，但是殘差中依然存在一定的連續(xù)結(jié)構(gòu)信息；文獻(xiàn)［12］和文獻(xiàn)［13］圖像質(zhì)量提升較好，但是殘差中存有少量的連續(xù)紋理以及“大顆?！卑郀钤肼?，說明去噪過度；BM3D 方法和本文方法去噪結(jié)果中，幾乎不含連續(xù)的紋理結(jié)構(gòu)，去噪殘差中的噪聲顆粒分布均勻大小適中，表明去噪結(jié)果好，對紋理細(xì)節(jié)保護(hù)好。

由于實(shí)際B 超圖像幾乎沒有不含噪聲的圖像，因此常規(guī)算法很難評(píng)價(jià)，針對上述情況，本文采用完全無參考算法（Natural Image Quality Evaluator，NIQE）對去噪結(jié)果進(jìn)行評(píng)價(jià)，同時(shí)統(tǒng)計(jì)不同結(jié)果的專家評(píng)估值，則B 超圖像處理前后的評(píng)價(jià)值如表2所示。

表2 實(shí)際B超去噪結(jié)果的評(píng)價(jià)值

由圖7、圖8 和表2 可知，本文方法和BM3D 方法的去噪效果最好，二者結(jié)果的去噪評(píng)價(jià)值最高。然而BM3D 算法屬于混合變換域去噪方法，復(fù)雜度高，其處理時(shí)間約為本文算法的3 倍左右，因此很難滿足圖像處理實(shí)時(shí)計(jì)算的需求，很大程度上限制BM3D算法的使用范圍。

5 結(jié)語

本文提出一種基于邊緣約束的改進(jìn)NLM 算法，在去除噪聲的前提下能夠保持邊緣/紋理細(xì)節(jié)。經(jīng)一系列的仿真實(shí)驗(yàn)證實(shí)，本文方法相對于常規(guī)去噪算法，在邊緣/紋理細(xì)節(jié)保持和噪聲壓制等方面具有更好的優(yōu)勢。和BM3D 算法相比，在具有相近去噪能力情況下，具有更高的計(jì)算效率。