肿瘤防治研究  2018, Vol. 45 Issue (12): 1009-1013
本刊由国家卫生和计划生育委员会主管,湖北省卫生厅、中国抗癌协会、湖北省肿瘤医院主办。
0

文章信息

基于Tirm序列联合直方图及纹理特征分析鉴别乳腺良恶性病变
Differential Diagnosis of Benign and Malignant Breast Lesions by Texture Feature Combined with Histogram Parameter Based on Tirm Sequence
肿瘤防治研究, 2018, 45(12): 1009-1013
Cancer Research on Prevention and Treatment, 2018, 45(12): 1009-1013
http://www.zlfzyj.com/CN/10.3971/j.issn.1000-8578.2018.18.0815
收稿日期: 2018-06-15
修回日期: 2018-09-25
基于Tirm序列联合直方图及纹理特征分析鉴别乳腺良恶性病变
李晓峰 ,    刘彩云 ,    石岩 ,    罗小虎 ,    王刚 ,    姚标     
221004 徐州,徐州市肿瘤医院放射科
摘要: 目的 分析乳腺良恶性病变在T2反转恢复压脂序列上的直方图及纹理特征,并评价诊断乳腺良恶性病变的效能。方法 回顾性分析100例行MRI检查并经病理证实的乳腺肿瘤资料。纹理特征提取采用灰度共生矩阵法得到病灶能量值、熵值、惯量值、自相关性、逆差矩及群显著性;直方图分析得到平均值、偏度、峰度。分别评价不同参数值在乳腺良恶性病变中的诊断效能。构建受试者工作特征曲线(ROC)结合临床实际评价乳腺良恶性组间的差异,确定诊断临界值。结果 各参数在良恶性组间比较,自相关性、群显著性、平均值及峰度值差异均有统计学意义(均P < 0.001);能量值、熵值、惯量值、逆差距及偏度值差异无统计学意义(P > 0.05)。构建ROC将自相关性、群显著性、平均值及峰度值分别用于诊断乳腺良恶性病变,均有统计学意义(P < 0.001),联合四组参数诊断良恶性病变,曲线下面积为0.868,敏感度为90.57%,特异性为72.34%。结论 基于Tirm序列的直方图及纹理分析可用于乳腺良恶性病变的鉴别诊断,自相关性、群显著性、平均值及峰度值在Tirm序列上对诊断乳腺良恶性病变有一定意义,联合诊断可提高诊断效能。
关键词: 磁共振成像     乳腺肿瘤     直方图分析     纹理分析     鉴别诊断    
Differential Diagnosis of Benign and Malignant Breast Lesions by Texture Feature Combined with Histogram Parameter Based on Tirm Sequence
LI Xiaofeng , LIU Caiyun , SHI Yan , LUO Xiaohu , WANG Gang , YAO Biao     
Department of Radiology, Xuzhou Cancer Hospital, Xuzhou 221004, China
Corresponding author: SHI Yan, E-mail:xzshiyan@163.com.
Abstract: Objective To explore histogram parameters and texture features of benign and malignant breast lesion on turbo inversion recovery magnitude(Tirm) sequence, and evaluate which parameter could help best differentiate benign from malignant breast lesion. Methods This retrospective study included 100 breast cancer patients who underwent conventional MRI and confirmed pathologically. Texture features were derived from the gray level co-occurrence matrix(GLCM), and entropy, energy, correlation, inertia, inverse difference moment, cluster prominence and mean value, skewness, kurtosis of histogram parameters were calculated. Then we assessed the diagnosis efficacy with these parameters among the variety kinds of benign and malignant breast lesions respectively, and establish the receiver-operating characteristic curve(ROC). We assessed the differences between benign and malignant groups by the Yoden index combined with clinic for the cut-off values. Results The differences of correlation, cluster prominence, mean value and kurtosis were statistically significant between the benign and malignant groups(all P < 0.001); the differences of energy, entropy, inertia, inverse difference moment and skewness were not statistically significant(all P > 0.05). The results of ROC with correlation, cluster prominence, mean value and kurtosis on the diagnosis of benign and malignant breast lesions were statistically significant, respectively(P < 0.001). Combined the four parameters on diagnosis benign and malignant lesions, the AUC was 0.868, the sensitivity was 90.57% and the specificity was 72.34%. Conclusion The histogram analysis and texture analysis based on Tirm sequence could be used for the differential diagnosis of benign and malignant breast lesions. Correlation, cluster prominence, mean value and kurtosis have certain significance in the diagnosis. The combined diagnosis could improve the differential ability of benign and malignant breast lesions.
Key words: MRI     Breast tumor     Histogram analysis     Texture analysis     Differential diagnosis    
0 引言

