基于定量影像组学的肺肿瘤良恶性预测方法

张利文 刘侠 汪俊 董迪 宋江典 臧亚丽 田捷

基于定量影像组学的肺肿瘤良恶性预测方法

张利文1,2刘侠1汪俊1董迪2宋江典3臧亚丽2田捷2

肺癌是世界范围内致死率最高的癌症之一,肺肿瘤的良恶性诊断对于治疗方式选择意义重大.本文借助影像组学(Radiomics)方法利用LIDC(Lung imaging database consortium)肺癌公开数据库中619例病人的肺癌计算机断层(Computed tomography,CT)影像数据,分割出病变区域,并结合肿瘤医学特性和临床认知,提取反映肿瘤形状大小、强度和纹理特性的60个定量影像特征,然后利用支持向量机(Support vector machine,SVM)构建诊断肺肿瘤良恶性的预测模型,筛选出对诊断肺肿瘤良恶性有价值的20个影像组学特征.为肺肿瘤良恶性预测提供了一种非入侵的检测手段.随着CT影像在肺癌临床诊断中的广泛使用,应用样本量的不断增加,本文方法有望成为一种辅助诊断工具,有效提高临床肺肿瘤良恶性诊断准确率.

影像组学,肺癌,图像分割,特征提取,支持向量机

目前,医学影像已经成为肿瘤诊断和指导治疗的主要手段之一,活体检查是诊断肿瘤良恶性的金标准,它经侵入式手术取出小部分肿瘤组织进行分析[1].但是由于肿瘤在空间和时间上的基因异质性,活检并不能全面评估肿瘤的信息,而且病人很难承受连续多次的活检[2].医学影像可以无创地提供整个肿瘤的信息,并且借助图像分析监测疾病的发生、发展及对治疗的反应[3],在指导治疗方面有巨大的潜力.肺癌是发病率和死亡率增长最快、对人类健康和生命威胁最大的恶性肿瘤之一.根据世界卫生组织的报告,到2030年全球将有1000万人死于肺癌[4].在肺癌诊断中最广泛使用的成像模态是X射线计算机断层(Computed tomography,CT)成像.肺癌在CT图像中呈现很强的对比度,包括图像中肿瘤灰度值强度差异、肿瘤内纹理差异和肿瘤形状差异.早期诊断对肺癌治疗方案的选择和病人生存期的延长有很大帮助[5].然而,结合肿瘤的医学特性和临床认知进行良恶性判断上并没有一个量化的标准.医生通过肿瘤图像的毛刺特征、肿瘤大小、是否有分叶、边缘清晰度等特征往往只能主观地做出初步判断,难于给出非常明确的诊断,必须对照活检病理信息才能给出最终的检查报告.往往不仅没有通过影像检查得到便捷的检查结果,还因为进一步的病理活检而延长检查周期.因此目前在临床实践上还缺乏一种具有高精度、高效率和非侵入式的诊断肺肿瘤良恶性的方法.

影像组学是一个新兴领域[6],其目标是从医学图像(CT扫描、正电子发射扫描或磁共振成像等方式获取的医学图像)中提取和分析大量具有高通量的定量图像特征,并构建描述肿瘤和预测临床表型的模型,建立图像特征和临床表型或基因分子标志的关联,进而进行肿瘤的诊断和临床表型预测[7].本文采用影像组学的方法构建肺肿瘤良恶性预测模型,为临床医生诊断肺癌提供帮助.

1 数据和方法

1.1 数据的选择

本文应用公开数据集LIDC(Lung imaging database consortium)[8]的619例肺癌病人的CT影像数据.其中426例用于构建预测模型的训练集,193例用于测试集.根据LIDC数据集提供的结节恶性度信息,分为良性和恶性两类.训练集包括213例恶性肿瘤和213例良性肿瘤.测试集中恶性肿瘤有122例,良性肿瘤有71例.

1.2 图像分割

图像分割是影像组学的关键步骤,肺结节病变区域的精准分割对于随后的特征提取和模型的构建十分重要.目前,放射科医生的手动分割通常被视为病灶分割的金标准.这种手动分割有很高的评估者个体差异,且非常耗时,效率低,重复性不强,不适合大数据集上的影像组学分析.影像组学分析需要一种自动、重复性高和精确的图像分割算法.

