林帆， 胡若凡， 梁超， 余娟， 劉俠靜， 雷益

·乳腺影像學(xué)·

磁共振紋理動(dòng)態(tài)特征鑒別乳腺良惡性腫塊

林帆， 胡若凡， 梁超， 余娟， 劉俠靜， 雷益

目的提取乳腺病灶的時(shí)空變化特征作為新的DCE-MRI標(biāo)記(稱為紋理動(dòng)態(tài)特征)并證明其鑒別良惡性腫塊的能力。方法回顧性分析52個(gè)乳腺腫塊，其中惡性腫瘤30個(gè)，良性腫塊22個(gè)，提取并對(duì)動(dòng)態(tài)特征信號(hào)強(qiáng)度特征、紋理特征、形態(tài)特征、邊緣特征進(jìn)行分組。為了更好評(píng)估這些特征，采用不同的特征類建立分組模型，計(jì)算正確率，敏感度，特異性及曲線下面積(AUC)。結(jié)果結(jié)合紋理動(dòng)態(tài)特征所建立的良惡性腫瘤分類器具有最大的AUC=0.94，準(zhǔn)確率90%，敏感度92%，特異性85%，優(yōu)于其他各組分類器，與信號(hào)強(qiáng)度特征所建立的模型差異有統(tǒng)計(jì)學(xué)意義(AUC=0.80,P<0.05)。結(jié)論磁共振紋理動(dòng)態(tài)特征有助于鑒別良惡性腫塊，甚至優(yōu)于臨床上最流行的DCE-MRI標(biāo)記信號(hào)強(qiáng)度動(dòng)態(tài)特征。

磁共振成像； 紋理； 乳腺腫瘤； 人工智能

乳腺動(dòng)態(tài)磁共振增強(qiáng)掃描(dynamic contrast-enhanced magnetic resonance imaging，DCE)已成為臨床常規(guī)檢查手段，通過測(cè)量信號(hào)強(qiáng)度的動(dòng)態(tài)變化(如時(shí)間-信號(hào)強(qiáng)度曲線)來分析病灶的血流動(dòng)力學(xué)[1-3]。然而根據(jù)乳腺BI-BADS分級(jí)的要求和實(shí)際工作經(jīng)驗(yàn)，除了信號(hào)強(qiáng)度以外，病灶內(nèi)部紋理特征、邊緣特征等隨時(shí)間的變化同樣重要[4]。動(dòng)態(tài)增強(qiáng)圖像去分析這些特征，耗費(fèi)時(shí)間和精力[5-10]，并容易出錯(cuò)。

因此，本文擬采取計(jì)算機(jī)輔助的方法解決問題，通過圖像處理算法提取乳腺病灶內(nèi)部及邊緣等的動(dòng)態(tài)變化特征，然后基于機(jī)器學(xué)習(xí)的方法進(jìn)行良惡性樣本訓(xùn)練，篩選出診斷良惡性病變最有價(jià)值的特征，從而成為DCE-MRI新的標(biāo)記。

材料與方法

1.研究對(duì)象

回顧分析2013年-2017年的乳腺動(dòng)態(tài)增強(qiáng)MRI病例52例。納入標(biāo)準(zhǔn)：乳腺M(fèi)RI上顯示腫塊病灶，經(jīng)過手術(shù)或者穿刺病理證實(shí)惡性腫瘤30個(gè)(均為浸潤性導(dǎo)管癌)，良性腫塊22個(gè)(纖維腺瘤19例，腺病3例)。

2.掃描參數(shù)

掃描機(jī)器為西門子Avanto 1.5T磁共振，使用乳腺專用breast線圈掃描。行T1WI,T2WI平掃后，使用3D GRE序列行動(dòng)態(tài)增強(qiáng)掃描：TR 4.5 ms，TE 1.2 ms， 翻轉(zhuǎn)角10°，層厚1.2 mm，層間隔0 mm，視野280 mm×280 mm，矩陣384×257，激勵(lì)次數(shù)1.0，掃描次數(shù)7次，第一次為蒙片，最后再掃描一次矢狀面的增強(qiáng)及其他序列圖像。乳腺增強(qiáng)對(duì)比劑是對(duì)比劑釓噴酸葡胺通過肘靜脈以0.1 mol/kg的劑量、2.0 mL/s的流率注入，注射完畢后加注20 mL生理鹽水。

3.病灶分割

