方玲玲，王欣

遼寧師范大學(xué)計(jì)算機(jī)與信息技術(shù)學(xué)院,大連 116000

0 引 言

顱內(nèi)腦腫瘤分為原發(fā)性腫瘤和繼發(fā)性腫瘤(單翔 等，2018；賴小波 等，2019)。原發(fā)性腫瘤是一種由大腦異常增生而引起的神經(jīng)性疾病，經(jīng)常出現(xiàn)在腦組織、腦神經(jīng)和腦膜位置。由身體其他部位(如肝臟、肺部等)的惡性腫瘤轉(zhuǎn)移至腦部而形成的稱之為繼發(fā)性腦腫瘤。據(jù)統(tǒng)計(jì)，顱內(nèi)腫瘤的患病率高達(dá)0.003 2%。多數(shù)情況下，腦腫瘤易發(fā)生于腦部幕上的額葉、顳葉等，通常由病毒感染或化學(xué)誘因等高危因素造成，其患病比例約為75%。具體來說，幕下腦腫瘤發(fā)病概率略低于幕上腫瘤，一般在小腦半球、橋腦和側(cè)腦室等位置發(fā)病率較高，常見的腫瘤為星形細(xì)胞瘤、髓母細(xì)胞瘤等。其中，最常見的腦部疾病為腦膠質(zhì)細(xì)胞瘤，其患病比例接近50%(葛婷 等，2017；Bagherimofidi等，2019)。磁共振成像(magnetic resonance imaging, MRI)具有對(duì)腫瘤位置和軟組織結(jié)構(gòu)反應(yīng)靈敏的特點(diǎn)，可在對(duì)人體損害最小的情況下獲得多角度、多模態(tài)成像，良好反映人體結(jié)構(gòu)和功能信息(Galv?o等，2020；Leblond和Kauffmann，2016)。醫(yī)師通過觀測(cè)MRI腦部醫(yī)學(xué)圖像中腫瘤的位置及大小，進(jìn)而判斷其良惡性狀況，實(shí)現(xiàn)對(duì)病灶區(qū)域進(jìn)行精確檢測(cè)。但是，由于醫(yī)學(xué)設(shè)備引起的場(chǎng)偏移和患者自身的活動(dòng)導(dǎo)致所獲得的MRI腦部醫(yī)學(xué)圖像中存在點(diǎn)灰度不均、信噪比差異等問題，需要一種普適算法輔助醫(yī)師對(duì)腦腫瘤進(jìn)行精確分割。超像素的概念最早由Ren和 Malik(2003)提出，該算法一經(jīng)提出便普遍應(yīng)用于圖像分割領(lǐng)域。經(jīng)過十余年的發(fā)展在醫(yī)學(xué)圖像分割方面使用得愈加廣泛，具體包括輔助診療(李孟歆 等，2020；姜濤 等，2020；劉靜和黃玉清，2018)、圖像配準(zhǔn)(劉晨 等，2019；Ji等，2014；Vupputuri等，2020)、多器官分割(Bechar等，2018；Khosravanian等，2021)、3維重建(吳連耀和萬旺根，2020；張志忠和李曉明，2017；Zhao等，2017)、軟組織提取(王海鷗 等，2019)和病灶良惡性判斷(Angulakshmi和 Lakshmi Priya，2019；Soltaninejad 等，2018)等。

基于腦部MRI醫(yī)學(xué)圖像的超像素/體素分割算法(Huang 等，2018；Randrianasoa 等，2018；祁琪，2018；Saha等，2018)通過設(shè)定權(quán)值或初始種子點(diǎn)的個(gè)數(shù)和位置，將2維圖像或3維圖像的點(diǎn)按照顏色、紋理和頻譜等特征，劃分成可表示圖像局部特征或功能子結(jié)構(gòu)的區(qū)域塊。要求得到的超像素/體素能緊湊整齊地獲得具有較強(qiáng)依附性，且最大程度保留腫瘤特征的邊界信息。同時(shí)，需要保證算法對(duì)圖像中的噪聲具有更強(qiáng)的魯棒性，有效減少計(jì)算冗余和大規(guī)模對(duì)象的運(yùn)行時(shí)間，從而生成形狀規(guī)則、邊界精準(zhǔn)且具有較低計(jì)算復(fù)雜度的圖像分割算法。

由于超像素/體素分割算法可平滑準(zhǔn)確地描述圖像邊緣信息，將具有相同結(jié)構(gòu)的點(diǎn)劃分至同一子區(qū)域，從而大幅度減少數(shù)據(jù)信息的規(guī)模計(jì)算，可得到腦腫瘤解剖區(qū)域內(nèi)部結(jié)構(gòu)或功能信息，其反映出腫瘤與周圍組織間的動(dòng)態(tài)關(guān)系，輔助醫(yī)師獲得腫瘤的位置與大小等信息，有效提高醫(yī)師疾病診斷和確定后續(xù)治療方案的效率。

本文使用超像素/體素的最新研究算法對(duì)MRI腦腫瘤醫(yī)學(xué)圖像進(jìn)行分割，通過邊緣召回率、欠分割錯(cuò)誤率、緊密度評(píng)測(cè)和可達(dá)分割準(zhǔn)確率4項(xiàng)指標(biāo)分析各算法性能。此外，針對(duì)超像素/體素算法在醫(yī)學(xué)圖像分割領(lǐng)域中的應(yīng)用現(xiàn)狀來闡述其未來發(fā)展趨勢(shì)和局限性，具有一定的應(yīng)用價(jià)值。

1 MRI醫(yī)學(xué)圖像特征

MRI圖像(Kong等，2019；Irving等，2016)利用氫原子核在外加磁場(chǎng)中不同的組織結(jié)構(gòu)對(duì)能量的吸收狀況不同這一特性，生成具有較高空間分辨率和時(shí)間分辨率的全方位成像。MRI圖像最大的優(yōu)點(diǎn)是在不損傷患者身體健康的情況下，獲得反映不同結(jié)構(gòu)信息的多模態(tài)圖像，因此常常應(yīng)用于腦疾病檢測(cè)。下面針對(duì)MRI腦部醫(yī)學(xué)圖像的不同特性進(jìn)行分析。

1.1 MRI腦部醫(yī)學(xué)圖像的多角度特性

MRI腦部醫(yī)學(xué)圖像可以多平面成像，主要通過水平面、矢狀面和冠狀面3個(gè)角度進(jìn)行觀測(cè)，分別將人體橫斷為上下兩面、縱切為左右兩部以及縱斷為前后兩部，如圖1所示。圖1(a)(b)分別展示水平序列為128和133的腦腫瘤多角度解剖圖像(來自開放存取系列影像研究數(shù)據(jù)集(Open Access Series of Imaging Studies, OASIS)中序號(hào)為OAS2_0099_MR2的多角度成像)。

圖1 多角度下的MRI腦腫瘤醫(yī)學(xué)圖像

