趙秀蘭，劉印文

1. 茌平縣樂平鎮(zhèn)郝集中心衛(wèi)生院 CT室，山東 聊城 252126；2. 聊城市東昌府區(qū)中醫(yī)院 CT室，山東 聊城 252003

引言

腹部動(dòng)脈為多個(gè)器官供血，精確定位動(dòng)脈血管有助于制定外科手術(shù)計(jì)劃和術(shù)中導(dǎo)航[1-3]。腹部動(dòng)脈血管位置、長(zhǎng)度和分支特征個(gè)體差異較大，常借助CT成像獲得腹部器官和血管的三維立體信息，但如何從腹部組織中精確分離動(dòng)脈血管是一個(gè)難題，圖像分割是一種能對(duì)血管定位并提取的方法。

傳統(tǒng)的腹主動(dòng)脈血管分割常采用水平集方法、區(qū)域增長(zhǎng)法、高斯混合模型等算法[4-5]。劉海華等[6]采用高斯混合模型進(jìn)行腹主動(dòng)脈CT圖像分割，分割性能優(yōu)于k均值聚類等算法，但存在不能分割分支血管、分割精度不足等缺點(diǎn)。近些年，深度學(xué)習(xí)方法開始應(yīng)用于醫(yī)學(xué)圖像分割[7-8]，主要通過訓(xùn)練集從CT圖像空間灰度分布中產(chǎn)生一個(gè)預(yù)測(cè)值，然后采用損失函數(shù)獲得測(cè)試標(biāo)簽值與真實(shí)分割圖像標(biāo)簽值的最小差異，完成圖像分割?；谏疃葘W(xué)習(xí)的分割算法已在肝臟、腎臟等大器官分割方面顯示出顯著優(yōu)勢(shì)，但在微小組織如血管分割上精度不足，主要難度在于血管和周圍組織體積顯著不平衡，如腹部動(dòng)脈血管僅占腹部體積的0.2%。許友彬等[9]提出一種基于卷積神經(jīng)網(wǎng)絡(luò)和圖論的肺部血管CT圖像分割算法，能分割出絕大多數(shù)肺部大血管，但不能分割小血管。張睿等[10]提出一種基于全卷積神經(jīng)網(wǎng)絡(luò)（Fully Convolutional Network，F(xiàn)CN）和塊分割的肝臟血管CT圖像分割方法，其中塊由多層圖像裁剪而成，但是交叉熵作為損失函數(shù)不能有效分割小血管。目前基于FCN算法的血管分割多是基于單一斷面進(jìn)行分割、重建，獲得血管三維結(jié)構(gòu)，而單一斷面不能提取足夠數(shù)量的血管特征，血管分割精度難以保證，且FCN算法分割精細(xì)度不高，空間一致性不夠精確。

本研究提出一種基于U型全卷積神經(jīng)網(wǎng)絡(luò)（U-Net）的分割算法，并應(yīng)用于腹部動(dòng)脈CT圖像分割，其中動(dòng)脈血管特征由橫斷面、冠狀面和矢狀面圖像聯(lián)合生成，并采用區(qū)域不平衡塊生成法獲得相同數(shù)量的動(dòng)脈血管和非血管區(qū)域。

1 材料與方法

1.1 材料

選取120例在本院放射科行腹部血管增強(qiáng)CT檢查的患者圖像作為研究對(duì)象，橫斷面CT圖像尺寸為512×512，冠狀面和矢狀面CT圖像尺寸根據(jù)患者體型而定，CT圖像層數(shù)為753-1051，層厚1 mm。全卷積神經(jīng)網(wǎng)絡(luò)訓(xùn)練階段選用90例患者CT圖像，CT三維圖像產(chǎn)生300000個(gè)小塊用于訓(xùn)練，剩余30例患者CT圖像用于分割測(cè)試。

1.2 訓(xùn)練部分

