一种基于CT影像学特征与非小细胞肺癌患者预后情况的相关性分析方法技术

本发明专利技术提供了一种基于CT影像学特征与非小细胞肺癌患者预后生存情况的相关性分析方法，属于计算机辅助医学领域。本发明专利技术首先对肺部肿瘤进行分割和特征提取、优化；然后，通过机器学习的方法将肿瘤影像学特征与患者生存状况结合，构建非小细胞肺癌的预后评估模型来对患者的预后生存情况进行预测，同时设计实验对预后分析模型的性能进行评估。本发明专利技术为探寻影像学特征和病人预后生存情况的关系提供了一种新的方案。

A correlation analysis method based on CT imaging features and prognosis of patients with non-small cell lung cancer

The invention provides a correlation analysis method based on CT imaging characteristics and prognosis survival of patients with non-small cell lung cancer, belonging to the field of computer-aided medicine. Firstly, the lung tumors are segmented, extracted and optimized; then, the prognostic evaluation model of non-small cell lung cancer is constructed to predict the prognosis and survival of patients by combining the imaging features of tumors with the survival status of patients through machine learning method, and the performance of the prognostic analysis model is evaluated by designing experiments. The invention provides a new scheme for exploring the relationship between imaging characteristics and prognosis and survival of patients.