乳腺癌作为女性常见恶性肿瘤之一,其良恶性鉴别诊断已成为当下研究热点。利用动态增强MRI(dynamic contrast-enhanced MRI, DCE-MRI)定量分析等研究已广泛开展。近几年,影像组学正受到越来越多学者的关注,其中直方图及纹理分析作为一种重要的研究方法,在临床工作中的作用及科研价值也越来越大。本试验旨在利用直方图及纹理特征分析在T2反转恢复压脂(turbo inversion recovery magnitude, Tirm)序列上鉴别诊断乳腺良恶性病变。

1 资料与方法 1.1 研究对象

回顾性分析徐州市肿瘤医院2016年2月—2018年1月间经病理确诊的术前行MRI扫描的乳腺病变患者100例(年龄26~79岁,中位年龄52岁,良性:恶性=47:53)。病变组纳入标准:(1)良性病变组包括囊肿、纤维腺瘤、腺病;恶性病变组包括黏液癌、浸润癌及导管内癌;(2)全部病灶最大截面直径 > 1 cm;(3)检查前未经手术、放疗、化疗及靶向药物治疗者;(4)图像无明显运动伪影。排除标准:(1)年龄在25岁以下或大于80岁的患者;(2)怀孕、哺乳期妇女;(3)患有严重心、肝、肾等器质性病变;(4)MR检查禁忌证。本研究已通过徐州市肿瘤医院伦理委员会审查,所有患者检查前均已签署知情同意书。

1.2 MRI扫描方法

使用Siemens MAGNETOM Avanto 1.5T超导磁共振仪,8通道乳腺专用矩阵线圈进行扫描,患者取俯卧位。将双乳自然悬垂于乳腺线圈内;并嘱患者在整个检查过程中保持平静呼吸。Tirm序列具体参数:TR 3700 ms,TE 85 ms,TI 160 ms,层厚3.0 mm,矩阵256×256;FOV 220 mm×220 mm,距离因子10%,层数13,平均次数2,总扫描时间为174 s。

1.3 病理结果

所有患者均于检查后3天内手术。手术切除病变标本后,经中性福尔马林液固定,石蜡包埋、HE染色(部分实验结果见图 1~3)。经病理会诊后得到病理诊断结论。

图 1 64岁女性左乳浸润性导管癌MRI图像(A)及病灶直方图(B)、显微镜下表现(C) Figure 1 MRI images(A), histogram(B) and microscopic findings(C) of invasive ductal carcinoma in left breast of a 64-yearold female patient

图 2 47岁女性左乳纤维腺瘤MRI图像(A)及病灶直方图(B)、显微镜下表现(C) Figure 2 MRI images(A), histogram(B) and microscopic findings(C) of fibroadenoma in left breast of a 47-year-old female patient

图 3 29岁女性左乳黏液癌MRI图像(A)及病灶直方图(B)、显微镜下表现(C) Figure 3 MRI images(A), histogram(B) and microscopic findings(C) of mucinous breast carcinoma in left breast of a 29-yearold female patient
1.4 MRI图像分析

