張夢(mèng)夢(mèng), 張涇周, 周三平, 張永濤

(1. 西北工業(yè)大學(xué)自動(dòng)化學(xué)院, 陜西 西安 710072;2. 西安交通大學(xué)人工智能與機(jī)器人研究所, 陜西 西安 710049)

?

基于局部熵的邊界與區(qū)域水平集圖像分割模型

張夢(mèng)夢(mèng)1, 張涇周1, 周三平2, 張永濤1

(1. 西北工業(yè)大學(xué)自動(dòng)化學(xué)院, 陜西 西安 710072;2. 西安交通大學(xué)人工智能與機(jī)器人研究所, 陜西 西安 710049)

基于圖像局部熵提出了一種改進(jìn)的結(jié)合邊界信息和區(qū)域信息的水平集圖像分割模型。利用局部熵構(gòu)造自適應(yīng)權(quán)重系數(shù),使其能夠根據(jù)圖像性質(zhì)自適應(yīng)的決定演化方向,準(zhǔn)確引導(dǎo)演化曲線向目標(biāo)方向移動(dòng);然后,根據(jù)自適應(yīng)權(quán)重系數(shù)定義新的邊界指示函數(shù),提高了模型檢測(cè)弱邊界能力,加快了曲線的演化速度;引入Chan-Vese (C-V)模型作為外部能量項(xiàng),提高了模型的抗噪性,增強(qiáng)了模型分割灰度不均勻圖像的能力。通過(guò)圖像分割實(shí)驗(yàn),驗(yàn)證模型對(duì)初始輪廓以及噪聲的魯棒性、分割灰度不均勻圖像的能力,并采用客觀數(shù)值指標(biāo),將所提模型與另外三種模型在分割效率和分割準(zhǔn)確性方面進(jìn)行比較。結(jié)果表明,提出的模型增強(qiáng)了對(duì)噪聲的魯棒性,提高了分割弱邊界圖像的能力,而且分割灰度不均勻的圖像時(shí)也取得了比較滿意的效果。

圖像分割; 水平集方法; 局部熵; 邊界信息; 區(qū)域信息

0 引 言

圖像分割是一個(gè)根據(jù)區(qū)域內(nèi)相似性以及區(qū)域間的相異性而把圖像分割成若干區(qū)域的過(guò)程[1]。圖像分割技術(shù)是圖像處理、圖像分析、圖像理解、圖像識(shí)別和計(jì)算機(jī)視覺(jué)領(lǐng)域一項(xiàng)基本而又關(guān)鍵的技術(shù)[2-3]。水平集方法的演化實(shí)質(zhì)就是偏微分方程的求解,由圖像信息定義能量泛函,根據(jù)梯度下降流方法得到一類演化方程,然后對(duì)該方程進(jìn)行求解。對(duì)于拓?fù)渥兓瘡?fù)雜的圖像有較強(qiáng)的適應(yīng)性并有穩(wěn)定的數(shù)值解。根據(jù)約束條件的不同,我們將水平集分割方法劃分為兩大類:基于邊界信息[4-6]和基于區(qū)域信息[7-14]的分割模型。

基于邊界信息的模型是根據(jù)圖梯度信息來(lái)約束演化曲線,將梯度變化大的區(qū)域看作圖像邊界。這類模型對(duì)于邊界明確的圖像有很好的分割效果,但對(duì)于弱邊界圖像分割效果不理想。例如,測(cè)地活動(dòng)輪廓(geodesic active contour, GAC)模型[4]以及文獻(xiàn)[5]提出的無(wú)需重新初始化的活動(dòng)輪廓模型。在文獻(xiàn)[5]模型的基礎(chǔ)上,文獻(xiàn)[6]利用圖像梯度信息建立了自適應(yīng)權(quán)重系數(shù)模型。由于圖像梯度受噪聲的影響較大,基于邊界信息的模型對(duì)噪聲有很強(qiáng)的敏感性?；趨^(qū)域信息的模型利用圖像灰度的統(tǒng)計(jì)信息來(lái)檢測(cè)目標(biāo)邊界,這類模型不但對(duì)噪聲具有很強(qiáng)的魯棒性,而且可以分割弱邊界的目標(biāo)。例如,文獻(xiàn)[7]根據(jù)圖像的全局灰度信息提出了著名的C-V模型,分割弱邊界和強(qiáng)噪聲圖像時(shí)取得了很好的效果。為了分割灰度不均勻圖像,文獻(xiàn)[5]利用圖像局部灰度信息建立了局部二值擬合(local binary fitting, LBF)模型[8],該模型對(duì)于分割灰度不均勻的圖像取得了極大的成功。然而,局部灰度信息的引入,不但使得模型對(duì)輪廓初始化比較敏感,而且極易導(dǎo)致演化曲線陷入局部最優(yōu)解。

