湯衛(wèi)霞，盛美紅，龔沈初，季娟，張玉嬌

近年來，隨著全球女性乳腺癌發(fā)病率的逐年增高，早發(fā)現(xiàn)、早診斷、早治療可明顯提高乳腺癌患者的生存率并改善預(yù)后。磁共振對軟組織的分辨率有著獨(dú)特的優(yōu)勢，檢測乳腺癌敏感度較高，磁共振成像(magnetic resonance imaging，MRI)在臨床應(yīng)用較廣，但傳統(tǒng)MRI序列如表觀擴(kuò)散系數(shù)(apparent diffusion coefficient ，ADC)、動(dòng)態(tài)對比增強(qiáng)MRI(dynamic contrast enhancement MRI，DCE-MRI)圖像中隱含許多肉眼無法直接觀測的信息，這些圖像紋理信息與腫瘤組織的顯微結(jié)構(gòu)、生物學(xué)特性等可能存在一定的關(guān)系，所以紋理分析作為一種高通量提取特征的無創(chuàng)性新技術(shù)應(yīng)運(yùn)而生，通過提取圖像紋理特征，為病變鑒別診斷提供肉眼無法觀察到的腫瘤內(nèi)部異質(zhì)性信息，有助于明確病灶性質(zhì)[1]。本文就紋理分析在乳腺癌診斷、治療后反應(yīng)及預(yù)后、預(yù)測有無腋窩淋巴結(jié)轉(zhuǎn)移評估中的應(yīng)用進(jìn)行綜述，并分析基于紋理分析的無創(chuàng)性預(yù)測性標(biāo)記物的臨床應(yīng)用價(jià)值。

常用的紋理特征提取方法

紋理分析的目的是提取興趣區(qū)(region of interest，ROI)內(nèi)所表征的圖像紋理的高通量信息。圖像分割是提取紋理特征的前提和關(guān)鍵，常用方法有手動(dòng)、半自動(dòng)和全自動(dòng)三種，前兩種方法應(yīng)用較多，目前最常用的紋理分析平臺有MaZda軟件、GE Omni-Kinetics軟件及基于相關(guān)算法的自主研發(fā)軟件等。病灶ROI的勾畫方法包括面積法和體積法，面積法是選擇病灶最大截面勾畫ROI，體積法則是勾畫病灶整體。另外有研究表明瘤周組織的紋理特征也與疾病的總體特征相關(guān)聯(lián)[2]。

最常用的紋理特征可以通過圖像中編碼的體素信息的統(tǒng)計(jì)順序來分層。具體的說，灰度級在圖像中的像素分配通過影像圖像的紋理特征的方式來定義，這些紋理特征使用一階、二階或更高階的統(tǒng)計(jì)方法進(jìn)行數(shù)學(xué)提取。一階紋理特征基于灰度直方圖導(dǎo)出的常見統(tǒng)計(jì)度量，例如平均值、中值和偏度。二階紋理特征由灰度共生矩陣(gray-level co-occurrence matrix，GLCM，1973年Haralick和Shanmugam研究定義)[3]和灰度游程長度[4]導(dǎo)出，而高階紋理特征為編碼結(jié)構(gòu)和基于頻率的紋理信息。

1.一階：灰度直方圖特征

通過計(jì)算每個(gè)灰度強(qiáng)度值的體素?cái)?shù)量的頻率計(jì)數(shù)來生成灰度直方圖，其中離散灰度的總數(shù)通常是用戶選擇的參數(shù)。從所得的直方圖可以得出一階統(tǒng)計(jì)特征，包括均值、中位數(shù)和方差。還可以從直方圖中提取更高的矩特征，包括偏度，直方圖分布對稱性的度量、峰度、直方圖分布形狀的度量。雖然可以從直方圖中提取許多描述符，但它們通常只是粗略洞察底層紋理，不考慮圖像內(nèi)的灰度級強(qiáng)度空間關(guān)系，并也取決于用戶選擇的參數(shù)。

2.二階：灰度共生矩陣和灰度游程長度