水平面角度可直觀反映腦腫瘤大小和占位效應(yīng)，且根據(jù)密度可獲得腦腫瘤的水腫區(qū)域或瘤轉(zhuǎn)移導(dǎo)致的鈣化區(qū)域。矢狀面和冠狀面角度進(jìn)一步直觀反映腦腫瘤的位置以及由于腦瘤轉(zhuǎn)移導(dǎo)致壓迫神經(jīng)或組織器官的情況，有效輔助醫(yī)師判斷腫瘤的良惡性狀況并確定后續(xù)治療方案。

1.2 MRI腦部醫(yī)學(xué)圖像的多模態(tài)特性

多模態(tài)MRI腦部醫(yī)學(xué)圖像具體表現(xiàn)為T1加權(quán)圖像(T1 weighted)、T2加權(quán)圖像(T2 weighted)、T1c圖像(T1 contrast,在T1加權(quán)圖像中使用造影劑)和流體衰減反轉(zhuǎn)恢復(fù)圖像(fluid attenuated inversion recovery, FLAIR)。不同模態(tài)的MRI圖像對(duì)顱內(nèi)腫瘤顯示的側(cè)重點(diǎn)也有所不同，如在T1加權(quán)圖像中，腦組織間的對(duì)比度較明顯，可有效分離腦部結(jié)構(gòu)和水腫區(qū)域，鑒別有明顯水腫占位效應(yīng)的高級(jí)別膠質(zhì)瘤和轉(zhuǎn)移瘤等疾病。T2加權(quán)圖像可有效分離人腦健康組織和顯著水腫區(qū)域，可有效鑒定由于細(xì)胞毒性擴(kuò)散或大面積水腫導(dǎo)致的中樞神經(jīng)系統(tǒng)原發(fā)性淋巴瘤等疾病。T1c圖像可有效劃分腦部結(jié)構(gòu)的囊變區(qū)域和腫瘤區(qū)域，針對(duì)腦腫瘤中心壞死部分彌散或腫瘤內(nèi)血管充盈情況進(jìn)行腦出血情況判斷。FLAIR圖像對(duì)蛛網(wǎng)膜下腔和腦實(shí)質(zhì)內(nèi)的病灶較為敏感，當(dāng)靠近腦脊液與健康腦組織的邊界處發(fā)生病變時(shí)，F(xiàn)LAIR圖像中的弛豫時(shí)間會(huì)隨之發(fā)生改變，病灶區(qū)域?qū)?huì)顯現(xiàn)為高信號(hào)。因此，F(xiàn)LAIR成像可體現(xiàn)腦腫瘤區(qū)域與大腦其他組織之間的明顯差異，有效分離腦部灰質(zhì)、白質(zhì)和腦脊液等結(jié)構(gòu)。

圖2顯示了多模態(tài)下的MRI腦腫瘤醫(yī)學(xué)圖像(來自多模態(tài)腦腫瘤分割挑戰(zhàn)賽(Multimodal Brain Tumor Segmentation Challenge, BraTS)2018數(shù)據(jù)集序號(hào)為TCIA08_469_1和TCIA08_218_1的多模態(tài)圖像)，真實(shí)數(shù)據(jù)圖像(ground truth)中黃色區(qū)域表示整個(gè)腦腫瘤的占位面積，紅色區(qū)域表示腫瘤的核心位置，綠色區(qū)域表示增強(qiáng)腫瘤結(jié)構(gòu)。通過ground truth給定的腫瘤結(jié)構(gòu)劃分結(jié)果可判斷腦部壞死狀況和出血占位情況，且由于病灶在不同模態(tài)下的密度信號(hào)不同，通過多模態(tài)信息整合和對(duì)比可有效獲取高性能的分割結(jié)果。

圖2 多模態(tài)下的MRI腦腫瘤醫(yī)學(xué)圖像

2 超像素/體素分割算法

由于超像素/體素分割算法能夠有效提取圖像的中層特征，可大幅度降低算法對(duì)噪聲點(diǎn)的敏感度，具有較強(qiáng)的魯棒性。使用超像素/體素分割圖像可有效增強(qiáng)對(duì)象規(guī)模較大和后續(xù)處理復(fù)雜的醫(yī)學(xué)圖像計(jì)算效率，降低計(jì)算復(fù)雜度。且隨著圖像尺寸和維度增加，該算法能將結(jié)構(gòu)一致且表示相似局部信息的點(diǎn)劃分至同一區(qū)域，有效表達(dá)局部特征信息和功能子結(jié)構(gòu)。因此越來越多的超像素/體素分割算法應(yīng)用到圖像的預(yù)處理步驟。

根據(jù)算法原理，將超像素/體素分割算法分為兩類，一是基于圖論的分割算法，二是基于梯度上升的分割算法，具體可參考附表1。

附表1 超像素/體素分割算法對(duì)比

基于圖論的分割算法通過圖論知識(shí)，將圖像以G=〈V,E〉的形式表示成一張帶有權(quán)值ωi,j的無向圖(V表示頂點(diǎn)，E表示邊)，該值表示連接頂點(diǎn)Vi和Vj的邊Ei,j在特定特征上的差別。通過不同的準(zhǔn)則將分割任務(wù)轉(zhuǎn)換為求取能量方程最小值的代價(jià)問題，從而獲得超像素/體素圖像分割結(jié)果，如圖3(a)所示。

圖3 超像素分割算法

其經(jīng)典算法包括：基于圖論(graph-based)的算法(Felzenszwalb和 Huttenlocher，2004；張明慧 等，2016)、標(biāo)準(zhǔn)化分割(normalized cut)(Ghosh等，2018；Dutta等，2019；Wang 等，2019；王春瑤 等，2014)和基于熵率(Liu等，2014)的算法(entropy rate superpixel, ERS)等?；谔荻壬仙姆指钏惴ㄍㄟ^設(shè)定初始種子點(diǎn)個(gè)數(shù)和位置的方式，迭代遍歷點(diǎn)之間的差異性以及特征的相似程度，從而將具有相同特征的像素點(diǎn)/體素點(diǎn)劃分至同一子區(qū)域，如圖3(b)所示。其經(jīng)典算法有：Watershed算法(Vincent和 Soille，1991)、基于水平集理論的Turbopixels算法(Levinshtein等，2009；Xing等，2020)和簡單線性迭代聚類算法(simple linear iterative clustering algorithm, SLIC)(Achanta 等，2012)等。

2.1 基于圖論的分割算法

王亞杰等人(2014)提出了一種基于熵率法的超像素分割算法，其主要思想為通過計(jì)算像素點(diǎn)之間的熵對(duì)區(qū)域塊整體緊湊均勻的作用，進(jìn)而獲得每個(gè)區(qū)域的空間頻率，獲得較為平滑的圖像邊界信息，劃分出的超像素符合圖像的功能信息和結(jié)構(gòu)信息。但是該算法對(duì)處于邊緣附近的噪聲點(diǎn)靈敏，導(dǎo)致其分割準(zhǔn)確性稍顯不足。在王亞杰等人(2014)方法的基礎(chǔ)上，Gao等人(2017)將圖像表示成以超像素或連接組件作為節(jié)點(diǎn)的加權(quán)圖，用梯度和空間距離線索對(duì)每個(gè)節(jié)點(diǎn)進(jìn)行排序，從而得到超像素和基于連通分量的圖像線索。該算法可有效克服王亞杰等人(2014)方法中的圖像噪聲和血管偽影等問題，具有較強(qiáng)的魯棒性。但該算法在處理與周圍組織器官存在黏連現(xiàn)象的具有較強(qiáng)依附性的醫(yī)學(xué)圖像時(shí)，不能較好地獲得病灶位置和分割信息。