将得到的所有MRI图像经工作站处理并导入血流动力学定量分析软件(GE Omni-Kinetics Version V2.20软件)进行分析,纹理分析基于灰度共生矩阵(gray level co-occurrence matrix, GLCM),得到能量值、熵值、惯量值、逆差矩及自相关性、群显著性;直方图特征参数包括平均值、偏度、峰度。每个病灶选取病灶最大截面层面及邻近上下各一层共3层图像,建立3D-ROI数据。

1.5 统计学方法

应用SPSS25.0统计学软件对数据进行分析。正态分布数据描述以频数、均数±标准差表示,非正态分布以中位数+四分位间距表示。良恶性组间各参数比较采用两独立样本的非参数检验方法;此外,构建受试者工作特征曲线(receiver operating characteristic curve, ROC)并计算曲线下面积,根据约登指数结合临床实际评价乳腺良、恶性组间的差异,确定诊断临界值。所有统计分析均选取α=0.05作为假设检验水准进行双侧检验。

2 结果 2.1 形态学特征

100例病灶中,良性病变中:纤维腺瘤25例、囊肿17例、腺病5例;45例病变境界清晰、2例腺病表现为边缘及境界欠清晰。恶性病变中:浸润性癌27例、导管内癌8例、黏液癌18例;53例均表现为边缘模糊,境界不清晰,并见毛刺形成。在良性组、恶性组中,能量值、熵值、惯量值、自相关性、逆差距、群显著性、平均值、偏度及峰度的中位数、四分位间距以及两组间的Z值、P值见表 1

表 1 直方图参数及纹理参数在良恶性组间比较 Table 1 Comparison of histogram parameters and texture parameters between benign and malignant groups
2.2 自相关性、群显著性、平均值及峰度值分别鉴别乳腺良恶性病变的诊断价值

ROC分析结果显示:自相关性的最佳诊断界值为4.935×10-5,此时敏感度为92.45%,特异性为55.32%,约登指数为0.478。峰度值的最佳诊断界值为3.000,此时敏感度为75.47%,特异性为70.21%,约登指数为0.457。群显著性的最佳诊断界值为1.063×1010,此时敏感度为77.36%,特异性为72.34%,约登指数为0.497。平均值的最佳诊断界值为1.753×102,此时敏感度为58.49%,特异性为100%,约登指数为0.585,见图 4表 2

图 4 自相关性、群显著性及平均值、峰度值在良恶性乳腺病变鉴别诊断的ROC Figure 4 ROC of correlation, cluster prominence, mean value and kurtosis in differential diagnosis of benign and malignant breast lesions

表 2 自相关性、群显著性及平均值、峰度值的ROC曲线下面积 Table 2 Area under the curve(AUC) of correlation, cluster prominence, mean value and kurtosis
2.3 联合自相关性、峰度值、群显著性及平均值4个指标鉴别乳腺良恶性病变

ROC分析结果显示:将自相关性、峰度值、群显著性及平均值作为联合指标预测乳腺良恶性病变,则约登指数为0.629,敏感度为90.57%,特异性为72.34%,见图 5表 3

图 5 联合自相关性、群显著性及平均值、峰度值指标在良恶性乳腺病变鉴别诊断的ROC Figure 5 ROC of correlation, cluster prominence, mean value and kurtosis in differential diagnosis of benign and malignant breast lesions

表 3 联合自相关性、群显著性及平均值、峰度值指标的ROC曲线下面积 Table 3 AUC of combination of correlation, cluster prominence, mean value and kurtosis parameters
3 讨论