GLCM最常提取的紋理特征，對兩個(gè)特定灰度強(qiáng)度的體素在指定的圖像方向上彼此隔開指定距離的頻率進(jìn)行編碼[3]。GLCM紋理特征通常在0°、45°、90°和135°四個(gè)對角線方向上進(jìn)行量化，然后可以從共現(xiàn)矩陣中提取二階紋理特征。這樣的特征示例包括對比度，由距離參數(shù)確定的像素與其相鄰像素之間的強(qiáng)度對比度的描述符、相關(guān)性，線性灰度相關(guān)性的描述符以及同質(zhì)性，像素中分布的緊密性的描述符，共現(xiàn)矩陣到矩陣對角線。還可以提取其他二階特征，例如能量，確定灰度共生的描述符以及簇陰影(灰度值不對稱的描述符)。另還有熵或GLCM的隨機(jī)性，通常表示圖像異質(zhì)性?；叶扔纬涕L度特征是指定線性方向上圖像的粗糙度[4-6]。

3.高階：基于結(jié)構(gòu)和變換

結(jié)構(gòu)紋理特征捕獲了中央體素與周圍相鄰體素之間的強(qiáng)度變化[7]?；诩y理的擴(kuò)展將紋理分析擴(kuò)展到灰度強(qiáng)度值的空間方向之外，可以捕獲在不同空間(例如頻率空間)中編碼的紋理信息?；谧儞Q的方法包括傅里葉變換，Gabor變換和小波變換。由于小波變換能夠捕獲不同圖像比例和多個(gè)指定方向的MR圖像頻率內(nèi)容[8-10]，所以較為常用。

紋理分析在乳腺癌中的應(yīng)用

1.鑒別乳腺病變良惡性

Gibbs等[11]是早期使用紋理分析來鑒別乳腺良惡性病變的研究者之一，分析了79名女性(其中45名被診斷為乳腺癌)患者DCE-MR圖像，共提取14個(gè)GLCM紋理特征，通過logistic回歸分析發(fā)現(xiàn)方差、熵、和熵的紋理特征在鑒別乳腺良惡性病變方面具有重要意義，說明量化圖像紋理特征可以作為病變描述的有用工具。這一初步研究結(jié)論給隨后的學(xué)者帶來啟發(fā)，他們利用更復(fù)雜的圖像紋理測量方法來鑒別良惡性。張竹偉等[12]回顧性分析了69例乳腺病變患者，在乳腺平掃T1WI、T2WI及延時(shí)增強(qiáng)T1WI圖像上選取病灶最大層面圖像，在MaZda軟件中手動(dòng)勾畫病灶ROI，提取病變的紋理特征參數(shù)六類共287個(gè)紋理特征參數(shù)值，使用MaZda軟件提供的紋理特征選擇方法及分類分析統(tǒng)計(jì)方法分析，最后得出T2WI序列的紋理特征在鑒別乳腺良惡性病變時(shí)的誤判率最小，紋理特征參數(shù)選擇方法中，使用FPM[交互信息(MI)、Fisher系數(shù)、分類錯(cuò)誤概率聯(lián)合平均相關(guān)系數(shù)(POE+ACC)三者聯(lián)合]選擇的紋理特征參數(shù)鑒別病變的誤判率最低，紋理特征分類分析方法中，NDA(非線性分類分析)區(qū)分兩種病變的誤判率具有最優(yōu)的鑒別診斷效能。鄧義等[13]同樣使用MaZda軟件，對早期DCE-MR圖像進(jìn)行紋理分析，提取整個(gè)病變的紋理參數(shù)，采用NDA的誤判率最低，其中POE+ACC聯(lián)合NDA的誤判率最低，誤判率為5.88%。在乳腺病變的30個(gè)最優(yōu)紋理參數(shù)中，有10個(gè)紋理參數(shù)在良惡性組間的差異具有統(tǒng)計(jì)學(xué)意義。Nagarajan等[14]發(fā)現(xiàn)基于GLCM的紋理特征量化可以顯著提高小病灶的診斷效能。雖然這些研究中有不少在具有高度鑒別能力的特定紋理特征方面有相似的結(jié)論，但是缺乏文獻(xiàn)對同一數(shù)據(jù)集進(jìn)行比較研究，并為提取的特征定義參數(shù)。