本文采用实验室前期研发的基于雪橇自动生长分割算法(Toboggan based growing automatic segmentation approach,TBGA)进行图像分割,该算法自动初始化种子点,不需要任何人工交互,且具有较高的精准度[9],为后续的影像组学过程提供了重要的数据基础,分割的结果如图1所示.

1.3 特征提取和特征选择

从肺结节医学特性和临床认知的角度来讲,医生主要从以下几点来判断肿瘤的病变程度:1)肺结节的大小:小于5.0mm的结节,良性可能性大:5.0∼10.0mm需要密切随访;10.0∼30.0mm的结节中恶性肿瘤占50∼80%以上.结节越大其为恶性的可能性越高,也就需要手术治疗;2)形态:有浅分叶、短细毛刺、胸膜凹陷特征的提示为肺癌;边缘光滑,有粗、长、稀疏毛刺的提示为良性.但有时很难判别;3)密度:密度不均匀,结节内有小泡、小管,或者有小结节堆积,有偏心小空腔的结节多为恶性；密度均一者多为良性;4)周围:周围肺组织呈现炎症、气肿、不张、侵润和支气管、血管束增粗的多为恶性结节;而有卫星灶、“晕环”的等多提示良性;5)强化:增强CT显示不均匀强化或是中度强化多为恶性,不强化或高度强化则多见于良性结节.以上特征只是一种经验化的结果,并没有量化成为固定的标准以给出确切的检查结果.本文结合以上肿瘤的医学特性和临床认知,从CT图像分割出来的肺部病变区域中总共提取了60个定量图像特征描述肿瘤的表型特征.特征可分为三组:1)形状大小特征,描述了肿瘤区域的几何形状及体积;2)影像强度特征,采用一阶统计量化的肿瘤强度特性,从所有肿瘤的像素强度直方图计算而得;3)纹理特征,可以量化肿瘤内异质差异.

图1 肺部病变区域分割Fig.1 Segmentation result of lung lesion

特征之间的相关性和冗余性会降低分类的准确率,同时医学图像通常属于小样本学习,特征过多将会提高分类器的复杂度,造成过拟合,降低分类器的泛化能力,因此需要对特征集合进行选择和优化[10].

本文采用最小冗余最大相关(Minimum redundancy maximum correlation,mRMR)方法进行特征选择[11].目的是在特征空间中找到与目标类别具有最大相关性且相互之间冗余度最小的包含m个特征的特征子空间,最小冗余最大相关的定义为

其中,S为特征集合,C为分类类别,式(1)选出与分类结果互信息均值最大的特征子集D(S,C),式(2)选出特征之间互信息均值最小的特征子集R(S).用特征i与分类类别C之间的互信息I(xi;C)度量特征i与分类结果的相关度,用集合中的特征i与特征j之间的互信息I(xi,xj)度量特征之间的相关性.

两个随机变量x和y的互信息I(x,y),根据其概率密度函数p(x),p(y)和联合分布概率密度函数P(x,y)计算为

可以得到最小冗余最大相关的特征选择标准为

式(4)表示应该选择与分类结果最大相关而与其他候选特征之间冗余度最小的特征.假定已确定一个拥有m个特征的特征集Sm,下一步需要从剩余特征集{S−Sm}中选择使得式(4)最大化的第m+1个特征为

特征选择的具体步骤:

步骤1.从全部候选特征空间中选出最小冗余最大相关前50个特征,形成50个特征子集,子集中的特征个数从1依次递增到50.

步骤2.计算每个特征子集在十折交叉验证下预测错误率.

步骤3.选出预测错误率最小时对应特征子集.

1.4 建立良恶性诊断模型

使用最小冗余最大相关方法选取特征,将训练集特征空间作为输入,然后使用遗传算法对支持向量机(Support vector machine,SVM)的参数进行优化,最后利用SVM 训练出肺部肿瘤良恶性预测模型,并在训练集上对模型的预测能力进行验证.要求在具有良好预测精确度的同时,模型的复杂度尽可能低,模型构建和验证的具体步骤如图2所示.

图2 肺肿瘤良恶性预测模型的生成和验证示意图Fig.2 Diagram of generation and validation of benign and malignant discrimination model of lung tumor

