基于活動(dòng)輪廓模型的腦梗死圖像分割

李智 陳業(yè)航 馮寶 張紹榮 李昌林 陳相猛 劉壯盛 龍晚生

(1.桂林電子科技大學(xué) 電子工程與自動(dòng)化學(xué)院，廣西 桂林 541004；2.桂林航天工業(yè)學(xué)院 電子信息與自動(dòng)化學(xué)院， 廣西 桂林 541004；3.中山大學(xué)附屬江門(mén)市中心醫(yī)院 醫(yī)學(xué)影像研究所，廣東 江門(mén) 529000)

腦梗死是目前世界范圍內(nèi)致殘率較高的疾病之一，大約有70%～80%的腦梗死患者在接受治療后仍然遺留不同程度的運(yùn)動(dòng)功能障礙[1]。腦梗死患者的早期發(fā)現(xiàn)和早期干預(yù)能夠明顯地改善患者預(yù)后。彌散加權(quán)成像(DWI)是一種檢測(cè)人體組織內(nèi)水分子布朗運(yùn)動(dòng)情況的核磁共振成像(MRI)序列[2]，它能夠反映早期腦梗死患者腦組織內(nèi)水分子變化情況，對(duì)腦梗死的早期臨床診斷具有很高的敏感性，是當(dāng)前腦梗死患者影像診斷的主要技術(shù)之一。

