王彩玲，孫昌霞

(1.河南警察學(xué)院網(wǎng)絡(luò)安全系，河南 鄭州 450046;2.河南農(nóng)業(yè)大學(xué)信管學(xué)院計(jì)算機(jī)系，河南 鄭州 450046)

0 引言

由于圖像在獲取、傳輸和數(shù)字化過(guò)程中會(huì)遇到不同光照的影響，傳感器溫度發(fā)生變化等情況，獲取的圖像中會(huì)引入隨機(jī)噪聲。圖像去噪的方法很多，常分為頻域?yàn)V波和空間域?yàn)V波。常用的頻域?yàn)V波方法有小波去噪[1]、高通濾波、維納濾波[2]等；空域?yàn)V波的方法很多，如偏微分方程、變分法[3]、統(tǒng)計(jì)方法[4]等，這些方法均有較多的應(yīng)用。經(jīng)典的偏微分方程法是一種各向同性的濾波算法，能很好地去除噪聲，但這類方法在對(duì)噪聲成分進(jìn)行平均的同時(shí)也對(duì)圖像信息進(jìn)行了平均，導(dǎo)致處理后的圖像變得較為模糊。Yaroslavsky根據(jù)像素灰度的差異來(lái)定義像素點(diǎn)之間的相似性，得到了Yaroslavsky濾波[5]。Tomasi在Yaroslavsky濾波的基礎(chǔ)上，結(jié)合像素點(diǎn)間的灰度差異和空間距離來(lái)定義了像素間的相似性，得到了雙邊濾波[6]。近年來(lái)，對(duì)雙邊濾波也出現(xiàn)了一些改進(jìn)算法[7-11]，主要集中于濾波參數(shù)的選擇、算法效率的提升等方面，其中，文獻(xiàn)[12]提出了距離核的概念，將其應(yīng)用于像素相似性權(quán)值的定義中，在細(xì)節(jié)豐富的區(qū)域中，算法的去噪性能有很大的提升。此外，文獻(xiàn)[13]提出了一種基于方向波域混合高斯模型的合成孔徑雷達(dá)圖像去噪方法。與小波域不同，文獻(xiàn)[14]對(duì)合成孔徑雷達(dá)圖像進(jìn)行復(fù)Shearlet變換，并對(duì)系數(shù)進(jìn)行高斯混合模型建模，有效提升去噪圖像的峰值信噪比。

針對(duì)圖像像素間相似性衡量不準(zhǔn)確而導(dǎo)致去噪性能下降的問(wèn)題，本文利用高斯混合模型對(duì)像素周圍鄰域內(nèi)的像素信息進(jìn)行建模，結(jié)合高斯混合模型間的L2范數(shù)與空間距離來(lái)定義像素點(diǎn)之間的相似性權(quán)值，提出基于L2范數(shù)的高斯混合模型空域去噪方法。

1 基于高斯混合模型的圖像去噪原理

基于高斯混合模型的圖像去噪的基本思想：首先利用每一個(gè)像素點(diǎn)周圍鄰域的像素信息估計(jì)出高斯混合模型的參數(shù)，將像素信息建模為高斯混合模型，得到預(yù)濾波圖像；然后計(jì)算兩像素點(diǎn)對(duì)應(yīng)高斯混合模型間的L2范數(shù)，結(jié)合L2范數(shù)值和兩像素點(diǎn)間的空間位置距離來(lái)定義兩像素點(diǎn)間的相似性權(quán)值；最后對(duì)每一個(gè)像素點(diǎn)進(jìn)行加權(quán)平滑濾波，得到去噪后的圖像。

與經(jīng)典的空域去噪方法不同，基于高斯混合模型的圖像去噪方法是根據(jù)像素點(diǎn)周圍信息的統(tǒng)計(jì)差異性大小來(lái)定義像素點(diǎn)間的相似性權(quán)值，其具體步驟如下：

1) 對(duì)某一個(gè)像素點(diǎn)，利用其周圍鄰域內(nèi)的像素信息估計(jì)出高斯混合模型的參數(shù)，得到一個(gè)高斯混合模型；

2) 對(duì)圖像中的每一個(gè)像素均估計(jì)出一個(gè)高斯混合模型，所有的高斯混合模型構(gòu)成了一個(gè)預(yù)濾波圖像；

3) 計(jì)算兩像素點(diǎn)對(duì)應(yīng)高斯混合模型間的L2范數(shù)，并結(jié)合兩像素點(diǎn)間的空間距離來(lái)定義像素點(diǎn)間的相似性；

4) 對(duì)圖像中的每一個(gè)像素均進(jìn)行加權(quán)平滑濾波，得到濾波后的圖像。