塊是由橫斷面、冠狀面和矢狀面裁剪所得的2D圖像[11]，若按照同等間隔裁剪圖像，則包含血管區(qū)域的塊數(shù)量遠(yuǎn)少于非血管區(qū)域，難以提取足夠多的血管特征用來分割圖像。為了克服此局限性，本研究采用區(qū)域不平衡生成法產(chǎn)生同等數(shù)目的血管和非血管塊[12]，其中血管塊的中心像素為血管像素，非血管塊的中心像素是非血管像素，采用密集取樣獲取血管塊，稀疏采樣獲取非血管塊，均進(jìn)入U(xiǎn)-Net網(wǎng)絡(luò)訓(xùn)練。將從CT測(cè)試血管圖像產(chǎn)生的塊稱為CT塊，塊中灰度值歸一化為[0,1]。將從真實(shí)血管圖像產(chǎn)生的塊稱為標(biāo)簽塊，塊中血管的灰度值為1，非血管區(qū)域灰度值為0。

首先采用3×3×3中值濾波器平滑處理CT血管圖像，然后采用區(qū)域不平衡塊產(chǎn)生方法生成橫斷面、冠狀面和矢狀面像素大小為s×s的CT塊，標(biāo)簽塊取自兩位臨床經(jīng)驗(yàn)豐富的放射科醫(yī)生標(biāo)注的平均結(jié)果作為訓(xùn)練集的，由橫斷面產(chǎn)生的塊集合表示為Pa={Ia，La}，Ia是CT塊的集合，La是標(biāo)簽塊的集合，冠狀面和矢狀面產(chǎn)生的塊集合為Pc和Ps。

接著采用U-Net訓(xùn)練橫斷面、冠狀面和矢狀面產(chǎn)生的塊集合，U-Net最優(yōu)化處理選用自適應(yīng)矩估計(jì)，優(yōu)化過程中損失函數(shù)采用二進(jìn)制交叉熵，最小塊規(guī)模設(shè)置為100，U-Net訓(xùn)練次數(shù)為25，此時(shí)分割圖像體素與真實(shí)圖像體素匹配率最高。U-Net是FCN算法的延伸與改進(jìn)，在CT、MR、顯微鏡等不同模態(tài)圖像分割方面均表現(xiàn)優(yōu)越[13-14]。U-Net結(jié)構(gòu)見圖1所示，U-Net包含一個(gè)左側(cè)的收縮路徑和右側(cè)對(duì)稱的擴(kuò)張路徑，前者用于捕捉上下文特征信息，后者用于目標(biāo)精確定位，二者通過裁剪和復(fù)制將高低層特征圖拼接起來，實(shí)現(xiàn)淺層位置紋理特征與深層語(yǔ)義信息的結(jié)合[15]。具體地，收縮路徑每一層先采用2個(gè)卷積核大小為3×3的卷積，而后通過線性整流函數(shù)（ReLU）和批歸一化操作進(jìn)行激活，接著采用1個(gè)2×2大小的最大池化層進(jìn)行下采樣操作，實(shí)現(xiàn)特征通道數(shù)量翻倍。擴(kuò)張路徑使用大小為2×2的轉(zhuǎn)置卷積層進(jìn)行上采樣，實(shí)現(xiàn)特征通道數(shù)減半。使用跳躍連接操作拼接收縮路徑特征圖與當(dāng)前特征圖，而后采用2個(gè)卷積層和激活函數(shù)提取圖像特征，最后采用大小1×1的卷積核將特征圖通道數(shù)縮減到與分割類別一致。

圖1 U-Net網(wǎng)絡(luò)結(jié)構(gòu)

1.3 分割部分

分割處理分別在橫斷面、冠狀面和矢狀面生成的塊上進(jìn)行，以橫斷面分割為例，首先采用中值濾波器平滑CT圖像，然后沿著坐標(biāo)（x，y）以步長(zhǎng)（s/2，s/2）裁剪CT圖像獲得CT塊，大小為s×s，x，y是沿著橫斷面相互垂直的坐標(biāo)系，CT塊像素灰度值同樣被歸一化為[0,1]，裁剪過程見圖2所示。

圖2 基于二維測(cè)試標(biāo)簽塊的圖像裁剪過程

