尘肺病等级智能判别方法技术

本发明专利技术涉及医学影像分析技术领域，公开了尘肺病等级智能判别方法，其包括如下步骤：将包含有多张DR胸片样本图像分为训练样本一和训练样本二。所述DR胸片样本图像为具有不同大小尘肺纤维阴影像素的各个肺区图像。基于卷积神经网络模型输入训练样本一，训练肺部纤维化阴影的特征的模型参数，以生成肺部纤维化智能判定模型。本发明专利技术的尘肺病等级智能判别方法中，输入的图像由训练用的图像块替换为整张DR胸片，特征学习模型设计为将最后三层全连接层替换为卷积层，这样设计的目的是为了实现输入图像和输出也是图像的端到端检测，同时可以直接检测出不同肺区的病灶区域，减少训练时间，同时提高了对肺纤维化的识别率。同时提高了对肺纤维化的识别率。同时提高了对肺纤维化的识别率。

【技术实现步骤摘要】
尘肺病等级智能判别方法


[0001]本专利技术涉及医学影像分析
，尤其涉及尘肺病等级智能判别方法。

技术介绍

[0002]肺纤维化是各种肺部疾病的共同结局，以肺组织疤痕化为主要表现，如果受累范围广泛则导致肺容积缩小，肺功能明显下降，严重影响患者生存质量。尤其是特发性间质性肺炎是最典型的代表，其病理和/或影像表现为寻常性间质性肺炎的慢性进展性肺部疾病。IPF病因不明，预后极差，确诊后平均生存期仅为3
‑
5年。IPF目前认为患者个体生存期差异很大，有些长期稳定存活多年，有些病情缓慢进展，有些则发生急性加重导致肺功能急速恶化，短时间内死亡。如何根据患者的情况做出较为精准的病情严重度评估和预后判断，目前尚缺乏被广泛接受的评估方法。
[0003]当前主要采用的人工评估方法，是通过选择DR肺窗的4个代表层面：主动脉弓层面、气管分叉层面、下叶上基底段气管分叉层面和右侧肺隔顶下层面；每个层面有左右两个肺野，全肺共分为8个肺野，将8个肺野总共划分为100个小片；通过对每个小片上蜂窝状改变的面积是否大于该小片面积的二分之一为判断阳性还是因阴性的标准，将阳性的片相加后得到全肺蜂窝状改变的占比。
[0004]当前临床识别肺纤维化病灶及评估其严重程度主要依赖于临床医生对肺部dr图片进行肉眼判断并粗略评估，评估准确性完全依赖于医生个人的经验。而且肺部DR影像数据呈现爆炸式增长，极大地增强了医生的工作量，随之在疾病检测过程中容易造成误诊和漏诊。
[0005]随着大数据计算机视觉的发展，计算机辅助诊断技术也被用帮助完成基于医学影像的肺纤维化疾病诊断，以减轻医生的工作负担。但现有的方法对肺纤维化识别率较低，且尘肺病等级判定不准确。

技术实现思路