乳腺癌早期诊断、分期及疗效评估作为乳腺癌重要的研究方向,一直以来备受关注[1]。目前,利用直方图及纹理分析等方法对乳腺肿瘤鉴别诊断及疗效评估已成为影像组学研究的重要组成部分[2-4]。Ko等[5]已发现乳腺癌肿块异质性与乳腺癌MR图像纹理特征具有一定的相关性。国内亦有文献报道[6]运用DCE-MRI纹理分析对乳腺癌新辅助化疗的患者进行预测评价时发现,能量值、熵值、惯量值和自相关性对于乳腺癌新辅助化疗后早期疗效评估具有一定的价值。Prevos等[7]的研究也证实了这一点。近几年,DCE-MRI作为乳腺癌MRI研究的热点,其最大优势在于能够完整描述乳腺病变整体的血流动力学特点,因而也成为当前研究乳腺癌的主要方法[8]。然而乳腺病变患者在DCE-MRI图像采集过程中时间分辨率要求较高[9],Tirm序列能够作为其重要补充,我们认为其优势首先在于Tirm序列相较于DCE-MRI图像质量能够得到保证,且对于病变的细节显示较常规T1WI及T2WI好;其次,Tirm序列由于脂肪抑制技术的应用使富含脂肪成分的乳腺组织显示对比度较高,也在一定程度上提高了病变的检出率。然而尽管如此,Tirm序列也只是作为常规MRI扫描的一部分,因此提供的信息没有DCE-MRI丰富,这也成为本次研究的不足。在今后的研究中将不同序列进行比较从而评价哪一个参数的诊断价值更高也是重要的研究方向。同时纹理分析及直方图分析作为一种定量分析方法,能够更加准确描述病变的生理信息及血流动力学特征,从而对乳腺良恶性病变的术前鉴别诊断提供更多有价值的信息。张竹伟等[10]利用T2WI序列进行纹理分析得到了较为理想的结果。

本研究中,结合乳腺良恶性组数据特点我们可以推测对于自相关性及峰度值这两个指标采用较大的检验结果表示对乳腺癌有一定的阳性诊断价值:即当患者的自相关性及峰度值大于诊断临界值时更倾向于诊断乳腺癌;对群显著性和平均值小于诊断临界值时更倾向于诊断为乳腺癌。此外,本组病例中能量值、熵值、惯量值及偏度值在乳腺良恶性病变间的差异无统计学意义,说明这四个参数对乳腺良恶性病变Tirm序列上的鉴别无明显提示作用。推测可能原因为乳腺良性病变中纤维腺瘤、腺病与乳腺癌在Tirm序列上的表现存在部分重叠。本研究中乳腺黏液癌与浸润性导管癌的直方图分布不一致,纤维腺瘤直方图形态与黏液癌的直方图形态也存在一定的差别,说明直方图分析在鉴别乳腺良恶性病变及乳腺癌不同病理类型中可能具有一定的价值。类似的结果在冯红梅等[11]研究中也得到了验证。本研究中,乳腺良恶性病变的偏度均为负,其中乳腺癌的偏度较小,可能与乳腺癌中存在微小钙化有关,这与Ko等[5]在研究中发现的熵值在T2WI及增强的T1WI上诊断浸润性乳腺癌具有一定的价值也是一致的。除此之外,我们发现2例腺病病灶境界不清晰,这可能与腺病引起腺体变形、纤维化,并且在MRI上表现与浸润性癌类似有关。

本研究数据经后处理软件分析可得到57个参数,我们仅选取部分较常用参数进行比较分析,且该研究中尚未对纹理特征及直方图参数进行选择以排除对ROI选择敏感的参数,这将在以后的研究中逐步完善。再者,本研究的诊断阈值仍需要大样本量数据的验证。此外,该研究仅对Tirm序列上乳腺良恶性病变的直方图及纹理的特征进行分析,未做其他序列的对照研究,故也具有一定的局限性。综上所述,常规MRI图像纹理分析可为鉴别乳腺良恶性病变提供可靠的信息,作为乳腺DCE-MRI序列的补充,Tirm序列对于鉴别乳腺良恶性病变的纹理特征也能够提供有价值的信息。