支持向量机基于结构风险最小化原理(Structural risk minimization,SRM),为了获得更好的泛化能力,需要控制经验风险和置信范围值两个因素[12].结构风险最小化在训练误差约束下,使置信范围最小化.它最终求解的是一个线性约束条件下的凸二次规划问题,所以支持向量机的解是全局最优的唯一解[13].

为了获取更高的预测准确率,本文对SVM的参数进行了优化.在参数优化前,对训练集和验证集数据做归一化处理.考虑到基本的线性归一化处理没有改变变量在不同样本间取值的离散程度,同时,因为这是一种对数据的无量纲化、标准化操作,可以使模型拟合的过程中在精度不变的情况下参数更易收敛[14].归一化公式如下:

SVM的分类性能受到多种因素影响,其中两个关键的因素为:1)误差惩罚参数C;2)核函数形式及其参数g.误差惩罚参数通过调节特征子空间中置信范围和经验风险使学习机器的泛化能力达到最好.在确定的特征子空间中,C取值代表了对经验误差惩罚程度,C的值小说明对误差惩罚程度小,也表明预测模型的复杂度小而经验风险值较大,出现欠拟合现象;反之则会出现过拟合现象.在每个特征子空间中,都存在使SVM 获得最好的泛化能力的C值.随着C的增大,SVM的复杂度不断增大直至数据子空间允许的最大复杂度,而经验风险和泛化能力几乎不再变化.不同核函数对分类性能有影响[15],相同的核函数其参数不同也会造成分类性能的不同.核函数及其参数的选择直接影响到SVM分类性能的好坏.

SVM常用的核函数有:1)多项式核函数:

2)径向基核函数

3)Sigmoid核函数(两层神经网络)

径向基核函数具有非线性、参数少、能够把原始特征映射到无限维,本文选择径向基核函数作为支持向量机的核函数.

本文采用遗传算法(Genetic algorithm,GA)优化误差惩罚参数C和核函数参数g.GA基于进化论的适者生存、优胜劣汰理论,使用遗传学的基本操作在种群中反复求解,从而不断地产生新的种群,根据个体适应度不断进行优胜劣汰,使种群得到进化,最终求得最优解[16].由于遗传算法的操作对象是编码后的结构,解除了求解过程中要求可求导以及函数必须连续性等的限制.如果采用概率化的寻优方式,则寻优过程不需要确定的规则,搜索方向自适应调整,搜索策略能够实现自动获取搜索空间.

在遗传算法中,每个个体将要优化问题的所有参数进行编码,形成一个有限长度的字符串.每个个体是优化问题的一个可行解.种群是遗传算法的搜索空间,由多个个体组成.开始搜索时先随机选取一定数量个体的初始种群.然后从初始种群开始运用优胜劣汰的准则进行个体的选择,就像自然界进化的过程.计算个体的适应度,适应度高的个体得以生存,适应度低的个体将被淘汰,所有幸存的个体组成新的种群,最后利用交叉、变异等遗传操作将留下来的种群特性遗传给下一代.如此不断重复“选择–交叉–变异–再选择”的过程,下一代不断积累上代遗传下来的好的特性,个体的适应度和种群的平均适应度不断增加,直到满足迭代结束条件[17].

遗传算法优化支持向量机参数的主要流程如图3所示,将支持向量机在十折交叉验证下的预测准确率作为目标函数值,即个体的适应度,个体对应的准确率越高,则个体的适应度越大,K为进化代数.

图3 遗传算法流程图Fig.3 Flow chart of genetic algorithm

2 结果

本文构建了基于定量影像组学的肺肿瘤良恶性预测方法,采用LIDC数据集上的619例肺癌病人的CT影像数据进行分割,对分割后的每例数据分别提取三类,共60个定量图像特征、形状大小、影像强度和纹理特征.使用最小冗余最大相关选取不同的特征空间,确定一个具有最高预测准确度的优化特征空间,该特征空间包括20个量化图像特征.然后采用“最小–最大标准化”对选择后的特征数据进行归一化处理.在相同数据集上分别使用归一化后的和原始数据进行预测,比较两者差异,发现归一化后数据求解速度明显比未归一化数据快很多.为了使SVM分类准确率更高,在训练之前利用GA算法对SVM参数优化,得到最优C为1.0052,g为2.1480,GA算法的适应度曲线见图4.利用SVM训练出相应的肺肿瘤良恶性预测模型,对模型预测的灵敏度、特异性、阳性预测值、阴性预测值和总体准确度进行计算.其中模型在训练集上的总体准确度为82.4%(351/426),灵敏度为82.2%(175/213),特异性为82.6%(176/213),阳性预测值为82.5%(175/212),阴性预测值82.2%(176/214).在测试集上的总体准确度为77.7%(150/193),灵敏度为75.4%(92/122),特异性为81.6%(58/71),阳性预测值为87.6%(92/105),阴性预测值65.9%(58/88),详细数据均列在表1中.

