針對弱邊緣信息的左心室圖像分割算法

王 南，許道云

貴州大學(xué) 計算機(jī)科學(xué)與技術(shù)學(xué)院，貴陽 550025

1 引言

左心室MR 圖像分割是運(yùn)用圖像分割的方法對左心室內(nèi)外膜進(jìn)行準(zhǔn)確分割，其結(jié)果可以為醫(yī)生定量地分析左心室收縮功能，提供參數(shù)指標(biāo)。目前，關(guān)于左心室的分割方法主要分為圖像驅(qū)動模型、基于圖譜的方法[1]、基于深度學(xué)習(xí)的方法[2-5]以及主動輪廓模型等。圖像驅(qū)動模型包括基于閾值、基于區(qū)域信息以及像素相似性質(zhì)分類等方法，該類方法通常具有分割速度快、實現(xiàn)較為簡單的特點，但其不能與左心室的先驗信息有效結(jié)合，對左心室邊界弱邊界、對比度不高的圖像分割能力有限?；趫D譜的方法在分割給定的圖像時，需要與大型訓(xùn)練集匹配分割圖像，該方法的分割結(jié)果受配準(zhǔn)圖像數(shù)據(jù)集的大小所影響[6]?；谏疃葘W(xué)習(xí)的方法，在處理左心室圖像時，具有魯棒提取圖像特征以及實現(xiàn)左心室全自動分割的優(yōu)點，但完成該方法的神經(jīng)網(wǎng)絡(luò)訓(xùn)練，需要大量圖像數(shù)據(jù)，并且訓(xùn)練時間較長。主動輪廓模型相較于其他方法，其優(yōu)點是在高噪聲的情況下，也能得到較好的分割效果，并且容易與其他幾種方法進(jìn)行結(jié)合，分割速度較快。

主動輪廓模型分為參數(shù)輪廓模型和幾何輪廓模型。參數(shù)輪廓模型是Kass 等提出的傳統(tǒng)Snake 模型及其改進(jìn)模型[7-9]，其主要是將輪廓曲線向檢測物邊緣演化的過程與求解能量泛函最小值的過程相關(guān)聯(lián)。該方法在給定適當(dāng)?shù)某跏驾喞螅梢宰詣邮諗康侥芰孔钚顟B(tài)，具有對噪聲和對比度不敏感的優(yōu)點，但其對初始輪廓的位置具有很強(qiáng)的依賴性，在圖像的弱邊界處容易陷入局部最優(yōu)。

