朱云龍，翁桂榮

?

面積項(xiàng)能量加強(qiáng)的距離規(guī)則水平集演化模型

朱云龍，翁桂榮

(蘇州大學(xué)機(jī)電工程學(xué)院，江蘇 蘇州 215000)

針對(duì)距離規(guī)則水平集演化(DRLSE)模型存在易陷入虛假邊界、對(duì)噪聲敏感、收斂速度慢以及容易從弱邊緣處泄露等不穩(wěn)定問題，提出了面積項(xiàng)能量加強(qiáng)的水平集演化函數(shù)對(duì)水平集方法進(jìn)行改進(jìn)。首先提出了一個(gè)自適應(yīng)邊緣指示函數(shù)，其根據(jù)圖像信息來調(diào)整函數(shù)參數(shù)，從而控制演化速度以及對(duì)噪聲敏感度，使水平集演化更加快速穩(wěn)定。同時(shí)結(jié)合區(qū)域生長(zhǎng)方法，將圖像處理成一個(gè)二值矩陣，并據(jù)此矩陣增加一加強(qiáng)項(xiàng)，使得面積項(xiàng)能量得到加強(qiáng)，令水平集函數(shù)隨著距離目標(biāo)遠(yuǎn)近而自動(dòng)調(diào)整能量大小，降低計(jì)算成本，有效解決對(duì)噪聲敏感、易陷入虛假邊界等問題。為驗(yàn)證模型的有效性，采用多張實(shí)圖進(jìn)行分割實(shí)驗(yàn)并與DRLSE等模型進(jìn)行對(duì)比，實(shí)驗(yàn)結(jié)果表明，提出的模型能有效解決存在問題，有更高的計(jì)算效率和準(zhǔn)確率。

圖像分割；水平集；主動(dòng)輪廓模型；區(qū)域生長(zhǎng)

近二十年來，基于水平集方法的主動(dòng)輪廓模型已成為圖像分割最有效的方法之一[1-3]，并被廣泛應(yīng)用于圖像分割中。主動(dòng)輪廓模型的圖像分割算法是一種結(jié)合上層知識(shí)的圖像處理技術(shù)，其在算法優(yōu)化過程中可以引入圖像灰度信息、邊界信息[4-6]，并允許引入各種先驗(yàn)知識(shí)進(jìn)行魯棒性圖像分割[7-10]，同時(shí)提供光滑的封閉輪廓作為分割結(jié)果，具有形式多樣、結(jié)構(gòu)靈活等優(yōu)越性能。

1988年，OSHER和SETHIAN[11]在分析流體力學(xué)時(shí)，創(chuàng)造性地提出了水平集方法。熱力學(xué)方程中火苗的外形具有高動(dòng)態(tài)性以及隨意的拓?fù)浣Y(jié)構(gòu)變化，用參數(shù)活動(dòng)輪廓模型中的Snake方法描述火苗的外形變化較為復(fù)雜，而在高一維空間用水平集就能夠簡(jiǎn)單地描述。

文獻(xiàn)[12]和[13]分別于1993年及1995年提出了幾何活動(dòng)輪廓(geometric active contour, GAC)模型。此模型成功地將水平集方法應(yīng)用到圖像處理和計(jì)算機(jī)視覺領(lǐng)域。其主要思想是二維平面上的閉合曲線可用高一維的函數(shù)與水平面的交集表達(dá)，并稱之為零水平集，通過求解這個(gè)零水平集就可得到目標(biāo)的邊界。其缺點(diǎn)是初始輪廓曲線必須設(shè)置在目標(biāo)邊界的內(nèi)部或外部，并且在演化的過程中需要周期性地重新初始化水平集函數(shù)[14-15]。盡管重新初始化能夠獲得規(guī)則的水平集函數(shù)，但未必能使零水平集函數(shù)運(yùn)動(dòng)至目標(biāo)邊界，且在實(shí)際操作中何時(shí)重新初始化也是個(gè)關(guān)鍵問題。