針對(duì)急性期(發(fā)病一周內(nèi)為急性期)的腦梗死患者進(jìn)行影像評(píng)估時(shí)，病灶區(qū)的位置和大小與腦梗死患者的功能損傷程度密切相關(guān)，是設(shè)計(jì)臨床治療方案的重要參考指標(biāo)[3-4]。準(zhǔn)確分割腦梗死病灶是診斷患者情況的一個(gè)重要的預(yù)處理步驟。然而，對(duì)于急性期腦梗死患者，由于梗死周?chē)[情況復(fù)雜多變[5]，使得整個(gè)腦梗死病灶在DWI影像上具有邊界模糊、形狀不規(guī)則且內(nèi)部亮度不均勻的特點(diǎn)，給病灶的準(zhǔn)確分割帶來(lái)了挑戰(zhàn)。因此，如何在DWI影像上準(zhǔn)確分割急性期患者的腦梗死病灶，對(duì)于影像輔助評(píng)估具有重要的臨床意義。

活動(dòng)輪廓模型(ACM)是目前醫(yī)療影像分割的常用方法之一。它將圖像分割問(wèn)題轉(zhuǎn)化為水平集函數(shù)求解問(wèn)題，并可進(jìn)行拓?fù)渥兓艿搅藦V泛的關(guān)注[6-8]。張明慧等[9]提出了基于多圖譜的局部區(qū)域活動(dòng)輪廓模型，該方法將多圖譜的形狀先驗(yàn)信息引入活動(dòng)輪廓模型中實(shí)現(xiàn)了對(duì)腦結(jié)構(gòu)的自動(dòng)精確分割。欒紅霞等[10]利用目標(biāo)輪廓的先驗(yàn)知識(shí)對(duì)目標(biāo)進(jìn)行邊緣約束，然后將邊緣約束后的邊緣吸引力場(chǎng)進(jìn)行正則化并融入邊界活動(dòng)輪廓模型中，可以實(shí)現(xiàn)邊界模糊的腦MRI圖像中腦組織的分割。上述兩種方法都是基于目標(biāo)先驗(yàn)形狀信息的活動(dòng)輪廓模型，但是腦梗死的發(fā)病部位和形狀都是隨機(jī)的，很難找到合適的先驗(yàn)形狀信息。Sandhya等[11]提出一種自組織映射網(wǎng)絡(luò)與區(qū)域活動(dòng)輪廓模型相結(jié)合的分割方法，通過(guò)自組織映射網(wǎng)絡(luò)訓(xùn)練神經(jīng)元對(duì)圖像灰度分布進(jìn)行全局建模，進(jìn)而獲得圖像灰度分布的拓?fù)渥兓?，然后將?xùn)練好的神經(jīng)元引入活動(dòng)輪廓模型中，驅(qū)動(dòng)輪廓曲線的演化，但自組織映射的參數(shù)(如學(xué)習(xí)率、學(xué)習(xí)終止)需要人為控制，且神經(jīng)元的數(shù)目難以確定。馬超等[12]提出了一種整合級(jí)聯(lián)隨機(jī)森林與活動(dòng)輪廓模型的腦MR圖像分割方法，該方法通過(guò)級(jí)聯(lián)隨機(jī)森林對(duì)演化曲線進(jìn)行初始化，并利用活動(dòng)輪廓模型進(jìn)一步演化曲線。Ma等[13]提出基于隨機(jī)森林的活動(dòng)輪廓模型，首先對(duì)隨機(jī)森林進(jìn)行訓(xùn)練，預(yù)測(cè)出待分割圖像的腦腫瘤結(jié)構(gòu)，然后將其作為局部活動(dòng)輪廓模型的初始輪廓和空間先驗(yàn)信息進(jìn)行細(xì)分割。但由于腦梗死的DWI圖像噪聲較大，且發(fā)生腦梗死后，其水腫區(qū)域會(huì)以小時(shí)為單位不斷變化，隨機(jī)森林模型容易陷入過(guò)擬合。王醒策等[14]提出了一種結(jié)合圖像全局信息和局部信息的活動(dòng)輪廓模型，使得該模型可以有效分割亮度不均勻的腦血管圖像。Charoensuk等[15]利用局部區(qū)域活動(dòng)輪廓模型在DWI圖像上對(duì)腦梗死病灶進(jìn)行分割。唐文杰等[16]提出一種雙水平集方法對(duì)腦部MR圖像進(jìn)行分割。上述文獻(xiàn)[14]至文獻(xiàn)[16]都是基于局部區(qū)域的活動(dòng)輪廓模型，原理是根據(jù)圖像的局部灰度統(tǒng)計(jì)信息驅(qū)動(dòng)輪廓曲線的演變，可較好的分割亮度不均勻圖像，但由于腦梗死核心周?chē)兴[，致使DWI圖像中梗死區(qū)域出現(xiàn)邊緣模糊的特點(diǎn)，基于局部灰度統(tǒng)計(jì)信息的活動(dòng)輪廓模型無(wú)法準(zhǔn)確收斂至目標(biāo)邊界。近年來(lái)，一些學(xué)者也將深度學(xué)習(xí)應(yīng)用到腦MRI圖像分割。Moeskops等[17]提出了一種多尺度的分段卷積神經(jīng)網(wǎng)絡(luò)來(lái)分割嬰兒和年輕人的大腦組織。Dou等[18]提出一種3D卷積神經(jīng)網(wǎng)絡(luò)來(lái)分割腦MRI圖像中腦出血病灶。但由于醫(yī)學(xué)領(lǐng)域的敏感性和特殊性，較少的訓(xùn)練數(shù)據(jù)限制了深度學(xué)習(xí)方法的性能，深度學(xué)習(xí)模型龐大的模型規(guī)模所具有的優(yōu)勢(shì)此時(shí)并不明顯。

文中針對(duì)急性期腦梗死病灶在DWI圖像中存在的邊界模糊、形狀不規(guī)則和亮度不均勻等問(wèn)題，提出一種模糊速度函數(shù)驅(qū)動(dòng)下的活動(dòng)輪廓模型，來(lái)完成腦梗死病灶的準(zhǔn)確分割：(1)在模型初始化方面，利用小波能量圖下的貝葉斯概率獲取模型初始輪廓，使得初始輪廓快速定位于腦梗死病灶的真實(shí)邊界附近，增強(qiáng)模型的魯棒性和準(zhǔn)確性。(2)將圖像局部熵引入活動(dòng)輪廓模型中。圖像局部熵反映了該局部圖像亮度信息的有序性，可表征腦梗死病灶和正常腦組織之間水分子分布的差異性，在一定程度上解決圖像的亮度不均勻性和噪聲的干擾問(wèn)題。(3)根據(jù)腦梗死病灶區(qū)的邊界模糊特性，提出結(jié)合圖像局部熵和灰度的模糊聚類(lèi)算法計(jì)算模糊隸屬度，進(jìn)一步加強(qiáng)腦梗死病灶與正常組織的區(qū)分。(4)將基于模糊隸屬度的模糊速度函數(shù)引入到模型中，使輪廓曲線在腦梗死病灶的模糊邊界處停止演變，完成腦梗死病灶的分割。(5)最后利用臨床數(shù)據(jù)驗(yàn)證本文方法的有效性。

1 模糊速度函數(shù)驅(qū)動(dòng)下的活動(dòng)輪廓模型

如圖1所示，文中提出的分割算法主要包括4個(gè)步驟：1)在小波變換域下計(jì)算圖像的貝葉斯概率，進(jìn)行活動(dòng)輪廓模型的輪廓曲線初始化，并得到基于初始輪廓的局部區(qū)域；2)結(jié)合圖像局部熵和灰度的模糊聚類(lèi)算法計(jì)算模糊隸屬度；3)將基于模糊隸屬度的模糊速度函數(shù)引入到活動(dòng)輪廓模型中，構(gòu)建基于模糊速度函數(shù)的活動(dòng)輪廓模型；4)活動(dòng)輪廓模型能量泛函演化。

