張 靜，仇清濤，段敬豪，姜慶軍，孫 鋼，鞏貫忠，李登旺，尹 勇*

(1.山東師范大學(xué)物理與電子科學(xué)學(xué)院，山東 濟(jì)南 250358；2.山東省醫(yī)學(xué)物理圖像處理技術(shù)重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250358；3.山東大學(xué)附屬山東省腫瘤醫(yī)院放療科，山東 濟(jì)南 250117；4.濟(jì)南軍區(qū)總醫(yī)院醫(yī)學(xué)影像科，山東 濟(jì)南 250031)

DOI：10.13929/j.1003-3289.201709043

肝硬化是多種肝臟疾病發(fā)展的終末期，以廣泛纖維瘢痕、殘留組織增生為特征，3%～5%的肝硬化可發(fā)展為肝細(xì)胞癌[1]。近年來(lái)，研究[2]發(fā)現(xiàn)基于醫(yī)學(xué)圖像提出的影像組學(xué)——紋理分析技術(shù)有望成為新的腫瘤影像學(xué)生物標(biāo)記物。DWI紋理分析技術(shù)已應(yīng)用于腦膠質(zhì)瘤[3]、鼻咽癌[4]等的鑒別、分級(jí)以及腹部惡性腫瘤的診斷、定性、分級(jí)和預(yù)后評(píng)估，但因b值算法標(biāo)準(zhǔn)尚未統(tǒng)一，其可比性有待進(jìn)一步探索。Becker等[5]發(fā)現(xiàn)b值對(duì)正常人不同器官的DWI紋理特征存在顯著影響。本研究分析不同b值對(duì)提取肝硬化DWI紋理特征的影響，并探討不同b值與紋理特征變化的相對(duì)趨勢(shì)。

1 資料與方法

1.1 一般資料 收集2016年1月—2017年2月于我院接受MR檢查的肝硬化患者53例。納入標(biāo)準(zhǔn)：①肝臟解剖結(jié)構(gòu)完整；②經(jīng)3名影像學(xué)診斷醫(yī)師確診為肝硬化；③無(wú)介入及手術(shù)治療史；④無(wú)最大徑>2 cm的占位性病變及彌漫性病變。排除標(biāo)準(zhǔn)：①圖像存在偽影；②大量腹腔積液；③圖像和數(shù)據(jù)缺失。最終納入39例患者(肝硬化組)，男29例，女10例，年齡33～81歲，平均(55.7±11.1)歲，均伴不同程度門(mén)靜脈高壓及脾大，其中食管、胃底靜脈曲張10例，伴少量腹腔積液10例，膽囊炎12例，肝癌8例、糖尿病2例。另選取同期27名肝臟正常、并接受肝臟DWI檢查的受檢者為對(duì)照組，其中男12名，女15名，年齡30～78歲，平均(57.5±15.1)歲。2組間性別差異有統(tǒng)計(jì)學(xué)意義(χ2=6.068，P=0.02)，年齡差異無(wú)統(tǒng)計(jì)學(xué)意義(t=0.563，P=0.576)。本研究經(jīng)我院倫理委員會(huì)批準(zhǔn)，所有受檢者均簽署知情同意書(shū)。

1.2 儀器與方法 采用GE Discovery 750 3.0T MR掃描儀，8通道相控陣腹部線圈。掃描序列及參數(shù)：T1WI，TR 3.7 ms，TE 1.7 ms，F(xiàn)OV 38 cm×34.2 cm，層厚5.2 mm，間隔-2.6 mm，NEX 1；體素內(nèi)不相干運(yùn)動(dòng)-擴(kuò)散加權(quán)平面回波成像(intravoxel incoherent motion-diffusion weighted imaging echo planar imaging, IVIM-DWI EPI)序列，TR 10 000 ms，TE 66.1 ms，層厚7 mm，間隔1 mm，NEX 1，矩陣256×256，b值分別為0、20、50、100、200、400、800、1 000、1 200、1 500 s/mm2。