表1 肺肿瘤良恶性预测模型的诊断准确率Table 1 Diagnostic accuracy of benign and malignant discrimination model of lung tumor

图4 遗传算法适应度曲线Fig.4 Fitness curve of genetic algorithm

3 总结

实验结果显示本文方法在训练集上的分类准确度达到82.4%,在测试集上的分类准确度为77.7%.本文构建的肺癌影像组学预测模型应用于肺癌良恶性的定量预判,辅助临床医生进行诊断.本文方法有两个优点:1)应用非侵入CT肺部图像从整个肿瘤层面预测肺肿瘤的良恶性;2)肺部的CT图像容易获取,利用肺癌影像组学预测模型进行良恶性评估易于临床应用.

本文方法主要针对肺癌良恶性的预判,也可应用于其他癌症的病理分型和预后预测[18].在后续工作中,将更紧密结合肿瘤医学和临床等相关知识,对初步开发的肿瘤良恶性预测软件进一步完善,在各省市周边医院推广来辅助当地医生对患者完成高效、方便和低费用的检查.由于国内现阶段病例数很难收集,训练数据限于公共数据集的426个,以后需要补充更多的多中心数据进行模型训练,提高预测模型的分类精度.

1 Lambin P,Rios-Velazquez E,Leijenaar R,Carvalho S,van Stiphout R G P M,Granton P,Zegers C M L,Gillies R,Boellard R,Dekker A,Aerts H J W L.Radiomics:extracting more information from medical images using advanced feature analysis.European Journal of Cancer,2012,48(4):441−446

2 Aerts H J W L,Velazquez E R,Leijenaar R T H,Parmar C,Grossmann P,Carvalho S,Bussink J,Monshouwer R,Haibe-Kains B,Rietveld D,Hoebers F,Rietbergen M M,Leemans C R,Dekker A,Quackenbush J,Gillies R J,Lambin P.Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach.Nature Communications,2014,5:Article No.4006

3 Huang Y Q,Liang C H,He L,Tian J,Liang C S,Chen X,Ma Z L,Liu Z Y.Development and validation of a radiomics nomogram for preoperative prediction of lymph node metastasis in colorectal cancer.Journal of Clinical Oncology,2016,34(18):2157−2164

4 World Health Organization.Description of the global burden of NCDs,their risk factors and determinants.Global Status Report on Noncommunicable Diseases 2010.World Health Organization,2011.1−176

5 Parmar C,Leijenaar R,Grossmann P,Velazquez E R,Bussink J,Rietveld D,Rietbergen M M,Haibe-Kains B,Lambin P,Aerts H J W L.Radiomic feature clusters and prognostic signatures speci fi c for lung and head and neck cancer.Scienti fi c Reports,2015,5:Article No.11044

6 Kumar V,Gu Y H,Basu S,Berglund A,Eschrich S A,Schabath M B,Forster K,Aerts H J W L,Dekker A,Fenstermacher D,Goldgof D B,Hall L O,Lambin P,Balagurunathan Y,Gatenby R A,Gillies R J.Radiomics:the process and the challenges.Magnetic Resonance Imaging,2012,30(9):1234−1248

7 Gillies B,Kinahan P E,Hricak H.Radiomics:images are more than pictures,they are data.Radiology,2016,278(2):563−577