LI等[16]針對(duì)水平集演化過程中的重新初始化問題，提出避免水平集重新初始化的距離正則項(xiàng)的概念，進(jìn)而提出了距離規(guī)則水平集演化(distance regularized level set evolution，DRLSE)模型。該模型在原模型的內(nèi)部能量項(xiàng)中添加了距離規(guī)則項(xiàng)，用于控制曲線的演化過程，消除水平集函數(shù)與符號(hào)距離函數(shù)之間的偏差，以避免水平集函數(shù)在演化過程中的震蕩問題，最終達(dá)到在演化過程中無需重新初始化的目標(biāo)[15, 17-19]。這個(gè)模型使得水平集函數(shù)的構(gòu)造變得更加容易，可簡(jiǎn)單地初始化為一個(gè)二值函數(shù)，是傳統(tǒng)水平集方法的重大突破。但該方法的缺陷有對(duì)噪聲敏感、收斂速度慢、容易陷入虛假邊界中以及從弱邊緣處泄露等問題。

針對(duì)DRLSE模型存在的上述缺點(diǎn)，本文提出了一種改進(jìn)的距離規(guī)則水平集模型。以DRLSE模型為基礎(chǔ)，選取一個(gè)自適應(yīng)邊緣指示函數(shù)，通過圖像的標(biāo)準(zhǔn)差來表達(dá)目標(biāo)邊界的清晰程度，并根據(jù)清晰程度自動(dòng)調(diào)整邊緣指示函數(shù)，既有利于加快演化速度，又可提升穩(wěn)定性。同時(shí)，運(yùn)用區(qū)域生長(zhǎng)算法處理圖像，得到一個(gè)二值矩陣，并將此矩陣添加到面積項(xiàng)，使得面積項(xiàng)能量加強(qiáng)，由此令水平集函數(shù)隨著距離目標(biāo)遠(yuǎn)近而自動(dòng)調(diào)整能量大小，有效解決了對(duì)噪聲敏感、易陷入虛假邊界等問題。最后，本文模型與DRLSE模型，自適應(yīng)水平集演化模型(adaptive level set evolution，ALSE)[20]和結(jié)合Log邊緣檢測(cè)的區(qū)域可收縮擬合模型(region-scalable fitting and optimized laplacian of gaussian energy, RSFOL)模型[21]進(jìn)行了一系列對(duì)比實(shí)驗(yàn)，驗(yàn)證改進(jìn)模型的有效性。

1 距離規(guī)則水平集介紹

1.1 模型簡(jiǎn)介

為了解決傳統(tǒng)主動(dòng)輪廓模型在演化過程中需要不斷周期性地重復(fù)初始化的問題[15-19]，DRLSE模型引入了一個(gè)新的距離規(guī)則項(xiàng)用于糾正水平集函數(shù)與符號(hào)距離函數(shù)之間的偏差，其做法是在傳統(tǒng)模型中的能量泛函上再添加一項(xiàng)內(nèi)部能量泛函。DRLSE模型的能量泛函定義為

其中，?為梯度算子；為水平集函數(shù)；(?)為勢(shì)阱函數(shù)；(?)為一維Dirac函數(shù)；(?)是Heaviside函數(shù)；(?)為邊緣指示函數(shù)；為常數(shù)。

其能量泛函由兩部分組成：①內(nèi)部能量項(xiàng)，由水平集正則項(xiàng)組成(即())，主要作用是矯正水平集與符號(hào)距離間的誤差，避免曲線在演化過程中重復(fù)初始化；②外部能量項(xiàng)，由長(zhǎng)度項(xiàng)與面積項(xiàng)構(gòu)成。其中L()為長(zhǎng)度項(xiàng)，是用于驅(qū)使零水平集輪廓朝目標(biāo)邊界演化；而A()為面積項(xiàng)，用于加快零水平集的演化速率。

DRSLE采用的雙勢(shì)阱函數(shù)(?)可表達(dá)為