本文提出了一種基于邊界信息與區(qū)域信息的水平集圖像分割模型。首先利用圖像信息熵建立了權(quán)重系數(shù)模型與邊界指示函數(shù)模型,使得模型能夠自適應(yīng)地引導(dǎo)曲線演化方向,并且準(zhǔn)確地停止在目標(biāo)邊界上。其次,引入C-V模型作為外部能量項(xiàng),增強(qiáng)了模型對(duì)弱邊界圖像的處理能力。最后,實(shí)驗(yàn)結(jié)果表明本文模型不但對(duì)噪聲具有很強(qiáng)的魯棒性,而且分割弱邊界與灰度不均勻的圖像都能取得了令人滿意的效果。

1 相關(guān)背景

1.1 自適應(yīng)距離保持水平集演化模型

在文獻(xiàn)[5]提出的模型中,v是一個(gè)常數(shù),缺乏方向和大小的自適應(yīng)性。為解決這個(gè)問(wèn)題, 文獻(xiàn)[6]提出了一種自適應(yīng)的演化模型,該模型的能量泛函可以表示為

(1)

(2)

式中,g為邊界指示函數(shù)；H(z),δ(z)分別為Heaviside函數(shù)和Dirac函數(shù);λ和u分別為長(zhǎng)度項(xiàng)和距離懲罰項(xiàng)的系數(shù)。長(zhǎng)度項(xiàng)和面積項(xiàng)能約束演化曲線向目標(biāo)邊界運(yùn)動(dòng),距離懲罰項(xiàng)能夠自動(dòng)補(bǔ)償水平集函數(shù)與符號(hào)距離函數(shù)之間的偏差,避免了迭代過(guò)程中的重新初始化。v(I)為自適應(yīng)權(quán)重系數(shù),Iσ表示經(jīng)標(biāo)準(zhǔn)差為σ的高斯濾波器平滑后的圖像,由圖像二階導(dǎo)數(shù)的符號(hào)決定演化方向。

雖然該模型能夠自適應(yīng)的調(diào)整演化方向,但是,自適應(yīng)參數(shù)的大小僅依賴于圖像梯度,這仍然存在一些缺點(diǎn):

(1) 非邊界噪聲點(diǎn)處圖像梯度也很大,甚至超過(guò)邊界處梯度的大小,因此僅依據(jù)梯度信息判斷圖像邊界,不能處理噪聲圖像；

(2) 對(duì)灰度不均勻,以及弱邊界圖像處理效果不好。

1.2 C-V模型

Chan和Vese在Mumford-Shah (M-S)模型的基礎(chǔ)上提出了C-V模型,能量泛函定義如下:

λ1∫Ω(I(x,y)-c1)2H(φ)dxdy+

λ2∫Ω(I(x,y)-c2)2(1-H(φ))dxdy

(3)

式中,H(z)，δ(z)分別為Heaviside函數(shù)和Dirac函數(shù);u為長(zhǎng)度項(xiàng)系數(shù);v為面積項(xiàng)系數(shù);λ1,λ2為權(quán)重系數(shù),通常取λ1=λ2=1;c1,c2分別為演化曲線內(nèi)外的灰度均值。當(dāng)且僅當(dāng)演化曲線到達(dá)目標(biāo)邊界時(shí),能量泛函有最小值。C-V模型對(duì)噪聲具有魯棒性,但不能分割灰度不均勻的圖像。

1.3 局部熵的定義

信息熵是圖像信息論中用來(lái)度量信息量的一個(gè)概念,Pun提出圖像熵的定義如下:

(4)

式中,對(duì)數(shù)函數(shù)一般取2作為底數(shù),也可以取其他數(shù)底。