由1名放射科12年年資的副主任醫(yī)生通過計(jì)算機(jī)進(jìn)行，先在動(dòng)態(tài)增強(qiáng)序列上選擇病灶邊界最明顯的時(shí)相，分別在3個(gè)方位均取病灶徑線最大的層面，在自行編寫的半自動(dòng)的程序下勾畫病灶的邊界。

4.提取病灶特征

通過圖像處理算法提取病灶內(nèi)部各種特征[11-12，17]。對(duì)不同時(shí)相(t0, t1，t2，t3，t4，t5，t6)的病灶提取紋理特征。

①形態(tài)學(xué)特征：a.病灶與最小外接圓的面積比；b.平均半徑比；c.半徑標(biāo)準(zhǔn)差；d.光滑性；e.緊湊性。

②病灶內(nèi)部特征：a.一階統(tǒng)計(jì)量(平均值，標(biāo)①準(zhǔn)差，變異系數(shù) )；b.Haralick紋理(4個(gè)角度)；c.Gabor濾波參數(shù)(8個(gè)角度，5個(gè)波長);d.Wavelet小波，以symlet wavelet濾波把每副圖像分解4次，每次分解生成4組小波系數(shù)，分別是低頻(LL)，水平 (HL)，垂直(LH)和對(duì)角(HH)。因此共有16個(gè)系數(shù)(表1)。Gabor函數(shù)公式如下：

X′=Xcosθ+Ysinθ

Y′=-Xsinθ+Ycosθ

其中λ波長，θ方向，φ相位，σ高斯標(biāo)準(zhǔn)差，γ spatial aspect ratio

方向θ={0,π/8,2π/8,3π/8,4π/8,5π/8,6π/8,7π/8}

波長λ={4π/5,π,4π/3,2π,4π}

③邊緣分析(擬合sigmoid曲線，以曲線參數(shù)表示邊緣變化程度)，先沿病灶輪廓插值重建以病灶為中線的圖，再沿中線平均各數(shù)值，以垂直中線的方向?yàn)閤軸，以平均值為y軸，再擬合sigmoid函數(shù)。sigmoid函數(shù)公式如下：

其中a值的絕對(duì)值越大，表示曲線越陡，反映病灶邊緣的變化程度。

④增強(qiáng)動(dòng)態(tài)特征(時(shí)間特征量曲線各參數(shù)),形態(tài)學(xué)特征，邊緣特征以及內(nèi)部構(gòu)造特征等均采用以下方法提取動(dòng)態(tài)增強(qiáng)特征。a.多項(xiàng)式擬合的系數(shù)；b.次序。根據(jù)不同時(shí)間點(diǎn)特征變量數(shù)值的大小排序，有7×6×5×4×3×2=5040種可能性，每種可能性對(duì)照一個(gè)號(hào)作為特征向量。

5.機(jī)器學(xué)習(xí)

由于病灶特征數(shù)相對(duì)于病例數(shù)較大(總的提取特征1818個(gè)特征 )，需對(duì)特征進(jìn)行分類和篩選：首先將特征分為4組，分別是形態(tài)學(xué)組、內(nèi)部紋理組、病灶邊緣組、病灶平均信號(hào)強(qiáng)度。然后在分類器上進(jìn)行篩選變量及分類。本文采用的分類器為隨機(jī)森林，它是一種集成統(tǒng)計(jì)學(xué)習(xí)方法[13-15]，使用多個(gè)決策樹進(jìn)行組合，從而獲得較高的預(yù)測(cè)精度。特征向量按不同組別輸入良惡性二分類模型(隨機(jī)森林參數(shù)為mtry=4，ntrees=500)，再用交叉檢驗(yàn)的方法評(píng)估學(xué)習(xí)效果，分別計(jì)算出各組特征的敏感度和特異性、受試者工作特征曲線 (receiver operating characteristic curve，ROC)曲線下面積(area under curve，AUC)等來評(píng)價(jià)該組特征的診斷能力，最后比較各組的AUC差異。統(tǒng)計(jì)過程在軟件R[16]上進(jìn)行。

表1 具體提取的特征向量

結(jié) 果