其定量地表示了水平集函數(shù)偏離符號(hào)距離函數(shù)的程度，用來糾正水平集函數(shù)與符號(hào)距離函數(shù)的偏差。

一維Dirac函數(shù)δ(?)和Heaviside函數(shù)H(?)分別表達(dá)為(為參數(shù))

其中，(?)為邊緣指示函數(shù)，可定義為

當(dāng)曲線演化到目標(biāo)邊界時(shí)，式(1)得到最小值。所以，曲線的演化可以轉(zhuǎn)化為求解式(1)能量泛函的最小化。因此，梯度下降流方程的穩(wěn)定狀態(tài)解，可表達(dá)為

在式(5)和(6)中，是標(biāo)準(zhǔn)偏差為的二維高斯濾波算子；*表示卷積。div(?)為向量的散度算子；d(|?|)是一雙勢(shì)阱函數(shù)，即

采用有限差分法對(duì)式(6)進(jìn)行網(wǎng)格離散化，引入時(shí)間步長(zhǎng)Δ，可得離散化的水平集函數(shù)方程

傳統(tǒng)的幾何活動(dòng)輪廓模型中，水平集必須反復(fù)初始化為符號(hào)距離函數(shù)。而在DRLSE模型中避免了反復(fù)初始化的情況，且水平集正則項(xiàng)在迭代演化中驅(qū)使變?yōu)榉?hào)距離函數(shù)。為了進(jìn)一步提高計(jì)算效率和更好地保持演化的規(guī)則性，可以定義水平集初始值0為

其中，0為大于零的常數(shù)；0是圖像區(qū)域中的某一部分區(qū)域。

1.2 DRLSE模型分析

分析DRLSE模型，可以發(fā)現(xiàn)2個(gè)問題：

(1) 外部能量項(xiàng)過多且依賴邊緣指示函數(shù)，而邊緣指示函數(shù)依賴的是梯度信息，梯度越大，其值越接近于零。而圖像經(jīng)過平滑濾波后，會(huì)使邊界變得模糊并變寬，因此目標(biāo)邊界附近即邊緣指示函數(shù)趨于零的區(qū)域也同樣變寬。當(dāng)二個(gè)目標(biāo)靠近時(shí)，目標(biāo)邊界可能會(huì)因?yàn)V波算子的平滑而相互連接導(dǎo)致分割失敗。圖1(a)中每個(gè)像素點(diǎn)代表原圖像中對(duì)應(yīng)像素點(diǎn)處的梯度值經(jīng)過邊緣指示函數(shù)處理之后的值。圖中在趨于零的區(qū)域(黑色部分)很粗，且靠近的兩個(gè)物體不易分割，出現(xiàn)相連的情況。圖1(b)中，運(yùn)用了DRLSE模型進(jìn)行分割，得到了失敗的結(jié)果。

圖1 DRLSE模型對(duì)靠近目標(biāo)處理結(jié)果圖

(2) 式(6)中是人為設(shè)定的常數(shù)，其大小嚴(yán)重影響分割的結(jié)果。過小則降低了演化速度，演化容易陷入虛假邊界中；適度增大可增加演化速度，使演化越過虛假邊界得到正確結(jié)果；但是過大也易出現(xiàn)邊界泄漏現(xiàn)象，導(dǎo)致分割失敗。如圖2所示，對(duì)于邊界較為模糊的情況，圖2(a)因偏小而出現(xiàn)了陷入虛假邊界的情況；圖2(c)因偏大出現(xiàn)了邊界泄露的情況。而對(duì)于圖3，目標(biāo)邊界較為清晰，圖3(a)~(c)均分割正確，但分割速度因的變化而改變，其中代表迭代次數(shù)。由于的敏感性大大降低了模型的實(shí)用性。

圖2 DRLSE模型分割效果與a的關(guān)系

圖3 DRLSE模型分割速度與a的關(guān)系(N指迭代次數(shù))

2 改進(jìn)模型

