影像組學(xué)在腦膠質(zhì)瘤分級(jí)中的應(yīng)用研究

楊志煒，郭天，謝海濱，楊光

作者單位：華東師范大學(xué)物理系，上海市磁共振重點(diǎn)實(shí)驗(yàn)室，上海 200062

隨著生物醫(yī)學(xué)的進(jìn)步，人們對(duì)不同疾病的理解逐漸深入，但癌癥仍然是尚未攻克的難題[1]。腦膠質(zhì)瘤的臨床治療方案一般由治療過(guò)程中腫瘤的分級(jí)情況決定。目前，腦膠質(zhì)瘤分級(jí)標(biāo)準(zhǔn)嚴(yán)格基于病理學(xué)檢查，這種方法有一定的局限性，包括侵入性、采樣誤差、解釋的變化性等[2]，而基于影像學(xué)檢查的分級(jí)方法[3]沒(méi)有侵入性，逐步成為研究的熱點(diǎn)。但是，早期醫(yī)生往往只能通過(guò)肉眼來(lái)獲取圖像信息，耗時(shí)較長(zhǎng)，自動(dòng)化的分級(jí)方法則能有效縮短診斷與治療時(shí)間。磁共振成像常被用于腦膠質(zhì)瘤檢查，它能夠提供大量的腦部軟組織解剖、功能、代謝等信息，有些與疾病相關(guān)的特征，肉眼很難發(fā)現(xiàn)，這也需要人們開(kāi)發(fā)自動(dòng)的工具或算法[4-5]。因此越來(lái)越多的研究者參與到腦膠質(zhì)瘤磁共振圖像的自動(dòng)化分級(jí)中。

影像組學(xué)是近年來(lái)興起的一種圖像分析方法，它從醫(yī)學(xué)圖像中自動(dòng)地高通量地提取大量定量特征。假設(shè)對(duì)于給定模態(tài)的醫(yī)學(xué)圖像數(shù)據(jù)，通過(guò)自動(dòng)的或半自動(dòng)的軟件對(duì)它們進(jìn)行定量分析，可能可以從圖像中發(fā)掘出更多更好的信息[6]，就能將影像組學(xué)與常規(guī)磁共振圖像相結(jié)合，使其在腦膠質(zhì)瘤分級(jí)的應(yīng)用中發(fā)揮巨大作用。本研究用影像組學(xué)的方法對(duì)基于常規(guī)磁共振圖像的腦膠質(zhì)瘤進(jìn)行了分級(jí)。

1 材料與方法

1.1 數(shù)據(jù)

1.1.1 腦膠質(zhì)瘤分級(jí)

根據(jù)2016年世界衛(wèi)生組織(World Health Organization，WHO)中樞神經(jīng)系統(tǒng)腫瘤分類(lèi)標(biāo)準(zhǔn)，腦膠質(zhì)瘤分類(lèi)包括[7]星形細(xì)胞腫瘤、少突膠質(zhì)細(xì)胞腫瘤、室管膜腫瘤、脈絡(luò)叢腫瘤、神經(jīng)元及混合神經(jīng)元神經(jīng)膠質(zhì)腫瘤等。腦膠質(zhì)瘤可以分級(jí)為低評(píng)級(jí)(low grade glioma，LGG)和高評(píng)級(jí)(high grade glioma，HGG)。

1.1.2 數(shù)據(jù)集介紹

本研究使用了2017腦膠質(zhì)瘤分割競(jìng)賽數(shù)據(jù)集(brain tumor segmentation challenge 2017，BRATS2017)[8-9]，該數(shù)據(jù)集提供了專(zhuān)家手工標(biāo)記的腫瘤區(qū)域以及腫瘤的分級(jí)情況，其中提供腫瘤分級(jí)信息的有285例數(shù)據(jù)，LGG有75例，HGG有210例，數(shù)據(jù)模態(tài)為T(mén)1加權(quán)像、增強(qiáng)T1加權(quán)像、T2加權(quán)像和液體衰減反轉(zhuǎn)恢復(fù)序列(fluid attenuated inversion recovery，F(xiàn)LAIR)像。數(shù)據(jù)來(lái)源于19個(gè)不同機(jī)構(gòu)，大多采集自3.0 T磁共振成像系統(tǒng)。數(shù)據(jù)全部經(jīng)過(guò)預(yù)處理的，包括配準(zhǔn)到相同的解剖模板、插值到相同分辨率等。