所有病灶在自主編寫的半自動(dòng)分割程序下完成分割，醫(yī)師檢查分割結(jié)果，分割區(qū)域基本反映病灶輪廓范圍(圖1)。鑒別良惡性分類器顯示：形態(tài)學(xué)組的診斷ROC曲線下面積為0.90，診斷正確率為83%，敏感度90%，特異度75%。紋理組的ROC曲線下面積0.92，準(zhǔn)確率86%，敏感度90%，特異度84%。邊緣特征組的ROC曲線下面積:0.75，準(zhǔn)確率76%，敏感度80%，特異度65%。而單獨(dú)以信號(hào)強(qiáng)度為特征的ROC曲線下面積0.80，準(zhǔn)確率81%，敏感度92%，特異度62%。綜合各組診斷所得的ROC曲線下面積0.94，準(zhǔn)確率90%，敏感度92%，特異度85%(圖2)。綜合診斷組ROC曲線下面積要高于包括內(nèi)部紋理特征組在內(nèi)的所有特征類別，但與內(nèi)部紋理特征和形態(tài)學(xué)特征差異無統(tǒng)計(jì)學(xué)意義(P>0.05)，而內(nèi)部紋理特征組的ROC曲線下面積要高于信號(hào)強(qiáng)度特征及邊緣特征，差異有統(tǒng)計(jì)學(xué)意義(P<0.05，表2)。

圖1 a) 乳腺腺瘤病灶，邊緣光整，輕度分葉; b) 乳腺腺瘤病灶分割的結(jié)果，紅色區(qū)為分割的部分，基本能反映病灶的實(shí)際范圍； c) 乳腺浸潤性導(dǎo)管癌病灶呈邊緣毛糙，輪廓呈小波浪狀； d) 乳腺導(dǎo)管癌病灶的分割結(jié)果，紅色的分割區(qū)域基本反映病灶的實(shí)際范圍。

圖2 各類動(dòng)態(tài)特征ROC曲線(Signal Intensity：信號(hào)強(qiáng)度；Texture:內(nèi)部紋理；Morphology：形狀或形態(tài)；Margin Feature：邊緣特征；All Feature:所有特征)

表2 不同特征組的AUC的P值比較

注：*差異有統(tǒng)計(jì)學(xué)意義。

討 論

乳腺動(dòng)態(tài)增強(qiáng)掃描有助于醫(yī)生了解病灶的血流動(dòng)力學(xué)改變，然而，由此產(chǎn)生的大量的圖像會(huì)耗費(fèi)醫(yī)生大量的時(shí)間和精力進(jìn)行后處理和分析。目前臨床上多用時(shí)間-信號(hào)曲線來簡(jiǎn)化原始動(dòng)態(tài)圖像來達(dá)到降維分析的效果，這種曲線雖然簡(jiǎn)便可行，但它只能反映平均信號(hào)強(qiáng)化的變化，然而，目前，無論是BI-RAS分類還是平時(shí)的經(jīng)驗(yàn)，都提示病灶的紋理特征及其隨著動(dòng)態(tài)增強(qiáng)的改變更為重要，然而，如此大量的圖像信息，對(duì)人腦是很具挑戰(zhàn)性的。

因此，本文解決該問題的關(guān)鍵是通過計(jì)算機(jī)輔助的方法來實(shí)現(xiàn)，具體來說是通過圖像處理的方法，首先要解決的問題是在特征選擇方面，本文比以往國內(nèi)外的研究采用更多的特征量，同時(shí)兼顧了病灶的形態(tài)特征、內(nèi)部特征和邊緣特征；病灶的形態(tài)特征包括病灶與外接圓的面積比、半徑比、光滑度和緊湊度等。而病灶的內(nèi)部特征有一階統(tǒng)計(jì)紋理征， Haralick紋理特征(包括角二階矩、熵、對(duì)比度、逆差分矩、相關(guān)性)，Gabor濾波，Wavelet小波；邊緣特征：利用Sigmoid函數(shù)擬合邊緣特征[17]，進(jìn)而提取參數(shù)，這項(xiàng)在國外文獻(xiàn)鮮有使用。

由于本文采用較多的特征個(gè)數(shù)，同時(shí)，每個(gè)特征都有5個(gè)動(dòng)態(tài)增強(qiáng)的時(shí)間點(diǎn)，輸入的特征量相對(duì)樣本數(shù)來說較為巨大，增大了噪聲并容易出現(xiàn)過擬合。因此，本文將特征降維成主要的幾種特征類，分別為代表形態(tài)學(xué)的特征，邊緣的特征，病灶內(nèi)部紋理的特征還有信號(hào)強(qiáng)度的特征，而特征的動(dòng)態(tài)變化信息蘊(yùn)藏于各組特征之中。

