自適應(yīng)擴(kuò)散系數(shù)優(yōu)化的圖像降噪算法

鄭香香，劉艷莉

(1. 北京市健宮醫(yī)院呼吸與危重癥醫(yī)學(xué)科，北京100054；2. 中北大學(xué)信息與通信工程學(xué)院，山西 太原 030051)

1 引言

數(shù)字圖像處理技術(shù)在當(dāng)今的信息獲取中起著舉足輕重的地位，而圖像降噪又是數(shù)字圖像處理重要的根基。圖像在生出或傳輸?shù)倪^程中都會受到外界環(huán)境或設(shè)備等的影響，使采集到的圖像一些重要信息被掩蓋，對圖像的后續(xù)處理帶來很大困擾，尤其對醫(yī)療疾病診斷行業(yè)的影響巨大，因此，如何有效去除噪聲，并更好地保持圖像原有的紋理細(xì)節(jié)等重要信息是圖像降噪的難題。

Perona和Malik于1990年提出了各向異性擴(kuò)散(Anisotropic Diffusion)模型之后(即PM模型)，國內(nèi)外研究者開始從各個(gè)角度研究如何提高PM模型的降噪性能[1-3]。為了有效權(quán)衡邊緣細(xì)節(jié)和噪聲間的矛盾，研究人員引入局部方差，信息熵，差分曲率、非局部均值等信息[4-8]，構(gòu)造不同的擴(kuò)散系數(shù)。Gilboa 等人[9]提出了一種保持斜坡的復(fù)擴(kuò)散方法，主要利用圖像的虛部作為拉普拉斯算子的近似來控制擴(kuò)散，Barbu 等[10]討論了在極小增長條件下基于梯度凸函數(shù)極小化的變分模型。郭業(yè)才等人[11]提出了基于脈沖耦合神經(jīng)網(wǎng)絡(luò)和圖像熵的各向異性擴(kuò)散模型，將噪聲圖像獲取的熵序列和圖像梯度共同構(gòu)成檢測因子，可以更有效保留圖像紋理等微弱信息?；贖essian矩陣特征值的結(jié)構(gòu)檢測算子[12-13]，能根據(jù)圖像特征調(diào)整擴(kuò)散系數(shù)的梯度閾值，以均衡噪聲與細(xì)節(jié)的擴(kuò)散程度。wang等人[14]提出了一個(gè)新的二階和四階各向異性方程來有效地去噪。Yu等[15]提出了一種核各向異性擴(kuò)散的圖像去噪方法；Chao[16]等將局部灰度方差引入到PM模型的擴(kuò)散系數(shù)中。這些方法的擴(kuò)散系數(shù)都依賴于圖像梯度值，因此也削弱了圖像的紋理和細(xì)節(jié)。

本文提出一種新的基于PM模型的圖像降噪算法。首先利用經(jīng)典的PM模型處理圖像，獲得殘差能量，主要包含圖像的噪聲和紋理信息；然后提取圖像殘差信息中的紋理部分，在二次PM迭代中，將其進(jìn)一步應(yīng)用到原圖像中；此時(shí)的紋理與噪聲同時(shí)存在，利用絕對差值排序算子很好地區(qū)分紋理與噪聲的特性，進(jìn)一步處理圖像；將圖像梯度信息，殘差能量及絕對差值排序算子融合到模型中，使得降噪后的圖像可以較完整的保留紋理細(xì)節(jié)，實(shí)現(xiàn)整體均衡降噪，更利于人類視覺觀察。

2 相關(guān)理論知識

Perona和Malik提出的經(jīng)典的各向異性擴(kuò)散模型，其表達(dá)式為

(1)

式中：I(x，y，t)為輸入圖像；div為散度算子；?為梯度；c(|?I|)稱為擴(kuò)散系數(shù)函數(shù)，其表達(dá)式為

(2)