所有圖像根據(jù)相同的標(biāo)準(zhǔn)由1～4個(gè)人進(jìn)行分割，并由經(jīng)驗(yàn)豐富的神經(jīng)放射學(xué)家驗(yàn)證?？偣卜譃?個(gè)區(qū)域，即：增強(qiáng)的核心(enhanced core)、水腫(edema)、非增強(qiáng)的實(shí)性核心(non-enhancing solid core)和壞死/囊性核心(necrotic/cystic core)，區(qū)域之間沒(méi)有重疊。數(shù)據(jù)集中將壞死/囊性核心與非增強(qiáng)的實(shí)性核心合并為一個(gè)區(qū)域，簡(jiǎn)稱(chēng)為NET區(qū)域，以感興趣區(qū)(region of interest，ROI)NET表示。增強(qiáng)的核心區(qū)域和水腫區(qū)域分別簡(jiǎn)稱(chēng)為ET區(qū)域和ED區(qū)域，以ROIET、ROIED表示。圖1為BRATS2017數(shù)據(jù)集中一例患者的膠質(zhì)瘤示意圖，藍(lán)色為T(mén)2加權(quán)像上顯示的整個(gè)腫瘤區(qū)域，即NET區(qū)域、ET區(qū)域、ED區(qū)域的總和；紅色為增強(qiáng)T1加權(quán)像上顯示的ET區(qū)域，綠色為FLAIR像上顯示的水腫區(qū)域，即ET區(qū)域與ED區(qū)域的總和。

1.2 特征提取

本研究自動(dòng)提取了117個(gè)特征，特征類(lèi)型包括全局特征、形態(tài)特征、灰度特征和Gabor特征。以下簡(jiǎn)要介紹各類(lèi)提取的特征。

1.2.1 全局特征

不同機(jī)器、環(huán)境、患者掃描得到的磁共振圖像具有不同的灰度分布。有些研究[2，10]利用直方圖匹配方法將不同圖像匹配到同一灰度分布下，但是這種方法會(huì)改變腫瘤區(qū)域的灰度分布，進(jìn)而丟失腫瘤的灰度信息。筆者直接將患者大腦的特征作為腫瘤分級(jí)所用的特征的一部分，這些特征稱(chēng)為全局特征。

全局特征包括腦脊液、白質(zhì)、灰質(zhì)在T1、T1增強(qiáng)、T2和FLAIR像上的灰度平均值。要獲取各組織的灰度值需要對(duì)大腦組織進(jìn)行分割，現(xiàn)有的最精確的分割方法是基于深度學(xué)習(xí)的[11]。另外，Kumar等[12]比較了不同的圖像分割算法，指出K-means算法具有較好的分割表現(xiàn)和較低的計(jì)算開(kāi)銷(xiāo)，與其他算法結(jié)合能夠提升分割效能，考慮到特征提取步驟需要對(duì)精確性與效率進(jìn)行權(quán)衡，因此，針對(duì)該文的應(yīng)用場(chǎng)景，采用K-means算法。

1.2.2 形態(tài)特征

形態(tài)特征包括腫瘤總體積V、腫瘤表面積A、腫瘤似圓度SPH[13]、腫瘤表面徑向分布的熵值RE、ET區(qū)域體積VET、ED區(qū)域體積VED以及NET區(qū)域體積VNET。另外，本研究還提取了VET與V的比值rET、VNET與V的比值rNET以及VED與V的比值rED。