眾所周知，方差可度量數(shù)據(jù)的離散程度。對(duì)于一幅圖像，用方差可以在一定程度上顯示圖像的光滑性與邊界清晰程度。而由上述DRLSE模型分析中可知，邊緣指示函數(shù)是控制演化速度的重要因素，演化速度快有利于避免陷入虛假邊界，但容易產(chǎn)生邊界泄露現(xiàn)象；演化速度慢時(shí)情況正好相反，有利于避免邊界泄露現(xiàn)象而易產(chǎn)生陷入虛假邊界的情況。因此，將方差作為參數(shù)放入邊緣指示函數(shù)，用于針對(duì)不同的圖像而自動(dòng)調(diào)節(jié)演化速度的變化。同時(shí)，對(duì)于實(shí)際圖像的背景區(qū)域，均會(huì)有或多或少的噪聲干擾，為了減少干擾，可以先用區(qū)域生長(zhǎng)法對(duì)圖像進(jìn)行粗分割，將圖像分為目標(biāo)區(qū)域與背景區(qū)域。而式(6)中的右側(cè)第三項(xiàng)控制著演化速度，因此可以根據(jù)粗分割結(jié)果調(diào)整此項(xiàng)，使得演化曲線在背景區(qū)域時(shí)增大演化速度，防止陷入虛假邊界，而進(jìn)入目標(biāo)區(qū)域時(shí)減小演化速度，避免邊界泄露。根據(jù)以上思想，設(shè)計(jì)了一個(gè)關(guān)于演化曲線的能量泛函，對(duì)能量泛函進(jìn)行最小化，得到其梯度下降流方程

其中，為經(jīng)過區(qū)域生長(zhǎng)算法所得到的二值矩陣；為參數(shù)。

與DRLSE模型相比，此模型改進(jìn)了邊緣指示函數(shù)；同時(shí)基于區(qū)域生長(zhǎng)算法增加了一個(gè)加強(qiáng)項(xiàng)()，與原來的面積項(xiàng)結(jié)合，得到了加強(qiáng)的面積項(xiàng)(+)δ()。

2.1 自適應(yīng)邊緣指示函數(shù)

其中，是經(jīng)過二維高斯濾波算子所濾波過的圖像；為標(biāo)準(zhǔn)差；1為邊界閾值函數(shù)，其表達(dá)式為

其中，()圖像的方差，表示圖像的離散程度。方差小，表明數(shù)據(jù)聚集好；方差大，表明數(shù)據(jù)離散度大。1的大小表示目標(biāo)相距平均值的遠(yuǎn)近關(guān)系。當(dāng)圖像目標(biāo)相距平均值較遠(yuǎn)時(shí)，邊界較為清晰，即1較大，(|?|)自動(dòng)增大，可加快演化速度，避免陷入虛假邊界中；當(dāng)圖像目標(biāo)相距平均值較近時(shí)，邊界則模糊，1較小，(|?|)自動(dòng)減小，演化速度降低，不會(huì)出現(xiàn)邊界泄漏現(xiàn)象。圖4為二種邊緣指示函數(shù)的示意圖，其中軸表示|?|，軸表示，藍(lán)色曲線表示DRLSE模型中的邊緣指示函數(shù)；而紅色曲線表示自適應(yīng)邊緣指示函數(shù)。根據(jù)式(13)，不同的圖像有不同的方差，1因此會(huì)隨著不同的圖像有所變化；根據(jù)式(12)，函數(shù)受1的影響，因此，自適應(yīng)函數(shù)也會(huì)隨著不同的圖而自動(dòng)調(diào)整。

圖4 邊緣指示函數(shù)曲線