以DWI(b=0)為參考圖像，避開(kāi)下腔靜脈、門(mén)靜脈及匯管區(qū)，采用MIM 6.7.6(www.mimsoftware.com)軟件，根據(jù)Couinaud肝段劃分法，于肝臟Ⅶ段、Ⅱ/Ⅲ段和Ⅴ/Ⅵ段選取3個(gè)直徑20 mm的圓形ROI(圖1)，并使用勾畫(huà)克隆功能將其映射至其余9個(gè)b值的DWI圖像。然后采用IBEX(imaging biomarker explorer, IBEX)軟件對(duì)每個(gè)ROI分別行灰度統(tǒng)計(jì)，獲得灰度共生矩陣(gray-level co-occurrence matrix, GLCM)、灰度游程矩陣(gray-level run-length matrix, GLRLM)和鄰域灰度差矩陣(neighbor gray-tone difference matrix, NGTDM)，對(duì)這三類(lèi)矩陣共提取37個(gè)紋理特征(表1)；計(jì)算不同b值下的平均紋理特征值，并以百分比變異系數(shù)(percent coefficient of variation, %COV)≤15%評(píng)價(jià)2組平均紋理特征[6]的穩(wěn)定性，公式為%COV=SD/Mean×100%，其中SD為標(biāo)準(zhǔn)差，Mean為所有受檢者3個(gè)ROI的不同b值下紋理特征均值。%COV>15%為紋理特征不穩(wěn)定。

表1 37個(gè)紋理特征名稱及其分類(lèi)

1.3 統(tǒng)計(jì)學(xué)分析 采用Matlab 2014a統(tǒng)計(jì)分析軟件，采用指數(shù)擬合方法分析紋理特征值與b值的擬合程度，R2為擬合系數(shù)，R2≥0.8為擬合程度較好，0.5≤R2<0.8為擬合程度中等，R2<0.5為擬合程度較差。采用兩獨(dú)立樣本t檢驗(yàn)比較2組%COV值和R2值，P<0.05為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 2組平均紋理特征值%COV結(jié)果 37個(gè)紋理特征中，1個(gè)(1/37，2.70%)紋理特征(聚類(lèi)陰影)因數(shù)值出現(xiàn)正負(fù)交叉現(xiàn)象、%COV>100%而被剔除。20個(gè)(20/37，54.05%)紋理特征值%COV均>15%，為不穩(wěn)定紋理特征，其中GLCM中8個(gè)(8/21，38.09%)、GLRLM中7個(gè)(7/11，63.63%)、NGTDM中5個(gè)(5/5，100%)；余16個(gè)(16/37，43.24%)紋理特征值%COV均≤15%，為穩(wěn)定紋理特征。

2組間相關(guān)性(t=-4.352，P<0.001)、能量(t=-2.743，P=0.008)、同質(zhì)性1(t=-2.285，P=0.026)、同質(zhì)性2(t=-2.498，P=0.015)、逆方差(t=-2.268，P=0.027)、粗糙度(t=-2.791，P=0.007)及對(duì)比度(t=-2.312，P=0.024)的%COV值差異均有統(tǒng)計(jì)學(xué)意義，其余紋理特征%COV值差異無(wú)統(tǒng)計(jì)學(xué)意義(P均>0.05)，見(jiàn)圖2。

2.2 2組不穩(wěn)定平均紋理特征值指數(shù)擬合模型分析結(jié)果 2組中自相關(guān)、聚類(lèi)突、聚類(lèi)趨勢(shì)、對(duì)比度、和方差、方差、高灰度游程因子、短游程高灰度因子、復(fù)雜度及對(duì)比度隨b值增大呈指數(shù)上升趨勢(shì)；而相關(guān)性、能量、長(zhǎng)游程因子及粗糙度隨b值增大呈指數(shù)下降趨勢(shì)；無(wú)法定義長(zhǎng)游程高灰度因子、長(zhǎng)游程低灰度因子、低灰度游程因子、短游程低灰度因子、頻度及紋理強(qiáng)度與b值的相對(duì)變化趨勢(shì)，見(jiàn)圖3。

2組不穩(wěn)定紋理特征的R2值比較見(jiàn)表2。2組間聚類(lèi)突、聚類(lèi)趨勢(shì)、相關(guān)性、方差、長(zhǎng)游程低灰度因子、粗糙度、紋理強(qiáng)度R2值差異有統(tǒng)計(jì)學(xué)意義(P均<0.05)，其中僅紋理強(qiáng)度R2值肝硬化組優(yōu)于對(duì)照組(0.71 vs 0.03，P<0.001)。

3 討論