定理 1 對(duì)于圖像的一個(gè)局部窗口,當(dāng)且僅當(dāng)局部窗口內(nèi)的所有灰度值均相等(不為0)時(shí),窗口的熵值最大。

由以上定理可知,在圖像的平坦區(qū)域局部熵值較大,圖像的邊界區(qū)域局部熵值較小。

2 本文模型

本文提出的模型將邊界信息與區(qū)域信息結(jié)合起來(lái),能量泛函為

E(φ)=wELi+(1-w)ECV(φ)+ER(φ)

(5)

式中,第一項(xiàng)為邊界能量項(xiàng);第二項(xiàng)為CV外部能量項(xiàng);第三項(xiàng)為正則項(xiàng)。正則項(xiàng)又包括長(zhǎng)度約束項(xiàng)以及距離懲罰項(xiàng)。

2.1 自適應(yīng)權(quán)重系數(shù)

本文定義的自適應(yīng)權(quán)重系數(shù)為

(6)

式中,Iσ表示經(jīng)標(biāo)準(zhǔn)差為σ的高斯濾波器平滑后的圖像;m>0為局部熵的權(quán)系數(shù),當(dāng)圖像邊界弱,噪聲比較大時(shí),可適當(dāng)?shù)奶岣適的值,一般取值在(6～15)之間。E為圖像的局部熵,增強(qiáng)了邊界處響應(yīng)。

v(I)的選擇主要基于以下兩方面的考慮。

(1)v(I)的方向:圖像的二階導(dǎo)數(shù)在目標(biāo)邊界處異號(hào),演化曲線在目標(biāo)邊界外部時(shí),v>0,曲線向內(nèi)朝目標(biāo)邊界運(yùn)動(dòng);演化曲線在目標(biāo)邊界內(nèi)部時(shí)v<0,與之前方向相反,曲線向外朝目標(biāo)方向運(yùn)動(dòng),達(dá)到自適應(yīng)的效果。

(2)v(I)的大小是由圖像局部熵和梯度共同決定的。在圖像平坦區(qū)域,局部熵大,梯度值小,v的值較小,加快曲線演化速度;在目標(biāo)邊界處,局部熵小,梯度值大,v的值較大,增強(qiáng)了零水平集檢測(cè)弱邊界及多層輪廓的能力;在圖像噪聲點(diǎn)處,局部熵大,梯度值大,v的值較小,曲線繼續(xù)演化,極大地提高了模型的抗噪性。

2.2 邊界指示函數(shù)

根據(jù)圖像局部熵信息,定義如下的邊界指示函數(shù):

(7)

式中,m為局部熵的權(quán)系數(shù);β(v)是v的函數(shù),可以定義為

(8)

將v中大于其平均值的數(shù)記作1,表示邊界;將v中小于平均值的數(shù)記作0,各向同性區(qū)域。在各向同性區(qū)域,邊界指示函數(shù)取值大于0,繼續(xù)演化;在目標(biāo)邊界區(qū)域,邊界指示函數(shù)取值為0,停止演化。

2.3 水平集方法

將自適應(yīng)權(quán)重系數(shù)模型(6)和邊界指示函數(shù)模型(7)帶入到能量函數(shù)(5)中,本文提出的模型總的能量泛函為

E(φ)=wELi+(1-w)ECV(φ)+ER(φ)=

(1-w)[λ1∫Ω(I(x,y)-c1)2H(φ)dxdy+

λ2∫Ω(I(x,y)-c2)2(1-H(φ))dxdy]+

(9)

式中,w=0.5 為邊界能量項(xiàng)與區(qū)域能量項(xiàng)之間的權(quán)系數(shù)。根據(jù)梯度下降流方法,能量函數(shù)(9)對(duì)應(yīng)的演化方程可以表示為

(φ)+(1-w)

[δ(φ)(-λ1(I-c1)2+λ2(I-c2)2)]+

(10)

灰度均值c1,c2可以表示為

(11)

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