[0006]为解决上述
技术介绍
中提到的技术问题，本专利技术提供尘肺病等级智能判别方法。
[0007]本专利技术采用以下技术方案实现：尘肺病等级智能判别方法，其包括如下步骤：
[0008]S1、将包含有多张DR胸片样本图像分为训练样本一和训练样本二；所述DR胸片样本图像为具有不同大小尘肺纤维阴影像素的各个肺区图像；
[0009]S2、基于卷积神经网络模型输入训练样本一，训练肺部纤维化阴影的特征的模型参数，以生成肺部纤维化智能判定模型；
[0010]S3、基于肺部纤维化智能判定模型输入训练样本二，识别出训练样本二内的每张DR胸片图像中尘肺纤维阴影像素的阴影形状和相应坐标；
[0011]S4、分析预测每张DR胸片图像中尘肺纤维阴影形状，确定该每个阴影形状所属的阴影密度集，计算出该阴影形状的总体密度集，并结合该尘肺纤维阴影像素分布所在坐标的肺区数量，进行尘肺等级的判定训练，以训练出尘肺等级判定模型；
[0012]S5、通过尘肺病判定模型的各项参数判定尘肺等级。
[0013]作为上述方案的进一步改进，所述的步骤S1中，对于每张DR胸片样本图像，需要预先标定出每张DR胸片样本图像中各个肺区所包含的各个种类的尘肺纤维阴影像素，并将标定后的DR胸片样本图像作为样本标签。
[0014]作为上述方案的进一步改进，所述的步骤S2中，对于每一张DR胸片样本图像依序经过卷积神经网络模型的其中的几个卷积层卷积后，分别生成不同映射尺度的目标特征向量；
[0015]根据获得的各个所述目标特征向量和训练样本一中预先标定的样本标签，通过复合训练损失函数训练所述卷积神经网络模型的模型参数。
[0016]作为上述方案的更进一步改进，所述卷积神经网络模型包括依序相连的8个卷积层，各个卷积层均使用1
×
1的卷积核进行相应的特征提取；
[0017]当DR胸片样本图像经过卷积层conv4、conv5、conv6、conv8后，会分别生成多个用于与样本标签比较的目标特征向量。
[0018]作为上述方案的更进一步改进，所述卷积层conv1
‑
3的图像像素均为512
×
515，卷积层conv4、卷积层conv5、卷积层conv6、卷积层conv7、卷积层conv8的图像像素依次分别为256
×
256、128
×
128、64
×
64、32
×
32、16
×
16。
[0019]作为上述方案的更进一步改进，通过复合训练损失函数训练肺部纤维化阴影的特征的模型参数的方法包括如下步骤：
[0020]将目标分类误差函数定义为：
[0021]其中，W是卷积神经网络模型的模型参数，S＝{S1,S2,...,S
N
}是训练样本集合，N是卷积层的数量，α
n
是误差函数l
n
的权重，X
i
是标注好的样本图像块，Y
i
＝(y
i
,b
i
)是标注好的样本标签，y
i
∈{0,1,...,K}，真实目标框坐标，(x，y，w，h)分别表示目标框的中心点坐标和宽长值；
[0022]将坐标误差函数定义为：
[0023]其中，m是检测时的预测目标框所在卷积层的层数，i是检测的初始目标框坐标，j是真实目标框坐标，坐标误差函数采用smooth
L1
()求解预测目标框和真实目标框的坐标误差，则的计算方法为：
[0024][0025]其中，表示预测目标框的中心、长宽坐标，表示真实目标框的中心、长宽坐标，表示预测目标框和真实目标框的坐标差值；
[0026]将目标分类误差函数与坐标误差函数联合训练学习，融合每一个卷积层输出的目标特征向量，则误差函数表示为：
[0027][0028]其中，p(X)＝(p0(X),...,p
K
(X))是目标分类概率值，λ是平衡参数且通过交叉验证设置为1，分类误差函数采用交叉熵函数计算，L
obj
(p(X),y)＝
‑
logp
y
(X)。
[0029]作为上述方案的更进一步改进，所述初始目标框的宽和长计算公式分别表示为：
[0030]则每个初始目标框的第i行和第j列的中心点设置为m是m
×
m特征映射图的尺寸。
[0031]作为上述方案的更进一步改进，所述利用平铺初始检测框方式训练学习特殊尺度的不同目标，假设使用m个卷积层来预测，则每层的初始目标检测框的尺寸比例表示为：
[0032]其中，s
min
就是最低层的尺度系数为，s
max
最高层的尺度系数为，k为初始目标框的数量。
[0033]作为上述方案的进一步改进，所述的步骤S4中，计算出阴影形状的总体密度集的像素和该总体密度集的像素在其所在肺区像素的比值k，通过比值k来表示每个肺区小阴影密度等级的判定阀值表示为：
[0034][0035]其中，表示第i个肺区小阴影密度集属于j级目标检测框的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.尘肺病等级智能判别方法，其特征在于，其包括如下步骤：S1、将包含有多张DR胸片样本图像分为训练样本一和训练样本二；所述DR胸片样本图像为具有不同大小尘肺纤维阴影像素的各个肺区图像；S2、基于卷积神经网络模型输入训练样本一，训练肺部纤维化阴影的特征的模型参数，以生成肺部纤维化智能判定模型；S3、基于肺部纤维化智能判定模型输入训练样本二，识别出训练样本二内的每张DR胸片图像中尘肺纤维阴影像素的阴影形状和相应坐标；S4、分析预测每张DR胸片图像中尘肺纤维阴影形状，确定该每个阴影形状所属的阴影密度集，计算出该阴影形状的总体密度集，并结合该尘肺纤维阴影像素分布所在坐标的肺区数量，进行尘肺等级的判定训练，以训练出尘肺等级判定模型；S5、通过尘肺病判定模型的各项参数判定尘肺等级。2.如权利要求1所述的尘肺病等级智能判别方法，其特征在于，所述的步骤S1中，对于每张DR胸片样本图像，需要预先标定出每张DR胸片样本图像中各个肺区所包含的各个种类的尘肺纤维阴影像素，并将标定后的DR胸片样本图像作为样本标签。3.如权利要求2所述的尘肺病等级智能判别方法，其特征在于，所述的步骤S2中，对于每一张DR胸片样本图像依序经过卷积神经网络模型的其中的几个卷积层卷积后，分别生成不同映射尺度的目标特征向量；根据获得的各个所述目标特征向量和训练样本一中预先标定的样本标签，通过复合训练损失函数训练所述卷积神经网络模型的模型参数。4.如权利要求3所述的尘肺病等级智能判别方法，其特征在于，所述卷积神经网络模型包括依序相连的8个卷积层，各个卷积层均使用1
×
1的卷积核进行相应的特征提取；当DR胸片样本图像经过卷积层conv4、conv5、conv6、conv8后，会分别生成多个用于与样本标签比较的目标特征向量。5.如权利要求3所述的尘肺病等级智能判别方法，其特征在于，所述卷积层conv1
‑
3的图像像素均为512
×
515，卷积层conv4、卷积层conv5、卷积层conv6、卷积层conv7、卷积层conv8的图像像素依次分别为256
×
256、128
×
128、64
×
64、32
×
32、16
×
16。6.如权利要求3所述的尘肺病等级智能判别方法，其特征在于，通过复合训练损失函数训练肺部纤维化阴影的特征的模型参数的方法包括如下步骤：将目标分类误差函数定义为：其中，W是卷积神经网络模型的模型参数，S＝{S1,S2,...,S
N
}是训练样本集合，N是卷积层的数量，α
n
是误差函数l
n
的权重，X
i
是标注好的样本图像块，Y
i
＝(y
i
,b
i
)是标注好的样本标签，y
i
∈{0,1,...,K}，真实目标框坐标...

【专利技术属性】
技术研发人员：张源，江震，刘静怡，
申请(专利权)人：安徽医学高等专科学校，
类型：发明
国别省市：