1.2.3 灰度特征

本研究分別計(jì)算了T1、T1增強(qiáng)、T2、Flair圖像中ROINET、ROIET、ROIED的灰度平均值、標(biāo)準(zhǔn)差，其中最終的平均值為除以白質(zhì)的灰度平均值之后得到的平均值。由于腫瘤區(qū)域的灰度一般不滿(mǎn)足高斯分布，所以筆者還統(tǒng)計(jì)了上述3個(gè)區(qū)域?qū)?yīng)圖像灰度分布的峰度與偏度。峰度表征某一分布縱向偏離高斯分布的程度，偏度表征某一分布橫向偏離高斯分布的程度[14]。與上節(jié)相仿，以的形式表示灰度特征，其中A為NET、ET、ED，分別表示3類(lèi)區(qū)域，B為μ、σ、K和SK，分別表示平均值、標(biāo)準(zhǔn)差、峰度和偏度，C為T(mén)1、T2、T1CE、FLAIR。

1.2.4 Gabor特征

本研究對(duì)ROI所處的每一層都用二維蓋波濾波器進(jìn)行卷積，最后對(duì)ROI內(nèi)體素信號(hào)進(jìn)行平均，得到三維圖像的Gabor特征[15]。二維蓋波濾波器在位置x，y處的公式為：

λ為波長(zhǎng)，θ為方向角度，γ為縱橫比，它決定了卷積核的離心率[16]，本研究取1，σ為高斯函數(shù)的帶寬，它決定了像素不同領(lǐng)域范圍的權(quán)重。上式中的濾波器為復(fù)數(shù)形式，因此卷積得到的是復(fù)圖像，筆者取模圖來(lái)計(jì)算Gabor特征。θ取0到

之

本研究?jī)H計(jì)算T1增強(qiáng)像中腫瘤區(qū)域與Flair像中水腫區(qū)域的Gabor特征，以的形式表示，其中A為T(mén)UMOR、ED，分別表示腫瘤區(qū)域和水腫，B取1～5，表示5個(gè)方向角度，C為T(mén)1CE、FLAIR，D取1～5，表示5個(gè)波長(zhǎng)。

1.3 特征選擇

特征選擇方法主要包括：過(guò)濾(Filter)方法、包裝(Wrapper)方法和集成(Ensemble)方法[17]。本文結(jié)合了前兩種特征選擇方法進(jìn)行特征選擇。具體方法為：(1)使用過(guò)濾方法對(duì)特征進(jìn)行排序；(2)添加排名靠前的特征作為新的特征子集，運(yùn)用包裝方法進(jìn)行特征選擇。

1.3.1 過(guò)濾方法與t檢驗(yàn)

過(guò)濾方法通常從特征集合中剔除與預(yù)測(cè)變量最不相關(guān)的特征，將剩余的變量用于預(yù)測(cè)。它所用的判據(jù)往往是通用的與預(yù)測(cè)變量的相關(guān)性之類(lèi)的判據(jù)，與后續(xù)所使用的分類(lèi)器模型無(wú)關(guān)。過(guò)濾方法得到的特征子集常常會(huì)包含冗余的特征，因?yàn)樗玫呐袚?jù)忽略了特征間的相關(guān)性。由于過(guò)濾方法與特定的分類(lèi)器無(wú)關(guān)，所以計(jì)算開(kāi)銷(xiāo)小，泛化能力強(qiáng)。對(duì)于大數(shù)據(jù)集和高維度特征，過(guò)濾方法應(yīng)用十分廣泛。

本研究使用雙總體t檢驗(yàn)值的大小作為特征排序規(guī)則。該類(lèi)t檢驗(yàn)?zāi)軌蚝饬績(jī)蓚€(gè)分布的差異，因此也能用來(lái)檢驗(yàn)每個(gè)特征的分類(lèi)能力。假設(shè)特征在不同分類(lèi)下均服從高斯分布，并且它們的方差相等，本研究利用t檢驗(yàn)選擇在兩個(gè)分類(lèi)間分布差異最大的30個(gè)特征。