幾何活動輪廓模型是基于曲線演化理論與水平集的方法，將平面曲線隱含的表達(dá)為三維曲面函數(shù)，通過曲面的演化來隱含地求解曲線演化[10]。C-V模型[11]是幾何輪廓模型中一個具有代表性的基于區(qū)域信息進(jìn)行圖像分割的算法，其具有易于處理拓?fù)浣Y(jié)構(gòu)變化、可分割較為復(fù)雜的對象、容易實現(xiàn)低維到高維的擴(kuò)展、對初始輪廓不敏感等優(yōu)勢，但同時具有一定局限性，該模型對強(qiáng)度不均勻圖像會產(chǎn)生不良分割結(jié)果，需要反復(fù)初始水平集，造成計算成本高等限制，眾多學(xué)者在該算法的基礎(chǔ)上進(jìn)行改進(jìn)與創(chuàng)新[12-17]。在文獻(xiàn)[13]中，盧小鵬等人將目標(biāo)區(qū)域與背景區(qū)域的灰度平均值之間的差值添加到C-V模型面積項中，使得C-V模型對亮度不均勻圖像的分割能力有所提升。在文獻(xiàn)[14]中，Wang等人提出一種將全局和局部統(tǒng)計信息結(jié)合的水平集新方法，克服了C-V模型不能成功分割強(qiáng)度不均勻圖像的缺點。Murtaza等人在文獻(xiàn)[14]基礎(chǔ)上，使用變異系數(shù)代替方差，提出新的圖像變分模型PSM[15]，在對圖像分割實驗中取得更好的分割結(jié)果。Li 等人在文獻(xiàn)[16]中提出一種對原圖像進(jìn)行擬合的方法。文獻(xiàn)[17]中，盧振泰等人提出了一種基于熵和局部鄰域信息的高斯約束C-V 模型對左心室進(jìn)行分割，該模型具有算法執(zhí)行速度快、避免重新水平集初始化、提高分割精度。文獻(xiàn)[18]中，Li 等人提出距離正則化水平集演化（DRLSE）方法在能量函數(shù)中增加了一個正則項，以避免在水平集演化過程中出現(xiàn)形狀不規(guī)則和不穩(wěn)定性等情況，并且在該模型不需要重新初始化水平集，但該方法對初始輪廓敏感且在對比度不高的圖片中演化過程容易陷入局部最小值。在運(yùn)用DLRSE算法對左心室MR圖像內(nèi)膜進(jìn)行分割時，會因為受到乳突肌與小梁組織的干擾而造成曲線陷入局部最優(yōu)的結(jié)果，為了避免這種情況的出現(xiàn)，文獻(xiàn)[19]中，根據(jù)左心室內(nèi)膜的先驗知識，劉宇提出一種保持凸性的水平集演化模型，該方法能保證曲線演化過程中具有良好的凸性，最終演化到左心室內(nèi)膜邊緣，但該算法存在對初始輪廓敏感，對圖像弱邊緣分割不精確的缺點。

綜上所述，在圖像分割中DRLSE 算法具有對初始輪廓敏感，演化過程僅依賴圖像梯度信息的缺陷。為降低DRLSE 算法對初始輪廓的敏感程度，提升分割弱邊緣、強(qiáng)度不均勻圖像的能力，提高算法對左心室MR 圖像內(nèi)外的分割精度，提出一種適用于弱邊緣信息左心室MR圖像的DRLSE分割算法，實現(xiàn)對左心室內(nèi)膜與外膜分割。在本文的內(nèi)膜算法中，首先，使用Li提出的圖像擬合方法，擬合得到新圖像，計算出原圖像與新圖像的圖像差，替換原PSM模型中的圖像差，形成PSM新的局部項；其次，將文獻(xiàn)[13]中C-V 面積項與上一步得到的PSM 新的局部項添加進(jìn) DRLSE 中；再次，將 CPLSE 模型的保持曲線凸性的曲率因子添加到DRLSE 模型中；最終實現(xiàn)對左心室MR 圖像的內(nèi)膜分割。在對左心室外膜進(jìn)行分割時，首先，運(yùn)用內(nèi)膜分割結(jié)果，結(jié)合左心室心肌各方向厚度基本一致的形狀先驗，構(gòu)造外膜形狀約束項，形成左心室外膜形狀約束項；其次，將外膜形狀約束項與內(nèi)膜分割算法中的PSM新的局部項、文獻(xiàn)[13]中的C-V 模型面積項嵌入到DRLSE 算法中；最終實現(xiàn)對左心室MR圖像外膜的分割。最后通過實驗，對本文提出的模型進(jìn)行測試，對比實驗結(jié)果進(jìn)行分析。

2 基礎(chǔ)知識

2.1 C-V模型

在C-V 模型中，將Ω?R2作為圖像范圍，u:Ω→R作為給出的灰度圖像，定義Ωin是在圖像u上被閉合曲線C所劃分的內(nèi)部區(qū)域，Ωout是在圖像u上被閉合曲線C所劃分的外部區(qū)域。C-V模型的能量泛函定義如下：

其中，μ≥0 ，v≥0 ，λ1＞0 ，λ2＞0 是各個項的系數(shù)，Length(C)是曲線長度項有平滑曲線的作用，Area(inside(C))是Ωin面積項，一般設(shè)其系數(shù)v=0，后兩項為保真項是用于驅(qū)動曲線C向真實邊界演化。