然后采用U-Net進(jìn)行圖像分割，CT塊輸入U(xiǎn)-Net，產(chǎn)生測(cè)試標(biāo)簽塊。三維CT血管圖像由測(cè)試標(biāo)簽塊重建獲得，測(cè)試標(biāo)簽塊像素值被拷貝到對(duì)應(yīng)CT重建圖像上，且在重建過程中被重疊4次（見圖2），計(jì)算重疊塊像素的平均值作為重建后圖像的體素值，以此基于橫斷面塊的三維血管分割結(jié)果。同理可得基于冠狀面和矢狀面塊的三維血管分割結(jié)果。

選取三斷面所得三維血管圖像中同一位置體素最大值作為后處理圖像體素，然后采用閾值法分割血管，采用膨脹和腐蝕等形態(tài)學(xué)操作連通相鄰區(qū)域，并清除體積小于υm3的孤立區(qū)域，獲得最終血管分割圖像，分割閾值設(shè)置為t，形態(tài)學(xué)操作微球直徑為rmm。圖像后處理參數(shù)設(shè)置為t=0.75，r=2 mm，υ=200 mm3。

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

圖像分割效果采用精確率（Precision rate）、召回率（Recall rate）和Dice系數(shù)評(píng)估[16-18]，三個(gè)指標(biāo)越大，表明分割性能越佳，見公式(1)～(3)，其中TP、FP和FN分別代表分割結(jié)果中真陽(yáng)性、假陽(yáng)性和假陰性的體素?cái)?shù)量。

2 結(jié)果與討論

2.1 基于不同塊尺寸的U-Net算法分割結(jié)果

基于不同塊尺寸的FCN算法血管分割結(jié)果和評(píng)價(jià)指標(biāo)見圖3和表1。塊特征圖像尺寸大小直接影響U-Net算法性能，塊尺寸越小，U-Net訓(xùn)練時(shí)間越短，內(nèi)存消耗越小。由圖3看出，基于不同塊尺寸的U-Net算法均能完整分割腹主動(dòng)脈等較粗壯血管，且三種塊尺寸分割血管視覺差異較小，但均不能有效分割部分細(xì)血管。由表1可知，基于不同塊尺寸的U-Net算法血管分割所得平均Dice系數(shù)隨著塊尺寸增大而減小。塊尺寸最小時(shí)（s=32），平均Dice系數(shù)取得最大值（87.2 %），較塊尺寸s=48（86.8 %）和s=64（86.5 %）的精確率分別高0.4 %和0.7 %。血管分割平均精確率隨著塊尺寸增加先上升后下降，且變化幅度微小，塊尺寸s=48時(shí)，血管分割所得平均精確率最高（86.2 %），較塊尺寸s=32的精確率（85.9 %）高0.3 %，較塊尺寸s=64的精確率高2.3 %。血管分割平均召回率隨著塊尺寸增加先下降后上升，變化范圍微小，塊尺寸s=64時(shí)，本所得平均召回率最高（89.1 %），僅較塊尺寸s=32（88.5 %）高0.6 %，塊尺寸s=48（87.7 %）時(shí)最小。此外，三種尺寸基于不同塊尺寸的U-Net算法所得評(píng)價(jià)指標(biāo)最大值與最小值相差較大，主要是由患者腹主動(dòng)脈形態(tài)位置差異所致。綜上可知，基于不同塊尺寸的U-Net算法血管分割準(zhǔn)確性大致相等，且分割精度較高，因此本研究U-Net算法選取塊尺寸s=32進(jìn)行血管分割。

圖3 基于不同塊尺寸的U-Net 算法所得血管分割結(jié)果

表1 基于不同塊尺寸的U-Net算法血管分割結(jié)果評(píng)價(jià)指標(biāo)比較（%）

2.2 基于不同平面數(shù)量的U-Net算法血管分割結(jié)果