参考文献
[1] Gallego-Ortiz C, Martel AL. Using quantitative features extracted from T2-weighted MRI to improve breast MRI computer-aided diagnosis (CAD)[J]. PLoS One, 2017, 12(11): e0187501. DOI:10.1371/journal.pone.0187501
[2] Imbriaco M, Cuocolo R. Does Texture Analysis of MR Images of Breast Tumors Help Predict Response to Treatment?[J]. Radiology, 2018, 286(2): 421–3. DOI:10.1148/radiol.2017172454
[3] Wu J, Gong G, Cui Y, et al. Intratumor partitioning and texture analysis of dynamic contrast-enhanced (DCE)-MRI identifies relevant tumor subregions to predict pathological response of breast cancer to neoadjuvant chemotherapy[J]. J Magn Reson Imaging, 2016, 44(5): 1107–15. DOI:10.1002/jmri.25279
[4] Pickles MD, Lowry M, Gibbs P. Pretreatment Prognostic Value of Dynamic Contrast-Enhanced Magnetic Resonance Imaging Vascular, Texture, Shape, and Size Parameters Compared With Traditional Survival Indicators Obtained From Locally Advanced Breast Cancer Patients[J]. Invest Radiol, 2016, 51(3): 177–85. DOI:10.1097/RLI.0000000000000222
[5] Ko ES, Kim JH, Lim Y, et al. Assessment of Invasive Breast Cancer Heterogeneity Using Whole-Tumor Magnetic Resonance Imaging Texture Analysis[J]. Medicine(Baltimore), 2016, 95(3): e2453.
[6] 孙赛花, 周纯武, 赵莉芸, 等. 动态增强磁共振成像纹理分析预测乳腺癌新辅助化疗疗效[J]. 中华肿瘤杂志, 2017, 39(5): 344–9. [ Sun SH, Zhou CW, Zhao LY, et al. Texture analysis based on contrast-enhanced MRI can predict treatment response to neoadjuvant chemotherapy of breast cancer[J]. Zhonghua Zhong Liu Za Zhi, 2017, 39(5): 344–9. DOI:10.3760/cma.j.issn.0253-3766.2017.05.005 ]
[7] Prevos R, Smidt ML, Tjan-Heijnen VC, et al. Pre-treatment differences and early response monitoring of neoadjuvant chemotherapy in breast cancer patients using magnetic resonance imaging: a systematic review[J]. Eur Radiol, 2012, 22(12): 2607–16. DOI:10.1007/s00330-012-2653-5
[8] Waugh SA, Purdie CA, Jordan LB, et al. Magnetic resonance imaging texture analysis classification of primary breast cancer[J]. Eur Radiol, 2016, 26(2): 322–30.
[9] Salem A, O'Connor JPB. Assessment of Tumor Angiogenesis: Dynamic Contrast- enhanced MR Imaging and Beyond[J]. Magn Reson Imaging Clin N Am, 2016, 24(1): 45–56. DOI:10.1016/j.mric.2015.08.010
[10] 张竹伟, 华婷, 徐婷婷, 等. 常规MRI纹理分析鉴别乳腺良、恶性病变的价值初探[J]. 中华放射学杂志, 2017, 51(8): 588–91. [ Zhang ZW, Hua T, Xu TT, et al. Differentiation of benign and malignant breast lesions using texture analysis of conventional MRI:a preliminary study[J]. Zhonghua Fang She Xue Za Zhi, 2017, 51(8): 588–91. DOI:10.3760/cma.j.issn.1005-1201.2017.08.006 ]
[11] 冯红梅, 郭彩平, 徐志锋, 等. 乳腺X线摄影和MRI直方图在鉴别乳腺纤维腺瘤和浸润性导管癌中的价值[J]. 医学影像学杂志, 2017, 27(1): 75–8. [ Feng HM, Guo CP, Xu ZF, et al. The value of X-ray photography and MRI histogram in distinguishing of breast fibroadenoma and infiltrating ductal car-cinoma[J]. Yi Xue Ying Xiang Xue Za Zhi, 2017, 27(1): 75–8. ]