水平集定義表示如下：

E(C)表示為水平集形式：

其中，c1、c2分別為Ωin、Ωout的灰度平均值，其計算公式如下：

δ是Dirac函數(shù)，H是Heaviside函數(shù)。

當(dāng)曲線演化到目標(biāo)邊緣的時，泛函式（3）將達(dá)到最小，其最小化可以通過以下的梯度下降流來求解：

C-V模型可以很好地對圖像進(jìn)行分割，但對灰度不均勻、均勻強(qiáng)度重疊的圖像分割效果很差。2010 年Wang等人提出LCV模型，該模型對圖像灰度不均勻圖像有很好的分割效果，但在圖像均勻強(qiáng)度一致部分所產(chǎn)生的分割結(jié)果較差。

2.2 PSM模型

PSM 模型是 2011 年由 Murtaza 等人在 LCV 模型基礎(chǔ)上提出，該模型具體如下：

其中，u*=gk?u-u，gk是高斯核，u*:Ω*→R作為給出的新灰度圖像，是曲線C將u*內(nèi)外兩部分區(qū)域，m1、m2分別表示區(qū)域灰度的平均值。水平集下的PSM公式如下：

其中：

PSM模型的最小化求解過程可以通過以下方程：

PSM 模型在灰度非均勻圖像的分割結(jié)果要優(yōu)于LCV、C-V模型。

2.3 DLRSE模型

DRLSE 模型作為一種基于區(qū)域的水平集算法，是由Li等人于2010年提出的模型。

定義水平集函數(shù)：

DRLSE模型的能量泛函定義為：

其中，μ≥0 ，λ ＞0 ，α∈R，Rg(φ)是距離正則項，Lg(φ)是加權(quán)長度項可以驅(qū)使邊緣向零水平集向著目標(biāo)邊緣方向演化，Ag(φ)是加權(quán)面積項可以改變水平集的演化速度，它們分別被定義為：

其中，P為勢函數(shù)，是邊緣指示函數(shù)，Gσ是高斯核；? 是梯度算子。在DRLSE方法中，選?。?/p>

能量泛函式（14）的最小化可以通過以下梯度下降流來求解：

DRLSE模型允許有更靈活、更有效的初始化，在曲線演化過程中顯著地減少了迭代次數(shù)與計算時間，但該方法對初始輪廓敏感，僅依靠梯度信息驅(qū)動曲線演化，模型容易陷入局部最優(yōu)。

3 弱邊緣信息左心室圖像分割DRLSE模型

對左心室圖像分割時，DRLSE 方法存在對初始位置敏感且無法有效分割弱邊界圖像的問題，為改善以上情況，本文運(yùn)用Li提出的圖像擬合方法得到新圖像，將新圖像與原圖像之間的圖像差，替換新的PSM 模型中局部項的圖像差，將新的局部統(tǒng)計信息與梯度信息結(jié)合改進(jìn)DRLSE模型。為降低新局部項受到圖像灰度不均勻情況的干擾，將對分割灰度不均勻圖像有一定貢獻(xiàn)的文獻(xiàn)[13]中C-V的面積項加入到DRLSE模型中，得到初步改進(jìn)DRLSE，其水平集下公式為：

其中，Hε1(x)是選取與DRLSE中相同的Hε(x) ，Hε2(x)選取與PSM中相同的Hε(x)，u*=f-u，m1、m2是曲線劃分u*產(chǎn)生區(qū)域內(nèi)外的灰度平均值，擬合圖像f的公式為：

其中：

gε是卷積核函數(shù)。

運(yùn)用|?φ|來消除傳統(tǒng)C-V模型中的Dirac函數(shù)對檢測遠(yuǎn)邊緣的抑制，并將邊停函數(shù)g加入到演化函數(shù)中?；诖耍倪M(jìn)的DRLSE 能量泛函方程的最小化求解對應(yīng)的梯度下降流方程：

