一种基于参数与特征迁移的胃肿瘤CT图像分类方法技术

本发明专利技术属于CT图像分类方法技术领域，且公开了一种基于参数与特征迁移的胃肿瘤CT图像分类方法，包括如下步骤：S1：基于Densenet121网络进行卷积特征提取；S2：基于MMD对提取的特征进行特征分析；S3：通过正则化极限学习机分类算法对分析后的特征进行分类，生成诊断模型。本发明专利技术通过将迁移学习思想和ELM算法相结合，提出基于最大均值差异的迁移极限学习机的深度特征分类算法，为减少源域数据和目标域数据特征的分布差异，获取更有效的特征，本发明专利技术采用最大均值差异方法衡量动脉期和静脉期深度学习特征的概率分布，剔除冗余信息，融合两期图像的有效信息，得到泛化能力强、鲁棒性高的胃CT影像识别模型。的胃CT影像识别模型。的胃CT影像识别模型。

【技术实现步骤摘要】
一种基于参数与特征迁移的胃肿瘤CT图像分类方法


[0001]本专利技术属于CT图像分类方法
，具体涉及一种基于参数与特征迁移的胃肿瘤CT图像分类方法。

技术介绍

[0002]胃癌与原发性胃淋巴瘤为胃最常见的两种恶性肿瘤，在影像学上表现极其相似，但两种恶性肿瘤的治疗方法及预后迥异，因此，术前准确诊断胃癌或胃淋巴瘤至关重要，是制定正确治疗计划的前提，计算机断层扫描即CT，可以根据人体不同组织对X射线的吸收与透过率的不同进行组织图像重建，是胃肿瘤的主要检查手段，由于原发性胃淋巴瘤病灶常位于粘膜下层，而医生主要依据CT主观征象进行胃肿瘤的初步诊断，其诊断存在主观差异，敏感度及特异度相对不高，具有一定的局限性，迫切需要引入新方法对胃肿瘤进行辅助化评估，近年来，随着医学影像技术和计算机技术的发展，计算机辅助诊断备受关注，胃癌的辅助诊断主要分为传统的放射组学和深度学习两大类，传统放射组学方面，是利用逻辑回归方法构建放射学标签预测胃癌的复发概率，在内外测试集中诊断性能分别为0.840和0.730，另外一种是基于荧光高光谱成像技术结合支持向量机对早期胃癌患者进行诊断，还有一种是利用极限学习机进行胃癌和卵巢癌的诊断，在胃癌和卵巢癌影像上显示出较好的聚类性能，但需指出，传统方法采用人工定义的较低阶的放射学特征构建模型，需要医生精确勾画病灶，工作量大且存在主观性，导致构建的模型诊断精度和泛化性较差，相比于上述方法，深度学习因其不需要医生精确勾画病灶，通过对打框的病灶图像进行卷积操作来提取特征，能自动学习到与任务更相关的高阶特征的优势，不断有学者展开研究，基于深度学习对胃癌淋巴结转移问题展开研究，开发胃癌术前CT图像分类系统，显示出较好的预测性能(准确率为0.876)，还有一种是利用VGG16卷积神经网络对胃癌的浸润程度进行预测，在测试集上准确率达到0.851，基于常见的深度学习框架进行胃癌诊断，研究表明基于Densenet121网络框架的分类效果最佳，然而，深度学习方法的应用基于两个前提：第一个前提为训练数据和测试数据服从同分布；第二个前提为目标诊断数据量充足，但在临床实践中，受限于疾病发生率，仪器扫描参数不统一等问题，导致医学影像有效数据集小，训练数据和测试数据分布存在差异。因此，有必要寻找一种有效的数据整合算法，解决医学图像数据量不足以及特征差异问题，提高小数据集下模型的性能。

技术实现思路