為了說明改進(jìn)的邊緣指數(shù)函數(shù)的有效性，根據(jù)DRLSE模型，通過改變函數(shù)進(jìn)行驗(yàn)證。圖5為3幅目標(biāo)物體都靠得較近的細(xì)胞圖。從圖5可以看出，DRLSE模型無法準(zhǔn)確地分割出多個(gè)靠近的物體，并易陷入虛假邊界。而使用本文提出的邊緣指示函數(shù)后，分割精度大大提高。實(shí)驗(yàn)參數(shù)為：=0.2，=5，=1.5，=1.5，=2，Δ=5，0=2。

圖5 二種邊緣指示函數(shù)方法對(duì)比分割圖

2.2 基于區(qū)域生長(zhǎng)方法的加強(qiáng)面積項(xiàng)

針對(duì)DRLSE模型的分割效果易受參數(shù)敏感性影響的問題，改進(jìn)的基本思想是通過區(qū)域生長(zhǎng)方法對(duì)圖像進(jìn)行處理，得到一個(gè)二值矩陣。根據(jù)此矩陣增加一個(gè)加強(qiáng)項(xiàng)使得面積項(xiàng)能量得到優(yōu)化，讓其可以根據(jù)演化曲線位置自動(dòng)調(diào)整演化速率的大小，減小對(duì)固定值的敏感性。

區(qū)域生長(zhǎng)的基本思想是將具有相似性質(zhì)的像素集合起來構(gòu)成區(qū)域。首先對(duì)于目標(biāo)區(qū)域，需要找到一個(gè)種子像素作為生長(zhǎng)的起點(diǎn)，隨后根據(jù)事先確定的生長(zhǎng)或相似準(zhǔn)則，搜尋種子點(diǎn)周圍的像素中是否有相同或相似性質(zhì)的像素。若有，則將其包括進(jìn)來，并以此點(diǎn)為新的種子點(diǎn)繼續(xù)向四周搜尋，直到再無滿足條件的像素可以包括進(jìn)來為止，一個(gè)區(qū)域就生長(zhǎng)而成[22]。具體步驟如下：

(3) 在目標(biāo)區(qū)域進(jìn)行區(qū)域生長(zhǎng)。首先對(duì)目標(biāo)區(qū)域進(jìn)行收縮，若目標(biāo)區(qū)域無空洞，則收縮成一個(gè)點(diǎn)，若有空洞則收縮成一個(gè)連通環(huán)，分別將收縮后得到的這些點(diǎn)作為種子點(diǎn)。生長(zhǎng)準(zhǔn)則為相鄰像素的梯度差小于3，以此生成新的二值矩陣；

(4) 以為基礎(chǔ)增加一加強(qiáng)項(xiàng)δ()。

圖6(c)為二值矩陣。當(dāng)演化曲線上的點(diǎn)位于圖中白色區(qū)域時(shí)，矩陣對(duì)應(yīng)的像素點(diǎn)的值為1，則面積項(xiàng)為(+)δ()，其值較大，有利于加快演化速率；當(dāng)演化曲線上的點(diǎn)位于圖中黑色區(qū)域時(shí)，矩陣對(duì)應(yīng)的像素點(diǎn)的值為0，相應(yīng)加強(qiáng)項(xiàng)δ()也同樣為0，此時(shí)面積項(xiàng)為gδ()，面積項(xiàng)的值較小，演化速率變慢。也就是說，當(dāng)水平集距離目標(biāo)區(qū)域較遠(yuǎn)時(shí)，面積項(xiàng)值較大，演化速率提升，有效避免了陷入虛假邊界；當(dāng)水平集演化到距離目標(biāo)區(qū)域較近時(shí)，面積項(xiàng)值變小，演化速率降低，有效避免出現(xiàn)邊界泄漏現(xiàn)象。

圖6以3組加強(qiáng)面積能量的水平集模型的過程分析為示例。通過第2列DRLSE模型分割圖與第5列僅僅加強(qiáng)了面積項(xiàng)之后的分割圖(邊緣指示函數(shù)未變換，依然為原DRLSE模型中定義的)，結(jié)果說明改進(jìn)的加強(qiáng)面積項(xiàng)對(duì)分割有著更好的效果。加強(qiáng)面積能量的模型參數(shù)：=0.2，=5，=1.5，=1.5，=1.5，Δ=5，0=2，=2。

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

