張大偉，王海寶，方 杰，李之豪，徐麗艷，余長亮，王萬勤，伏樹奇，李曉舒

甲狀腺微小結節(jié)(直徑≤1.0 cm)在臨床上基本無明顯癥狀，容易造成漏診或誤診。如何鑒別甲狀腺微小結節(jié)，特別是對甲狀腺微小乳頭狀癌(papillary microcarcinoma of thyroid,PTMC)的早診早治，對改善患者的預后有重要臨床價值。目前影像學檢查在診斷甲狀腺微小結節(jié)時，主要通過觀察結節(jié)形態(tài)學如：形態(tài)、大小、鈣化、邊緣及有無頸部和遠處淋巴結腫大等基本特征進行主觀判斷結節(jié)的良惡性。影像組學是通過計算機數(shù)學檢測方法從圖像中提取圖像數(shù)字化紋理特征參數(shù)，并進行定性分析診斷。既往研究報道，CT平掃的影像組學應用于良惡性甲狀腺大結節(jié)(直徑>1.0 cm)具有良好的鑒別意義。目前，對于甲狀腺微小結節(jié)(直徑≤1.0 cm)的良惡性鑒別診斷CT影像組學的研究報道較少。該研究探討CT平掃結合增強的影像組學在甲狀腺微小結節(jié)良惡性鑒別診斷中的價值。

1 材料與方法

1.1 病例資料

回顧性收集2016年1月—2019年10月在安徽醫(yī)科大學第一附屬醫(yī)院及安慶市第一人民醫(yī)院術前行甲狀腺CT平掃及增強檢查甲狀腺結節(jié)患者203例，所有病例均經(jīng)手術或穿刺病理學證實，其中甲狀腺微小乳頭狀癌(PTMC)148例 ，男37例，女111例，年齡20～70(44.5±10.5)歲; 微小結節(jié)性甲狀腺腫(micro nodular goiter, MNG)55例，男11例，女44例，年齡20～77(53.9±11.0) 歲。納入標準: ① 甲狀腺微小結節(jié)(直徑≤1.0 cm)；② 同時有高質(zhì)量清晰 CT圖像、病理結果；③ 完整的原始圖像數(shù)據(jù)(DICOM)。排除標準：① 甲狀腺直徑>1.0 cm結節(jié)；② 亞急性甲狀腺炎；③ 偽影較多、質(zhì)量較差的圖像。

1.2 主要儀器與試劑

GE Discovery CT 750 HD 64排CT購自美國通用電氣醫(yī)療系統(tǒng)有限公司；西門子SOMATOM Definition AS+64排CT掃描儀器購自德國西門子醫(yī)療系統(tǒng)有限公司；碘海醇(濃度按碘計50 ml：17.5 g)購自江蘇揚子江藥業(yè)集團有限公司。

1.3 邏輯回歸等機器學習模型的建立

Logistic是一個分類算法并非回歸算法。通常是利用已知的自變量來預測一個離散型因變量的值，它的輸出值在0到1之間。隨機森林是基于機器學習中常見的決策樹概念，該方法對一組決策樹進行訓練，引入了兩個層次的隨機化。所有的決策樹都被訓練來預測所有的單個樹，并選擇最高頻率類別作為最終結果。支持向量機(SVM)：是一種常用的機器判讀方法。在機器學習領域，它是一種有監(jiān)督的學習模型，通常用于模式識別、分類和回歸分析。KNN(k-Nearest Neighbor)算法，可根據(jù)最近的一個或多個樣本的類別來確定待分類樣本的類別。決策樹(Tree)是基于樹的模型中最基礎的概念，它可用于解決分類或回歸問題,該模型可用于擬合人們的決策行為，因變量既可以是分類變量，也可以是連續(xù)變量，該研究應用分類樹進行預測。貝葉斯模型(Bayes)是運用貝葉斯統(tǒng)計進行的一種預測，該研究運用樸素貝葉斯統(tǒng)計模型進行預測。該研究對CT平掃及增強影像組學數(shù)據(jù)特征降維后，隨機分組(訓練組和測試組)并構建6種不同深度學習模型，包括Forest模型、 SVM模型、KNN模型、Tree模型、 Bayes模型及 Logistic模型，進一步運用模型對兩組特征數(shù)據(jù)進行統(tǒng)計分析和預測，以評估鑒別甲狀腺微小癌(PTMC)和微小結節(jié)性甲狀腺腫(MNG)的各影像組學模型效能，并獲得準確度、特異性、敏感度及AUC(曲線下面積ROC)。

1.4 影像組學特征提取、降維及建模、降維及組成模型分析

將選定患者CT平掃及增強DICOM 圖像數(shù)據(jù)通過影像組學智能分析軟件(AK軟件)重建(1.0 mm×1.0 mm×1.0 mm)并標準化處理，然后運用圖像處理軟件(ITK-SNAP)，由2名像診斷醫(yī)師對圖像進行共同處理，其中1名5年以上工作經(jīng)驗診斷醫(yī)師對甲狀腺微小結節(jié)進行診斷，并對感興趣區(qū)( region of interest，ROI)進行手動勾畫，運用三維模式勾畫病灶作為興趣區(qū)(ROI)，另1名更高年資診斷醫(yī)師復核完善。在勾畫ROI前認真閱片、觀察圖像質(zhì)量、剔除偽影圖像，平掃圖像獲得217個ROI，增強圖像獲得96個ROI，通過AK軟件共提取402個影像組學特征，包括一階特征和二階特征6類參數(shù)：直方圖參數(shù)、形態(tài)學矩陣、Haralick、灰度共生矩陣、灰度區(qū)域大小矩陣、游程矩陣參數(shù)，導入IPMS軟件進行統(tǒng)計分析，通過2次特征性降維，首先通過Univariate-Logistic方法分別降到207個、157個特征參數(shù)，通過Lasso L1 Regularized Logistic Regression方法再次特征降維，分別剩余20、19個特征參數(shù)。隨機分組為訓練組：選擇70%病例，測試組：剩余30%病例。并用這些影像組學特征組構建6種鑒別甲狀腺微小癌(PTMC)和微小結節(jié)性甲狀腺腫(MNG)的影像組學模型，分別為Forest模型、SVM模型、KNN模型、Tree模型、Bayes模型及Logistic模型并進行分類，進一步統(tǒng)計分析獲得準確度、特異性、敏感度及AUC。