[0003]本专利技术的目的在于提供一种基于参数与特征迁移的胃肿瘤CT图像分类方法，以解决上述
技术介绍
中提出的问题。
[0004]为实现上述目的，本专利技术提供如下技术方案：一种基于参数与特征迁移的胃肿瘤CT图像分类方法，包括如下步骤：
[0005]S1：基于Densenet121网络进行卷积特征提取；
[0006]S2：基于MMD对提取的特征进行特征分析；
[0007]S3：通过正则化极限学习机分类算法对分析后的特征进行分类，生成诊断模型。
[0008]优选地，所述基于Densenet121网络进行卷积特征提取步骤具体包括：
[0009]S11：将ImageNet预训练的密集连接卷积神经网络作为公共特征提取器，基于公共特征提取器分别对胃部CT动静脉期病灶图像进行特征提取。
[0010]优选地，所述将ImageNet预训练的密集连接卷积神经网络作为公共特征提取器，基于公共特征提取器分别对胃部CT动静脉期病灶图像进行特征提取步骤之后包括：
[0011]S12：假设前l
‑
1层的输出特征为[F0,F1,...,F
l
‑1]，则第l层的特征输出F
l
为：
[0012]F
l
＝Z
l
([F0,F1,...,F
l
‑1])；其中，为Z
l
(*)复合函数，包括3
×
3和1
×
1的卷积、池化、批量归一化以及非线性激活函数四种运算；
[0013]S13：通过KDE对源域和目标域特征数据的概率密度函数差异进行估计。
[0014]优选地，所述通过正则化极限学习机分类算法对分析后的特征进行分类，生成诊断模型的步骤具体包括：
[0015]S31：假设动脉期影像是满足P分布的源领域静脉期影像是满足Q分布的目标领域令H为RKHS，那么X
(s)
和X
(t)
在RKHS中的MMD为：
[0016][0017]其中sup表示求上界，E
P
和E
Q
表示期望，f(*)表示映射函数，H表示RKHS，在此空间内，概率分布P的平均嵌入表示为μ(P)，则μ(P)是空间H中的唯一元素，则有：
[0018]E
P
[f(X
(s)
)]＝<f(X
(s)
,μ(P)>H；
[0019]E
Q
[f(X
(t)
)]＝<f(X
(t)
,μ(Q)>H；
[0020]S32：令d(P,Q)表示数据分布P和Q之间的MMD，则此距离的平方等价于两个分布的平均嵌入距离，即：
[0021][0022]由于μ(P)和μ(Q)无法直接获取，这里采用均值作为其无偏估计的结果，则公式转换为：
[0023][0024]上式中，n表示源域样本数，m表示目标域数据的样本数，f(*)：X
→
H表示原始特征空间映射到RKHS的映射函数，满足＜f(x),f(y)＞＝k(x,y)，该文选取k(x,y)为径向基核函
数，即：
[0025]k(x,y)＝exp(
‑
||x
‑
y||2/δ2)；
[0026]其中δ表示高斯核大小，使用多个径向基核函数加权平均得到最终结果；
[0027]S33：假设d
MMD
＝[d2(X
(s)
,X
(t)
)1,d2(X
(s)
,X
(t)
)2,...,d2(X
(s)
,X
(t)
)
M
]表示动静脉期影像特征MMD值的集合，其中M表示特征总维数，对应M＝11264，然后选取所有的特征的MMD的均值为阈值σ，即：
[0028][0029]其中i表示第i维特征；
[0030]S34：基于上述阈值，MMD小于σ的特征为相似度更高的更有效的特征，分别选出对应的动脉期更有效的特征和静脉期更有效的特征。
[0031]优选地，所述基于上述阈值，MMD小于σ的特征为相似度更高的更有效的特征，分别选出对应的动脉期更有效的特征和静脉期更有效的特征步骤之后包括：
[0032]S35：对于N个病例U＝(x
i
,t
i
)∈[R
m<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于参数与特征迁移的胃肿瘤CT图像分类方法，其特征在于：包括如下步骤：S1：基于Densenet121网络进行卷积特征提取；S2：基于MMD对提取的特征进行特征分析；S3：通过正则化极限学习机分类算法对分析后的特征进行分类，生成诊断模型。2.根据权利要求1所述的一种基于参数与特征迁移的胃肿瘤CT图像分类方法，其特征在于：所述的基于Densenet121网络进行卷积特征提取步骤具体包括：S11：将ImageNet预训练的密集连接卷积神经网络作为公共特征提取器，基于公共特征提取器分别对胃部CT动静脉期病灶图像进行特征提取。3.根据权利要求2所述的一种基于参数与特征迁移的胃肿瘤CT图像分类方法，其特征在于：所述的将ImageNet预训练的密集连接卷积神经网络作为公共特征提取器，基于公共特征提取器分别对胃部CT动静脉期病灶图像进行特征提取步骤之后包括：S12：假设前l
‑
1层的输出特征为[F0,F1,...,F
l
‑1]，则第l层的特征输出F
l
为：F
l
＝Z
l
([F0,F1,...,F
l
‑1])；其中，为Z
l
(*)复合函数，包括3
×
3和1
×
1的卷积、池化、批量归一化以及非线性激活函数四种运算；S13：通过KDE对源域和目标域特征数据的概率密度函数差异进行估计。4.根据权利要求1所述的一种基于参数与特征迁移的胃肿瘤CT图像分类方法，其特征在于：所述的通过正则化极限学习机分类算法对分析后的特征进行分类，生成诊断模型的步骤具体包括：S31：假设动脉期影像是满足P分布的源领域静脉期影像是满足Q分布的目标领域令H为RKHS，那么X
(s)
和X
(t)
在RKHS中的MMD为：其中sup表示求上界，E
P
和E
Q
表示期望，f(*)表示映射函数，H表示RKHS，在此空间内，概率分布P的平均嵌入表示为μ(P)，则μ(P)是空间H中的唯一元素，则有：E
P
[f(X
(s)
)]＝&lt;f(X
(s)
,μ(P)&gt;
H
；E
Q
[f(X
(t)
)]＝&lt;f(X
(t)
,μ(Q)&gt;
H
；S32：令d(P,Q)表示数据分布P和Q之间的MMD，则此距离的平方等价于两个分布的平均嵌入距离，即：由于μ(P)和μ(Q)无法直接获取，这里采用均值作为其无偏估计的结果，则公式转换为：
上式中，n表示源域样本数，m表示目标域数据的样本数，f(*)：X
→
H表示原始特征空间映射到RKHS的映射函数，满足＜f(x),f(y)＞＝k(x,y)，选取k(x,y)为径向基核函数，即：k(x,y)＝exp(
‑
||x
‑
y||2/δ2)；其中δ表示高斯核大小，使用多个径向基核函数加权平均得到最终结果；S33：假设d
MMD
＝[d2(X
(s)
,X
(t)
)1,d2(X...

【专利技术属性】
技术研发人员：刘昱，匡梦，韦松林，
申请(专利权)人：桂林航天工业学院，
类型：发明
国别省市：