Bhooshan[15]提出了一種將T2WI和T1WI序列組合多參數(shù)DCE-MRI計(jì)算機(jī)輔助診斷方法，采用基于模糊c均值聚類的方法對DCE-MR圖像自動(dòng)分割乳腺腫瘤，提取GLCM 3D紋理特征，從DCE-MR圖像中提取血流動(dòng)力學(xué)特征，建立多參數(shù)特征向量。研究發(fā)現(xiàn)，當(dāng)區(qū)分良惡性病變時(shí)，多參數(shù)特征向量優(yōu)于從單個(gè)成像序列中提取的特征，這表明紋理分析應(yīng)用于T2WI圖像時(shí)，提供了除T1WI圖像外的其他腫瘤異質(zhì)性信息。Yun等[16]在IDL平臺上開發(fā)了一種用于逐像素獲取ADC映射的自制軟件包，在顯示腫瘤最大徑層面繪制ROI，發(fā)現(xiàn)其紋理分析能更精確的提供關(guān)于腫瘤細(xì)胞增殖能力信息，并可以對腫瘤的異質(zhì)性進(jìn)行定量測量。據(jù)此可猜測DWI的紋理分析將成為非侵入性的檢測乳腺惡性腫瘤治療后異質(zhì)性的方法。另外，傳統(tǒng)的紋理分析在乳腺診斷中的應(yīng)用主要集中在從圖像中提取二維紋理特征。Chen等[17]擴(kuò)展這一思想從3D乳腺病變分割中提取3D-GLCM紋理特征，發(fā)現(xiàn)在區(qū)分乳腺良惡性病變時(shí)，3D-GLCM特征比從2D-ROI提取的2D-GLCM特征具有更高的診斷符合率。

2.鑒別乳腺癌組織病理及分子分型

MRI紋理分析有廣泛的應(yīng)用前景，利用紋理特征或許可以對乳腺良惡性病變的影像特征之間的結(jié)構(gòu)性差異進(jìn)行量化。近來不少學(xué)者提出假設(shè)，潛在的腫瘤生物學(xué)差異是否可以通過MRI和紋理分析來表征，并對這一設(shè)想進(jìn)行了拓展研究。有學(xué)者[18]在區(qū)分浸潤性小葉癌和浸潤性導(dǎo)管癌研究中提取了277個(gè)直方圖、GLCM、游程和小波特征，研究發(fā)現(xiàn)乳腺癌與正常乳腺組織不同的紋理特征，使用減影后的圖像進(jìn)行分類效果最好，在鑒別導(dǎo)管癌和小葉浸潤癌方面也取得了令人滿意的結(jié)果。2016年，Waugh等[19]的研究也得出相似的結(jié)論，在小葉性病變和導(dǎo)管性病變之間熵特征有顯著差異，MR增強(qiáng)圖像的紋理差異可反映潛在的病變亞型。同樣紋理分析有助于表征乳腺癌分子亞型。乳腺癌常用的生物學(xué)標(biāo)記物有雌激素和孕激素受體(ER，PR)、人表皮生長因子受體2型(HER2)和Ki67?；谶@些標(biāo)志物狀態(tài)將乳腺癌分為不同亞型，激素受體陽性乳腺癌通常分為管腔A型和管腔B型(有或無HER2過表達(dá)，乳腺癌管腔A亞型表達(dá)激素受體水平高，預(yù)后好，管腔B亞型預(yù)后差)，HER2陰性型、 HER2陽性型、基底樣型、三陰性型及其他特殊類型。Holli-Helenius等[20]對27例經(jīng)病理證實(shí)的浸潤性乳腺導(dǎo)管癌患者進(jìn)行初步研究，在非壓脂T1WI增強(qiáng)前MR圖像上，采用MaZda軟件提取每個(gè)腫瘤基于共生矩陣的紋理特征，結(jié)果顯示和熵、和方差紋理參數(shù)在區(qū)分luminal A和luminal B亞型中差異有統(tǒng)計(jì)學(xué)意義，和熵及和方差與Ki-67指數(shù)呈顯著正相關(guān)，測量隨機(jī)性、異質(zhì)性、平滑性和均勻性的紋理特征可以直接或間接反映乳腺腫瘤的潛在生長模式。近來，F(xiàn)an等[2]采用leave one-out交叉驗(yàn)證(LOOCV)方法對隨機(jī)森林模型進(jìn)行訓(xùn)練和測試，以預(yù)測乳腺癌亞型，也得出相似結(jié)論，另外研究者們還發(fā)現(xiàn)在進(jìn)行紋理分析研究時(shí)擴(kuò)大范圍，即擴(kuò)增瘤周區(qū)域，這一方法可更準(zhǔn)確的識別腫瘤動(dòng)力學(xué)特征，從而提高乳腺癌亞型的預(yù)測水平。