從上述步驟可知，基于高斯混合模型的圖像去噪是根據(jù)兩像素點(diǎn)x與y間的空間距離和其對(duì)應(yīng)高斯混合模型間的L2范數(shù)來(lái)定義相似性權(quán)值w(x,y)，并進(jìn)行加權(quán)平滑濾波得到濾波后的圖像，濾波公式如下：

(1)

(2)

式(2)中，G(x)和G(y)分別表示像素點(diǎn)x與y處的高斯混合模型，d2(x,y)=|x-y|2表示像素點(diǎn)x與y間的空間距離，參數(shù)r是濾波控制系數(shù)。r越大，圖像的平滑效果越好，圖像的細(xì)節(jié)信息損失越多；r越小，圖像的平滑效果越差，圖像細(xì)節(jié)信息損失越小。因此，在濾波過(guò)程中需要選擇合適的濾波控制系數(shù)r?；诟咚够旌夏Ｐ偷膱D像去噪方法的基本思想可表示如圖1所示。

圖1 基于高斯混合模型的圖像去噪Fig.1 Image denoisingbased on Gaussian mixture model

2 基于L2范數(shù)的高斯混合模型去噪

對(duì)每一個(gè)像素點(diǎn)，利用其周圍鄰域內(nèi)的像素信息估計(jì)出高斯混合模型的參數(shù)，得到一高斯混合模型。像素點(diǎn)對(duì)應(yīng)的高斯混合模型表示該像素周圍鄰域內(nèi)像素信息的局部灰度強(qiáng)度和細(xì)節(jié)豐富程度，利用兩高斯混合模型間的L2范數(shù)，并結(jié)合像素點(diǎn)間的空間距離來(lái)定義像素間的相似性權(quán)值。依據(jù)該相似性權(quán)值對(duì)圖像中的每一個(gè)像素進(jìn)行加權(quán)平滑濾波，下面主要介紹像素信息的高斯混合模型建模以及兩高斯混合模型間的L2范數(shù)計(jì)算。

2.1 像素信息的高斯混合模型估計(jì)

對(duì)某一像素點(diǎn)x，其周圍鄰域窗口Ψ(x)的大小為M×M，鄰域內(nèi)包含M2個(gè)像素點(diǎn)，位置(i,j)處的像素灰度值為I(i,j)，利用這些像素信息估計(jì)出一高斯混合模型，公式如下：

(3)

(4)

式(4)中，z表示類別的標(biāo)簽，P(·)表示概率大小。式(4)給出γim,i=1,…,n2,m=1,…,K的計(jì)算公式，進(jìn)一步，依據(jù)式(4)可計(jì)算出其余的參數(shù)：

(5)

根據(jù)式(4)和式(5)可以利用EM算法迭代估計(jì)出高斯混合模型的參數(shù)，參數(shù)估計(jì)的具體步驟如下：

1) 對(duì)參數(shù)θ取初值；

4) 重復(fù)步驟2)和步驟3)，直至收斂。

采用上述步驟可以估計(jì)出每個(gè)像素點(diǎn)的高斯混合模型的參數(shù)，所有像素點(diǎn)對(duì)應(yīng)的高斯混合模型即構(gòu)成了預(yù)濾波圖像，在預(yù)濾波圖像中，結(jié)合兩像素點(diǎn)對(duì)應(yīng)高斯混合模型間的L2范數(shù)和空間距離來(lái)計(jì)算兩像素點(diǎn)之間的相似性權(quán)值。為了分析利用EM算法進(jìn)行高斯混合模型參數(shù)估計(jì)的有效性，以某次參數(shù)估計(jì)為例，對(duì)于給定7×7的圖像區(qū)域的像素信息，利用此49個(gè)像素進(jìn)行高斯混合模型的參數(shù)估計(jì)，選取模型的成份數(shù)為5，其參數(shù)估計(jì)的中間結(jié)果如表1所示。

表1 高斯混合模型的參數(shù)估計(jì)的中間結(jié)果Tab.1 Intermediate results of parameter estimation

從表1的結(jié)果可以看出，利用EM算法進(jìn)行高斯混合模型的參數(shù)估計(jì)，在第7步時(shí)算法收斂，并且不同高斯成分對(duì)應(yīng)的權(quán)值以及均值和方差均不相同。

2.2 L2范數(shù)的計(jì)算

假設(shè)每個(gè)像素點(diǎn)周圍的像素信息建模為混合成份數(shù)均為K的高斯混合模型，對(duì)于兩像素點(diǎn)x與y對(duì)應(yīng)的高斯混合模型為：

(6)