其中，|?I|為梯度幅值，k為梯度閾值。PM模型主要利用圖像梯度模值，通過設(shè)定擴(kuò)散系數(shù)來調(diào)控方程的擴(kuò)散效果，使擴(kuò)散主要作用于圖像的非邊緣區(qū)。盡管該思想保留了圖像的邊緣，但其在去噪過程中仍會將一些紋理和細(xì)節(jié)濾除。降噪方法都是假設(shè)噪聲是振蕩的，圖像是光滑的或分段光滑的，因此，試圖將光滑部分和震蕩部分分開，即將噪聲與圖像紋理細(xì)節(jié)分開。但是圖像的紋理和細(xì)節(jié)具有類似于噪聲的振蕩特性，因此在消除噪聲的同時(shí)，不可避免的損害圖像紋理細(xì)節(jié)信息。而紋理細(xì)節(jié)體現(xiàn)了圖像的豐富信息，其在圖像檢測和圖像內(nèi)容理解上起著舉足輕重的作用，尤其對于醫(yī)學(xué)影像疾病診斷、工業(yè)缺陷識別的作用影響巨大，因此基于PM的降噪算法對紋理部分的處理顯得尤為重要。

2.1 基于PM的紋理檢測算子

傳統(tǒng)的梯度算子對紋理部分的提取顯得比較困難，文獻(xiàn)[17]提出了基于PM算法本身的紋理檢測算子，利用紋理區(qū)含有較大的殘留局部能量，可以很好的區(qū)分圖像邊緣、紋理和細(xì)節(jié)等信息。其基本思想是：在PM模型中，如果圖像的某些紋理信息被錯(cuò)誤的去除，那么它將出現(xiàn)在殘留的局部能量中，可通過計(jì)算局部殘差能量得到圖像紋理算子。此外，它僅檢測PM算法消除的紋理區(qū)域，而不檢測PM算法保留下的紋理區(qū)域。PM局部殘差I(lǐng)r定義

Ir(t)=I0-I(t)=Is+In

(3)

其中，I0為含噪原圖像，I(t)為經(jīng)典PM濾波后圖像，Is包含一些信息成分，如紋理和細(xì)節(jié)，In為噪聲成分。圖像信息和噪聲是不相關(guān)的。因此局部殘差能量可表示為

Pr=Ps+Pn

其中，Pr，Ps，Pn分別為Ir，Is，In的局部能量(局部方差)，由于Pn趨于恒定，因此可以將其視為圖像Pr的強(qiáng)度偏移。在紋理區(qū)域中，Ps的取值比其它區(qū)域大，即紋理區(qū)域的Pr值大于其它區(qū)域的Pr值。局部能量(紋理檢測算子)Pr的計(jì)算步驟如下：

1)用式(1)的PM算法分離出PM殘差，即噪聲和相應(yīng)的紋理細(xì)節(jié)信息，直到殘差能量大于設(shè)定值時(shí)停止迭代，即滿足下述公式

(4)

其中，μ(·)是均值算子，I0為含噪原圖像，I(t)為PM迭代t次濾波后圖像。迭代次數(shù)取決于殘差信息。

2)計(jì)算此時(shí)的局部殘差能量，即紋理檢測算子

(5)

其中，ωx，y(x，y)是歸一化高斯窗，Ωs表示窗口區(qū)域，局部殘差能量Pr(x，y)體現(xiàn)了紋理信息。

文獻(xiàn)[17]的擴(kuò)散系數(shù)是將圖像的邊緣檢測算子(圖像梯度幅值)和上一步計(jì)算出的紋理檢測算子(殘留局部能量)之和構(gòu)成的遞減函數(shù)，共同控制擴(kuò)散程度。即：將下式應(yīng)用于含噪圖像進(jìn)行二次PM迭代，達(dá)到降噪效果：

(6)

下圖顯示了PM紋理檢測算子Pr(x，y)的有效性，如下圖1、圖2所示，lena含噪圖像和腹部低劑量含噪圖像紋理區(qū)的亮度均高于其它區(qū)域，即紋理細(xì)節(jié)處的局部殘差能量較大，因此，基于PM的殘差局部能量可以作為紋理檢測器。該算法可以從經(jīng)典的PM算法濾波后獲得的殘差信息中提取圖像的紋理信息，并在PM的第二次運(yùn)行過程中可以將殘差中的紋理信息提取并保留在原圖。如下圖1、圖2中的(c)圖，即殘差信息圖，只有噪聲，沒有圖像的結(jié)構(gòu)信息，可見算法的有效性，但同時(shí)，從殘差圖像中也發(fā)現(xiàn)在紋理細(xì)節(jié)較多的地方，圖像的噪聲去除程度也減弱，即對紋理區(qū)的降噪效果不佳，整體去噪不夠均勻。