1.3.2 包裝方法與遺傳算法

包裝方法對(duì)特征子集的性能進(jìn)行評(píng)價(jià)，它可以發(fā)現(xiàn)特征間的關(guān)系，選出最優(yōu)子集。包裝方法分類(lèi)能力很強(qiáng)，但計(jì)算開(kāi)銷(xiāo)大，容易過(guò)擬合，限制了該方法在高維問(wèn)題上的應(yīng)用。本研究采用了遺傳算法來(lái)搜索特征子集空間。

遺傳算法(genetic algorithm，GA)[18]是模擬生物在自然環(huán)境中的遺傳和進(jìn)化過(guò)程而形成的一種自適應(yīng)全局優(yōu)化概率搜索算法。對(duì)于函數(shù)最優(yōu)化的問(wèn)題，解的空間一般非常大，在很多情況下幾乎求不出解析解。遺傳算法能夠利用其隨機(jī)搜索的特性，求出其近似最優(yōu)解。

進(jìn)行優(yōu)化時(shí)，首先設(shè)定初始群體，對(duì)群體中的所有個(gè)體進(jìn)行評(píng)價(jià)，即計(jì)算適應(yīng)值。然后將一系列遺傳算子作用到群體上，判斷是否達(dá)到算法終止條件，若是，將群體中具有最優(yōu)適應(yīng)值的個(gè)體作為問(wèn)題的最優(yōu)解，若否，則繼續(xù)迭代。遺傳算法的典型流程如圖2所示。

遺傳算子主要包含選擇算子、交叉算子和變異算子[19]。選擇運(yùn)算將群體中適應(yīng)度函數(shù)值較優(yōu)的個(gè)體以某種規(guī)則遺傳到下一代中。交叉運(yùn)算以一定概率交換群體兩個(gè)個(gè)體間的部分染色體，產(chǎn)生新的個(gè)體。變異運(yùn)算以一定概率改變?nèi)后w中某些個(gè)體的一個(gè)或者多個(gè)等位基因。

本研究最終選取群體人口為50，最大進(jìn)化的代數(shù)為100，選擇算法為錦標(biāo)賽算法，選擇個(gè)數(shù)為2，即每次選取2個(gè)個(gè)體中最好的那個(gè)，共選取50次。交叉算法為分散交叉，交叉率為0.8，變異算法為均勻變異，變異率為0.1。

1.4 分類(lèi)

本研究使用過(guò)濾特征選擇方法，對(duì)非線性的SVM[20]、K近鄰法(K nearest neighbour，KNN)[21]和決策樹(shù)(decision tree，DTREES)[16]進(jìn)行了比較，步驟如下：從最具有區(qū)分度的特征開(kāi)始，逐步添加次一級(jí)區(qū)分度的特征，找到分類(lèi)性能最好的特征組合。

由于數(shù)據(jù)是不平衡的，需要對(duì)數(shù)據(jù)進(jìn)行平衡。常用的方法有對(duì)數(shù)量多的分類(lèi)樣本下采樣或?qū)?shù)量小的分類(lèi)樣本上采樣，前者不能充分利用數(shù)據(jù)，后者可能會(huì)增加噪聲點(diǎn)。根據(jù)每種分類(lèi)算法的特性來(lái)分別解決這個(gè)問(wèn)題。支持向量機(jī)中一個(gè)重要的參數(shù)是懲罰因子C。較大的C值會(huì)給誤分類(lèi)一個(gè)較大的懲罰。本研究還采用帶權(quán)重的支持向量機(jī)[22]來(lái)加大對(duì)小樣本分類(lèi)錯(cuò)誤的懲罰。在K近鄰法中，使用帶距離權(quán)重的分類(lèi)判決方法，這使得距離較近的點(diǎn)對(duì)分類(lèi)具有較大貢獻(xiàn)。對(duì)于決策樹(shù)，則使用了對(duì)數(shù)據(jù)不平衡不敏感的CART算法[23]。1.4.1 SVM