其間的差異性度量計(jì)算方法較多，如KL散度、蒙特卡洛方法等，但KL散度并沒(méi)有解析表達(dá)式，只能通過(guò)近似的方式來(lái)計(jì)算，蒙特卡洛方法的計(jì)算量太大，均不適合。為此，本文采用L2范數(shù)來(lái)計(jì)算兩高斯混合模型間的差異性，高斯混合模型間的L2范數(shù)具有解析表達(dá)式，便于算法的實(shí)施。

對(duì)于高斯混合模型G(x)與G(y)，其間的L2范數(shù)可計(jì)算如下：

(7)

上式展開(kāi)可得，

(8)

將式(6)代入式(8)，化簡(jiǎn)可得，

(9)

根據(jù)高斯分布的積分性質(zhì)可知，

(10)

將式(10)代入式(9)，化簡(jiǎn)可得，

(11)

根據(jù)式(11)，可計(jì)算兩高斯混合模型間的L2范數(shù)距離。

The acoustic impedance of the perforated plate is related to the acoustic resistance and mass reactance of the orifice on the assumption that the interaction between the orifices,if they are sparsely distributed,can be neglected.

3 實(shí)驗(yàn)驗(yàn)證

3.1 評(píng)價(jià)指標(biāo)

對(duì)去噪算法的評(píng)價(jià)須從兩個(gè)指標(biāo)進(jìn)行評(píng)價(jià)，即去除噪聲的多少和圖像細(xì)節(jié)信息損失多少。噪聲去除的越多，表明去噪算法的去噪能力越強(qiáng)；噪聲去除的越少，表明算法的去噪能力越差。同時(shí)，圖像細(xì)節(jié)信息保留越多，表明圖像的質(zhì)量越好；圖像的細(xì)節(jié)信息保留越少，表明圖像的質(zhì)量越差。通常在去除圖像噪聲的同時(shí)，盡可能多地保留圖像的細(xì)節(jié)信息。為此，本文采用峰值信噪比(peak signal to noise ratio，PSNR)來(lái)衡量算法的去噪能力，同時(shí)，采用結(jié)構(gòu)相似度因子(structural similarity index，SSIM)[15]來(lái)評(píng)價(jià)去噪后圖像的質(zhì)量好壞，也即圖像的細(xì)節(jié)信息保留的能力。PSNR的值越大，SSIM的值越接近1，說(shuō)明算法的去噪能力越強(qiáng)，且圖像的細(xì)節(jié)信息保留越多。

3.1.1PSNR

在圖像去噪的算法評(píng)價(jià)中，PSNR常用于圖像去噪能力的評(píng)價(jià)指標(biāo)，PSNR可計(jì)算如下：

(12)

式(12)中，vmax表示圖像中的最大像素值，在0～255的圖像中，vmax=255；v(i)和u(i)分別表示噪聲圖像和無(wú)噪聲圖像的灰度值，mean{[v(i)-u(i)]2}是指噪聲的平均功率。

3.1.2SSIM

結(jié)構(gòu)相似度因子是根據(jù)人類視覺(jué)機(jī)制來(lái)定義的圖像去噪的評(píng)價(jià)指標(biāo)，從視覺(jué)上來(lái)看圖像，圖像的細(xì)節(jié)信息越多，則圖像越清晰，反之，圖像越模糊。SSIM是結(jié)合圖像的對(duì)比度、邊緣結(jié)構(gòu)和圖像亮度等三方面的特性，通過(guò)加權(quán)乘積得到的評(píng)價(jià)指標(biāo)，SSIM可由如下公式計(jì)算：

(13)

在本文中，c1=(k1L)2，c2=(k2L)2，L表示像素值的取值范圍，k1和k2是權(quán)系數(shù)，一般情況下，k1=0.01，k2=0.03。一幅圖像的SSIM是不同窗口內(nèi)SSIM的均值mSSIM。

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

為了驗(yàn)證本文算法(GMMD)的有效性，首先需要分析算法的處理過(guò)程，包括相似性權(quán)值的分析，待濾波圖像的求解，然后分析去噪效果。算法在聯(lián)想筆記本上運(yùn)行，其環(huán)境為Windows 10系統(tǒng)，處理軟件為matlab 2019a，CPU采用英特爾i7-8850H，內(nèi)存16 G。

對(duì)于一幅“Lena”圖，從中選擇A、B、C、D 4個(gè)區(qū)域，依據(jù)相似性權(quán)值的計(jì)算公式分別計(jì)算A與B、C、D之間的相似性。結(jié)果如圖2所示：A與B之間的像素信息差異較大，其相似性較低，只有0.31；A與C之間的像素信息比較接近，其相似性高達(dá)0.92；A與D之間的顏色信息有一定的差異，其相似性適中，有0.75。上述相似性的結(jié)果能較好地反映出不同區(qū)域間的差異性大小，為加權(quán)平均濾波奠定了基礎(chǔ)。