本文提出的模型將通過實(shí)圖進(jìn)行驗(yàn)證。實(shí)驗(yàn)采用Matlab R2015b軟件，電腦處理器參數(shù)為：2.6 GHz Inter(R) Core(TM) i7-4720HQ。

將自適應(yīng)邊緣指示函數(shù)與基于區(qū)域增長(zhǎng)方法的加強(qiáng)面積項(xiàng)通過改進(jìn)結(jié)合為一個(gè)新模型，如式(11)。模型算法如下：

(1) 設(shè)置參數(shù)：=0.2，=5，=1.5，=1.5，=1.5，Δ=5，0=2，=0.2。內(nèi)迭代次數(shù)為20，每計(jì)算演化20次記錄一次分割結(jié)果。外迭代次數(shù)為50，即最多允許50次外迭代；

圖6 加強(qiáng)面積能量的水平集模型分割過程圖

(2) 求經(jīng)過高斯平滑后圖像的像素平均值，并以此平均值為閾值對(duì)圖像進(jìn)行分割；

(3) 對(duì)于閾值分割后的圖像進(jìn)行區(qū)域生長(zhǎng)，得到新的二值矩陣；

(4) 根據(jù)圖像信息計(jì)算邊界閾值函數(shù)1()，邊緣指示函數(shù)(|?|)；

(5) 根據(jù)式(11)進(jìn)行內(nèi)迭代計(jì)算；

(6) 若當(dāng)前與上次記錄的分割結(jié)果的面積相對(duì)誤差小于1%，則停止演化，輸出分割結(jié)果；否則返回到步驟(5)；

(7) 若一直不滿足步驟(6)中的條件，當(dāng)?shù)螖?shù)達(dá)到設(shè)置的最大值時(shí)停止，輸出分割結(jié)果。

圖7為DRLSE模型與本文模型的分割對(duì)比圖，其中藍(lán)色曲線表示初始框，綠色曲線表示為迭代過程，紅色曲線為分割結(jié)果。第1~4行為背景噪聲較為明顯的單目標(biāo)分割圖，邊界較為模糊，DRLSE模型則易陷入虛假邊界，使演化停止而本文模型則能夠根據(jù)圖像調(diào)整參數(shù)從而達(dá)到良好的效果。第5~8行為靠近的多目標(biāo)分割圖，對(duì)于背景噪聲較大的，易陷于外部噪聲處，而對(duì)于背景噪聲較小的，演化曲線易陷入兩目標(biāo)之間的凹陷處，而本文模型依然得到良好的分割效果。

為了對(duì)分割效果進(jìn)行定量分析，采用Dice相似性系數(shù)(dice similarity coefficient，DSC)來評(píng) 判[21]，即

其中，1和2分別代表得到的分割和真實(shí)的分割區(qū)域的面積。

為了進(jìn)一步說明本文算法的有效性與實(shí)用性，從Weizmann數(shù)據(jù)庫[23]收集了如下6組實(shí)景圖(圖8)對(duì)本文模型繼續(xù)進(jìn)行驗(yàn)證，并將其分別與DRLSE模型、ALSE模型、RSFOL模型進(jìn)行對(duì)比。初始曲線均距邊框6個(gè)像素。表1為4種模型的DSC值與迭代時(shí)間比較。較DRLSE模型、RSFOL模型，本文模型不僅在時(shí)間上更具優(yōu)勢(shì)(除了(c)列第4行的DRLSE模型很快進(jìn)入虛假邊界，用時(shí)較短的情況)，而且有著更高的分割精度。與ALSE模型相比，本文模型在分割速度上不具優(yōu)勢(shì)，但分割效果更佳。其中本文模型對(duì)于第5行DSC值較低的原因是因?yàn)槟繕?biāo)較小，使式(14)分母較小，小的偏差會(huì)對(duì)其值影響較大。