【技术实现步骤摘要】
一种基于CT影像学特征与非小细胞肺癌患者预后情况的相关性分析方法
本专利技术涉及计算机辅助医学领域，涉及一种基于CT影像学特征与非小细胞肺癌患者预后情况的相关性分析方法。
技术介绍
世界卫生组织(WHO)国际癌症研究机构(IARC)近日发布最新报告称，肺癌是全球范围内发病率和死亡率增长最快的恶性肿瘤，预计2018年将造成180万人死亡，占预计癌症死亡总人数的18.4％。根据组织学类型分类，肺癌分为非小细胞肺癌与小细胞肺癌，其中，非小细胞肺癌(non-smallcelllungcancer，NSCLC)占肺癌患病总人数的80％～85％，包括鳞状细胞癌(鳞癌)、腺癌、大细胞癌。非小细胞肺癌相较于小细胞肺癌，生长分裂较慢，扩散转移较晚，致死性也相对较弱，但由于不同非小细胞肺癌患者的病灶之间存在大量个体间差异，使得不同患者患癌预后病情的发展速度存在很大的差别。流行病学统计表明，大量的非小细胞肺癌患者由于未得到准确的病情发展预测而未能及时接受合适的治疗，使得该类肺癌患者的死亡率高达75％。因此，迫切需要有效的患者生存时间预测模型来对治疗和复查方案的选择进行辅助，以提高非小细胞肺癌的治疗效果，进而提高患者的治愈率和存活率。放射组学是医学中的一个新兴领域，该技术的诞生与放射基因组学在疾病研究中的良好性能以及医学影像在疾病的诊断和治疗中的巨大潜力密不可分。放射组学通过特征提取，将肿瘤区域的信息映射到高维特征空间，接着通过机器学习的方法构建疾病的预后模型来对疾病的未来发展进行预测，从而对疾病的治疗和复查方案的选择进行指导。CT影像数据拥有采集便易和结果易于比对等特点，作为放射组学数据库中的重要模态之一，在放射组学研究中得到广泛的应用。近年来，国内外针对计算机辅助诊断(computeraideddiagnosis，CAD)技术和精准医疗(PrecisionMedicine)的研究越来越火热。CAD技术和精准医疗均需借助影像学手段，通过提取大量影像学特征对肿瘤进行量化分析，从而达到辅助临床诊断的目的。而利用提取的影像学特征对患者预后生存情况进行相关性分析，构建预后评估模型预测患者的预后生存情况，从而更好的指导医生对患者的治疗和复查方法进行选择。从目前的国内外研究现状来看，非小细胞癌预后分析的研究普遍从临床的角度出发，首先，根据病例样本的直观的临床特征(临床分期、吸烟史、有无脑转移、肿瘤标记物、医学征象等)作为病例样本的量化指标；然后，通过传统的统计学方法对临床特征与预后关系进行单因素生存分析，得到与患者预后相关的临床特征；最后，将单因素分析中与患者预后相关的临床特征代入COX回归模型进行多因素分析，得到非小细胞肺癌的预后因素，帮助医生根据预后因素对非小细胞肺癌患者进行更加精准的预后评估，设计出更好的治疗与复查方案，以延长患者生存时间。而这类方法也存在着局限性，所能利用的临床信息种类较少，且医学征象类的特征只表现出肿瘤区域的部分形态特性，而从影像学角度出发，可以得到数量更多且种类更为丰富的影像学特征，从而反映肿瘤更多的隐含信息,可以有效解决肿瘤异质性难以定量评估的问题。针对目前研究存在不足，本专利技术设计了新的肺癌预后研究方法，对非小细胞肺癌的预后生存状况进行分析探究，得出基于CT影像特征的非小细胞肺癌预后分析模型对患者的预后生存时间进行预测；同时，设计实验对实验方法进行验证，进而完善目前存在不足，获得更好的预后分析效果。
技术实现思路
本专利技术的目的是针对现有研究中的不足，提出一种基于半自动分割提取的非小细胞肺癌CT影像特征与非小细胞肺癌患者预后情况的相关性分析方法。为实现上述目的，本专利技术采用的技术方案包括以下步骤：CT影像的处理、特征数据处理、影像特征与生存情况的关联分析和实验结果的验证分析。每部分的具体过程描述如下：一.CT影像的处理：首先，通过交互式医学影像控制系统RadiAntViewer软件从CT影像中找到肿瘤并框选出其大致区域，得到肿瘤序列图像。然后，运用半自动分割方法分割出肿瘤，对不同类型的肿瘤采取不同的分割方案：(1)孤立型肿瘤：取灰度阈值分割算法和区域生长算法分割结果的交集；(2)粘连肺壁型肿瘤：首先，基于链码法与“滚球法”对肺区边缘进行修补。然后，参照孤立型肿瘤分割方法进一步分割，最后，检查上述分割结果，对过分割或欠分割情况进行人工修补；(3)粘连血管型肿瘤：首先，运用孤立型肿瘤分割方法得到图像。然后，采用“吹球法”与模糊C均值聚类的方法，并基于阈值得到二值图像，最终的分割结果取两者的交集。最后，检查上述分割结果，对过分割或仍有未去除血管的情况进行人工修正。在肿瘤区域分割的基础上，提取了包括灰度特征、形态特征、纹理特征、医学征象在内的258种定量特征来描述肿瘤；二.特征数据处理：首先，由于机器学习所使用的数据正负样本不平衡会影响训练的模型性能，为了解决这一问题，本专利技术采用SMOTE(SyntheticMinorityOversamplingTechnique)算法平衡正负样本；然后，为了探究CT影像学特征与患者预后生存情况之间的相关性，必须尽可能地多尺度多类别的提取肿瘤的定量指标，以避免造成信息的遗漏，导致分析结果不准确。而提取的特征维数太高经常会导致特征匹配过于复杂，消耗系统资源，而且绝大多数同类特征描述的是不同种类肿瘤之间的共性，还有一部分特征之间具有相关性，使用此类特征会造成大量冗余的计算。为得到性能良好的分类模型，需要对表征冗余信息的特征数据进行过滤，同时将具有类特异性的特征进行保留，因此，为得到性能良好的分类模型，需要对表征冗余信息的特征数据进行过滤，同时将具有类特异性的特征进行保留。所以，采用PCA主成分分析(PrincipalComponentAnalysis)的方法对提取到的特征数据进行降维、优化，减少冗余信息所造成的误差,提高分类模型的性能；最后，对特征数据进行归一化处理，将特征的取值区间缩放到[0,1]的范围；三.影像特征与预后情况的关联分析：利用相关性分析算法分析影像特征与预后情况之间的相关性，以P<0.05为指标筛选出统计学上显著关联的影像特征和预后情况；四.实验结果的验证分析：采取多种分类器建立影像学特征与病人预后生存情况的相关性模型，为了尽可能地构建一个性能最优的分类模型，本文通过大量的研究和测试，最终挑选了以下7种泛化能力强、可用于小样本数据集训练的代表性分类器，分别为决策树(Decisiontree，DT)、鉴别分析分类器(DiscriminantAnalysisClassifiers，DAC)、逻辑回归分类器(LogisticRegression，LR)、支持向量机(SupportVectorMachine，SVM)、K近邻分类器(k-NearestNeighbor,KNN)、集合分类器(Ensembleclassifiers，EC)，随机森林分类器(RandomForest，RF)对数据进行分类和测试。以具有临床意义的3年生存期为分界进行二分类，从而预测出病人的预后情况，并以准确率(ACC)、敏感性(SE)、特异性(SP)、ROC曲线下面积(AUC)作为指标，使用十折交叉验证方法来验证预后分析模型的性能。附图说明图1为基于CT影像学特征与非小细胞肺癌患者预后情况本文档来自技高网...