1.1 基于貝葉斯概率的活動(dòng)輪廓模型的初始化

模型輪廓曲線的初始位置對(duì)于模型的分割精度有很大的影響。為了解決對(duì)初始輪廓敏感的問(wèn)題，在小波域下構(gòu)建貝葉斯概率對(duì)模型進(jìn)行初始化。首先，計(jì)算圖像的小波能量，小波能量可以增強(qiáng)腦梗死病灶區(qū)與正常組織的區(qū)分。然后，在小波能量圖上計(jì)算每個(gè)像素的貝葉斯概率，獲取活動(dòng)輪廓模型的初始輪廓，使得初始輪廓位于腦梗死病灶的真實(shí)邊界附近，提高分割的準(zhǔn)確性。

對(duì)DWI圖像的像素(x，y)進(jìn)行二維離散小波變換后，小波能量可以定義為[19]

W(x，y)=EA(x，y)+EH(x，y)+EV(x，y)+

ED(x，y)

(1)

式中：EA(x，y)、EH(x，y)、EV(x，y)和ED(x，y)分別為低頻小波能量和3個(gè)高頻小波能量，定義如下：

圖1 文中算法流程圖

(2)

式中：DA(x，y)為低頻小波系數(shù)；DH(x，y)、DV(x，y)和DD(x，y)分別為水平、垂直、對(duì)角線方向的高頻小波系數(shù)；Q是整數(shù)的集合；r和s分別為鄰域Q內(nèi)的像素的橫坐標(biāo)和縱坐標(biāo)；K(·)為高斯核函數(shù)。采用小波框架方法，將小波頻域下的圖像大小擴(kuò)展到原始圖像的大小。

假設(shè)小波能量圖符合混合高斯模型[22]，則將小波能量圖像分為腦梗死區(qū)域(r=1)與背景區(qū)域(r=2)兩部分?；旌细咚鼓Ｐ涂杀硎緸?/p>

(3)

L(θ)=logP(w|θ)

(4)

(5)

則模型的初始輪廓可定義為

(6)

圖2顯示了小波域下采用貝葉斯概率獲取的初始輪廓。圖2(a)是腦梗死病灶區(qū)的DWI原圖，圖2(b)是小波能量分布圖，圖2(c)是小波域下基于貝葉斯概率獲得的初始輪廓，可看出該初始輪廓位于梗死病灶真實(shí)邊界附近。圖2(d)是水平集函數(shù)的初始化圖，中間白色曲線為水平集初始化曲線對(duì)應(yīng)的輪廓曲線。