Zhang等人(2017)提出了一種基于圖論的圖像分割算法，其主要思想是從預(yù)分割核開始定義核的粗略中心和大小，通過笛卡兒坐標(biāo)系中的橢圓狀邊界映射到極坐標(biāo)系統(tǒng)中的線，并利用動(dòng)態(tài)規(guī)劃和確保最佳閉合輪廓的迭代方法來識(shí)別所構(gòu)建圖中的全局最佳路徑，具有良好的圖搜索性能。但是該算法生成的超像素邊界極不規(guī)則，像素塊在區(qū)域密度大的地方分布較密，區(qū)域密度小的地方分布較稀疏。該算法同樣具有Gao等人(2017)方法中對(duì)有明顯占位效應(yīng)的腦腫瘤、受側(cè)腦室擠壓導(dǎo)致腫瘤周圍水腫和瘤內(nèi)出血的弱邊界分割準(zhǔn)確率略有不足的問題。為有效克服算法對(duì)弱邊界信息提取不足的情況，崔強(qiáng)(2018)提出了一種基于圖切(graph cut)和隨機(jī)游走(random walk)的交互式超體素分割算法，該算法的主要思想是通過處理代表體素點(diǎn)所處位置的權(quán)值圖譜，表示該點(diǎn)行至另一體素點(diǎn)的概率。并通過相鄰已知MRI切片圖像的形狀特征初始化未知圖像的權(quán)值和圖譜概率，從而劃分超體素區(qū)域。該算法有效克服了傳統(tǒng)算法需要大量人工標(biāo)記的分割圖譜導(dǎo)致時(shí)間成本過高的問題,解決了Zhang等人(2017)和Gao等人(2017)方法中存在的與血管黏連的病灶弱邊界分割準(zhǔn)確率不高的問題，有效提高超體素分割算法的準(zhǔn)確性。但是該算法忽略了切片之間的空間特征，僅從2維角度預(yù)測(cè)腦腫瘤生長規(guī)律，不能準(zhǔn)確地描述實(shí)際腦腫瘤的生長趨勢(shì)。且使用已知切片作為未知切片的初始量導(dǎo)致誤差迭代，不具有較強(qiáng)的魯棒性。

Ghosh等人(2018)提出了一種在迭代k-means聚類框架中結(jié)合normalized cut超像素的分割算法，其主要思想是通過譜聚類將normalized cut算法與k-means聚類框架結(jié)合，通過迭代簡單的k-means聚類算法來優(yōu)化歸一化割的代價(jià)函數(shù)，可有效測(cè)量像素之間的相似性，且由于算法的計(jì)算復(fù)雜度是線性的，可減少迭代計(jì)算量，提高算法效率。該算法可生成整齊平整的超像素邊界，但是分割的準(zhǔn)確率較低，且得到的像素塊大小基本相同，導(dǎo)致無法獲得細(xì)小邊界信息。

2.2 基于梯度上升的分割算法

Xiang等人(2010)提出了基于Turbopixels的渦輪像素分割算法，該算法對(duì)于每個(gè)局部窗口中的像素都引入線性變換，將顏色向量映射為聚類指示向量。該算法通過將線性變換下的估計(jì)誤差相加以獲得目標(biāo)函數(shù)，并從該目標(biāo)函數(shù)中獲得全局最優(yōu)分割結(jié)果。該方法效率較高，但對(duì)于紋理圖像而言，生成的超像素中會(huì)包含不同的對(duì)象，使得病灶的弱邊界分割效果較差。Li和Chen(2015)在Xiang 等人(2010)方法的基礎(chǔ)上進(jìn)行改進(jìn)，提出了一種線性譜聚類(linear spectral clustering, LSC)的超像素分割算法，該算法的主要思想是通過核函數(shù)將像素點(diǎn)的位置與大小信息映射到高維的特征空間中，從而得到像素點(diǎn)的相似程度，獲得超像素邊界?？朔iang 等人(2010)方法對(duì)紋理圖像提取特征略有不足的問題，可得到較為規(guī)整，且分割像素塊大小基本一致的超像素邊界。但該算法只有在像素較為平坦的區(qū)域才能獲得光滑整齊的超像素，在邊界區(qū)域(非均勻區(qū)域)獲得的超像素邊界較為雜亂。

為實(shí)現(xiàn)圖像弱邊界的準(zhǔn)確分割，Shen 等人(2016)提出了一種基于密度的聚類算法。該算法克服了經(jīng)典基于密度的聚類算法(density-based spatial clustering of applications with noise, DBSCAN)中計(jì)算種子點(diǎn)到每個(gè)像素點(diǎn)距離導(dǎo)致運(yùn)行速度過慢的問題。其主要思想是通過計(jì)算鄰近距離項(xiàng)和種子距離項(xiàng)判斷像素點(diǎn)的相似特征來決定點(diǎn)歸屬情況，獲得初始的超像素分割結(jié)果。然后，計(jì)算顏色特征和空間特征將小簇的像素塊進(jìn)行合并，從而得到具有較高邊界粘附性的超像素邊界。但是，該算法為每個(gè)像素塊設(shè)置一個(gè)最大區(qū)域值，若達(dá)到該值，即使像素點(diǎn)與該塊具有相同特征也不進(jìn)行合并。然后通過計(jì)算特征距離進(jìn)行合并操作，從而造成計(jì)算的冗余，導(dǎo)致算法運(yùn)行時(shí)間較長。且生成的超像素邊界較為雜亂，無法按照?qǐng)D像的結(jié)構(gòu)子功能進(jìn)行劃分。在Shen 等人(2016)方法的基礎(chǔ)上，Wu 等人(2017)提出了一種基于超像素和分水嶺的多光子顯微圖像分割方法，其主要思想是用超像素代替像素來分析多光子顯微圖像，將空間、CIE Lab 顏色空間和相位一致性特征的超像素結(jié)合，從而將圖像分割成保持邊界細(xì)節(jié)的小塊。接著，利用標(biāo)記控制分水嶺從重建圖像中得到細(xì)胞邊緣信息，使算法可有效獲得弱邊界信息。該算法的分割準(zhǔn)確率較高，可有效克服Xiang 等人(2010)方法中弱邊界分割結(jié)果略有不足的問題。但該算法無法有效克服噪聲，當(dāng)噪聲點(diǎn)位于圖像邊界附近時(shí)，將對(duì)最終的超像素分割結(jié)果產(chǎn)生極大影響，不具備較強(qiáng)的魯棒性。