3.預(yù)測乳腺癌預(yù)后

最近不少學(xué)者在探索乳腺病灶紋理特征與有無復(fù)發(fā)風(fēng)險(xiǎn)之間的關(guān)系，試圖把紋理特征作為無創(chuàng)性生物學(xué)標(biāo)記物。Kim等[21]探討了原發(fā)性乳腺癌患者基于MR圖像紋理分析評價(jià)腫瘤異質(zhì)性與生存率的關(guān)系，對203例診斷為原發(fā)性乳腺腫瘤患者進(jìn)行回顧性分析，采用T2WI及增強(qiáng)后T1WI減影MR圖，從這些序列圖像中提取基于直方圖的均勻性和熵，確定紋理參數(shù)與無復(fù)發(fā)生存率(relapse free survival，RFS)之間的關(guān)系，結(jié)果顯示在T2WI圖像較高熵值的乳腺癌患者的RFS較低，這項(xiàng)研究表明用紋理量化腫瘤的異質(zhì)性可以作為一個(gè)獨(dú)立的預(yù)后指標(biāo)。Park等[22]也得出了類似的結(jié)論，他們基于形態(tài)學(xué)、直方圖紋理和GLCM紋理特征生成了一個(gè)多變量特征向量，從中特別選擇了GLCM聚類趨勢、GLCM方差和GLCM和方差組成模型，用于對有復(fù)發(fā)風(fēng)險(xiǎn)的患者進(jìn)行分層。這些研究都已經(jīng)證明了MRI提取的紋理特征有潛力成為無創(chuàng)性預(yù)測標(biāo)記物。

4.乳腺癌治療療效評估

新輔助化療(neoadjuvant chemotherapy，NAC)是晚期乳腺癌患者的標(biāo)準(zhǔn)治療方法，可以縮小腫瘤體積，降低臨床分期，從而提高保乳手術(shù)的概率，在臨床上得到廣泛認(rèn)可和應(yīng)用。早期并且準(zhǔn)確預(yù)測腫瘤對化療的反應(yīng)能夠?qū)εR床提供幫助，有利于及時(shí)調(diào)整治療方案。為了預(yù)測腫瘤患者對NAC治療的反應(yīng)，孫賽花等[23]探討了磁共振紋理分析早期預(yù)測NAC療效的可行性，使用GE Omni-Kinetics軟件，在增強(qiáng)第3期病灶強(qiáng)化明顯圖像上選擇病灶最大層面，手動(dòng)勾畫ROI，共收集47例患者，分別測量NAC前和2個(gè)化療周期后的MRI紋理參數(shù)(能量、熵、慣量、相關(guān)和逆差距)，比較兩次MRI紋理參數(shù)的變化量之間的差異，結(jié)果表明這些參數(shù)的差異均具有統(tǒng)計(jì)學(xué)意義，其中熵值變化量的AUC最大，預(yù)測腫瘤治療的pCR敏感度較高。另外還有不少類似的研究[24-25]，其中Eun等[24]對136名乳腺癌患者分別在治療前、中后(第3或第4周期)進(jìn)行MRI檢查，采用T2WI、增強(qiáng)后T1WI、DWI及ADC圖像進(jìn)行紋理參數(shù)提取，應(yīng)用隨機(jī)森林法建立基于紋理參數(shù)的pCR分類預(yù)測模型，并與其他6種機(jī)器學(xué)習(xí)分類器(自適應(yīng)boosting、決策樹、k近鄰、線性支持向量機(jī)、樸素貝葉斯和線性判別分析)進(jìn)行比較，最后得出治療中期應(yīng)用MRI增強(qiáng)圖像隨機(jī)森林模型提取的乳腺癌相關(guān)紋理參數(shù)較有價(jià)值，并與乳腺癌pCR相關(guān)聯(lián)。因此定量MRI紋理特征有望成為基于圖像特征評估癌癥復(fù)發(fā)風(fēng)險(xiǎn)的方法，并能與基因組數(shù)據(jù)相結(jié)合改善生存預(yù)測指標(biāo)。