SVM分類(lèi)算法的主要思想是找到一個(gè)獨(dú)一無(wú)二的超平面，這個(gè)超平面能最大化兩個(gè)分類(lèi)之間的間隔。給出l組訓(xùn)練數(shù)據(jù)，x為輸入樣本，y為標(biāo)簽。

軟間隔的SVM方法需要得到下列優(yōu)化問(wèn)題的解：

約束條件為：

訓(xùn)練集的特征矢量被映射到高維空間中，C是懲罰因子，用來(lái)平衡間隔最大化的要求與分類(lèi)錯(cuò)誤之間的矛盾，取C=0.5。若對(duì)不同標(biāo)簽的樣本進(jìn)行不同懲罰，則稱(chēng)為帶權(quán)重的，這里取正類(lèi)負(fù)類(lèi)的權(quán)重比為1∶1.5。該優(yōu)化問(wèn)題的一種簡(jiǎn)便解法是解它的對(duì)偶問(wèn)題，即：

約束條件為：

最后給出一個(gè)數(shù)據(jù)x，得到該數(shù)據(jù)的分類(lèi)為，令，則函數(shù)K稱(chēng)為核函數(shù)，這里取高斯核函數(shù)為，，n為特征數(shù)量。

1.4.2 K近鄰法

KNN算法常用來(lái)解決分類(lèi)問(wèn)題，它沒(méi)有顯式的學(xué)習(xí)過(guò)程。典型的KNN算法步驟如下：給定一個(gè)訓(xùn)練數(shù)據(jù)集，對(duì)新的輸入實(shí)例，在訓(xùn)練數(shù)據(jù)集中找到與該實(shí)例距離最相近的K個(gè)實(shí)例，根據(jù)分類(lèi)決策規(guī)則將該實(shí)例分為某個(gè)類(lèi)。該算法有3個(gè)基本要素[24]：距離度量、k值選擇和分類(lèi)決策規(guī)則。常用的分類(lèi)決策規(guī)則有多數(shù)表決，即K個(gè)實(shí)例中較多的那一類(lèi)為最終分類(lèi)。這里采用歐氏距離度量，這個(gè)距離還能帶上權(quán)重，取權(quán)重為。較小的K值會(huì)使模型變復(fù)雜容易過(guò)擬合，而較大的K值會(huì)增大分類(lèi)誤差，這里取K為3。

1.4.3 決策樹(shù)

決策樹(shù)模型是定義在特征空間與類(lèi)空間上的條件概率分布，主要優(yōu)點(diǎn)有：模型具有可讀性，分類(lèi)速度快[24]。決策樹(shù)算法學(xué)習(xí)一般包含3個(gè)步驟：特征選擇、決策樹(shù)生成和決策樹(shù)修剪。本文決策樹(shù)是基于分類(lèi)回歸樹(shù)(classification and regression tree，CART)算法的，該算法通過(guò)遞歸地構(gòu)建二叉樹(shù)來(lái)生成決策樹(shù)。它使用基尼(Gini)指數(shù)選擇最優(yōu)特征，同時(shí)決定該特征的最優(yōu)二值切分點(diǎn)。對(duì)于給定樣本集合D，基尼指數(shù)定義為：

K為標(biāo)簽個(gè)數(shù)，當(dāng)K=2時(shí)即為二分類(lèi)問(wèn)題， 標(biāo)簽為k的樣本個(gè)數(shù)。根據(jù)特征A的取值分割成D1、D2兩部分，則在該條件下基尼指數(shù)定義為：