圖1 lena圖像的紋理檢測及殘差結(jié)果

圖2 腹部低劑量含噪圖的紋理檢測及殘差結(jié)果

2.2 絕對差值排序算子

為進(jìn)一步解決紋理區(qū)的噪聲問題，本文引入一種局部圖像統(tǒng)計(jì)算子：絕對差值排序算子(rank-ordered absolute differences，ROAD)[18]，ROAD可以很好地檢測圖像噪聲與細(xì)節(jié)，原文實(shí)驗(yàn)證明，圖像噪聲點(diǎn)像素的ROAD值較大，而沒有被噪聲污染的像素點(diǎn)的ROAD值較小。ROAD的基本原理如下：設(shè)中心像素點(diǎn)的坐標(biāo)為a=(a1，a2)，則：

Ωa(N)：={a+(i，j)：-N≤i，j≤N}

(7)

da，b=|Ia-Ib|

(8)

然后，將局部da，b中值按升序排列，并定義

(9)

r(a)=sort(da，b)

(10)

ROAD值的計(jì)算過程如圖3所示。

圖3 ROAD的計(jì)算過程

較小的四個(gè)絕對差值排序分別為：r1=10，r2=21，r3=33，r4=40。

ROAD的計(jì)算結(jié)果：

ROAD主要通過中心像素與其周圍四個(gè)像素絕對差值的大小，來衡量它們之間的相似性。一般而言，圖像的內(nèi)部及邊緣均具有較好的連續(xù)性，因此在中心像素的八個(gè)鄰域中，至少應(yīng)有四個(gè)像素間是相似的，即具有較小的ROAD值，而噪聲一般會使中心像素與周邊大多數(shù)鄰域像素具有較大的強(qiáng)度差，即具有較大的ROAD值，因此，ROAD值可很好地區(qū)分圖像的噪聲與細(xì)節(jié)。

3 本文算法

PM的梯度幅值算子可以很好的檢測邊緣，保留邊緣信息，但會除掉一些紋理和細(xì)節(jié)信息，這主要是因?yàn)镻DE的梯度相關(guān)算法只能在邊緣(梯度變化較大的區(qū)域)和其它區(qū)域之間產(chǎn)生差異，紋理和細(xì)節(jié)具有類似于噪聲的振蕩特性，因此其在去除噪聲的過程中也會被消除一部分。因此文獻(xiàn)[17]在擴(kuò)散系數(shù)中，除了考慮梯度以外，還結(jié)合紋理檢測算子，有效地提取出紋理信息，使紋理區(qū)的重要信息都被保留，但其在紋理區(qū)的降噪效果明顯弱于邊緣及背景區(qū)。主要是PM濾波過程中濾除了大部分紋理、邊緣及噪聲，沒有很好地區(qū)分紋理區(qū)中噪聲的干擾，使得PM新算法對紋理區(qū)的保護(hù)過度，降噪效果又不夠理想。受絕對差值排序檢測算子可以很好地區(qū)分噪聲和紋理信息的啟發(fā)，本文算法結(jié)合ROAD算子，將絕對差值排序算子引入到紋理檢測算子獲得的紋理區(qū)中，以進(jìn)一步區(qū)分紋理區(qū)的噪聲與細(xì)節(jié)。則基于PM的自適應(yīng)擴(kuò)散系數(shù)新模型降噪公式可重寫為如下

(11)

其中擴(kuò)散系數(shù)為

(12)

(13)

(14)

i=1，2，3，4，分別為上、下、左、右四個(gè)鄰域方向的梯度，在每次迭代過程中，擴(kuò)散系數(shù)都會有更新。即

(15)

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