8 Armato S G,McLennan G,Bidaut L,McNitt-Gray M F,Meyer C R,Reeves A P,Zhao B S,Aberle D R,Henschke C I,Hoあman E A,Kazerooni E A,MacMahon H,van Beek E J R,Yankelevitz D,Biancardi A M,Bland P H,Brown M S,Engelmann R M,Laderach G E,Max D,Pais R C,Qing D P Y,Roberts R Y,Smith A R,Starkey A,Batra P,Caligiuri P,Farooqi A,Gladish G W,Jude C M,Munden R F,Petkovska I,Quint L E,Schwartz L H,Sundaram B,Dodd L E,Fenimore C,Gur D,Petrick N,Freymann J,Kirby J,Hughes B,Vande C A,Gupte S,Sallam M,Heath M D,Kuhn M H,Dharaiya E,Burns R,Fryd D S,SalganicoあM,Anand V,Shreter U,Vastagh S,Croft B Y,Clarke L P.The lung image database consortium(LIDC)and image database resource initiative(IDRI):a completed reference database of lung nodules on CT scans.Medical Physics,2011,38(2):915−931

9 Song J,Yang C,Fan L,Wang K,Yang F,Liu S,Tian J.Lung lesion extraction using a toboggan based growing automatic segmentation approach.IEEE Transactions on Medical Imaging,2016,35(1):337−353

10 Dougherty E R,Hua J P,Sima C.Performance of feature selection methods.Current Genomics,2009,10(6):365−374

11 Peng H C,Long F H,Ding C.Feature selection based on mutual information criteria of max-dependency,max-relevance,and min-redundancy.IEEE Transactions on Pattern Analysis and Machine Intelligence,2005,27(8):1226−1238

12 Vapnik V.The Nature of Statistical Learning Theory.New York:Springer Science and Business Media,2013.123−167

13 Zhang Xue-Gong.Introduction to statistical learning theory and support vector machines.Acta Automatica Sinica,2000,26(1):32−42(张学工.关于统计学习理论与支持向量机.自动化学报,2000,26(1):32−42)

14 Chang C C,Lin C J.LIBSVM:a library for support vector machines.ACM Transactions on Intelligent Systems and Technology(TIST),2011,2(3):Article No.27

15 Zhang Kai-Jun,Liang Xun.An improved domain multiple kernel support vector machine.Acta Automatica Sinica,2014,40(10):2288−2294(张凯军,梁循.一种改进的显性多核支持向量机.自动化学报,2014,40(10):2288−2294)

16 Whitley D,Starkweather T,Bogart C.Genetic algorithms and neural networks:optimizing connections and connectivity.Parallel Computing,1990,14(3):347−361

17 Deb K,Pratap A,Agarwal S,Meyarivan T.A fast and elitist multiobjective genetic algorithm:NSGA-II.IEEE Transactions on Evolutionary Computation,2002,6(2):182−197

18 Song J D,Liu Z Y,Zhong W Z,Huang Y Q,Ma Z L,Dong D,Liang C H,Tian J.Non-small cell lung cancer:quantitative phenotypic analysis of CT images as a potential marker of prognosis.Scienti fi c Reports,2016,6:Article No.38282

Prediction of Malignant and Benign Lung Tumors Using a Quantitative Radiomic Method

ZHANG Li-Wen1,2LIU Xia1WANG Jun1DONG Di2SONG Jiang-Dian3ZANG Ya-Li2TIAN Jie2

Lung cancer is a leading cause of cancer mortality around the world.Accurate diagnosis of lung cancer is signi fi cant for treatment regimen selection.Radiomics refers to comprehensively quantifying the tumor phenotypes by applying a large number of quantitative image features.Here we analyze a computed tomography(CT)data set of 619 patients with lung cancer on the lung image database consortium(LIDC)by radiomic method.Combining with the medical character and clinical recognition of lung tumor,we present a radiomic analysis of 60 features.Then,we use SVM to build a prediction model and fi nd radiomic features which have predictive value for discrimination of malignant and benign lung tumors.Nowadays,as CT imaging is routinely used in lung cancer clinical diagnosis,there is an increase in data set size.We consider that our radiomic prediction model will be developed a valuable medical software and an auxiliary tool which can provide malignant and benign information of lung tumors eきciently.

Radiomics,lung cancer,image segmentation,feature extraction,support vector machine(SVM)

Zhang Li-Wen,Liu Xia,Wang Jun,Dong Di,Song Jiang-Dian,Zang Ya-Li,Tian Jie.Prediction of malignant and benign lung tumors using a quantitative radiomic method.Acta Automatica Sinica,2017,43(12):2109−2114