5.預(yù)測腋窩淋巴結(jié)狀態(tài)

在乳腺癌患者中，腋窩淋巴結(jié)轉(zhuǎn)移也是一個(gè)重要的預(yù)后因素，同時(shí)也是決定是否需要全身輔助化療或術(shù)后放療的關(guān)鍵[26-28]。目前臨床常用的是超聲引導(dǎo)下細(xì)針穿刺活檢淋巴結(jié)進(jìn)行確診，對于確診的腋窩淋巴結(jié)轉(zhuǎn)移患者進(jìn)行腋窩淋巴結(jié)清掃，以上這兩種有創(chuàng)方法不可避免帶來一些潛在的并發(fā)癥，如淋巴水腫、漿膜瘤、感染、血腫和手臂疼痛等[29]。紋理分析將有助于確定術(shù)前腋窩淋巴結(jié)狀態(tài)，減少不必要的淋巴結(jié)清掃，減輕患者痛苦?；谶@一設(shè)想，Chai等[30]回顧性分析了120例乳腺癌患者(59例腋窩淋巴結(jié)轉(zhuǎn)移)，從多個(gè)MRI序列中提取典型的腫瘤紋理學(xué)特征，這些特征的符合率均高于74%，又將增強(qiáng)后第2期圖像紋理特征與增強(qiáng)后動(dòng)力學(xué)特征組合成新模型，它的診斷性能得到了進(jìn)一步提高(符合率86%，AUC為0.91)，基于這一模型可以很好的預(yù)測腋窩淋巴結(jié)狀態(tài)。同樣，Tan等[31]研究發(fā)現(xiàn)基于脂肪抑制T2WI序列的紋理特征有助于準(zhǔn)確預(yù)測腋窩淋巴結(jié)轉(zhuǎn)移，這與前者研究結(jié)果相似。因此，將來有必要進(jìn)行大樣本的研究以進(jìn)一步驗(yàn)證這一發(fā)現(xiàn)。

總結(jié)與展望

紋理分析通過對醫(yī)學(xué)圖像進(jìn)行深度挖掘，揭示腫瘤異質(zhì)性，反映潛在的分子生物學(xué)變化[32]，在腫瘤的早期診斷、精準(zhǔn)治療和預(yù)后評估等方面表現(xiàn)出巨大的潛力，目前的研究主要集中在乳腺癌分子分型、預(yù)測乳腺癌治療后反應(yīng)及預(yù)后、以及預(yù)測腋窩淋巴結(jié)狀態(tài)等。但是其作為一個(gè)相對新興的領(lǐng)域，仍處于起步階段面臨很多挑戰(zhàn)，恰當(dāng)?shù)难芯吭O(shè)計(jì)和統(tǒng)計(jì)評價(jià)、協(xié)調(diào)好表征轉(zhuǎn)化紋理特征的量化和預(yù)測值的影響將成為未來人工智能領(lǐng)域進(jìn)行研究的重要課題。將紋理分析與傳統(tǒng)影像學(xué)指標(biāo)以及其他臨床特征相結(jié)合，將會(huì)在臨床診療中發(fā)揮重要作用，不斷提高乳腺疾病的診斷創(chuàng)新，推進(jìn)整個(gè)醫(yī)學(xué)界的精確性和個(gè)性化醫(yī)療。