決策樹(shù)修剪是通過(guò)剪枝算法完成的，算法從完全生長(zhǎng)的決策樹(shù)底端減去一些子樹(shù)，分為兩步[24]：首先從生成算法的決策樹(shù)底端開(kāi)始不斷剪枝，直到根節(jié)點(diǎn)，形成子數(shù)序列，然后通過(guò)交叉驗(yàn)證法在獨(dú)立的驗(yàn)證數(shù)據(jù)集上測(cè)試，選擇最優(yōu)子樹(shù)。

表1 不同方法分類(lèi)的準(zhǔn)確率、敏感性、特異性Tab.1 Classification accuracy, sensitivity, and specificity of various methods

表2 不同方法特征的選擇結(jié)果Tab.2 The feature selection results of various methods

圖1 BRATS2017膠質(zhì)瘤示意圖。藍(lán)色為腫瘤區(qū)域，紅色為增強(qiáng)的核心區(qū)域，綠色為水腫區(qū)域Fig. 1 The segmentation of BRATS2017 dataset. Blue: the whole tumor, Red: the enhanced core,Green: edema.

圖2 遺傳算法步驟示意圖Fig. 2 The procedure of the genetic algorithm.

2 結(jié)果

2.1 基于Filter方法的分類(lèi)

本研究用留一交叉驗(yàn)證法對(duì)BRATS2017數(shù)據(jù)集的包含分級(jí)信息的285例數(shù)據(jù)進(jìn)行分類(lèi)。首先，使用基于過(guò)濾方法的特征選擇方法，特征排序方法為t檢驗(yàn)法。選出前30個(gè)特征依次添加，并用支持向量機(jī)、K近鄰法和決策樹(shù)進(jìn)行分類(lèi)。本研究使用了LIBSVM工具箱[25]來(lái)進(jìn)行SVM算法分類(lèi)。圖1為支持向量機(jī)、K近鄰法、決策樹(shù)的分類(lèi)準(zhǔn)確率與特征數(shù)量的關(guān)系圖。

筆者確定3種方法達(dá)到最高準(zhǔn)確率時(shí)的特征子集，計(jì)算對(duì)應(yīng)的敏感性和特異性，見(jiàn)表1。由表1可知，支持向量機(jī)算法具有最高的準(zhǔn)確率與敏感性，所以選擇該算法作為下一步實(shí)驗(yàn)的分類(lèi)器。

2.2 基于遺傳算法特征選擇的分類(lèi)

本研究以過(guò)濾方法排序得到的前30個(gè)特征作為初值，進(jìn)一步使用基于遺傳算法的包裝方法進(jìn)行特征選擇。選擇上節(jié)中選擇的SVM作為分類(lèi)器。圖4為每一代適應(yīng)度函數(shù)的平均值、最大值與進(jìn)化代數(shù)的關(guān)系圖。圖5為支持向量機(jī)分類(lèi)的受試者工作特征曲線(receiver operating characteristic，ROC)圖，黑線為基于過(guò)濾方法的，紅線為基于遺傳算法的。

由圖4可以看出，算法迭代結(jié)束時(shí)最小函數(shù)值為-93.33，即達(dá)到的最大準(zhǔn)確率為93.33%，證明遺傳算法搜索得到的特征組合效果最好，優(yōu)于過(guò)濾方法的結(jié)果。

由圖5可以看出，GA-SVM方法ROC曲線下的面積(area under the curve，AUC)為0.94，大于Filter-SVM方法的0.90，證明具有較好的分類(lèi)表現(xiàn)，該方法在同等敏感性下有著更低的假陽(yáng)率，即達(dá)到相同效果所付出的代價(jià)更低。

2.3 特征選擇結(jié)果

表2為上兩節(jié)中4種方法選取得到的特征。由表2可知，基于過(guò)濾特征選擇方法的支持向量機(jī)與決策樹(shù)算法挑選出了同樣的特征，而K近鄰法選出了較多特征，但是特征之間有冗余。GA-SVM方法則在這些特征的基礎(chǔ)上又選擇出了Gabor特征和全局特征。