圖2 基于貝葉斯概率初始輪廓

1.2 計(jì)算圖像模糊隸屬度

文中根據(jù)腦梗死病灶邊界的模糊特性計(jì)算模糊速度函數(shù)，即采用結(jié)合局部熵和灰度的模糊聚類(lèi)算法，在初始輪廓確定的局部區(qū)域中，計(jì)算模糊隸屬度，從而進(jìn)一步加強(qiáng)腦梗死病灶和正常組織的區(qū)分。

由于急性腦梗死患者梗死中心區(qū)的神經(jīng)細(xì)胞已經(jīng)壞死，無(wú)法正常新陳代謝，且周?chē)殡S缺血半暗帶和水腫情況的出現(xiàn)，因此梗死區(qū)的水分子分布失去固定規(guī)律，與正常腦組織有顯著差異。而圖像的局部熵反映了該局部圖像亮度信息的有序性，可表征腦梗死病灶和正常腦組織之間水分子分布的差異性，在一定程度上解決圖像的亮度不均勻性和噪聲的干擾問(wèn)題，從而加強(qiáng)腦梗死病灶和正常腦組織的區(qū)分。對(duì)于一幅大小為M×N的圖像，設(shè)I(x，y)為圖像中像素點(diǎn)(x，y)處的灰度值，則局部熵可定義為[20]

(7)

其中，

(8)

式中，pxy為像素點(diǎn)(x，y)處的灰度分布概率，如果M×N是以(x，y)為中心的一個(gè)局部圖像窗口，則H為圖像的局部熵。文中的局部圖像窗口大小設(shè)置為9×9。從圖3(e)和圖3(f)可以看出正常組織的局部熵值小于腦梗死病灶局部熵值，這與腦梗死病灶的水分子分布失去固定規(guī)律，比正常組織具有更多信息的臨床背景相吻合。但從圖3(c)、圖3(d)、圖3(e)和圖3(f)可以看出，單一使用圖像的灰度和圖像局部熵很難將腦梗死病灶從正常組織中分離出來(lái)。所以考慮構(gòu)建結(jié)合局部熵特征和灰度特征的模糊聚類(lèi)算法，計(jì)算圖像的模糊隸屬度，進(jìn)一步加強(qiáng)腦梗死病灶和正常腦組織的區(qū)分。

模糊聚類(lèi)算法通過(guò)目標(biāo)函數(shù)迭代尋找最優(yōu)參數(shù)，目標(biāo)函數(shù)定義為

(9)

式中：n表示圖片像素點(diǎn)的總數(shù)；l表示第l個(gè)像素點(diǎn)，Ol為圖像局部熵和灰度構(gòu)成的特征向量；U為2×n的隸屬度矩陣；ulh是向量Ol隸屬h類(lèi)的程度；G為聚類(lèi)中心向量；m為一個(gè)模糊控制常數(shù)；d(Ol，Gh)為Ol和距離中心Gh之間的距離，該距離函數(shù)定義為

d(Ol，Gh)=‖Ol-Gh‖

(10)

重復(fù)計(jì)算ulh和Gh，當(dāng)所有的特征向量被正確分類(lèi)時(shí)，J(U，G)目標(biāo)函數(shù)達(dá)到最小值。ulh和Gh的計(jì)算公式如下：

(11)

(12)

根據(jù)ulh構(gòu)造隸屬腦梗死病灶的模糊隸屬度矩陣Z，則像素點(diǎn)(x，y)的模糊隸屬度可表示為Z(x，y)∈[0，1]。從圖3(g)和圖3(h)中可以看出，文中結(jié)合灰度值和局部熵的模糊聚類(lèi)算法，將基于灰度的圖像空間區(qū)域轉(zhuǎn)變到基于模糊隸屬度的圖像空間域，使得腦梗死病灶和正常組織的模糊隸屬度以0.5為界，從而加強(qiáng)了腦梗死病灶與正常組織的區(qū)分。