通過(guò)自然圖像與人工合成圖像的分割實(shí)驗(yàn),驗(yàn)證模型對(duì)初始輪廓以及噪聲的魯棒性、分割灰度不均勻圖像的能力、分割的效率和分割的準(zhǔn)確性。本文實(shí)驗(yàn)環(huán)境為Matlab2014Ra、Lenovo PC、Pentium IV處理器,主頻3.2 GHz,內(nèi)存2 GB。實(shí)驗(yàn)參數(shù)規(guī)定為:時(shí)間步長(zhǎng)Δt=5,懲罰項(xiàng)系數(shù)u=0.04,長(zhǎng)度項(xiàng)系數(shù)λ=0.2,權(quán)重w=0.5,λ1=λ2=1,高斯核函數(shù)的窗口大小為3×3,標(biāo)準(zhǔn)差σ=2.0,局部熵窗口大小為9×9。

最終算法的實(shí)現(xiàn)可按如下步驟進(jìn)行:

步驟 1 輸入待分割圖像,初始化演化曲線φ0；

步驟 2 如果|length(ci+1)-length(ci)|>η；

步驟 3 按照方程(11)計(jì)算c1(φn),c2(φn)；

步驟 4 分別按照式(6)和式(7)計(jì)算自適權(quán)系數(shù)v(I)以及邊界指示函數(shù)g；

步驟 5 由方程(10)更新水平集函數(shù)；

步驟 6 end；

步驟 7 輸出最終的演化結(jié)果。

實(shí)驗(yàn) 1 驗(yàn)證模型對(duì)初始輪廓的魯棒性

局部熵權(quán)系數(shù)m=6,初始輪廓的位置分別如圖1(a)～圖1(c)所示,圖1(d)～圖1(f)分別為經(jīng)過(guò)2次迭代得到的相應(yīng)分割結(jié)果?？梢钥闯?在初始輪廓的位置、大小都不相同的情況下,本文模型都能夠得到正確的分割結(jié)果。

圖1 模型對(duì)初始輪廓的魯棒性Fig.1 Robustness of our method to different initializations

實(shí)驗(yàn) 2 驗(yàn)證模型對(duì)噪聲的魯棒性

(1) 對(duì)不同噪聲的魯棒性

對(duì)未經(jīng)污染的圖像分別加入不同強(qiáng)度的高斯噪聲和椒鹽噪聲,驗(yàn)證模型對(duì)各種噪聲的抑制能力,實(shí)驗(yàn)結(jié)果如圖2所示。其中,圖2(a)為未受污染的原始圖像的分割結(jié)果,圖2(b)和圖2(c)分別為圖像添加方差0.01和0.1的高斯噪聲后的分割結(jié)果,圖2(c)和圖2(d)分別為圖像添加方差為0.005和0.01的椒鹽噪聲后的分割結(jié)果。由實(shí)驗(yàn)結(jié)果可以看出,本文模型對(duì)不同強(qiáng)度的高斯噪聲以及椒鹽噪聲都有較好的魯棒性。

圖2 對(duì)不同噪聲的魯棒性Fig.2 Robustness to different noisies

(2) 3種模型對(duì)噪聲魯棒性的比較

圖3中m=6,(a1)，(a2)，(a3)為初始輪廓在圖像中的位置；(b1),(b2),(b3)分別為本文模型迭代15次,3次,10次的結(jié)果;(c1),(c2),(c3)為C-V模型迭代100次,120次,65次的結(jié)果;(d1),(d2),(d3)為文獻(xiàn)[6]提出的模型迭代100次的結(jié)果。由圖3以及表1可以看出,本文模型對(duì)噪聲具有極好的魯棒性,而且效率很高,在很少的迭代次數(shù)內(nèi)就能得到準(zhǔn)確的分割效果。C-V模型具有一定的抗噪性,但效率不高,對(duì)于某些噪聲點(diǎn)并不能準(zhǔn)確判斷出來(lái)。而文獻(xiàn)[6]提出的模型,則對(duì)噪聲有很大的敏感性。因此,可以認(rèn)為本文模型在分割效率和噪聲的魯棒性方面優(yōu)于另兩種模型。這是因?yàn)?本文模型采用邊緣與區(qū)域相結(jié)合的方式,根據(jù)圖像邊緣信息可以快速的鎖定圖像邊緣,而區(qū)域信息又增加了模型的抗噪性。

圖3 3種模型對(duì)噪聲魯棒性比較Fig.3 Robustness to noisy comparisons between our model and two other models