Amami 等人(2019)克服了Wu 等人(2017)方法中為得到弱邊界準(zhǔn)確信息導(dǎo)致超像素形狀不規(guī)則的問題，提出了一種基于SLIC的醫(yī)學(xué)圖像自適應(yīng)超體素算法AdaSLIC。該算法主要通過求得子區(qū)域種子點(diǎn)到最近圖像輪廓的距離得到體素點(diǎn)的數(shù)值大小，進(jìn)而得到圖像的邊界信息。與Shen 等人(2016)算法相比，該算法具有更強(qiáng)的緊湊型和自適應(yīng)性，可得到較好的分割結(jié)果。但在計(jì)算過程中，其未將感興趣區(qū)域與不感興趣區(qū)域區(qū)分開，導(dǎo)致計(jì)算過程中出現(xiàn)大量冗余現(xiàn)象，致使計(jì)算效率較低。同年，Singh等人(2020)提出了一種基于SLIC超像素的圖像算法，該算法將簡單線性迭代聚類思想與自動(dòng)可調(diào)的基于高斯徑向基函數(shù)核的模糊c均值算法相結(jié)合。該算法可用于圖像分割以及深度學(xué)習(xí)網(wǎng)絡(luò)使用，具有較強(qiáng)的魯棒性，且提高了算法的運(yùn)行效率。但是該算法不能較好地獲得圖像的弱邊界信息，導(dǎo)致分割效率有所欠缺。

為有效增強(qiáng)算法的魯棒性和自適應(yīng)性，張志龍 等人(2020)提出了一種基于密度、圖論和聚類思想的采用密度峰值進(jìn)行聚類的超像素分割算法。該算法的主要思想是利用初始設(shè)置的圓形模板來獲取局部圖像的密度，利用處于邊緣的像素點(diǎn)局部密度較小的特性，將求得的局部密度最大點(diǎn)設(shè)置為種子點(diǎn)。然后用一個(gè)多叉的歸屬關(guān)系樹表示所計(jì)算的每個(gè)像素點(diǎn)之間的距離大小和所屬關(guān)系，從而求得超像素的分割結(jié)果。該算法的優(yōu)點(diǎn)是不需要進(jìn)行迭代優(yōu)化，計(jì)算效率較高，且獲得的超體素分割結(jié)果在初始設(shè)置不同大小的圓形模板或設(shè)置不同數(shù)目種子點(diǎn)的情況下，都能獲得具有較強(qiáng)依附性且貼近圖像邊緣的超像素分割邊界。但該算法獲得的超像素分割邊界不具有規(guī)則平滑的形狀，無法清晰地表達(dá)鄰域區(qū)域的信息內(nèi)容，不利于后續(xù)的相同特征提取操作。

2.3 超像素/體素分割算法評(píng)價(jià)指標(biāo)

基于MRI腦腫瘤醫(yī)學(xué)圖像的超像素/體素分割算法評(píng)價(jià)指標(biāo)(Stutz等，2018；Wang等，2017)主要關(guān)注：算法求得的病灶邊界信息與ground truth是否具有較高的邊界貼合度；另外，算法是否最大程度地保留腦腫瘤核心區(qū)域、增強(qiáng)區(qū)域和明顯水腫占位區(qū)域的邊界信息?，F(xiàn)有的評(píng)價(jià)指標(biāo)(羅學(xué)剛 等，2019；宋熙煜 等，2015)包括邊緣召回率(boundary recall, BR)、欠分割錯(cuò)誤率(under-segmentation error, UE)、緊密度評(píng)測(cè)(compactness measure, CM)、可達(dá)分割準(zhǔn)確率(achievable segmentation accuracy, ASA)以及其他標(biāo)準(zhǔn)，如圖4所示。其中，黑色實(shí)線表示超像素/體素算法分割的腦腫瘤邊界，黑色虛線表示對(duì)應(yīng)ground truth的病灶邊緣。

圖4 超像素/體素算法評(píng)價(jià)標(biāo)準(zhǔn)示意圖

2.3.1 邊緣召回率