在內(nèi)膜分割中，內(nèi)膜形狀類似于橢圓，而且DRLSE算法受到乳突肌或小梁的影響，內(nèi)膜分割結(jié)果可能是凹形曲線，出現(xiàn)欠分割情況，將保持曲線凸性的曲率因子加入到上述改進(jìn)的DRLSE 方法中，這里對文獻(xiàn)[19]中kmean公式作修改，定義：

最終內(nèi)膜分割算法水平集最小化演化方程如下：

在分割外膜時，外膜與肌肉組織的灰度接近，在使用DRLSE 模型對圖像外膜進(jìn)行分割時因受到弱邊緣、對比度低等因素的影響，造成過分割或欠分割等不好結(jié)果，觀察圖像中左心室的特性，發(fā)現(xiàn)左心室心肌各方向厚度基本一致[10]，本文設(shè)計的外膜形狀約束項是采用已分割好的內(nèi)膜與已知左心室最小厚度Rmin結(jié)合先驗信息[20]，結(jié)合先驗信息形成形狀約束，以實現(xiàn)對演化過程中曲線的約束，具體公式如下：

其中dis是外膜分割曲線c上的點與內(nèi)膜曲線的最短距離。

最終外膜分割算法水平集最小化演化方程如下：

其中，ω ＞0 ，γ ＞0 是兩項牽引力的系數(shù)。

本文方法的具體步驟如下：

（1）內(nèi)膜分割

①設(shè)定模型參數(shù)，對初始水平集φ0，即φt(t=0)，最大迭代次數(shù)T，0 ≤t≤T；

②根據(jù)當(dāng)前φt，結(jié)合公式（4），計算更新當(dāng)前區(qū)域內(nèi)外平均值c1、c2；

③ 根據(jù)公式（20），擬合出圖像ft，計算更新u*、m1、m2；

④ 據(jù)公式（22）、（23）、（24），計算更新k、s(k)in、kmean；

⑤根據(jù)公式（25），更新φt+1；

⑥ 判斷t ＞T，若是，則停止；否則判斷φt+1與φt輪廓是否相同，若是，則停止；否則，φt=φt+1，轉(zhuǎn)回②。

內(nèi)膜分割流程圖如圖1。

圖1 內(nèi)膜分割算法流程圖

（2）外膜分割

①設(shè)定模型參數(shù)，對初始水平集φ0，即φt(t=0)，最大迭代次數(shù)T，0 ≤t≤T；

②根據(jù)當(dāng)前φt，結(jié)合公式（4），計算更新當(dāng)前區(qū)域內(nèi)外平均值c1，c2；

③根據(jù)公式（20），擬合出圖像ft，計算更新u*、m1、m2；

④ 根據(jù)公式（26）、（27）、（28），計算更新dmean、s(k)out；

⑤根據(jù)公式（29），更新φt+1；

⑥ 判斷t ＞T，若是，則停止；否則判斷φt+1與φt輪廓是否相同，若是，則停止；否則，φt=φt+1，轉(zhuǎn)回②。

外膜分割流程圖如圖2。

4 實驗結(jié)果與分析

圖2 外膜分割算法流程圖