圖像名稱本文模型CPU時(shí)間/s迭代次數(shù)C?V模型CPU時(shí)間/s迭代次數(shù)文獻(xiàn)[6]模型CPU時(shí)間/s迭代次數(shù)圖3(a1)0．2340152．63641003．8688100圖3(a2)0．390032．77681206．3336100圖3(a3)0．5460102．8236656．9888100

實(shí)驗(yàn) 3 灰度不均勻圖像的分割

圖4為3種模型對(duì)灰度不均勻圖像項(xiàng)的分割比較,此時(shí)m=10,圖4(b)～圖4(d)分別為3種不同模型的分割結(jié)果。對(duì)分割結(jié)果進(jìn)行視覺(jué)對(duì)比可以看到,本文模型對(duì)灰度不均勻圖像的分割能力較強(qiáng)。

為了從數(shù)值對(duì)本文提出的模型進(jìn)行客觀的評(píng)價(jià),選擇4種圖像分割準(zhǔn)確性評(píng)價(jià)指標(biāo)進(jìn)行驗(yàn)證，分別為真正類(true positive, TP)、假正類(false positive, FP)、精確度(precision, P)、相似度(jaccard similarity, JS)。

圖4 3種模型對(duì)灰度不均勻圖像分割的比較Fig.4 Intensity inhomogeneity image segmentation comparisons between our model and two other models

其中,SG與ST分別為分割的目標(biāo)區(qū)域與實(shí)際的目標(biāo)區(qū)域的集合。實(shí)際目標(biāo)區(qū)域來(lái)源于標(biāo)準(zhǔn)庫(kù)中的分割結(jié)果。

由以上可知一個(gè)好的圖像分割算法應(yīng)該取得較高的TP、P、JS值以及較低的FP值。由表2的數(shù)據(jù)結(jié)果可以看出,本章模型的分割準(zhǔn)確性要高于其他兩種模型。原因在于,首先,本章模型中引入局部熵,在分割過(guò)程中考慮了圖像的局部信息。其次,定義了新的自適應(yīng)權(quán)重系數(shù),可以更加準(zhǔn)確的定位邊緣信息。

表2 3種模型分割準(zhǔn)確性比較

4 結(jié) 論

本文提出了一種基于局部熵的自適應(yīng)分割模型,引入了C-V模型的外部能量項(xiàng),使邊界信息與區(qū)域信息結(jié)合起來(lái),同時(shí)又可以將局部信息與全局信息結(jié)合起來(lái)。該模型可以分割噪聲圖像,對(duì)灰度不均勻圖像也有很好的分割效果?；诰植快氐淖赃m應(yīng)權(quán)重系數(shù),能夠根據(jù)圖像信息自動(dòng)的調(diào)節(jié)演化方向,即使是在噪聲的干擾下,也能夠準(zhǔn)確定位目標(biāo)邊界,極大的提高了模型分割弱邊界以及細(xì)小目標(biāo)的能力。此外,通過(guò)實(shí)驗(yàn)對(duì)模型的性能進(jìn)行了驗(yàn)證,人工合成圖像和自然圖像的分割實(shí)驗(yàn)結(jié)果表明本文提出的模型與文獻(xiàn)[6]提出的模型,以及基于局部灰度信息的活動(dòng)輪廓模型相比,除增強(qiáng)了對(duì)初始輪廓以及噪聲的魯棒性外,還具有較高的分割效率,對(duì)灰度不均勻圖像也有較好的分割效果。

[1] Pong F, Chen Z. A new level set method for inhomogenous image segmentation[J].ImageandVisionComputing,2013,(31):809-822.

[2] Zhang Y J. A survey on evaluation method for image segmentation[J].PatternRecognition, 1996, 29(8):1335-1346.

[3] Liu J W. Research of level set method for image segmentation and its applicationsin medical image[D].Hefei:University of Science and Technology of China,2009.(劉軍偉. 基于水平集的圖像分割方法研究及其在醫(yī)學(xué)圖像中的應(yīng)用[D].合肥:中國(guó)科學(xué)技術(shù)大學(xué),2009.)

[4] Kichenassamy S, Kumar A. Gradient flows and geomerric active contour models[C]∥Proc.oftheIEEEInternationalConferenceonComputerVision, 1995: 810-815.

[5] Li C M, Xu C Y, Gui C F, et al. Level set evolution without re-initialization: a new variational formulation[C]∥Proc.oftheIEEEComputerSocietyConferenceonComputerVisionandPatternRecognition, 2005: 430-436.