BR是超像素/體素算法分割MRI腦腫瘤邊界結(jié)果與ground truth相重疊的比率，用來計(jì)算位于病灶邊界周圍γ(一般取1或2，本文設(shè)置為2個(gè)單位點(diǎn)內(nèi)超像素/體素分割結(jié)果與對(duì)應(yīng)ground truth中點(diǎn)的個(gè)數(shù)比值，計(jì)算為(Liu等，2011)

(1)

式中，p和q分別代表ground truth和超像素/體素病灶分割結(jié)果中的目標(biāo)區(qū)域點(diǎn)，δG表示ground truth分割結(jié)果G中位于病灶邊緣γ范圍內(nèi)所有像素/體素點(diǎn)的集合，δR表示超像素/體素分割結(jié)果R中超像素/體素邊緣上所有點(diǎn)的集合，f(t)用來判別圖像上兩點(diǎn)間的歐氏距離是否小于γ，具體為

BR判斷分割結(jié)果好壞如圖4(a)所示，灰色區(qū)域表示病灶邊界允許存在γ范圍內(nèi)的誤差點(diǎn)，優(yōu)秀的超像素/體素分割算法虛線和實(shí)線基本覆蓋甚至完全覆蓋，表示超像素/體素分割的腦腫瘤邊緣與ground truth病灶的邊界貼合度高。

2.3.2 欠分割錯(cuò)誤率

UE用于計(jì)算超像素/體素腦腫瘤劃分區(qū)域超出ground truth對(duì)應(yīng)病灶的面積和比值，計(jì)算為(Buyssens等，2014)

(2)

式中，S表示超像素/體素腦腫瘤分割結(jié)果，G表示ground truth的病灶真值，gi表示G中的第i個(gè)病灶劃分區(qū)域，Sk表示S中第k個(gè)超像素/體素塊，|Sk-gi|表示符合真實(shí)病灶劃分范圍內(nèi)的超像素/體素中點(diǎn)的個(gè)數(shù)。Achanta等人(2012)克服了式(2)沒有考慮超像素/體素間可能存在獨(dú)立像素點(diǎn)的問題，對(duì)UE計(jì)算式進(jìn)行改進(jìn)，計(jì)算為

(3)

式中,N表示MRI醫(yī)學(xué)圖像中點(diǎn)的個(gè)數(shù)，Sk|Sk∩gi表示超像素/體素分割第k個(gè)超像素/體素與第i個(gè)ground truth病灶區(qū)域中不重合點(diǎn)的個(gè)數(shù)，B表示容錯(cuò)率，該實(shí)驗(yàn)將其設(shè)置為5%。UE判斷分割結(jié)果好壞如圖4(b)所示，灰色區(qū)域表示UE面積，UE區(qū)域越小，表示超像素/體素算法求得的腦腫瘤邊界越精準(zhǔn)，腦腫瘤邊緣與病灶邊界的貼合度越高。

2.3.3 緊密度評(píng)價(jià)

CM用于評(píng)價(jià)超像素/體素的形狀、大小的緊湊規(guī)則程度(Schick等，2013)，F(xiàn)u等人(2014)定義了CM的計(jì)算式為

(4)

式中，Len(Sk)和Area(Sk)分別表示第k個(gè)超像素/體素的周長和面積大小，n表示超像素/體素的分辨率，N表示MRI腦部醫(yī)學(xué)圖像中點(diǎn)的總數(shù)。CM值一定程度上取決于超像素/體素分辨率大小，其判斷分割結(jié)果好壞如圖4(c)所示，CM值越大，表明超像素/體素算法求得的腦腫瘤邊緣形狀越規(guī)則，越接近ground truth病灶邊界。

2.3.4 可達(dá)分割準(zhǔn)確率

ASA通過將超像素/體素腦腫瘤分割結(jié)果與ground truth病灶分割結(jié)果中所有區(qū)域進(jìn)行對(duì)比，為每一個(gè)超像素/體素匹配一個(gè)與之最相近的ground truth病灶區(qū)域塊，通過計(jì)算所有區(qū)域內(nèi)分配標(biāo)簽正確點(diǎn)的個(gè)數(shù)和MRI腦部醫(yī)學(xué)圖像點(diǎn)總數(shù)N的比值，計(jì)算為

(5)

ASA判斷算法分割結(jié)果好壞如圖4(d)所示，黑點(diǎn)表示正確分割點(diǎn)(即黑色實(shí)線內(nèi)正確劃分至黑色虛線內(nèi)的點(diǎn))，白點(diǎn)表示錯(cuò)分割點(diǎn)，灰點(diǎn)表示孤立點(diǎn)。優(yōu)秀的超像素/體素分割算法保證黑點(diǎn)與白點(diǎn)和灰點(diǎn)的比值最大，從而獲得更高準(zhǔn)確率的分割結(jié)果。

3 對(duì)比實(shí)驗(yàn)與分析

磁共振成像設(shè)備對(duì)腦腫瘤分割效果較好，為了便于腦腫瘤分割算法的研究，國際醫(yī)學(xué)圖像計(jì)算和計(jì)算機(jī)輔助干預(yù)會(huì)議(Medical Image Computing and Computer Assisted Interventions Conference, MICCAI)每年舉辦多模態(tài)腦腫瘤分割挑戰(zhàn)賽，期間提供了大量測(cè)試集(test set)和訓(xùn)練集(training set)，用以標(biāo)記腦腫瘤的位置與腦內(nèi)水腫、壞死和出血區(qū)域，ground truth劃分完整腫瘤(whole tumor)、增強(qiáng)腫瘤(enhance tumor)和腫瘤核心(tumor core)區(qū)域(附錄介紹了部分公開的MRI腦腫瘤醫(yī)學(xué)圖像數(shù)據(jù)集)。本文實(shí)驗(yàn)使用BraTS 2018數(shù)據(jù)集，包含 285個(gè)顱內(nèi)腫瘤病例，每個(gè)病例包含T1加權(quán)圖像、T2加權(quán)圖像、T1c圖像、FLAIR圖像和對(duì)應(yīng)的ground truth，圖像大小為240×240×155。圖5顯示了圖像中種子點(diǎn)個(gè)數(shù)設(shè)置為200的超像素算法分割表皮樣囊腫的結(jié)果(來自BraTS18數(shù)據(jù)集序號(hào)為TCIA08_469_1的圖像)。本文若無特殊說明，默認(rèn)實(shí)驗(yàn)設(shè)置超像素個(gè)數(shù)為200(分割代碼可參考https://pan.baidu.com/s/1wWEL72Tt0Ky18 Wg3 myuVzg)。

圖5 12種超像素分割MRI腦部醫(yī)學(xué)圖像的分割結(jié)果

為進(jìn)一步比較不同超像素算法分割腦腫瘤醫(yī)學(xué)圖像的性能，將graph-based算法(Abate等，2016)、normalized cut算法(Van Den Heuvel等，2008；王春瑤 等，2014；Cour等，2005)、基于熵率算法(Liu 等，2011)、拓?fù)浔３终齽t算法(Fu 等，2014；Tang 等，2012)、lazy random walk算法(Shen 等，2014)、turbopixels算法(Levinshtein等，2009)、DBSCAN算法(Shen 等，2016)、LSC算法(Li和Chen，2015)、SLIC算法(Angulakshmi 等，2019；Saha 等，2018)和SLICO(simple linear iterative cluster zero)算法等10種可直接對(duì)超像素種子點(diǎn)個(gè)數(shù)進(jìn)行設(shè)置的算法進(jìn)行性能對(duì)比。此外，分水嶺算法(Vincent和Soille，1991；Wu 等，2017)和superpixel lattice算法(Moore 等，2008)無法實(shí)現(xiàn)超像素?cái)?shù)目的精準(zhǔn)控制，因此本文對(duì)這兩種算法不進(jìn)行性能評(píng)測(cè)。實(shí)驗(yàn)對(duì)BraTS 2018數(shù)據(jù)集中15幅3維腦腫瘤圖像進(jìn)行超像素分割。BraTS 2018數(shù)據(jù)集提供將腦腫瘤劃分成完整腫瘤、核心腫瘤和增強(qiáng)腫瘤3個(gè)區(qū)域的ground truth，如圖6(a)所示(來自BraTS 2018數(shù)據(jù)集序號(hào)為TCIA08_218_1的圖像)。

圖6 腫瘤區(qū)域超像素分割圖像

圖6(b)顯示了圖5中12種超像素分割算法對(duì)應(yīng)的腫瘤結(jié)果圖像(灰色、黑色和白色分別代表增強(qiáng)腫瘤區(qū)域、核心腫瘤區(qū)域和完整腫瘤區(qū)域)。白線為相應(yīng)算法獲得的超像素邊界，若交界處沒有白線，則證明該算法不能有效劃分功能不同的區(qū)域邊界。由于醫(yī)學(xué)圖像分割更側(cè)重于對(duì)病灶的形狀、大小和位置等進(jìn)行測(cè)量和判斷，因此本文只對(duì)完整腦腫瘤區(qū)域進(jìn)行算法測(cè)評(píng)，圖7—圖10分別展示BR、UE、CM和ASA等4個(gè)評(píng)價(jià)指標(biāo)在超像素個(gè)數(shù)為200、400、600、800和1 000下分割完整腦腫瘤區(qū)域的定量分析曲線。其中，BR和UE指標(biāo)用來評(píng)測(cè)超像素分割算法得到的腦腫瘤邊界與ground truth中實(shí)際病灶的邊界依附性和最大程度保留腦部圖像中腫瘤邊界特征信息的性能。CM指標(biāo)用來測(cè)評(píng)生成的超像素分割邊界緊湊整齊的形狀規(guī)則程度及是否表示成具有較強(qiáng)圖像局部特征或功能子結(jié)構(gòu)的區(qū)域塊。ASA指標(biāo)用來測(cè)評(píng)分割邊界的精準(zhǔn)信息及是否對(duì)噪聲干擾具有較強(qiáng)魯棒性。