本文對提出的方法通過實驗進(jìn)行驗證，對其結(jié)果進(jìn)行定量分析，并與DRLSE、CPLSE方法進(jìn)行對比。近年來，基于深度學(xué)習(xí)的方法在圖像分割領(lǐng)域應(yīng)用成為主流方法，對左心室也有較好的分割效果，選擇U-Net 網(wǎng)絡(luò)加入到對比方法中。實驗進(jìn)行驗證所使用的數(shù)據(jù)集來源于http：//www.cse.yorku.ca/～mridataset/的Cardiac MRI dataset，該數(shù)據(jù)集是從33 個受試者獲得共7 980 張的心臟MR 圖像，并對圖像進(jìn)行手工分割獲得5 011 個分段MR 圖像和10 022 個輪廓。每個實驗受測對象由20 幀和8～15個切片組。對U-Net網(wǎng)絡(luò)進(jìn)行訓(xùn)練時，從數(shù)據(jù)集中選取4 200 張帶有標(biāo)簽的左心室圖像，將數(shù)據(jù)圖像分為兩部分，總數(shù)量10%作為測試數(shù)據(jù)集，90%作為訓(xùn)練數(shù)據(jù)集輸入到U-Net網(wǎng)絡(luò)之中進(jìn)行訓(xùn)練，采用二分類交叉熵作為損失函數(shù)，經(jīng)過訓(xùn)練100 次循環(huán)，得到訓(xùn)練好的U-Net網(wǎng)絡(luò)，其內(nèi)膜與外膜訓(xùn)練的損失函數(shù)值分別為0.001 與0.002。在數(shù)據(jù)集中選取不同受試者、不同幀、不同切片的左心室MR圖像，分別使用DRLSE、DRLSE-shape、CPLSE、訓(xùn)練好的U-Net 網(wǎng)絡(luò)和本文模型對MR圖像進(jìn)行分割。選擇的實驗環(huán)境為Windows 7 操作系統(tǒng)，CPU為AMD Athlon II X4 640，5.75 GB內(nèi)存，實驗平臺是Python3.5。

4.1 實驗過程

從數(shù)據(jù)集中選取不同受試者選取不同幀、不同層、像素為256×256的左心室MR圖像進(jìn)行實驗。分別賦予每幅圖各自不同的內(nèi)膜初始化位置，運(yùn)用本文左心室內(nèi)膜分割算法對MR圖像進(jìn)行分割，該算法涉及的參數(shù)設(shè)置為μ=0.04，λ1=2，α=-1.5，λ2=0.1，v=0.001，η=1。為每幅圖賦予各自不同的外膜初始化輪廓，運(yùn)用本文左心室外膜分割算法對MR 圖像進(jìn)行分割，其涉及的參數(shù)設(shè)置為μ=0.04，λ1=2，α=5，λ2=0.1，v=0.001，η=0.25，γ=5。

（1）給定一張左心室的MR圖像，分別運(yùn)用DRLSE、CPLSE及本文算法，在給與三種方法相同初始化位置的情況下，對MR 圖像進(jìn)行內(nèi)膜分割，其分割結(jié)果如圖3，其中每張圖的紅色虛線為手工輪廓，綠色曲線為算法分割結(jié)果。

圖3 同一MR圖像不同初始位置的分割結(jié)果

（2）運(yùn)用（1）中的三種方法對不同類型左心室MR圖片的內(nèi)膜進(jìn)行分割，結(jié)果如圖4。

圖4 正常的MR圖像分割結(jié)果

圖5～圖7 中，情形 1、2、3 依次代表圖像含有乳突肌、弱邊界、乳突肌和弱邊界的情況。紅色虛線為手工輪廓結(jié)果，綠色曲線為算法分割結(jié)果。

圖5 含情形1的MR圖像分割結(jié)果

圖6 含情形2的MR圖像分割結(jié)果

圖7 含情形3的MR圖像分割結(jié)果

（3）運(yùn)用DRLSE、DRLSE-shape、U-Net 網(wǎng)絡(luò)以及本文算法對六個不同左心室圖像進(jìn)行外膜分割，最終MR圖像內(nèi)外膜分割結(jié)果如圖8。

圖8 不同左心室內(nèi)外膜的分割結(jié)果

圖8（a）～（f）中的紅色曲線代表的是專家手工輪廓；白色虛線代表的是DRLSE算法對MR圖像內(nèi)外膜分割結(jié)果；紫色虛線代表的是DRSLE-shape對MR圖像外膜分割結(jié)果；藍(lán)色曲線代表的是U-Net 網(wǎng)絡(luò)對MR 圖像內(nèi)外膜的分割結(jié)果；綠色曲線代表的是本文算法對左心室內(nèi)與外膜的分割結(jié)果。