(a)原始圖像 (b)病灶區(qū)域放大圖

(c)正常組織的灰度值 (d)病灶組織的灰度值

(e)正常組織的局部熵值 (f)病灶組織的局部熵值

(g)正常組織的模糊隸屬度(h)病灶組織的模糊隸屬度

Fig.3 Gray value，local entropy and fuzzy membership of DWI image

1.3 構(gòu)建模糊速度函數(shù)

為了使模型在腦梗死病灶邊界處停止演變，根據(jù)模糊隸屬度矩陣，分別定義使輪廓曲線向內(nèi)演變的模糊速度函數(shù)V1、輪廓曲線向外演變的模糊速度函數(shù)V2和反映輪廓邊界的模糊速度函數(shù)V3，并有：1)當(dāng)輪廓曲線遠(yuǎn)離邊界時(shí)，V1、V2和V3越來(lái)越大，驅(qū)動(dòng)輪廓曲線快速演變；2)當(dāng)輪廓曲線趨近邊界時(shí)，V1、V2和V3逐漸變小，且在目標(biāo)邊界處時(shí)V1≈V2≈V3≈0，輪廓曲線停止演變。

文中所提出的模型中，V1、V2和V3分別定義為：

V1(x，y)=|e-t1[Z(x，y)-0.5]-1|

(13)

V2(x，y)=|e-t2[0.5-Z(x，y)]-1|

(14)

V3(x，y)=et3[Z(x，y)-0.5]2-1

(15)

式中，t1、t2和t3為整數(shù)，Z(x，y)為像素點(diǎn)(x，y)隸屬腦梗死病灶的模糊隸屬度。

圖4展示了一個(gè)腦梗死病灶的模糊速度函數(shù)。從圖4(b)和圖4(c)中可以看出，在腦梗死病灶邊界點(diǎn)A、點(diǎn)B、點(diǎn)C和點(diǎn)D處，模糊速度函數(shù)Vk(A)、Vk(B)、Vk(C)和Vk(D)，k∈(1，2，3)分別趨近0。

(a)原始圖像

(b)實(shí)型直線的模糊速度

(c)點(diǎn)型直線的模糊速度

1.4 構(gòu)建活動(dòng)輪廓模型的能量泛函

假設(shè)Ω是DWI圖像的定義域，輪廓曲線C將Ω分為輪廓曲線的內(nèi)部區(qū)域Ω1和外部區(qū)域Ω2。根據(jù)模糊速度函數(shù)V1和V2可定義活動(dòng)輪廓模型的區(qū)域項(xiàng)：

V2(x，y)|I(x，y)-c2|2(1-H(φ))]dxdy

(16)

式中，φ為水平集函數(shù)，I(x，y)為像素點(diǎn)(x，y)處的灰度值，c1和c2分別表示輪廓曲線內(nèi)部灰度平均值和輪廓曲線外部灰度平均值。H(φ)為正則化Heaviside函數(shù)，參數(shù)ε用以控制δ(φ)的有效寬度。

根據(jù)反映輪廓邊界的模糊速度函數(shù)V3可定義活動(dòng)輪廓模型的邊界檢測(cè)項(xiàng)：

(17)

結(jié)合式(16)和式(17)，活動(dòng)輪廓模型的能量泛函可定義為

V2(x，y)|I(x，y)-c2|2(1-H(φ))]dxdy+

(18)

式中：等號(hào)右邊第1項(xiàng)是控制輪廓曲線向目標(biāo)邊界演變的區(qū)域項(xiàng)；第2項(xiàng)是控制輪廓曲線在目標(biāo)邊界處停止演化的邊界檢測(cè)項(xiàng)；第3項(xiàng)是為了避免輪廓曲線在演變過(guò)程中對(duì)水平集函數(shù)φ進(jìn)行重新初始化[21]。μ、λ和γ分別是控制區(qū)域項(xiàng)、邊界檢測(cè)項(xiàng)和正則項(xiàng)的權(quán)重參數(shù)，在文中μ、λ、γ和t1、t2、t3的取值都為1[21]。