圖7 DRLSE模型與本文模型分割對(duì)比圖(1~4行為單目標(biāo)分割圖，5~8行為多目標(biāo)分割圖)

圖8 景物圖分割效果對(duì)比

表1 4種模型的DSC值與迭代時(shí)間比較(s)

4 結(jié) 論

提出了利用面積項(xiàng)能量加強(qiáng)的水平集演化函數(shù)對(duì)水平集方法進(jìn)行改進(jìn)。在DRLSE模型基礎(chǔ)上，提出了一個(gè)新穎的邊緣指示函數(shù)，并且結(jié)合區(qū)域生長(zhǎng)的方法，構(gòu)成一個(gè)加強(qiáng)面積項(xiàng)。利用該模型與DRLSE模型進(jìn)行了一系列實(shí)驗(yàn)對(duì)比。實(shí)驗(yàn)結(jié)果顯示，該模型無論在分割效果還是在分割速度上，較DRLSE模型都有較大的優(yōu)越性，對(duì)于噪聲較大、較為模糊的圖像，以及多目標(biāo)且較為靠近的圖像，也擁有良好地分割效果。較好地解決了DRLSE模型對(duì)噪聲敏感、易陷入虛假邊界、收斂速度慢以及容易從弱邊緣處泄露等不穩(wěn)定問題，魯棒性較強(qiáng)。同時(shí)，該模型也有一些待改進(jìn)的地方，例如對(duì)于初始邊界較為敏感，初始演化曲線只能完全包含目標(biāo)物體或完全在目標(biāo)曲線內(nèi)部，才有可能分割正確，一旦初始輪廓與目標(biāo)邊界相交，則分割失敗，這也是日后需要基于該模型繼續(xù)努力改進(jìn)的方向。

[1] 郭笑妍, 梅雪, 李振華, 等. 形狀約束下活動(dòng)輪廓模型冠脈血管圖像多尺度分割[J]. 中國(guó)圖象圖形學(xué)報(bào), 2016, 21(7): 923-932.

[2] 王琦, 潘振寬, 魏偉波, 等. 隱式曲面兩相圖像分割的變分水平集模型及對(duì)偶方法[J]. 電子學(xué)報(bào), 2011, 39(1): 207-212.

[3] 李亞峰. 基于圖像分解的稀疏正則化多區(qū)域圖像分割方法[J]. 電子學(xué)報(bào), 2015, 43(9): 1841-1849.

[4] DAI L Z, DING J D, YANG J. Inhomogeneity-embedded active contour fornatural image segmentation [J]. Pattern Recognition, 2015, 48(8): 2513-2529.

[5] 陳志彬, 邱天爽, SU R. 一種基于FCM和Level Set的MRI醫(yī)學(xué)圖像分割方法[J]. 電子學(xué)報(bào), 2008, 36(9): 1734-1736.

[6] CHEN Y M, TAGARE H D, THIRUVENKADAM S, et al. Using prior shapes in geometric active contours in a variational framework [J]. International Journal of Computer Vision, 2002, 50(3): 315-328.

[7] LI Q, DENG T Q, XIE W. Active contours driven by divergence of gradient vector flow [J]. Signal Processing, 2016, 120: 185-199.

[8] CHAN T F, VESE L A. Active contours without edges [J]. IEEE Transactions on Image Processing, 2001, 10(2): 266-277.

[9] 郝智慧, 郭滿才, 宋楊楊. 結(jié)合全局與局部信息活動(dòng)輪廓的非同質(zhì)圖像分割[J]. 中國(guó)圖象圖形學(xué)報(bào), 2016, 21(7): 886-892.

[10] 姜慧研, 馮銳杰. 基于改進(jìn)的變分水平集和區(qū)域生長(zhǎng)的圖像分割方法的研究[J]. 電子學(xué)報(bào), 2012, 40(8): 1659-1664.