DWI可從分子水平對(duì)生物體的組織結(jié)構(gòu)和功能狀態(tài)進(jìn)行無(wú)創(chuàng)性檢查，檢測(cè)組織水分子擴(kuò)散狀況[7]。DWI中b值影響局灶性病變檢出的敏感性，低b值對(duì)檢出最大徑1～3 cm的肝臟惡性腫瘤更具優(yōu)勢(shì)，高b值區(qū)分肝臟惡性腫瘤與肝纖維化和假性動(dòng)脈增強(qiáng)時(shí)的特異度較高[8]。紋理特征是通過(guò)醫(yī)學(xué)圖像處理技術(shù)將ROI的影像數(shù)據(jù)轉(zhuǎn)化為可發(fā)掘的高維特征數(shù)據(jù)，可為疾病診斷、治療監(jiān)測(cè)和預(yù)后評(píng)估提供參考[9]。目前肝纖維化、肝硬化、肝囊腫、肝血管瘤、肝癌及肝轉(zhuǎn)移瘤等均有T1WI、T2WI、DWI及DCE-MRI影像組學(xué)——紋理分析相關(guān)研究[10-12]。使用DWI進(jìn)行紋理分析時(shí)，有必要將不同的可變量b值作為潛在的混雜因素考慮在內(nèi)[5，13]。

表2 2組紋理特征R2比較

圖1 DWI中選取ROI示意圖(b=0 s/mm2) A.ROI位于肝Ⅶ段； B.ROI位于肝Ⅱ/Ⅲ段； C.ROI位于肝Ⅴ/Ⅵ段

圖2 2組平均紋理特征值%COV比較 (*：差異有統(tǒng)計(jì)學(xué)意義)

圖3 紋理特征與b值指數(shù)擬合曲線圖 A.自相關(guān)； B.相關(guān)性； C.低灰度游程因子； D.聚類(lèi)突； E.能量； F.紋理強(qiáng)度 (紅線：肝硬化組；藍(lán)線：對(duì)照組)

本組20個(gè)不穩(wěn)定紋理特征的平均紋理特征值%COV均>15%，其中GLCM 8個(gè)、GLRLM 7個(gè)、NGTDM 5個(gè)。GLCM包含灰度值信息，且能展示灰度值關(guān)于灰度空間的特征分布，其中差熵、熵、逆差鉅、逆差及和熵均是涉及圖像紋理復(fù)雜程度的特征。本研究中2組上述紋理特征%COV值差異均無(wú)統(tǒng)計(jì)學(xué)意義，可能因DWI雖然對(duì)組織中水分子擴(kuò)散敏感，但其組織空間分辨率較T1WI低[10]，故從不同b值DWI中提取的紋理特征對(duì)鑒別肝硬化組織中的纖維組織及再生結(jié)節(jié)不敏感。利用DWI進(jìn)行肝硬化及相關(guān)紋理分析研究時(shí)，需考慮b值差異的影響，進(jìn)一步優(yōu)化紋理特征組合，以輔助提高肝硬化及肝硬化背景下的肝癌及肝血管瘤等疾病的診斷及檢測(cè)能力。本研究2組間有7個(gè)紋理特征的%COV值差異有統(tǒng)計(jì)學(xué)意義(P均<0.05)，且“聚類(lèi)陰影”在各b值時(shí)出現(xiàn)正負(fù)交叉現(xiàn)象，可能是由于在不同病例間勾畫(huà)ROI時(shí)，因血管灌注、不等量纖維基質(zhì)成分及磁場(chǎng)均勻性等因素所致[14]。

肝硬化患者的DWI提取的紋理特征(包括GLCM、GLRLM、NGTDM)受不同b值影響。本研究2組中部分紋理特征隨b值增大而呈指數(shù)上升或指數(shù)下降趨勢(shì)，與Becker等[5]研究結(jié)果一致，可能因患者間的異質(zhì)性和不同的磁化衰減程度，使低b值時(shí)DWI圖像易受毛細(xì)血管灌注效應(yīng)影響，高b值時(shí)易受細(xì)胞和細(xì)胞外基質(zhì)的物理密度及其幾何排列影響[15]。另外，本研究中2組有6個(gè)紋理特征隨b值增長(zhǎng)變化趨勢(shì)雜亂無(wú)章，原因可能在于肝臟內(nèi)膽管、血管及肝實(shí)質(zhì)的比例變化同樣可引起肝臟紋理特征分析改變[16]。2組平均紋理特征的R2值比較提示對(duì)照組擬合效果普遍優(yōu)于肝硬化組，且7個(gè)紋理特征間的R2值比較差異有統(tǒng)計(jì)學(xué)意義(P均<0.05)。郭冬梅等[13]發(fā)現(xiàn)從門(mén)靜脈期提取的紋理特征分類(lèi)識(shí)別3組肝臟組織(正常肝臟、較早期肝硬化及中晚期肝硬化)的正確率優(yōu)于T1WI、動(dòng)脈期及T2WI；王波濤等[12]發(fā)現(xiàn)紋理特征中的對(duì)比度、自相關(guān)鑒別肝囊腫及肝血管瘤的診斷價(jià)值較高，提示采用MR對(duì)不同器官、疾病及病程進(jìn)行紋理分析時(shí)，需考慮序列選擇差異。本研究中肝硬化組聚類(lèi)突、聚類(lèi)趨勢(shì)、相關(guān)性、方差及粗糙度5個(gè)紋理特征受不同b值影響，與對(duì)照組變化趨勢(shì)一致，且2組間上述紋理特征R2值差異有統(tǒng)計(jì)學(xué)意義(P均<0.05)，提示上述紋理特征有助于提高肝硬化背景疾病分析的診斷、治療監(jiān)測(cè)和預(yù)后評(píng)估。