（4）將專家手工輪廓、本文算法分割好的內(nèi)外膜結(jié)果進(jìn)行二值化后得到結(jié)果如圖9。

4.2 實驗分析

為了對本文提出的模型進(jìn)行定量評價，采用相似系數(shù)（Dice Similarity Coefficient，DSC）、平均垂直距離（Average Perpendicular Distance，APD）。DSC是對兩個區(qū)域重疊度的測量，越接近1，兩幅圖的重疊度越高；APD 是自動分割所得輪廓曲線上的全部點與手動輪廓曲線的最短垂直距離的平均值，其值越小代表輪廓曲線重合越高；質(zhì)心距離（Centroid Distance，CD）是自動分割環(huán)形質(zhì)心與手工分割環(huán)形質(zhì)心的歐氏距離。兩者指標(biāo)定義如下：

圖9 不同左心室心肌分割結(jié)果的二值圖

其中，A代表專家已分割輪廓，Aor代表專家手工分割內(nèi)外膜輪廓之間區(qū)域，B代表經(jīng)過算法自動分割輪廓，Bre代表自動分割內(nèi)外膜輪廓之間區(qū)域，|Ω|代表計算Ω區(qū)域內(nèi)像素點的總數(shù)，例如，|A|+|B|代表兩個區(qū)域內(nèi)全部像素點個數(shù)，|A?B|代表兩個區(qū)域相交部分像素點個數(shù)，xi代表B中的像素點，(xor,yor)代表Aor的質(zhì)心，(xre,yre)代表Bre的質(zhì)心。

表1 給出的是運(yùn)用DRLSE、DRLSE-shape、CPLSE、U-Net 網(wǎng)絡(luò)與本文算法對左心室內(nèi)外膜分割所用平均時間。

表1 各分割算法的平均時耗s

表2給出的是運(yùn)用DRLSE、CPLSE和U-Net網(wǎng)絡(luò)與本文算法分割不同心臟左心室、不同層MR 圖像內(nèi)膜，所得分割結(jié)果的DSC與APD值。

表3給出的是運(yùn)用DRLSE、DRLSE-shape、U-Net網(wǎng)絡(luò)與本文算法分割不同心臟左心室、不同層MR圖像外膜，所得分割結(jié)果的DSC和APD值。

表4給出的是，針對不同圖像，分別采用本文所提算法對內(nèi)膜與外膜的分割結(jié)果、U-Net網(wǎng)絡(luò)對內(nèi)膜與外膜的分割結(jié)果、DRLSE對內(nèi)膜與外膜的分割結(jié)果、DRSLE對內(nèi)膜分割的結(jié)果與DRLSE-shape 對外膜分割的結(jié)果進(jìn)行組合得到左心室心肌結(jié)果，并將其結(jié)果與專家手工分割的心肌結(jié)果進(jìn)行比較，所得到的DSC和CD值。

表2 不同圖片內(nèi)膜分割結(jié)果的APD和DSC值

表3 不同圖片外膜分割結(jié)果的APD和DSC值

表4 不同左心室心肌分割結(jié)果的DSC和CD值

由圖3知DRLSE、CPLSE依據(jù)梯度信息進(jìn)行MR圖像分割時，存在對初始輪廓敏感的問題，本文算法同時依據(jù)梯度信息和圖像局部區(qū)域統(tǒng)計信息來驅(qū)動曲線演化，降低了算法對初始位置的依賴性。