圖3 分類(lèi)準(zhǔn)確率與添加特征數(shù)量關(guān)系圖。決策樹(shù)分類(lèi)算法隨著特征數(shù)增加準(zhǔn)確率趨于平穩(wěn)。K近鄰分類(lèi)算法在特征數(shù)很少時(shí)準(zhǔn)確率就達(dá)到最大，隨后趨于平穩(wěn)。支持向量機(jī)分類(lèi)算法能達(dá)到的準(zhǔn)確率最高，但是隨著特征數(shù)增加，準(zhǔn)確率急劇下降 圖4 適應(yīng)度函數(shù)平均值及最小值隨遺傳代數(shù)變化的情況 圖5 ROC曲線圖。Filter特征選擇方法的ROC曲線下面積(area under the curve，AUC)=0.90?；赪rapper與遺傳算法的特征選擇方法的AUC=0.94Fig. 3 Classification accuracy vs number of features used. Accuracy of decision tree is quite stable with the increase of the number of features. Accuracy of KNN reaches the maximum when the number of features is small. Accuracy of SVM is the highest among three classifiers, but it rapidly falls off when the number of features increases. Fig.4 The best or mean fitness function value vs generations. Fig.5 ROC curve. Filter feature selection method, AUC=0.90. Feature selection method based on Wrapper and genetic algorithm, AUC=0.94.

3 討論

3.1 腦膠質(zhì)瘤分級(jí)的自動(dòng)化

本文利用影像組學(xué)方法對(duì)腦膠質(zhì)瘤進(jìn)行了自動(dòng)分級(jí)。腦膠質(zhì)瘤的自動(dòng)分級(jí)一般有如下步驟：(1)圖像分割；(2)特征提取；(3)特征選擇；(4)分類(lèi)。最常用的特征提取方法有離散小波變換(discrete wavelet transform，DWT)和灰度共生矩陣(gray level co-occurrence matrix，GLCM)；在特征降維方面，主成分分析法(principal component analysis，PCA)和遺傳算法被廣泛使用[26]；隨機(jī)森林(random forests，RF)由于其算法特性，更適合處理高維問(wèn)題和多分類(lèi)問(wèn)題，因此該算法對(duì)于腦膠質(zhì)瘤分級(jí)最精確[8]。本文則提取了多種類(lèi)型的特征包括全局、形態(tài)、灰度、紋理各個(gè)方面；并結(jié)合兩種特征選擇方法，進(jìn)行了多步驟的特征降維，由于是單分類(lèi)問(wèn)題，本研究只比較了3種常見(jiàn)的分類(lèi)算法。筆者首先采用基于t檢驗(yàn)的過(guò)濾方法進(jìn)行特征選擇，發(fā)現(xiàn)支持向量機(jī)分類(lèi)算法能達(dá)到的準(zhǔn)確率和敏感性最高，但是隨著特征數(shù)增加，準(zhǔn)確率急劇下降。支持向量機(jī)算法具有較好性能的原因是：該算法的高斯核函數(shù)將特征映射到了高維空間，提高了特征的分類(lèi)能力。結(jié)果表明，支持向量機(jī)算法適合作為腦膠質(zhì)瘤分級(jí)問(wèn)題中的分類(lèi)器。

然后，采用基于遺傳算法的包裝方法進(jìn)行特征選擇，算法得到的準(zhǔn)確率和AUC值都超過(guò)了過(guò)濾方法，結(jié)果表明，遺傳算法挑選出了特定的特征組合，提升了腦膠質(zhì)瘤的分級(jí)診斷效能。包裝方法能夠搜索更大的特征空間，但是需要更長(zhǎng)的搜索時(shí)間，而遺傳算法是一種自適應(yīng)的全局優(yōu)化概率搜索算法，彌補(bǔ)了該方法的缺點(diǎn)，加快了搜索過(guò)程。

3.2 影像組學(xué)在常規(guī)磁共振圖像上的應(yīng)用