[6] He C J, Li M, Zhang Y,Adaptive distance preserving level set evolution for image segmentation[J].JournalofSoftware, 2008, 12:3161-3169.(何傳江,李夢(mèng),詹毅. 用于圖像分割的自適應(yīng)距離保持水平集演化[J]. 軟件學(xué)報(bào), 2008, 12:3161-3169.)

[7] Chan T F, Vese L A. Active contours without edges[J].IEEETrans.onImageProcessing, 2001, 10(2): 266-277.

[8] Li C M, Kao C Y. Minimization of region-scalable fitting energy for image segmentation[J].IEEETrans.onImageProcessing, 2008, 17(10): 1940-1949.

[9] Zhang K, Song H, Zhang L. Active contours driven by local image fitting energy[J].PatternRecognition, 2010, 43(4): 1199-1206.

[10] Li C M, Kao C Y. Minimization of region-scalable fitting energy for image segmentation[J].IEEETrans.onImageProces-sing, 2008, 17(10): 1940-1949.

[11] Zhang K, Zhang L, Song H, et al. Active contours with selective local or global segmentation: a new formulation and level set method [J].ImageandVisionComputing, 2010, 28(4): 668-676.

[12] Wang X F, Huang D S, Xu H. An efficient local Chan-Vese model for image segmentation[J].PatternRecognition, 2010, 43(3): 603-618.

[13] Zhou S, Wang J, Zhang S, et al. Active contour model based on local and global intensity information for medical image segmentation[J].Neurocomputing, 2016,186:107-118

[14] Wang L, Li C. Active contours driven by local and global intensity fitting energy with application to brain MR image segmentation[J].ComputerizedMedicalImagingandGraphics, 2009: 520-531.

Boundary and region level set method based on local entropy for image segmentation

ZHANG Meng-meng1, ZHANG Jing-zhou1, ZHOU San-ping2, ZHANG Yong-tao1

(1.SchoolofAutomation,NorthwesternPolytechnicalUniversity,Xi’an710072,China; 2.InstituteofArtificialIntelligenceandRobotics,Xi’anJiaotongUniversity,Xi’an710049,China)

An improved level set method for image segmentation based on image local entropy information is proposed, which combines the edge-based model and the region-based model into a joint framework. Using image local entropy information, an adaptive weighting function is built firstly, which enables the evolving curve choose the evolution direction and move to the object boundary, adaptively. A novel edge indicator function is proposed based on the weighting function, which improves the ability of detecting weak boundary and accelerates the speed of contour evolution. Finally, the Chan-Vese (C-V) model is introduced into the joint framework as an external energy, which enhances the model dealing images with intensity inhomogeneity. In the experiments, the robustness of the method is evaluated to initial contours and noises, and the ability of segmenting images with intensity inhomogeneity. The results show that the proposed method can not only enhance the robustness to noises and improve the ability of segmenting images with weak boundary, but also achieve the satisfying results in segmenting images with intensity inhomogeneity, as compared with the other three methods using objective numerical indicators.

image segmentation; intensity inhomogeneity; local entropy, level set method; active contour model

2016-05-24;

2016-08-03;網(wǎng)絡(luò)優(yōu)先出版日期:2016-09-30。

TN 957.52

A

10.3969/j.issn.1001-506X.2016.12.30

張夢(mèng)夢(mèng)(1992-),女,碩士研究生,主要研究方向?yàn)閳D形圖像處理、圖像分割。

E-mail:z_mm@mail.nwpu.edu.cn

張涇周(1960-),男,教授,碩士研究生導(dǎo)師,主要研究方向?yàn)樯镝t(yī)學(xué)成像。

E-mail:bme@nwpu.edu.cn

周三平(1988-),男,博士研究生,主要研究方向?yàn)閳D像分割、目標(biāo)檢測(cè)、多目標(biāo)跟蹤。

E-mail:sanpingzhou@stu.xjtu.edu.cn

張永濤(1992-),男,碩士研究生,主要研究方向?yàn)閳D像處理。

E-mail:252361378@qq.com

網(wǎng)絡(luò)優(yōu)先出版地址:http:∥www.cnki.net/kcms/detail/11.2422.TN.20160930.1313.030.html