2 結果

2.1 PTMC和MNG患者臨床基本信息

兩組患者年齡差異有統(tǒng)計學意義(

t

=5.559,

P

=0.000)，PTMC患者相對MNG患者明顯年輕；兩組性別差異無統(tǒng)計學意義，見表1。

表1 PTMC和MNG患者臨床基本信息

2.2 PTMC和MNG患者CT平掃及增強診斷結果

PTMC通常形態(tài)欠規(guī)則，邊緣中斷呈“咬餅征”，平掃呈低密度，增強呈不均勻強化(圖1A、B)；MNG病灶一般形態(tài)規(guī)則，邊緣光整，平掃多呈等密度，增強呈環(huán)形強化，也可不均勻強化(圖1C、D)。CT檢查正確診斷惡性結節(jié)112例，良性結節(jié) 30例，假惡性結節(jié)25例，假良性結節(jié)36例，診斷的準確度65.0%，特異性54.5%，敏感度81.8%。

圖1 甲狀腺結節(jié)A：甲狀腺左葉PTMC CT平掃；B：增強圖像(箭頭所示)；C：甲狀腺左葉MNG CT平掃；D：增強圖像(箭頭所示)

2.3 CT平掃影像組學分析結果

構建不同深度學習模型，進一步對CT平掃影像組學數(shù)據(jù)特征進行分類及統(tǒng)計分析評估，獲取不同組的準確度、特異性、敏感度及AUC(具體見方法部分)。Forest模型鑒別診斷效能最佳，見表2及圖2A。

表2 CT平掃影像組學不同模型鑒別診斷結果

2.4 CT增強影像組學分析結果

構建不同深度學習模型，進一步對CT增強影像組學數(shù)據(jù)特征進行分類及統(tǒng)計分析評估，獲取不同組的準確度、特異性、敏感度及AUC(具體見方法部分)。Forest模型鑒別診斷效能最佳。見表3及圖2B。

圖2 測試組Forest模型ROC曲線A：CT平掃 ；B：CT增強掃描

表3 CT增強影像組學不同模型鑒別診斷結果

3 討論

甲狀腺癌發(fā)病率從1975—2009年增加了3倍，其中甲狀腺微小乳頭狀癌(WHO定義直徑≤1.0 cm甲狀腺乳頭狀癌，不論有無淋巴結轉移)增加了49%。超聲結果顯示甲狀腺結節(jié)直徑低至2 mm，是目前首選診斷甲狀腺微小結節(jié)的影像檢查，但其受診斷醫(yī)師診斷水平等主觀因素影響較大。CT顯示結節(jié)內(nèi)砂礫樣鈣化及頸部淋巴結腫大高度提示惡性。MNG大多形態(tài)規(guī)則，邊緣光整，可有粗大鈣化，環(huán)形強化等特征，PTMC通常形態(tài)欠規(guī)則，邊緣中斷呈“咬餅征”，可有砂爍樣鈣化，中心強化。

荷蘭學者在2012年首次提出影像組學的概念。通過計算機技術從影像圖像中提取數(shù)字化紋理特征參數(shù)，量化評估病變圖像異質(zhì)性，已經(jīng)運用于腫瘤診斷、分期、治療和預后評估。既往研究表明甲狀腺惡性結節(jié)熵值比良性結節(jié)高，差異有統(tǒng)計學意義。Yu et al對甲狀腺結節(jié)超聲圖像紋理分析研究表明SVM及神經(jīng)網(wǎng)絡模型對甲狀腺良惡性結節(jié)鑒別診斷的準確率達99%～100%。杜丹丹 等運用CT影像組學聯(lián)合SVM模型研究表明，對大于1 cm的甲狀腺腺瘤及乳頭狀癌結節(jié)鑒別診斷的準確度、特異性及敏感度顯著提高。該研究對應用CT影像組學方法鑒別甲狀腺微小結節(jié)良惡性的可行性進行初步探索，構建Forest模型、 SVM模型 、KNN模型、Tree模型、Bayes模型及Logistic模型等6種影像組學模型用于鑒別甲狀腺微小乳頭狀癌和微小結節(jié)性甲狀腺腫，CT平掃及增強影像組學每個模型對甲狀腺良惡性微小結節(jié)鑒別診斷的準確度、特異性、敏感度均高于常規(guī)增強CT研究結果。CT 平掃及增強影像組學Forest模型對PTMC和MNG具有較高的鑒別價值。該研究CT平掃的準確度、特異性及敏感度與杜丹丹 等運用平掃CT平掃影像組學聯(lián)合SVM模型研究結果相比稍低，但較常規(guī)CT影像診斷的準確率有所提高?？赡茉颍篗NG是上皮濾泡增生與正常甲狀腺組織差異不大，同時MNG及PTMC病灶較小，少有囊變、壞死，密度分辨率?。籑NG以纖維、膠質(zhì)等成分為主，這些成分取代了大部分的毛細血管床，呈較低強化方式；而甲狀腺微小乳頭狀癌內(nèi)含纖維成分較多，但含相對較少的膠質(zhì)等成分，占據(jù)相應較少毛細血管床結構，強化程度高于MNG。