訓(xùn)練方法：本文選取較為經(jīng)典的隨機(jī)森林方法，將多個(gè)決策樹進(jìn)行組合以提高預(yù)測(cè)精度。隨機(jī)森林有以下特點(diǎn)：相對(duì)其它算法，準(zhǔn)確率很高；不容易形成過度擬合；速度快，能夠處理大數(shù)據(jù)；能處理很高維度的數(shù)據(jù)，自動(dòng)做特征選擇；由于抽樣的原因，會(huì)有一些未被抽中的樣本，這形成了所謂的“out of bag”數(shù)據(jù)，并據(jù)此來估計(jì)泛化誤差，比較適合本文的數(shù)據(jù)。

在結(jié)果中，動(dòng)態(tài)紋理特征顯示出其較高的診斷能力，其效果比其他幾種動(dòng)態(tài)特征要好，高于目前臨床常用的動(dòng)態(tài)信號(hào)強(qiáng)度特征(AUC 0.72～0.88)[18]，這一點(diǎn)與既往研究相一致[7,9]；另外，動(dòng)態(tài)形態(tài)學(xué)特征也比較重要，但筆者推測(cè)其中起診斷作用的因素與“動(dòng)態(tài)”無關(guān)，因?yàn)閯?dòng)態(tài)增強(qiáng)時(shí)病灶強(qiáng)化范圍的形態(tài)學(xué)一般變化不明顯，所以病灶的輪廓(如圓形、星芒狀、分葉等)可能更重要； 病灶的邊緣特征在肝臟病變的鑒別診斷中起重要作用[17]，而本研究中其診斷能力是4組特征中最弱的一組，筆者推測(cè)可能因?yàn)楦闻K腫瘤更容易出現(xiàn)有特征性的邊緣強(qiáng)化模式，比如肝血管瘤的邊緣結(jié)節(jié)狀強(qiáng)化和假包膜的出現(xiàn)，而乳腺腫瘤這種特征性的邊緣改變較少。最后，值得一提的是本文的結(jié)果顯示綜合所有特征所得診斷模型診斷效率最高，這也是符合預(yù)期的。

未來可以以增強(qiáng)時(shí)間點(diǎn)為橫軸，病灶特征為縱軸，建立“時(shí)間-病灶特征”曲線，像“時(shí)間-信號(hào)強(qiáng)度”曲線那樣應(yīng)用于臨床，令醫(yī)生更為直觀地理解和使用。本文結(jié)果已經(jīng)證明了內(nèi)部紋理的動(dòng)態(tài)特征比單純信號(hào)強(qiáng)度的變化更具優(yōu)勢(shì)。因此，筆者相信“時(shí)間-病灶特征”曲線比“時(shí)間-信號(hào)強(qiáng)度”曲線更有診斷價(jià)值，但此不在本文的研究范圍內(nèi)。另外，該研究與影像組學(xué)的研究一脈相承，但后者要獲得更佳的結(jié)果，需要更大量影像、臨床和分子病理資料等，有望在以后的研究中進(jìn)行。

本研究緊密結(jié)合機(jī)器學(xué)習(xí)、醫(yī)學(xué)圖像處理、醫(yī)學(xué)影像診斷學(xué)3個(gè)學(xué)科，應(yīng)用圖像處理的方法提取病灶特征，機(jī)器學(xué)習(xí)的方法選擇特征，以醫(yī)學(xué)影像診斷學(xué)的方法分析和利用這些特征，是目前計(jì)算機(jī)輔助影像甚至人工智能醫(yī)學(xué)的發(fā)展的方向。

[1] 許玲輝,顧雅佳.乳腺磁共振在乳腺癌診斷治療中的運(yùn)用[J].中國癌癥雜志，2013，23(8):613-617.

[2] 何東峰,馬大慶,靳二虎.動(dòng)態(tài)增強(qiáng)MRI在鑒別乳腺良、惡性腫瘤和提示惡性腫瘤病理分級(jí)中的作用[J].中華放射學(xué)雜志，2012，46(12)：1075-1078.

[3] Saslow D,Boetes C,Burke W,et al.American cancer society guidelines for breast screening with MRI as an adjunct to mammography[J].CA Cancer J Clin，2007，57(2):75-89.

[4] American College of Radiology (ACR):Breast imaging reporting and data system atlas (BIRADS Atlas).American College of Radiology,Reston,2003

[5] 陳文靜,牟瑋,張文馨,等.MR動(dòng)態(tài)增強(qiáng)圖像紋理分析判斷乳腺結(jié)節(jié)良惡性的價(jià)值[J].中國醫(yī)學(xué)影像技術(shù),2017,33(5):647-651.