圖2 圖像不同區(qū)域的相似性Fig.2 Similarity of different regions in image

為了進(jìn)一步分析算法的有效性，需要給出圖像的預(yù)濾波圖像。如圖3所示，分別給出了原圖像、噪聲方差為0.01的圖像、均值圖像以及方差圖像，其中均值和方差圖像為預(yù)濾波圖像，說(shuō)明了圖像在均值和方差域的信息。

從圖3的結(jié)果可以看出，均值圖能較好地反映出圖像的像素信息，而方差圖能反映出圖像的邊緣細(xì)節(jié)信息，利用高斯混合模型能較好地建模圖像區(qū)域像素信息的局部灰度強(qiáng)度和細(xì)節(jié)豐富程度。

從圖4—圖8的去噪結(jié)果可以看出，GMMD算法和KBF算法的去噪性能明顯優(yōu)于MF和BF的性能，GMMD算法能很好地去噪圖像噪聲，同時(shí)能較大程度上保留圖像的細(xì)節(jié)信息。為了從視覺(jué)上更好地看出不同方法的去噪性能，在圖像中選取某一區(qū)域，給出其去噪結(jié)果，以“Lena”圖為例，選取的圖像區(qū)域以及去噪結(jié)果如圖9所示。

圖3 圖像的預(yù)濾波結(jié)果Fig.3 Pre-filtering results

圖4 噪聲圖像Fig.4 Noise images

圖5 經(jīng)BF去噪后的圖像Fig.5 Denoised images by BF

圖6 經(jīng)MF去噪后的圖像Fig.6 Denoised images by MF

從圖9的結(jié)果可看出：圖像區(qū)域經(jīng)MF去噪后細(xì)節(jié)信息損失較多，圖像較模糊；利用BF和KBF進(jìn)行去噪后，得到的結(jié)果雖然細(xì)節(jié)保持較好，但噪聲較多，去噪效果較差；而本文方法得到的去噪結(jié)果不僅在圖像細(xì)節(jié)信息上保持較好，而且能去除噪聲，去噪效果較好。為了定量分析和比較GMMD算法的性能優(yōu)勢(shì)，在不同噪聲背景下，利用不同去噪算法進(jìn)行去噪實(shí)驗(yàn)，計(jì)算PSNR和mSSIM的值，所得結(jié)果如表2和表3所示。

圖7 經(jīng)KBF去噪后的圖像Fig.7 Denoised images by KBF

圖8 經(jīng)GMMD去噪后的圖像Fig.8 Denoised images by GMMDX

從表2和表3的結(jié)果可得，GMMD算法和KBF算法的去噪能力和細(xì)節(jié)信息保留能力均優(yōu)于MF和BF算法。在比較平坦的圖像“Lena”和“Einstein”去噪結(jié)果中，GMMD的去噪能力和細(xì)節(jié)信息保持能力均優(yōu)于KBF算法；然而，在圖像細(xì)節(jié)信息較豐富的圖像“Academy”和“Mandrill”去噪結(jié)果中，GMMD算法的去噪性能提升并沒(méi)有像平坦區(qū)域那么大。結(jié)合圖像像素信息的局部平均灰度強(qiáng)度和細(xì)節(jié)豐富程度來(lái)定義像素間的相似性權(quán)值，比傳統(tǒng)的方法更準(zhǔn)確。

圖9 局部區(qū)域的去噪結(jié)果Fig.9 Denoising results of local area

表2 不同去噪算法的PSNRTab.2 PSNR of different denoising algorithms dB

表3 不同算法的mSSIM(σ2=0.01)Tab.3 mSSIM of different denoising algorithms

4 結(jié)論

本文提出了基于L2范數(shù)的高斯混合模型空域去噪方法。該方法將像素點(diǎn)周圍像素的局部灰度信息的統(tǒng)計(jì)特性建模為高斯混合模型，結(jié)合高斯混合模型間的L2范數(shù)和像素位置的空間距離來(lái)定義像素點(diǎn)之間的相似性權(quán)值，以提高像素間相似性度量的精確度，利用加權(quán)平滑濾波來(lái)對(duì)圖像進(jìn)行去噪處理。實(shí)驗(yàn)結(jié)果表明，利用高斯混合模型可更加準(zhǔn)確地衡量像素間的相似性，同時(shí)，該去噪方法可提高圖像的去噪效果，并能較好地保留圖像的細(xì)節(jié)信息。如何在細(xì)節(jié)豐富的圖像區(qū)域中提升算法的去噪性能，是下一步研究的重點(diǎn)。