圖7 BR指標(biāo)變化曲線

圖8 UE指標(biāo)變化曲線

圖9 CM指標(biāo)變化曲線

圖10 ASA指標(biāo)變化曲線

通過圖5—圖10可看出graph-based算法可較為精準(zhǔn)地分割核心腫瘤，有效鑒別血管充盈的腦腫瘤區(qū)域。但對(duì)完整腫瘤區(qū)域和增強(qiáng)腫瘤區(qū)域分割準(zhǔn)確能力稍顯不足，易出現(xiàn)對(duì)灰度差異不明顯腫瘤的漏診情況。normalized cut算法在BR、UE、CM和ASA共4項(xiàng)指標(biāo)評(píng)測(cè)中性能較為突出，可獲得具有強(qiáng)邊界依附性的腦腫瘤邊界，極大程度保留腫瘤邊界的特征信息。但該算法將病灶區(qū)域與大腦結(jié)構(gòu)中腦部灰質(zhì)和白質(zhì)劃分至同一超像素內(nèi)，使得完整腫瘤區(qū)域劃分成多個(gè)不能有效表示腦部功能子結(jié)構(gòu)的區(qū)域塊，不利于后續(xù)腦腫瘤區(qū)域特征分區(qū)域提取。

superpixel lattice算法可較好地定位核心腫瘤位置，但分割的超像素邊界不具有較強(qiáng)的依附性，無法獲取精準(zhǔn)的腫瘤結(jié)構(gòu)信息?；陟芈仕惴軠?zhǔn)確獲取增強(qiáng)腫瘤邊界信息，獲得腦腫瘤區(qū)域與周圍腦腫瘤細(xì)胞有明顯密度差異的邊界信息，有效鑒別由不規(guī)則鈣化和骨化引起的砂礫體型腫瘤或由于顱骨受損導(dǎo)致硬化的病灶。但該算法生成的超像素邊界形狀不規(guī)則，容易出現(xiàn)將具有相同結(jié)構(gòu)的區(qū)域塊割裂的情況，不能有效實(shí)現(xiàn)鄰域信息表達(dá)。

拓?fù)浔３终齽t算法克服了基于熵率算法無法對(duì)像素邊界實(shí)現(xiàn)精確描述的問題，可生成緊湊整齊的超像素，在CM和ASA指標(biāo)下的分割性能良好。但該算法沒有解決由于初始超像素個(gè)數(shù)不同導(dǎo)致最終的分割結(jié)果質(zhì)量跨度較大問題，不能較好地劃分核心腫瘤和增強(qiáng)區(qū)域具有模糊線性的邊界。lazy random walk算法可生成較為規(guī)則的核心腫瘤超像素邊界，較為完整地保留增強(qiáng)腫瘤的結(jié)構(gòu)信息，有效診斷由血管增多導(dǎo)致的惡變腫瘤和頸內(nèi)動(dòng)脈供血形成的腦膜瘤等疾病。但該算法不能獲得增強(qiáng)腫瘤邊界信息，BR和UE指標(biāo)較低，未能最大程度保留腦部圖像中腫瘤邊界信息的特征。

分水嶺算法可獲得有明顯占位效應(yīng)的腦腫瘤或受側(cè)腦室擠壓導(dǎo)致的明顯腫瘤周圍水腫和瘤內(nèi)出血的弱邊界分割信息，可獲得較為完整的增強(qiáng)腫瘤邊界信息。但該算法得到的超像素不符合大腦功能子單元結(jié)構(gòu)，針對(duì)邊界較為模糊的核心腫瘤和完整腫瘤區(qū)域，其分割準(zhǔn)確率稍顯不足，會(huì)導(dǎo)致具有大腦功能塊與不具備腦部功能塊劃分至同一超像素內(nèi)，在計(jì)算聚類中心到超像素范圍內(nèi)點(diǎn)距離時(shí)增加腦部特征點(diǎn)到無用點(diǎn)的計(jì)算過程，導(dǎo)致計(jì)算效率過低。turbopixels算法有效克服初始設(shè)置的超像素個(gè)數(shù)不同導(dǎo)致分割結(jié)果準(zhǔn)確率差異過大問題，有效增強(qiáng)算法魯棒性，且可較為完整地保留增強(qiáng)腫瘤的結(jié)構(gòu)信息。該算法在4項(xiàng)指標(biāo)下的數(shù)值較為穩(wěn)定，隨著設(shè)置的種子點(diǎn)數(shù)目增加，獲得的分割性能較為穩(wěn)定。但該算法對(duì)整體灰度對(duì)比度較小，且對(duì)于腦腫瘤位置與周圍組織器官存在粘黏現(xiàn)象的醫(yī)學(xué)影像而言，其分割準(zhǔn)確率大大降低，尤其在不同的種子點(diǎn)數(shù)目下的ASA分割性能有待提高。DBSCAN算法可有效獲得精準(zhǔn)的完整腫瘤、核心腫瘤和增強(qiáng)腫瘤區(qū)域的邊界信息，根據(jù)圖像密度有效鑒定腫瘤實(shí)部、壞死和液化區(qū)域，提供腦部感染等并發(fā)癥的腫瘤信息，但該算法對(duì)噪聲點(diǎn)較為敏感，對(duì)處于邊緣信息的噪聲點(diǎn)不具有較強(qiáng)的魯棒性。

LSC算法可有效克服單模態(tài)梭形纖維增生性腦腫瘤和致密的脈管結(jié)構(gòu)性血管腫瘤與有絲分裂增加導(dǎo)致的非典型性腦腫瘤存在的邊界模糊情況和醫(yī)學(xué)成像設(shè)備本身存在的模糊性等問題。但該算法生成的超像素邊界將具有相同特征和功能子結(jié)構(gòu)的大腦區(qū)域劃分成多塊，無法正確反映腦腫瘤形狀、大小和外觀等形態(tài)和與周圍腦膜或血管的牽拉情況。從圖11看出，該算法在不同數(shù)目的種子點(diǎn)下分割性能相差較大，其中，CM和ASA指標(biāo)比相同條件下的分割算法數(shù)值偏低。SLIC算法具有較強(qiáng)的緊湊性，完整保留腦部醫(yī)學(xué)圖像紋理結(jié)構(gòu)和特征連續(xù)性，便于后續(xù)腦腫瘤特征提取和分割。但算法計(jì)算過程存在大量冗余，不利于大對(duì)象規(guī)模的分割操作。SLICO算法在SLIC算法上進(jìn)行改進(jìn)，具有較低的計(jì)算復(fù)雜度和自適應(yīng)性特點(diǎn)的高效率圖像分割算法，可有效降低醫(yī)師發(fā)生錯(cuò)診和漏診的情況，利用計(jì)算機(jī)輔助診療進(jìn)行高效精準(zhǔn)的病灶處理。由于高級(jí)別膠質(zhì)瘤通常由具有明顯核異型性的多形性星型細(xì)胞所構(gòu)成，在MRI醫(yī)學(xué)圖像中表現(xiàn)為具有明顯占位效應(yīng)的組織水腫和瘤內(nèi)出血導(dǎo)致增強(qiáng)腫瘤區(qū)域顯現(xiàn)明顯環(huán)形強(qiáng)化。且腫瘤實(shí)質(zhì)整體呈現(xiàn)高信號(hào)，具有較強(qiáng)的紋理特征和顏色特征。圖11展示了超像素個(gè)數(shù)設(shè)置為200和400分割高級(jí)膠質(zhì)瘤圖像的結(jié)果(來自BraTS18數(shù)據(jù)集序號(hào)為TCIA06_211_1的圖像)。