[6] Nie K,Chen JH,Yu HJ,et al.Quantitative analysis of lesion morphology and texture features for diagnostic prediction in breast MRI[J].Acad Radiol，2008，15(12):1513-1525.

[7] Agner SC,Soman S,Libfeld E,et al.Textural kinetics:a novel dynamic contrast-enhanced (DCE)-MRI feature for breast lesion classification[J].J Digit Imaging，2011，24(3):446-463.

[8] Gibbs P,Turnbull LW.Textural analysis of contrast-enhanced MR images of the breast[J].Magn Reson Med，2003，50(1):92-98.

[9] Zheng Y,Englander S,Baloch S，et al.STEP:spatiotemporal enhancement pattern for MR-based breast tumor diagnosis[J].Med Phys，2009，37(7):3192-3204.

[10] Honda E,Nakayama R,Koyama H,et al.Computer-aided diagnosis scheme for distinguishing between benign and malignant masses in breast DCE-MRI[J].J Digit Imaging,2016,29(3):388-393.

[11] Rakoczy M,McGaughey D,Korenberg M,et al.Feature selection in computer-aided breast cancer diagnosis via dynamic contrast-enhanced magnetic resonance images[J].J Digit Imaging.Springer-Verlag,2013,26(2):198-208.

[12] Abdel-Nasser M,Rashwan HA,Puig D,et al.Analysis of tissue abnormalityand breast density in mammographic images using a uniform local directional pattern[J].Expert Syst Appl，2015，42(24):9499-9511.

[13] Zhou Z.Machine Learning[M].Beijing:Tsinghua University Press,2016:229-232.

[14] Lu W,Li Z,Chu J.A novel computer-aided diagnosis system for breast MRI based on feature selection and ensemble learning[J].Comput Biol Med，2017，83(1):157-165.DOI:10.1016/j.compbiomed.2017.03.002.

[15] Strobl C,Boulesteix AL,Kneib T,et al.Conditional variable importance for random forests[J].BMC Bioinformatics,2008,9(307):471-2105.

[16] R Core Team.A language and environment for statistical computing.Foundation for Statistical Computing, Vienna, Austria. (Accessed December 10, 2012 at http://www.R-project. org/.)

[17] Fan L,Yi L.Image features related to margin and enhancement pattern improve the performance of computer-aided diagnosis for hepatic diseases using multi-phase computed tomography[J].Chinese Medical J,2014,127(19):3406-3417.

[18] El Khouli RH,Macura KJ,Barker PB,et al.The relationship of temporal resolution to diagnostic performance for dynamic contrast enhanced (DCE) MRI of the breast[J].JMRI,2009,30(5):999-1004.

DifferentiateddiagnosisbetweenmalignantandbenignbreastmassusingtexturalkineticMRIfeatures

LIN Fan,HU Ruo-fan,LIANG Chao,et al.

Department of Radiology,Shenzhen Second People's Hospital,Guangdong 518035,China

Objective:In this study,we proposed a new DCE-MRI descriptor (called textural kinetics) to capture the spatiotemporal changes of lesion texture for the better discrimination of benign and malignant lesions.MethodsWe retrospectively collected 52 breast mass lesions (30 malignant and 22 benign).Kinetic features were extracted from breast DCE-MRI and divided into four categories:signal intensity,texture,morphology,margin feature.In order to evaluate these features,we created classification model using different feature classes.Accuracy,sensitivity and specifity and area under curve (AUC) were calculated.ResultsWhen textural kinetic features were combined with other features,the classifier yielded 90% accuracy,92% sensitivity and 85% specificity,and 0.94 AUC,showing the best differentiation between benign and malignant lesions.The difference for AUC between combined classifier and signal intensity classifier is significant (AUC=0.80,P<0.05).ConclusionsTextural kinetic features in DCE-MRI is helpful in distinguishing malignant from benign lesions,which is even better than the most popular DCE-MRI descriptor of signal intensity kinetics feature in clinic.

Magnetic resonance imaging; Texture; Breast neoplasms; Artificial intelligence

R445.2; R737.9; R05

A

1000-0313(2017)10-1037-05

2017-09-19

2017-10-09)

518035 廣東，深圳市第二人民醫(yī)院 放射科

林帆(1981-)，男，廣東梅縣人，碩士，副主任醫(yī)師，主要從事計(jì)算機(jī)輔助診斷工作。

雷益，E-mail：13602658583@163.com

深圳市科技研究資金(JCYJ20150330102720117)

10.13609/j.cnki.1000-0313.2017.10.009