1.5 能量泛函演化

活動(dòng)輪廓模型輪廓曲線演化過(guò)程可歸結(jié)為最小化能量泛函E(φ)。首先，固定水平集函數(shù)φ，相對(duì)c1和c2最小化E(φ)，可得到：

(19)

為了簡(jiǎn)化表示，書(shū)寫(xiě)上省略變量(x，y)，E(φ)可簡(jiǎn)化為

F(φ，φx，φy)=λe1H(φ)+λe2(1-H(φ))+

(20)

然后，固定c1和c2，最小化E(φ)對(duì)應(yīng)于求解如下偏微分方程：

(21)

其中F對(duì)φ、φx、φy的偏導(dǎo)數(shù)可以表示為

(22)

(23)

(24)

進(jìn)一步對(duì)式(23)和(24)求導(dǎo)：

(25)

(26)

同時(shí)結(jié)合

(27)

將式(22)、(25)-(27)代入式(21)，可以得到F關(guān)于水平集函數(shù)φ的歐拉-拉格朗日方程：

(28)

由變分原理，可得φ的梯度下降流：

(29)

2 分割結(jié)果評(píng)價(jià)指標(biāo)

在評(píng)估算法有效性的過(guò)程中，分別定義陽(yáng)性率(RTP)、假陽(yáng)性率(RFP)和相似度(DS)3個(gè)指標(biāo)來(lái)評(píng)價(jià)算法的分割結(jié)果。

(30)

式中:Am為醫(yī)生手動(dòng)分割結(jié)果；Aa為分割算法結(jié)果；RTP的值越大，說(shuō)明Aa對(duì)Am的覆蓋率越好；RFP的值越小，說(shuō)明Aa包含越少的背景；DS的值越大，說(shuō)明Aa與Am越相似。

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

實(shí)驗(yàn)數(shù)據(jù)采用江門(mén)市中心醫(yī)院的臨床數(shù)據(jù)。由影像科醫(yī)生從臨床數(shù)據(jù)中隨機(jī)選擇11 個(gè)邊界模糊、亮度不均勻的腦梗死患者數(shù)據(jù)進(jìn)行分割結(jié)果的對(duì)比展示。將文中提出的方法與基于區(qū)域的ACM、基于邊界的ACM、結(jié)合二值水平集和形態(tài)學(xué)的ACM、基于局部高斯模型的ACM、基于模糊水平集的ACM和基于圖像局部熵的ACM[20]進(jìn)行分割結(jié)果的對(duì)比。

分割結(jié)果如圖5所示，藍(lán)色曲線包圍的區(qū)域?yàn)榉指罱Y(jié)果。圖5(d)和圖5(f)分別展示了基于區(qū)域的ACM和結(jié)合二值水平集和形態(tài)學(xué)的ACM的分割結(jié)果，可以看出在圖像邊界模糊、亮度不均勻的情況下這兩種方法的分割結(jié)果并不理想，如第K個(gè)數(shù)據(jù)分割結(jié)果出現(xiàn)了邊界泄漏。圖5(e)展示了基于邊界的ACM的分割結(jié)果，在模糊邊界處輪廓曲線無(wú)法收斂到目標(biāo)邊界且在凹陷邊界處都存在邊界泄漏的問(wèn)題，見(jiàn)第A、C、I和J個(gè)數(shù)據(jù)。圖5(g)是基于局部高斯模型的ACM的分割結(jié)果，第A、J、H、I和K個(gè)數(shù)據(jù)均發(fā)生了邊界泄漏問(wèn)題。圖5(h)展示了基于模糊水平集的ACM分割結(jié)果，可以看到在多病灶的情況下分割不完整(如第G、H和K個(gè)數(shù)據(jù))，且在邊界模糊和亮度不均勻的情況下存在著過(guò)分割(如第J個(gè)數(shù)據(jù))。圖5(i)是基于圖像局部熵的分割結(jié)果，當(dāng)腦梗死病灶存在邊界模糊的問(wèn)題時(shí)，該方法不能準(zhǔn)確的分割腦梗死病灶，如第A，F(xiàn)和G個(gè)數(shù)據(jù)。圖5(j)是文中提出方法的分割結(jié)果，較其他幾種分割算法更接近醫(yī)生手動(dòng)分割結(jié)果。