為了驗(yàn)證改進(jìn)算法的有效性，下面將從主觀和客觀兩個(gè)方面進(jìn)行分析。本文所提算法主要是針對高斯白噪聲在σ=50的情況下進(jìn)行降噪，同時(shí)對比算法的參數(shù)設(shè)置相同，客觀評價(jià)指標(biāo)主要用到信噪比(Signal to Noise Ratio， SNR)、均方根誤差(Root Mean Squared Error， RMSE)、通用質(zhì)量指標(biāo)[19](Universal quality index ，UQI)來評價(jià)處理后圖像的質(zhì)量。SNR值越大，RMSE值越小表示圖像降噪效果越好；通用質(zhì)量指標(biāo)(UQI)是一種客觀質(zhì)量指標(biāo)，它是通過對噪聲圖像進(jìn)行建模，將相關(guān)性，亮度和對比度失真損失綜合起來進(jìn)行設(shè)計(jì)的，UQI值越大，降噪效果越好。

(16)

(17)

本文將提出的新算法應(yīng)用于普通的lena圖像、barbara圖像和醫(yī)學(xué)胸腔體模的真實(shí)低劑量CT掃描圖像進(jìn)行有效性驗(yàn)證。分別采用PM模型、文獻(xiàn)[16]模型、文獻(xiàn)[17]模型和本文改進(jìn)模型對含噪圖像進(jìn)行濾波處理，圖4、5、6的(c)、(d)、(e)、(f)分別是不同算法去噪后的圖像， (g)、(h)、(i)、(j)分別是含噪圖像與對應(yīng)算法降噪后的差值圖，差值圖含有的噪聲越多，邊緣、紋理細(xì)節(jié)信息越少，則表明去噪效果越好。

分別將圖4 、圖5、圖6的(c)、(d)、(e)和(f)與圖4、圖5、圖6的(a)作對比，可以明顯看出：與原始圖像相比，PM模型處理后的圖像細(xì)節(jié)丟失較多，對應(yīng)的殘差圖中可明顯看到殘留的邊緣、紋理等結(jié)構(gòu)信息。文獻(xiàn)[16]、文獻(xiàn)[17]去噪后的圖像細(xì)節(jié)保留較多，改善了PM算法對紋理細(xì)節(jié)的濾除缺陷，但根據(jù)其差值圖像圖4、圖5、圖6對應(yīng)的(h)、(i)可看出，對應(yīng)的紋理區(qū)域的降噪效果較弱，整體降噪效果不均勻，如lena圖像在帽子的羽毛處，barbarat圖像圍脖、褲子等紋理部分，腹部低劑量圖的結(jié)構(gòu)部分這些紋理細(xì)節(jié)較多的區(qū)域，其降噪效果受到一定限制；本文算法處理后的圖像具有較好的視覺效果，圖像噪聲得到有效抑制，且邊緣、紋理等細(xì)節(jié)得到很好地保留，整體去噪效果較均衡，如圖(j)所示，各部分去除噪聲較為均勻，含有很少的結(jié)構(gòu)信息。從表1的客觀指標(biāo)上也可以看出，在同樣的噪聲條件下，本文模型得到的SNR(第一行), RMSE(第二行),UQI(第三行)優(yōu)于其它模型，因此，本文提出的模型具有更有效的去噪效果和紋理細(xì)節(jié)保護(hù)能力。

圖6 腹部低劑量圖像不同算法處理結(jié)果對比

表1 不同降噪模型處理后的質(zhì)量評估參數(shù)對比

5 結(jié)束語

本文以傳統(tǒng)的PM模型為基礎(chǔ)，結(jié)合局部殘差的紋理檢測算子和絕對差值排序算子，提出了新的基于PM模型的圖像降噪算法。基于局部殘差的紋理檢測算子可以很好地檢測到紋理信息，改善了PM算法本身對紋理信息錯(cuò)誤濾除的缺陷，同時(shí)結(jié)合絕對差值排序算子的噪聲與細(xì)節(jié)鑒別能力，很好地克服了PM算法的去噪和紋理保持的矛盾。所提算法很好地保留了紋理、邊緣、和弱細(xì)節(jié)等部分。為了更客觀地驗(yàn)證算法的有效性，文章的SNR、RMSE、UQI等參數(shù)指標(biāo)也有所提高。