圖11 不同超像素?cái)?shù)目對(duì)應(yīng)的超像素分割結(jié)果

由于超像素的數(shù)目與圖像的細(xì)節(jié)紋理有關(guān)，為有效測(cè)試在不同超像素尺寸下的腦部醫(yī)學(xué)圖像分割性能，本文進(jìn)行了200、400、600、800、1 000共5種不同的超像素種子點(diǎn)個(gè)數(shù)的對(duì)比實(shí)驗(yàn)。從圖7—圖11可以看出，當(dāng)設(shè)置的種子點(diǎn)數(shù)目較大時(shí)，算法獲得的腫瘤邊界信息更為精準(zhǔn)，BR、UE和ASA指標(biāo)會(huì)隨著超像素尺寸的減小而增加。但是，隨著超像素尺寸的減小，將會(huì)導(dǎo)致分割結(jié)果出現(xiàn)過擬合現(xiàn)象，使得位于同一功能塊的超像素劃分至不同的區(qū)域，從而無法正確地描述病灶與健康腦組織之間的鄰域關(guān)系。當(dāng)設(shè)置的種子點(diǎn)數(shù)目過小時(shí)，可獲得較為規(guī)則的超像素，但是獲得的邊界信息準(zhǔn)確度有待提高，無法精準(zhǔn)地描述病灶的占位效應(yīng)信息，CM分割指標(biāo)會(huì)隨著超像素尺寸的增加而提高。

從視覺效果來看，graph-based、normalized cut、lazy random walk、DBSCAN和LSC算法可較好地保留腫瘤邊界信息，具有較好的局部局灶信息特征。拓?fù)浔３终齽t、turbopixels、SLIC和SLICO算法分割的超像素形狀結(jié)構(gòu)上更加完整緊湊，對(duì)病灶邊界的特征描述較為平滑柔和，以此彌補(bǔ)算法對(duì)邊界描述的不足之處。

4 超像素/體素分割算法在醫(yī)學(xué)圖像上的發(fā)展前景

4.1 提高超像素/體素算法分割效率

隨著Hinton等人(2006)基于手寫數(shù)字識(shí)別的深度網(wǎng)絡(luò)算法的提出，解決了傳統(tǒng)機(jī)器學(xué)習(xí)中反向傳播(back propagation, BP)算法隨著網(wǎng)絡(luò)層數(shù)增加調(diào)參導(dǎo)致梯度爆炸或消失問題，實(shí)現(xiàn)將深度學(xué)習(xí)網(wǎng)絡(luò)退化成淺層網(wǎng)絡(luò)問題，從而開啟深度學(xué)習(xí)的熱潮。如今，深度學(xué)習(xí)在圖像分割、語義分析和語音識(shí)別等方面應(yīng)用愈加廣泛?；卺t(yī)學(xué)圖像的超像素/體素算法嘗試使用基于深度學(xué)習(xí)實(shí)現(xiàn)病灶精準(zhǔn)分割，便于后續(xù)對(duì)醫(yī)學(xué)圖像進(jìn)行處理，得到更為準(zhǔn)確的病灶分割結(jié)果。

金蘭依等人(2018)提出一種基于半監(jiān)督階梯網(wǎng)絡(luò)的醫(yī)學(xué)圖像分割算法，其主要思想為構(gòu)建一個(gè)內(nèi)嵌多層感知機(jī)與卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural networks, CNN)的半監(jiān)督機(jī)器學(xué)習(xí)架構(gòu)，為減少計(jì)算冗余，克服超像素分割邊界的不規(guī)則性導(dǎo)致信息丟失問題，算法通過構(gòu)建梯形網(wǎng)絡(luò)并實(shí)現(xiàn)網(wǎng)絡(luò)層之間的橫向連接加快網(wǎng)絡(luò)調(diào)參的收斂速度，解決多層網(wǎng)絡(luò)梯度消失的問題，且極大增強(qiáng)算法準(zhǔn)確率，實(shí)現(xiàn)病灶的精準(zhǔn)分割，較好地輔助醫(yī)師預(yù)測(cè)并診斷病灶良惡性狀況和確定后續(xù)治療方案。

呂朝暉(2018)提出一種基于超像素的多器官醫(yī)學(xué)圖像分割算法，其主要思想是通過卷積神經(jīng)網(wǎng)絡(luò)解決具有高灰度相似性的超像素醫(yī)學(xué)圖像無法準(zhǔn)確分割的問題，利用人工標(biāo)注的真實(shí)病灶圖像訓(xùn)練包含多個(gè)集成分類器的多器官分割機(jī)器學(xué)習(xí)模型，通過該模型可有效提高超像素特征提取的準(zhǔn)確率，獲取腫瘤生長的重要特征，實(shí)現(xiàn)對(duì)生長腫瘤進(jìn)行準(zhǔn)確分割，為后續(xù)腫瘤的治療提供幫助。

超像素/體素分割算法能有效捕捉醫(yī)學(xué)圖像病灶局部信息，并將其與圖像冗余信息分至不同超像素/體素內(nèi)，將具有相似紋理和顏色等特征的點(diǎn)劃分在一個(gè)區(qū)域塊中，且該塊能有效表達(dá)圖像功能子結(jié)構(gòu)的特征信息，因此常常用于圖像分割預(yù)處理階段，旨在獲得對(duì)圖像宏觀特征有更佳描述的最小區(qū)域塊分類。國內(nèi)外研究者目前致力于通過機(jī)器學(xué)習(xí)提高超像素/體素分割算法的準(zhǔn)確性以及特征信息的提取效率，在此基礎(chǔ)上實(shí)現(xiàn)病灶分割(Ji 等，2014；孫明杰，2017)、3維重建(胡軍，2010；李媛媛，2008)等操作，以滿足當(dāng)今大數(shù)據(jù)處理時(shí)代要求(周友龍等，2010；曾元一，2019)。

4.2 擴(kuò)展算法應(yīng)用維度