圖5 腦梗死病灶分割結(jié)果

對(duì)上述分割結(jié)果分別計(jì)算RTP、RFP和DS3個(gè)指標(biāo)，結(jié)果如表1所示。結(jié)合二值水平集和形態(tài)學(xué)的ACM雖然得到最好的DS值，但存在欠分割問(wèn)題；基于邊界的ACM雖然得到最好的RTP值，但存在過(guò)分割問(wèn)題；而文中方法DS為0.887，RTP為0.908，RFP為0.025，說(shuō)明文中方法的結(jié)果與醫(yī)生手動(dòng)分割的結(jié)果最為相近。

表1 11個(gè)臨床病例的分割結(jié)果

為了進(jìn)一步驗(yàn)證本文方法的準(zhǔn)確性，又隨機(jī)選擇58個(gè)病例進(jìn)行分割實(shí)驗(yàn)，結(jié)果如表2所示。綜合3個(gè)指標(biāo)，可以看出文中方法的準(zhǔn)確性更好。

從以上結(jié)果可知，由于腦梗死病灶存在模糊特性和對(duì)比度低等問(wèn)題，基于邊界的、區(qū)域的、結(jié)合二值水平集和形態(tài)學(xué)的、基于局部高斯模型的和基于模糊水平集的ACM無(wú)法準(zhǔn)確分割病灶。同樣的，基于圖像局部熵特征的ACM是一個(gè)結(jié)合全局信息和局部信息的ACM，該方法利用全局或局部的灰度統(tǒng)計(jì)信息驅(qū)動(dòng)活動(dòng)輪廓模型進(jìn)行演變，但由于腦梗死病灶具有邊緣模糊的特點(diǎn)，基于全局或局部的灰度統(tǒng)計(jì)信息無(wú)法有效處理邊界模糊的問(wèn)題。

表2 58個(gè)臨床病例的分割結(jié)果

文中方法利用結(jié)合局部熵特征和灰度特征的模糊聚類(lèi)算法計(jì)算圖像的模糊隸屬度，加強(qiáng)腦梗死病灶與正常組織的區(qū)分。然后分別計(jì)算3個(gè)基于模糊隸屬度的模糊速度函數(shù)構(gòu)建模型，使得模型在腦梗死病灶邊界處能夠停止演變。最后選取一個(gè)腦梗死病例的所有DWI圖像進(jìn)行分割，該病灶有8個(gè)連續(xù)層面。結(jié)果如圖6所示，對(duì)比圖6(b)和圖6(c)可以看出，文中提出方法的分割結(jié)果與醫(yī)生手動(dòng)分割結(jié)果相似。

圖6 文中提出的方法對(duì)腦梗死病灶DWI序列圖像分割

4 結(jié)論

為了準(zhǔn)確分割DWI圖像中的腦梗死病灶，文中提出用模糊速度函數(shù)驅(qū)動(dòng)下的活動(dòng)輪廓模型對(duì)其進(jìn)行分割。該模型一方面提出結(jié)合圖像局部熵特征和灰度特征的模糊聚類(lèi)算法計(jì)算模糊隸屬度，以加強(qiáng)腦梗死病灶與正常組織的區(qū)分，并利用基于模糊隸屬度的模糊速度函數(shù)構(gòu)建模型；另一方面利用小波變換域下的貝葉斯概率獲取模型初始輪廓，增強(qiáng)模型的魯棒性和準(zhǔn)確性。實(shí)驗(yàn)結(jié)果表明，本文的模型可有效分割DWI圖像中的腦梗死病灶。