基于影像組學(xué)的方法對(duì)腦膠質(zhì)瘤分級(jí)進(jìn)行的研究相對(duì)較少，值得注意的是，鄧慧媛等[27]利用Logistic多元回歸相關(guān)性分析篩選出了鑒別高低級(jí)別腦膠質(zhì)瘤效能較高的DWI序列特征，高于T1加權(quán)增強(qiáng)像、T2加權(quán)像的鑒別效能。另外，林坤等[28]的研究得出多模態(tài)MRI檢查技術(shù)的分級(jí)診斷效能高于常規(guī)增強(qiáng)MRI及DWI檢查。本文則結(jié)合公開(kāi)數(shù)據(jù)集的大數(shù)據(jù)量和影像組學(xué)發(fā)掘了常規(guī)磁共振圖像在腦膠質(zhì)瘤分級(jí)問(wèn)題上的巨大潛力。

在特征選擇方面，基于過(guò)濾特征選擇方法的支持向量機(jī)與決策樹(shù)挑選出了同樣的特征，但是它們的分類(lèi)性能不同；K近鄰法選出了較多特征，但是特征之間有冗余。這3種方法選出的特征基本與先驗(yàn)知識(shí)相符，即腫瘤越大、腫瘤區(qū)域灰度值越高，腫瘤惡性可能性越大?；谶z傳算法的特征選擇方法則在這些特征的基礎(chǔ)上又選擇出了Gabor特征和全局特征。這些特征增強(qiáng)了分類(lèi)性能，而且人的肉眼無(wú)法觀察到，這凸顯了影像組學(xué)方法與人相比可能具有的優(yōu)勢(shì)。本研究從T1、T1增強(qiáng)和FLAIR像中都提取到了有用的特征，這也表現(xiàn)了常規(guī)磁共振圖像在腦膠質(zhì)瘤分級(jí)中不可忽視的作用。

3.3 本文工作的不足與局限性

目前，影像學(xué)檢查經(jīng)常使用于腦膠質(zhì)瘤的分級(jí)診斷，常規(guī)磁共振圖像對(duì)于病變的解剖定位意義重大，但是不能同時(shí)顯示腫瘤的代謝、血液動(dòng)力學(xué)、分子及細(xì)胞結(jié)構(gòu)的信息[29]。這也提示需要克服常規(guī)磁共振圖像的局限性，將多種MR技術(shù)相結(jié)合，彼此取長(zhǎng)補(bǔ)短，可更好地對(duì)腦膠質(zhì)瘤進(jìn)行分級(jí)診斷。

另外，本研究所使用的公開(kāi)數(shù)據(jù)集中已經(jīng)提供腦膠質(zhì)瘤的分割結(jié)果。在腦膠質(zhì)瘤分級(jí)的實(shí)際應(yīng)用中，還需要在ROI區(qū)域?qū)δ[瘤進(jìn)行分割，在分割結(jié)果上進(jìn)行特征提取。因此，在之后的工作中，對(duì)腦膠質(zhì)瘤的分割方法進(jìn)行研究也是重點(diǎn)之一。

總之，本研究對(duì)基于常規(guī)磁共振圖像的腦膠質(zhì)瘤分級(jí)進(jìn)行了研究，實(shí)驗(yàn)結(jié)果表明，基于常規(guī)磁共振圖像，利用影像組學(xué)方法，選擇合適的特征組合，可以對(duì)腦膠質(zhì)瘤進(jìn)行自動(dòng)分級(jí)。通過(guò)對(duì)分類(lèi)方法的比較，發(fā)現(xiàn)非線性支持向量機(jī)對(duì)腦部腫瘤分級(jí)效果最好，在特征篩選方面，筆者提出先利用過(guò)濾方法篩選出部分特征，然后在過(guò)濾的基礎(chǔ)上進(jìn)一步使用包裝方法，用遺傳算法選擇了最后的特征子集。