基于單平面和三斷面U-Net算法分割結(jié)果比較見圖4和表2。目前對(duì)血管分割多基于橫斷面圖像進(jìn)行，導(dǎo)致特征提取不充分，三斷面同步提取可獲得更豐富的血管特征，對(duì)血管三維重建至關(guān)重要。本研究U-Net算法中塊特征由橫斷面、冠狀面和矢狀面聯(lián)合產(chǎn)生，塊尺寸采用s=32。由圖4可知，基于三斷面和橫斷面的U-Net算法均能有效分割腹主動(dòng)脈等較為粗壯血管，三斷面分割出的小血管數(shù)量明顯多于橫斷面。由表2可知，基于三斷面的U-Net算法的平均Dice系數(shù)（87.2%）高于橫斷面所得（83.2%），主要得益于三斷面可提供更多特定的組織特征，能辨別分支血管和周圍組織，而橫斷面只能提供單一平面組織特征，只能辨別相對(duì)粗壯血管?；谌龜嗝娴腢-Net算法的平均精確率（85.9%）略低于橫斷面（87.5%）所得，表明基于橫斷面的U-Net算法血管分割不完整，但包含更少的假陽(yáng)性體素（FP）?；谌龜嗝娴腢-Net算法的平均召回率（88.5%）高于橫斷面所得（75.5%），且提升了17.22%，表明三斷面分割顯著降低假陰性體素（FN），提高小血管分割精度。綜合所得，基于三斷面的U-Net算法優(yōu)于單一橫面，但需要訓(xùn)練更多的塊特征降低假陽(yáng)性。

圖4 基于不同斷面的U-Net 算法所得血管分割結(jié)果

2.3 不同血管分割算法比較

將本研究U-Net算法與淺層U-Net算法、基于殘差單元的U-Net算法、水平集方法、區(qū)域增長(zhǎng)法和高斯混合模型進(jìn)行比較，結(jié)果見表3。淺層U-Net包含3個(gè)池化層和3個(gè)上采樣操作[19]，較U-Net少1個(gè)卷積層和池化層，其余與U-Net無異?；跉埐顔卧腢-Net是一種更深層的U-Net[20]，在一定范圍內(nèi)，隨著卷積層數(shù)增加，網(wǎng)絡(luò)特征學(xué)習(xí)能力更強(qiáng)，但不斷加深網(wǎng)絡(luò)容易引起自損函數(shù)梯度下降，導(dǎo)致性能退化。由表3可知，當(dāng)塊尺寸s=32，塊特征由三斷面生成時(shí)，三種U-Net算法所得平均Dice系數(shù)、精確率和召回率差別微小，基于批量歸一化的淺層U-Net算法所得評(píng)價(jià)指標(biāo)均最低，基于殘差單元的N-Net算法的割性能略有提升，但均低于基于本研究批量歸一化的U-Net算法。三種U-Net分割算法所得平均Dice系數(shù)、精確率和召回率顯著高于水平集方法、區(qū)域增長(zhǎng)法和高斯混合模型，主要由于水平集算法無法辨別對(duì)邊界模糊區(qū)域，區(qū)域增長(zhǎng)法分割精度受限于初始種子點(diǎn)的選取，而高斯混合模型對(duì)圖像噪聲和灰度不均較為敏感，導(dǎo)致三種算法分割性能較弱，表明基于本研究批量歸一化的淺層U-Net算法較其他算法具有優(yōu)越性。

表2 基于不同平面的U-Net算法血管分割結(jié)果評(píng)價(jià)指標(biāo)比較（塊尺寸s=32）

表3 基于不同U-Net算法的血管分割結(jié)果評(píng)價(jià)指標(biāo)比較（%，塊尺寸s=32）

3 結(jié)論

本研究提出一種基于U型全卷積神經(jīng)網(wǎng)絡(luò)的腹部血管CT圖像分割算法。實(shí)驗(yàn)結(jié)果表明，當(dāng)塊尺寸s=32，塊特征由三斷面聯(lián)合生成時(shí)，本研究提出的U-Net算法能精確分割腹主動(dòng)脈血管和大部分分支血管，分割性能優(yōu)越，優(yōu)于其他兩種U-Net算法和傳統(tǒng)的血管分割算法，有助于臨床輔助診斷和外科手術(shù)引導(dǎo)。本研究的缺點(diǎn)是未能充分提取血管結(jié)構(gòu)特征，今后將引入全局血管結(jié)構(gòu)特征和血管3D塊特征，更加精確提取分支微細(xì)血管，且將本研究算法應(yīng)用于門靜脈等血管分割。