隨著醫(yī)學(xué)設(shè)備的逐漸完善和計(jì)算機(jī)存儲(chǔ)及計(jì)算能力的不斷升級(jí)，醫(yī)師診斷醫(yī)學(xué)影像、階段性診療報(bào)告和病歷數(shù)據(jù)信息呈爆炸式增長，傳統(tǒng)的超體素分割算法通過先分割出固定的3維面片，計(jì)算距離體素點(diǎn)最近的種子點(diǎn)標(biāo)簽值并進(jìn)行點(diǎn)的歸屬劃分，進(jìn)而求得圖像邊緣信息，該算法不能有效根據(jù)圖像自身的特征信息自適應(yīng)解決由于圖像初始超體素個(gè)數(shù)不同導(dǎo)致最終的分割結(jié)果質(zhì)量不同的問題，且該算法通過將3維醫(yī)學(xué)圖像分割成多幅2維醫(yī)學(xué)圖像進(jìn)行超像素分割，通過合成多幅2維圖像分割結(jié)果獲得最終的3維病灶邊界信息，導(dǎo)致病灶片與片之間的關(guān)系特征消失，忽略了病灶整體的空間特征，致使最終的分割效果差強(qiáng)人意。為增進(jìn)算法準(zhǔn)確性，劉斌(2017)將病灶片之間的特征信息考慮在超像素劃分范圍內(nèi)，通過將上一序列圖迭代求得的分割結(jié)果作為其相鄰切片的初始分割輪廓，以此增加相鄰病灶片之間的聯(lián)系，但該算法是基于2維圖像操作，并未真正實(shí)現(xiàn)3維圖像病灶分割。腫瘤細(xì)胞的生長受到多方面的因素影響，僅從2維角度預(yù)測(cè)腫瘤生長的規(guī)律和趨勢(shì)不夠全面，不能準(zhǔn)確描述實(shí)際腫瘤生長趨勢(shì)。

楊雨婷(2018)提出了一種基于3維大腦磁共振成像的超體素分割算法，該算法的主要思想為通過將全卷積神經(jīng)網(wǎng)絡(luò)和超體素分割思想結(jié)合形成可有效進(jìn)行并行計(jì)算提高算法效率的特征信息提取算法，通過對(duì)超體素提供的細(xì)節(jié)特征提取實(shí)現(xiàn)與卷積神經(jīng)網(wǎng)絡(luò)獲得的整體與局部信息結(jié)合，從而實(shí)現(xiàn)對(duì)腦部結(jié)構(gòu)中腦脊液、灰質(zhì)和白質(zhì)的有效分割。

基于3維的超體素算法實(shí)現(xiàn)對(duì)整體3維醫(yī)學(xué)圖像進(jìn)行病灶精準(zhǔn)分割操作，完整保留其圖像紋理結(jié)構(gòu)和特征連續(xù)性，避免由于將圖像拆分成多幅2維圖像導(dǎo)致圖像特征向量不連貫的問題，從而解決初始設(shè)置的體素種子點(diǎn)數(shù)目不同導(dǎo)致結(jié)果分割準(zhǔn)確性不同的狀況，有效增強(qiáng)算法的分割精度便于后續(xù)對(duì)3維醫(yī)學(xué)圖像進(jìn)行處理，得到更為準(zhǔn)確的病灶分割結(jié)果(Cao和Cao，2013；Zhang 等，2013)。

受計(jì)算機(jī)設(shè)備和病人階段性病理影響，僅憑一幅單模態(tài)3維醫(yī)學(xué)圖像無法提供足夠精準(zhǔn)的病例信息進(jìn)行診療，且對(duì)實(shí)現(xiàn)腫瘤模型重建，描述腫瘤位置、大小和特征信息顯然是不充分的。由于3維醫(yī)學(xué)影像分辨率較高，其包含的特征信息較為龐大，且腫瘤自身具有復(fù)雜性和多樣性，可以通過有效整合多模態(tài)圖像數(shù)據(jù)信息，降低數(shù)據(jù)存儲(chǔ)和處理的成本和時(shí)間，減少計(jì)算復(fù)雜度，提高數(shù)據(jù)處理的高效性(李鏘 等，2020；Lee 等，2014；Giraud等，2018；Vupputuri等，2018)。因此，如何針對(duì)單模態(tài)梭形纖維增生性腫瘤和致密的脈管結(jié)構(gòu)性血管腫瘤與有絲分裂增加導(dǎo)致的非典型性腫瘤存在的邊界模糊情況和醫(yī)學(xué)成像設(shè)備本身存在的模糊性等問題可通過多模態(tài)磁共振成像信息整合有效解決(劉宏 等，2014)。

5 結(jié) 語

超像素/體素算法可有效提高大對(duì)象規(guī)模和后續(xù)處理復(fù)雜的醫(yī)學(xué)圖像處理效率，隨著醫(yī)學(xué)圖像尺寸和維度增加，算法有利于結(jié)構(gòu)局部信息的表達(dá)，通過劃分特征相似的區(qū)域塊，為后續(xù)能夠獲得正確表達(dá)病灶結(jié)構(gòu)信息的目標(biāo)分割奠定基礎(chǔ)。基于醫(yī)學(xué)圖像的超像素/體素算法已成為當(dāng)前圖像分割領(lǐng)域熱點(diǎn)，受到國內(nèi)外學(xué)者的密切關(guān)注。通過介紹超像素/體素分割算法的評(píng)價(jià)指標(biāo)、國內(nèi)外最新發(fā)展動(dòng)態(tài)和對(duì)比實(shí)驗(yàn)與分析，可看出超像素/體素分割算法具有較高的分割性能，研究領(lǐng)域具有良好的發(fā)展前景。但該算法在醫(yī)學(xué)圖像上的應(yīng)用仍存在幾點(diǎn)不足：1)如何克服算法將具有大腦功能的區(qū)域塊與不具備腦部功能的區(qū)域塊劃分至同一子區(qū)域內(nèi)，且使得初始設(shè)置的種子點(diǎn)可最大程度接近病灶結(jié)構(gòu)以滿足當(dāng)今大數(shù)據(jù)和實(shí)時(shí)性要求，實(shí)現(xiàn)病灶結(jié)構(gòu)精確分割，提高算法精確度與性能為當(dāng)前研究一大難點(diǎn)。2)如何有效克服醫(yī)學(xué)圖像中存在于邊界附近的異常點(diǎn)和噪聲點(diǎn)造成算法過擬合分割或分割不足的情況，在實(shí)現(xiàn)分割區(qū)域具有強(qiáng)空間連通性并保證得到的分割邊界不跨越遮擋影像邊緣信息的條件下，實(shí)現(xiàn)病灶的3維重建為另一大難點(diǎn)。隨著機(jī)器學(xué)習(xí)的興起，如何有效實(shí)現(xiàn)2維至3維醫(yī)學(xué)圖像的轉(zhuǎn)變，避免圖像特征不連貫的問題，并將多模態(tài)病灶信息有效整合，降低數(shù)據(jù)存儲(chǔ)空間和處理成本，進(jìn)而高效實(shí)現(xiàn)病灶精準(zhǔn)分割，較好輔助醫(yī)師預(yù)測(cè)并診斷病灶良惡性狀況和確定后續(xù)治療方案是未來超像素/體素算法的重要發(fā)展方向。