【技术保护点】
1.一种基于CT影像学特征与非小细胞肺癌预后生存情况的相关性分析方法，其特征在于，包括以下步骤：CT影像的处理、特征数据处理、影像学特征与生存情况的关联分析和实验结果的验证分析。具体过程如下：1.CT影像的处理：首先，通过交互式医学影像控制系统RadiAntViewer软件从CT影像中找到肿瘤并框选出其大致区域，得到肿瘤序列图像。然后，运用半自动分割方法分割出肿瘤，对不同类型的肿瘤采取不同的分割方案：(1)孤立型肿瘤：取灰度阈值分割算法和区域生长算法分割结果的交集；(2)粘连肺壁型肿瘤：首先，基于链码法与“滚球法”对肺区边缘进行修补。然后，参照孤立型肿瘤分割方法进一步分割，最后，检查上述分割结果，对过分割或欠分割情况进行人工修补；(3)粘连血管型肿瘤：首先，运用孤立型肿瘤分割方法得到图像。然后，采用“吹球法”与模糊C均值聚类的方法，并基于阈值得到二值图像，最终的分割结果取两者的交集。最后，检查上述分割结果，对过分割或仍有未去除血管的情况进行人工修正。在肿瘤区域分割的基础上，提取了包括灰度特征、形态特征、纹理特征、医学征象在内的258个定量特征来描述肿瘤；2.特征数据处理：首先，平衡特征数据集正负样本数量。然后，对特征数据进行降维，优化；最后，归一化处理，将特征的取值区间缩放到[0,1]的范围；3.影像特征与预后情况的关联分析：利用Pearson相关性分析算法分析影像学特征与预后情况之间的相关性，以P<0.05为指标筛选出统计学上与预后生存分析模型显著关联的影像学特征；4.实验结果的验证分析：采取多种类分类器建立影像学特征与病人生存情况的相关性模型，选择性能更优的分类模型，以具有临床意义的3年生存期为界限进行分类，从而预测出病人的预后情况，并以准确率(ACC)、敏感性(SE)、特异性(SP)、ROC曲线下面积(AUC)作为指标，使用十折交叉验证方法来验证预后分析模型的性能。...

【技术特征摘要】
1.一种基于CT影像学特征与非小细胞肺癌预后生存情况的相关性分析方法，其特征在于，包括以下步骤：CT影像的处理、特征数据处理、影像学特征与生存情况的关联分析和实验结果的验证分析。具体过程如下：1.CT影像的处理：首先，通过交互式医学影像控制系统RadiAntViewer软件从CT影像中找到肿瘤并框选出其大致区域，得到肿瘤序列图像。然后，运用半自动分割方法分割出肿瘤，对不同类型的肿瘤采取不同的分割方案：(1)孤立型肿瘤：取灰度阈值分割算法和区域生长算法分割结果的交集；(2)粘连肺壁型肿瘤：首先，基于链码法与“滚球法”对肺区边缘进行修补。然后，参照孤立型肿瘤分割方法进一步分割，最后，检查上述分割结果，对过分割或欠分割情况进行人工修补；(3)粘连血管型肿瘤：首先，运用孤立型肿瘤分割方法得到图像。然后，采用“吹球法”与模糊C均值聚类的方法，并基于阈值得到二值图像，最终的分割结果取两者的交集。最后，检查上述分割结果，对过分割或仍有未去除血管的情况进行人工修正。在肿瘤区域分割的基础上，提取了包括灰度特征、形态特征、纹理特征、医学征象在内的258个定量特征来描述肿瘤；2.特征数据处理：首先，平衡特征数据集正负样本数量。然后，对特征数据进行降维，优化；最后，归一化处理，将特征的取值区间缩放到[0,1]的范围；3.影像特征与预后情况的关联分析：利用Pearson相关性分析算法分析影像学特征与预后情况之间的相关性，以P<0.05为指标筛选出统计学上与预后生存分析模型显著关联的影像学特征；4.实验结果的验证分析：采取多种类分类器建立影像学特征与病人生存情况的相关性模型，选择性能更优的分类模型，以具有临床意义的3年生存期为界限进行分类，从而预测出病人的预后情况，并以准确率(ACC)、敏感性(SE)、特异性(SP)、ROC曲线下面积(AUC)作为指标，使用十折交叉验证方法来验证预后分析模型的性能。2.根据权利要求1所述的基于CT影像学特征与非小细胞肺癌基因表达间的相关性分析方法，其特征在于，所述的灰度特征共14个，包括均值、标准差、极大值、极小值、最大值...

【专利技术属性】
技术研发人员：郑军，聂生东，郭小辉，王旭，顾海洋，徐江松，
申请(专利权)人：杭州英库医疗科技有限公司，上海理工大学，
类型：发明
国别省市：浙江,33