由表1可以直觀看出各算法的時間復(fù)雜度，運(yùn)用各種方法對左心室內(nèi)外膜分割，U-Net的平均消耗時間較短，但該方法對網(wǎng)絡(luò)進(jìn)行訓(xùn)練時，不僅需要大量圖片標(biāo)簽數(shù)據(jù)，還會消耗較長時間。結(jié)合表2、3 可知，相較CPLSE、DRLSE-shape 及本文方法，DRLSE 雖然平均時間最短但其分割效果最不理想。本文算法、CPLSE 及DRLSE-shape方法的運(yùn)行時間相近，本文算法在運(yùn)行速度方面表現(xiàn)并不突出，主要因為本文內(nèi)膜與外膜算法在DRLSE 基礎(chǔ)上，兩者除了都增加新局部信息及輪廓內(nèi)外區(qū)域灰度平均值的計算之外，內(nèi)膜算法添加了曲率因子的計算，外膜算法添加了形狀約束項的計算，從而影響了內(nèi)外膜算法的運(yùn)行速度。

由圖4～圖7 與表2 知，對比 DRLSE、CPLSE 模型、U-Net網(wǎng)絡(luò)，本文算法加入圖像局部信息，對邊緣處有加強(qiáng)的作用，在處理弱邊緣、對比度低且有乳突肌影響的左心室圖像時，可以得到更高的DSC平均值，本文算法所得輪廓包圍的內(nèi)部區(qū)域與專家手工分割所得輪廓包圍的內(nèi)部區(qū)域更加重合。分割結(jié)果APD平均值相對于對比中的其他模型分割所得APD 平均值更低，本文算法所得的輪廓更加靠近專家手工分割輪廓，有更高的分割精度。

由圖8 與表3 知，深度學(xué)習(xí)方法對左心室MR 圖像分割時，可能會無法對圖像做出分割。在對DRLSE 算法加入本文先驗形狀約束項后，提高了DRLSE 方法的分割能力。本文外膜分割算法基于本文內(nèi)膜分割效果的基礎(chǔ)上進(jìn)行的，內(nèi)膜分割結(jié)果優(yōu)于DRLSE、CPLSE模型、U-Net網(wǎng)絡(luò)分割結(jié)果，并且在外膜分割方法中依靠了圖像差的局部信息，所以本文外膜分割效果是對比方法中最好。

由圖9 與表4 知，在對左心室整體的分割結(jié)果對比中，U-Net 網(wǎng)絡(luò)對左心室心肌整體分割的質(zhì)心CD 平均值最低，更接近手工輪廓質(zhì)心，但其DSC均值在對比方法中低于本文算法的DSC均值，且本文的質(zhì)心CD值較傳統(tǒng)方法分割結(jié)果相較于傳統(tǒng)分割方法更接近手工輪廓的質(zhì)心，本文算法有更好的分割效果。

5 結(jié)論及下一步工作

針對運(yùn)用DRLSE 模型分割存在弱邊界、對比度低的左心室圖像時不能得到很好分割效果，為了提高左心室內(nèi)外膜的分割精度，本文提出改進(jìn)的DRSLE 算法內(nèi)膜分割算法和帶有形狀約束的外膜分割算法。通過區(qū)域擬合圖像與原圖像比較，提供局部區(qū)域信息，將改進(jìn)的C-V 面積項與局部區(qū)域信息與保持凸性的方法添加到DRLSE算法中，以此驅(qū)動曲線進(jìn)行演化。實驗表明，本文算法在對弱邊界、對比度低的左心室圖像進(jìn)行內(nèi)膜分割時，相對DRLSE、CPLSE算法、U-Net網(wǎng)絡(luò)有更高的分割精度；在對弱邊界外膜分割時，也取得相對穩(wěn)定、較高的分割精度。但本文算法，需要手動給出初始輪廓，以及調(diào)節(jié)部分參數(shù)，在分割過程中仍依賴人工經(jīng)驗。未來的工作將是，需求新的方法、減少變參數(shù)量、降低模型對人工依賴，以達(dá)到低成本人工成本同時達(dá)到更高的分割精度。