本研究的局限性：①樣本量較小；②未使用三維VOI；③未對(duì)肝硬化不同致病原因分類(lèi)；④未使用T1WI、ADC的特征作為參照，結(jié)果可能存在偏倚。

綜上所述， DWI的b值影響肝硬化背景下的紋理特征；部分紋理特征與b值存在相關(guān)性。對(duì)于依賴b值改變的紋理特征及相關(guān)指數(shù)模型能否成為潛在生理意義的生物參數(shù)尚需進(jìn)一步討論。

[

]

[1] Hemming AW, Berumen J, Mekeel K. Hepatitis B and hepatocellular carcinoma. Clin Liver Dis, 2016,20(4):703-720.

[2] Aerts HJ, Velazquez ER, Leijenaar RT, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach. Nat Commun, 2014,5:4006.

[3] Qin JB, Liu Z, Zhang H, et al. Grading of gliomas by using radiomic features on multiple magnetic resonance imaging (MRI) sequences. Med Sci Monit, 2017,23(1):2168-2178.

[4] Liu J, Mao Y, Li Z, et al. Use of texture analysis based on contrast-enhanced MRI to predict treatment response to chemoradiotherapy in nasopharyngeal carcinoma. J Magn Reson Imaging, 2016,44(2):445-455.

[5] Becker AS, Wagner MW, Wurnig MC, et al. Diffusion-weighted imaging of the abdomen: Impact of b-values on texture analysis features. NMR Biomed, 2017,30(1):1-10.

[6] Shafiq-Ul-Hassan M, Zhang GG, Latifi K, et al. Intrinsic dependencies of CT radiomic features on voxel size and number of gray levels. Med Phys, 2017,44(3):1050-1062.

[7] Koh DM, Collins DJ. Diffusion-weighted MRI in the body: Applications and challenges in oncology. AJR Am J Roentgenol, 2007,188(6):1622-1635.

[8] Vandecaveye V, De Keyzer F, Verslype C, et al. Diffusion-weighted MRI provides additional value to conventional dynamic contrast-enhanced MRI for detection of hepatocellular carcinoma. Eur Radiol, 2009,19(10):2456-2466.

[9] Lambin P, Rios-Velazquez E, Leijenaar R, et al. Radiomics: Extracting more information from medical images using advanced feature analysis. Eur J Cancer, 2012,48(4):441-446.

[10] Yu H, Touret AS, Li B, et al. Application of texture analysis on parametric T1 and T2 maps for detection of hepatic fibrosis. J Magn Reson Imaging, 2017,45(1):250-259.

[11] Zhou W, Zhang L, Wang K, et al. Malignancy characterization of hepatocellular carcinomas based on texture analysis of contrast-enhanced MR images. J Magn Reson Imaging, 2017,45(5):1476-1484.

[12] 王波濤，劉剛，樊文萍，等.紋理特征分析在肝囊腫及肝血管瘤磁共振成像鑒別診斷中的價(jià)值．中國(guó)醫(yī)學(xué)科學(xué)院學(xué)報(bào)，2017，39(2)：169-176．

[13] 郭冬梅，劉惠，邵瑩，等.基于紋理特征的多序列MRI的肝硬化識(shí)別研究．中國(guó)醫(yī)學(xué)影像技術(shù)，2014，30(7)：1105-1108．

[14] Efremidis SC, Hytiroglou P. The multistep process of hepatocarcinogenesis in cirrhosis with imaging correlation. Eur Radiol, 2002,12(4):753-764.

[15] Wurnig MC, Donati OF, Ulbrich E, et al. Systematic analysis of the intravoxel incoherent motion threshold separating perfusion and diffusion effects: Proposal of a standardized algorithm. Magn Reson Med, 2015,74(5):1414-1422.

[16] Jirák D, Dezortová M, Taimr P, et al. Texture analysis of human liver. J Magn Reson Imaging, 2002,15(1):68-74.