2016-03-11 录用日期2017-01-04

March 11,2016;accepted January 4,2017国家自然科学基金(81227901,81527805,61231004,81501616,813 01346,61672197),黑龙江省然科学基金(F201311,12541105),中国科学院科技服务网络计划(KFJ-SW-STS-160),中国科学院科研设备项目(YZ201502)资助

Supported by National Natural Science Foundation of China(81227901,81527805,61231004,81501616,81301346,61672197),Natural Science Foundation of Heilongjiang Province(F201311,12541105),Chinese Academy of Science Program of Scienti fi c Service Network(KFJ-SW-STS-160),and Chinese Academy of Science Program of Equipment(YZ201502)

本文责任编委朱朝喆

Recommended by Associate Editor ZHU Chao-Zhe

1.哈尔滨理工大学自动化学院哈尔滨150080 2.中国科学院自动化研究所北京100190 3.东北大学中荷生物医学与信息工程学院沈阳110819

1.School of Automation,Harbin University of Science and Technology,Harbin 150080 2.Institute of Automation,Chinese Academy of Sciences,Beijing 100190 3.School of Sino-Dutch Biomedical and Information Engineering,Northeastern University,Shengyang 110819

张利文,刘侠,汪俊,董迪,宋江典,臧亚丽,田捷.基于定量影像组学的肺肿瘤良恶性预测方法.自动化学报,2017,43(12):2109−2114

DOI10.16383/j.aas.2017.c160264

张利文 中国科学院自动化研究所和哈尔滨理工大学自动化学院联合培养硕士研究生.主要研究方向为医学图像处理.E-mail:zhangliwen2015@ia.ac.cn

(ZHANG Li-Wen Master student at the Institute of Automation,Chinese Academy of Sciences,and School of Automation,Harbin University of Science and Technology.His main research interest is medical imaging processing.)

刘 侠 哈尔滨理工大学自动化学院教授.2006年获得哈尔滨工程大学自动化学院博士学位.主要研究方向为模式识别与智能系统.

E-mail:liuxia@hrbust.edu.cn

(LIU Xia Professor at the School of Automation,Harbin University of Science and Technology.He received his Ph.D.degree from Harbin Engineering University in 2006.His research interest covers pattern recognition and intelligent systems.)

汪 俊 哈尔滨理工大学测控技术与通信工程学院硕士研究生.主要研究方向为图像处理和模式识别.

E-mail:wangjun.542@163.com

(WANG Jun Master student at the School of Measurement-Control Technology and Communications Engineering,Harbin University of Science and Technology.His research interest covers imaging processing and pattern recognition.)

董 迪 中国科学院自动化研究所副研究员.2013年获得中国科学院自动化研究所博士学位.主要研究方向为影像组学.E-mail:di.dong@ia.ac.cn

(DONG DiAssociate professor at the Institute of Automation,Chinese Academy of Sciences.He received his Ph.D.degree from Chinese Academy of Sciences in 2013.His research interest covers radiomics.)

宋江典 东北大学中荷生物医学与信息工程学院博士研究生,IEEE会员.主要研究方向为医学影像处理与分析.

E-mail:dr.j.song@ieee.org

(SONGJiang-Dian Ph.D.candidate at Biomedical Engineering,School of Sino-Dutch Biomedical and Information Engineering,Northeastern University,IEEE member.His research interest covers medical image processing and analysis.)

臧亚丽 中国科学院自动化研究所副研究员.2013年获得中国科学院自动化研究所博士学位.主要研究方向为影像组学.E-mail:yali.zang@ia.ac.cn

(ZANG Ya-LiAssociate professor at the Institute of Automation,Chinese Academy of Sciences.She received her Ph.D.degree from the Institute of Automation,Chinese Academy of Sciences in 2013.Her research interest covers radiomics.)

田 捷 中国科学院自动化研究所研究员.1992年获得中国科学院自动化研究所博士学位.主要研究方向为模式识别,医学图像处理和分子影像.本文通信作者.E-mail:tian@ieee.org

(TIAN Jie Professor at the Institute of Automation,Chinese Academy of Sciences.He received his Ph.D.degree from the Institute of Automation,Chinese Academy of Sciences in 1992.His research interest covers pattern recognition,medical image processing,and molecular imaging.Corresponding author of this paper.)