[11] OSHER S, SETHIAN J A. Fronts propagating with curvature dependent speed: algorithms based on Hamilton Jacobi formulations [J]. Journal of Computational Physics, 1988, 79(1): 12-49.

[12] CASELLES V, CATTé F, COLL T, et al. A geometric model for active contours in image processing [J]. Numerische Mathematik, 1993, 66(1): 1-31.

[13] MALLADI R, SETHIAN J A, VEMURI B C. Shape modeling with front propagation: A level set approach [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1995, 17(2): 158-175.

[14] 王艷. 圖像分割的自適應(yīng)快速水平集演化模型[J]. 計(jì)算機(jī)仿真, 2016, 33(2): 269-272.

[15] CASELLES V, KIMMEL R, SAPIRO G. Geodesic active contours [J]. International Journal of Computer Vision, 1997, 22(1): 61-79.

[16] LI C M, XU C, GUI C, et al. Level set evolution without re-initialization: a new variational formulation [C]//2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR05). New Youk: IEEE Prees, 2005: 430-436.

[17] 李夢(mèng). 圖像分割的變分水平集模型研究[D]. 重慶: 重慶大學(xué), 2011.

[18] 何傳江, 李夢(mèng), 詹毅. 用于圖像分割的自適應(yīng)距離保持水平集演化[J]. 軟件學(xué)報(bào), 2008, 19(12): 3161-3169.

[19] HAN B, WU Y Q. A novel active contour model based on modified symmetric cross entropy for remote sensing rover image segmentation [J]. Pattern Recognition, 2017, 67: 396-409.

[20] WANG Y, HE C J. An adaptive level set evolution equation for contour extraction [J]. Applied Mathematics and Computation, 2013, 219: 11420-11429.

[21] DING K Y, XIAO L F, WENG G R. Active contours driven by region-scalable ?tting and optimized Laplacian of Gaussian energy for image segmentation [J]. Signal Processing, 2017, 134: 224-233.

[22] GONZALEZ R C, WOODS R E, EDDINS S L. Digital image processing using MATLAM [M]. Beijing: Publishing House of Electronics Industry, 2005: 308-310.

[23] ALPERT R B S, GALUN M, BRANDT A, et al. Image segmentation by probabilistic bottom-up aggregation and cue integration [J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2012, 34(2): 315-327.

Improved Distance Regularized Level Set Evolution Model by Enhancing Energy of Area Term

ZHUYunlong, WENGGuirong

(College of Mechanical and Electronic Engineering, Soochow University, Suzhou Jiangsu 215000, China)

The distance regularized level set evolution (DRLSE) model has lots of shortcomings when applied in image segmentation. It is easy to stuck around false boundaries, is sensitive to noise, have slow convergence speed and may not detect weak edges. To solve these problems, this paper present a level set function whose area term is enhanced. Firstly, an adaptive edge indicator function is proposed, which can adjust some parameters automatically based on image information, to control evolution speed and the sensitivity to noise. Furthermore, combining the model with region growing method and processing the image into a binary matrix, a reinforcement term is added based on the binary matrix to intensify energy of the area term. So the energy functional could be automatically adjusted according to the distance between evolution curves and targets. This reduces the computational cost, makes the model insensitive to noise and isn’t easy to fall into false boundaries. Experiments on some synthetic and real images prove that the proposed model is not only robust, but also achieves higher segmentation accuracy and efficiency.

image segmentation; level set; active contour; region growing

TP 391

10.11996/JG.j.2095-302X.2018010012

A

2095-302X(2018)01-0012-09

2017-06-01；

2017-07-14

國(guó)家自然科學(xué)基金項(xiàng)目(61473201)

朱云龍(1992-)，男，江蘇如東人，碩士研究生。主要研究方向?yàn)閳D形圖像處理。E-mail：20154229024@stu.suda.edu.cn

翁桂榮(1963-)，男，江蘇蘇州人，教授，學(xué)士。主要研究方向?yàn)閳D形圖像處理、生物信息處理。E-mail：wgr@suda.edu.cn