CN113205153A - 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型 - Google Patents

一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型 Download PDF

Info

Publication number
CN113205153A
CN113205153A CN202110579827.2A CN202110579827A CN113205153A CN 113205153 A CN113205153 A CN 113205153A CN 202110579827 A CN202110579827 A CN 202110579827A CN 113205153 A CN113205153 A CN 113205153A
Authority
CN
China
Prior art keywords
feature
lbp
training
feature map
processing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110579827.2A
Other languages
English (en)
Other versions
CN113205153B (zh
Inventor
郑力新
阳潮鹏
严潭
苏秋玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huaqiao University
Original Assignee
Huaqiao University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huaqiao University filed Critical Huaqiao University
Priority to CN202110579827.2A priority Critical patent/CN113205153B/zh
Publication of CN113205153A publication Critical patent/CN113205153A/zh
Application granted granted Critical
Publication of CN113205153B publication Critical patent/CN113205153B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G06F17/153—Multidimensional correlation or convolution
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/25—Fusion techniques
    • G06F18/253—Fusion techniques of extracted features
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H30/00—ICT specially adapted for the handling or processing of medical images
    • G16H30/20—ICT specially adapted for the handling or processing of medical images for handling medical images, e.g. DICOM, HL7 or PACS
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Computation (AREA)
  • Evolutionary Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Epidemiology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Biomedical Technology (AREA)
  • Primary Health Care (AREA)
  • Databases & Information Systems (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • General Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Software Systems (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Image Analysis (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明提供一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型，训练方法是：获取儿科肺炎患者的医学图像及对应的医学诊断语句，医学图像处理为统一规格m×m像素大小后形成训练集；对所述训练集分别提取LBP特征图集合、多尺度特征图集合、以及多级别特征图集合；分别采用注意力模块进行注意力处理后进行特征融合，获得新的特征图集合；将所获得的特征图集合输入到分类器中进行分类，获得训练好的模型。本发明通过对现有儿科肺炎患者的医学图像进行训练，让训练所得的模型作为医生学习诊断的工具或为医生的临床诊断提供有效的参考意见。

Description

一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型
技术领域
本发明涉及医学计算机技术领域，特别涉及一种儿科肺炎辅助诊断模型及其训练方法。
背景技术
肺炎是一种影响肺部功能的急性呼吸道感染疾病，肺炎使得肺炎泡发炎导致呼吸道产生脓液，让人呼吸困难。根据世界卫生组织的报告，2015年肺炎估计造成世界约92万名5岁以下儿童死亡，占造成5岁以下儿童死亡人数的16％，直至2017年全球仍有80余万5岁以下儿童死于肺炎。目前诊断的方法主要是依靠胸部X光检查来诊断肺炎，肺炎在胸部X光图像中表现为不透明增加的区域。
部分偏远地区的医护人员由于经验不足等情况使得鉴别胸部X光图像是否为患有肺炎的图像成为一个具有挑战的性的工作，同时诊断过程中需要大量的时间和精力，而且还十分容易受到观察人员主观因素的影响导致误诊。
通过使用基于深度学习的儿科肺炎辅助诊断算法能够很好的辅助医务人员对肺炎的判别，从而减轻医生的负担、提高医院检测工作效率。近年来随着深度学习的研究不断深入，可以说深度学习模型是图像分类问题的一个强有力的工具。依靠深度学习来解决如图片分类等实际问题也越来越可靠。因此若能实现一种儿科肺炎辅助诊断模型，一方面可以为医生的诊断提供帮助，另一方面还可以让该模型作为医生学习诊断的工具，无疑会为医生的快速成长提供巨大的助力。
发明内容
本发明要解决的技术问题，在于提供一种儿科肺炎辅助诊断模型及其训练方法，通过对现有儿科肺炎患者的医学图像及对应的医学诊断语句进行训练，让训练所得的模型作为医生学习诊断的工具，帮助医生快速成长。
第一方面，本发明提供了一种儿科肺炎辅助诊断模型的训练方法，包括：
S1、获取儿科肺炎患者的医学图像及对应的医学诊断语句，将获得的医学图像处理为统一规格m×m像素大小的医疗图像集合，m由网络模型大小决定；根据医学诊断语句将医学图像分为有病和没病两类，并混合后形成训练集；
S2、对所述训练集分别提取LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw；
S3、将所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw分别采用注意力模块进行注意力处理；
S4、将注意力处理后的所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw进行特征融合，获得新的特征图集合Msum；
S5、将所获得的特征图集合Msum输入到分类器中进行分类，获得训练好的模型。
第二方面，本发明提供了一种儿科肺炎辅助诊断模型，其经由第一方面的训练方法训练获得。
本发明实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：在儿科肺炎患者的医学图像上提取出LBP特征后，得到特征图集合MLBP，还提取了多尺度特征图集合Mm和多级别特征图集合Mw，并与降维后的LBP特征相融合，作为新的特征再输入到分类器中进行训练，从而获得儿科肺炎辅助诊断模型。由于本发明捕捉了更多关于肺炎特征信息，从而提高了分类的准确度。另外，多尺度特征图集合Mm是通过采用不同的扩张因子进行空洞卷积获得多尺度的特征集合，从而达到通过多尺度窗口获取特征；同时多级别特征图集合是逐级提取，包括浅层特征信息和深层特征信息，从而获得的特征信息更多，并进行融合以达到获取全局上下文特征信息，进而能够从全局的角度推断肺炎目标与背景之间的关系，能够突出目标区域，淡化背景信息的干扰。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
下面参照附图结合实施例对本发明作进一步的说明。
图1为本发明用于训练儿科肺炎辅助诊断模型的方法的流程图；
具体实施方式
本申请实施例通过提供一种儿科肺炎辅助诊断模型及其训练方法，通过对现有儿科肺炎患者的医学图像及对应的医学诊断语句进行训练，让训练所得的模型作为医生学习诊断的工具，无疑会为医生的快速成长提供巨大的助力。
本申请实施例中的技术方案，总体思路如下：本发明除了在儿科肺炎患者的医学图像上提取出LBP特征，还提取了多尺度特征图集合Mm和多级别特征图集合Mw，并将三种特征进行融合，作为新的特征再输入到分类器中进行训练，从而获得儿科肺炎辅助诊断模型。由于捕捉了更多关于肺炎特征信息，从而提高了分类的准确度。
实施例一
如图1所示，本实施例提供一种儿科肺炎辅助诊断模型的训练方法，包括：
S1、获取5岁以下儿科肺炎患者的医学图像及对应的医学诊断语句，将获得的医学图像处理为m×m的像素大小，获得统一规格医疗图像集合，m由用于提取特征的网络模型的大小决定(如ResNet模型输入的图片大小是224x224，即m＝224)；根据医学诊断语句将医学图像分为有病和没病两类，再将所述医学图像集合按照64％、16％和20％的比例划分为训练集、验证集和测试集，其中训练集、验证集和和测试集均包含一定比例的有病和没病两类医学图像；比如将医学图像集合按照是否患有肺炎分类为0、1两个集合，0表示患有肺炎，1则表示健康，0集合中的图片命名都以0开头，1集合中的图片命名都以1开头。
训练集用来训练模型；通过验证集的反馈可以及时修改模型中的各种参数，避免过拟合的问题；测试集用来检测模型的效果。
S2、对所述训练集分别提取LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw；
S3、将所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw分别采用注意力模块进行注意力处理；
S4、将注意力处理后的所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw进行特征融合，获得新的特征图集合Msum；
S5、将所获得的特征图集合Msum输入到分类器中进行分类，获得训练好的模型。
其中，作为本实施例的一种更优或更为具体的实现方式，所述步骤S1中还对所述医学图像进行图像扩充处理，图像扩充处理具体是：
S11、将所述医学图像作为初始图像；
S12、对所述初始图像进行旋转处理，如旋转90°、180°和270°，得到旋转医学图像；在所述初始图像上添加噪声处理，得到噪声医学图像；
其中，所述旋转采用下述公式1计算：
Figure BDA0003085740230000041
公式1中，(i2,j2)是原图像F(i2,j2)的像素点的坐标，θ为旋转角度，(i1,j1)是将对应像素点(i2,j2)旋转之后的像素点的坐标；
所述添加噪声处理采用下述公式2计算：
(i1,j1)＝(i2,j2)+XMeans+sigma*G(d) 公式2；
公式2中，(i2,j2)是原图像F(i2,j2)的像素点的坐标，(i1,j1)是将对应像素点(i2,j2)加入高斯噪声之后的像素点的坐标，XMeans代表高斯分布的平均值，sigma代表高斯分布的标准方差，d为一个线性的随机数，G(d)是随机数的高斯分布随机值；
S13、将所述初始图像、所述旋转医学图像和所述噪声医学图像进行放缩处理，获得统一规格的医疗诊断图像，并作为所述训练图像集合。
所述缩放处理采用下述公式3计算：
Figure BDA0003085740230000051
公式2中，(i2,j2)是原图像某一像素点的坐标，
Figure BDA0003085740230000052
为水平方向缩放的比率，
Figure BDA0003085740230000053
为垂直方向缩放的比率，(i1,j1)是将对应像素点(i2,j2)缩放之后的像素点的坐标。
通过图片旋转、图片加噪声的方法，对肺部X光图片进行扩充。
所述步骤S2中，对所述训练集提取LBP特征图集合MLBP的具体过程是：
(11)将所述训练集中的医学图像划分成像素大小为n×n的子块，m和n为正整数且m能被n整除，如m＝224，n＝16，共可划分为196个子块，则该196个子块为一个检测窗口；
(12)取每一子块作为中心子块，在中心子块中选择一个像素点作为窗口中心像素点(xc,yc)，在邻域半径为R的圆形区域内选取p个相邻子块的像素点作为采样点(通常p＝8，包括中心子块的上、下、左、右、左上斜、左下斜、右上斜、右下斜8个相邻子块)，所述窗口中心像素点采用
Figure BDA0003085740230000054
的算子分别与p个采样点进行灰度值比较，若所述采样点的灰度值更大，则该窗口中心像素点的位置标记为1，否则标记为0，从而产生p位二进制数，得到该中心子块的窗口中心像素点(xc,yc)的LBP特征值；如公式4所示：
Figure BDA0003085740230000055
公式4中，ic是中心子块的窗口中心像素点(xc,yc)的灰度值，ip是采样点的灰度值，s是一个符号函数，如公式5所示：
Figure BDA0003085740230000056
对于所述窗口中心像素点(xc,yc)，它的采样点集合P中的某个采样点p(xp,yp)而言，p∈P用以下公式计算：
Figure BDA0003085740230000061
公式6中，如果计算的结果不在像素坐标上，则采用双线性插值进行近似处理，具体计算公式如下：
Figure BDA0003085740230000062
(13)计算每个子块的直方图，并对所述直方图进行归一化处理，从而将所述LBP特征值与图形的空间信息结合在一起，得到各子块的统计直方图LBPH；
(14)将各子块的统计直方图LBPH连接成为一张肺部X光图片的LBP特征集合a，所述LBP特征集合a中的LBP特征的维度为2×n2×m/n，若m＝224，n＝16，则LBP特征的维度是100352；
(15)通过PCA对所述LBP特征集合a中的LBP特征进行降维处理，得到h维的LBP特征，h为用户设定值，如h＝224。
PCA的做法就是通过原来的变量做线性投影形成新的变量，一般通过公式8来计算特征的主成分：
y＝UT(xi-x) 公式8；
公式8中，y表示主成分特征，x表示训练样本的特征均值，xi为需要降维的特征，UT为协方差矩阵计算公式，如公式9所示：
Figure BDA0003085740230000063
所述步骤S2中，所述多尺度特征图集合Mm的提取具体是：
(21)采用ResNet50模型通过次残差后提取特征图Mm1，
(21)使用不同的扩张因子对所述特征图Mm1进行空洞卷积处理；同时，采用自适应平均池化操作，获取相同分辨率、相同通道数以及不同尺度的上下文特征信息，将所述上下文特征信息进行级联，得到所述多尺度特征图集合Mm；所述空洞卷积处理的输出计算公式表示如下：
Mi＝Cov(Mm1,r),r∈[1,3,5,7]
Mag＝Avgpool(Mm1) 公式10；
Mm＝C(M1,M3,M5,M7,Mag)
公式10中，Mi为不同的扩张因子1，3，5，7卷积后的特征图，Cov(Mm1,r)为对特征图进行扩张率为r的空洞卷积，Mag为平均池化后的特征图，Avgpool(Mm1)为对特征图Mm1进行平均池化运算，C(M1，M3，M5，M7，Meg)为对所有的特征图即M1，M3，M5，M7和Meg进行级联。
所述步骤S2中，所述多级别特征图集合Mw的提取具体为：
(31)获得所述多尺度特征图集合Mm，并进行全局平均池化操作后得到全局上下文特征信息；
(32)校准所述全局上下文特征信息，生成具有全局信息的掩码；浅层特征经过卷积运算与所述掩码相乘，再经过卷积运算得到特征图的输出；
Mg＝Cov(Mlow)×Mask(Mglobal)
M1＝Cov(Mhigh)×Mask(Mlow)
Mh＝Mask(Mhigh)×Cov(Mlow)
Ma＝C(Mg,Ml,Mh) 公式11；
Mb＝C(Cov(Mglobal,Mlow,Mhigh))
Mw＝Add(Cov(Ma,Mb))
公式11中，Mlow为浅层特征信息，Mhigh为深层特征信息，Mglobal为全局上下文特征信息，Mg为全局上下文信息生成的掩码与浅层特征融合的结果，Ml为浅层特征生成的掩码与深层特征融合的结果，Mh为深层特征生成的掩码与浅层特征融合的结果，Cov(·)为卷积运算，Mask(·)为不同级别特征生成对应的掩码，Add(·)为特征图对应元素相加，Ma为Mg、Ml和Mh级联后的特征图,Mb为Mglobal、Mlow和Mhigh经过卷积后联的特征图,Mw为Ma和Mb相加后的特征图。
所述步骤S3中，将所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw分别采用注意力模块处理具体是：
S31、将上述特征图像集合Mm，Mw，MLBP分别进行两次卷案件操作，分别得到：第一次卷积后获得的多尺度特征集合Mmc1、多层次特征集合Mwc1和LBP特征集合MLBPc1，第二次卷积后获得的多尺度特征集合Mmc2、多层次特征集合Mwc2和LBP特征集合MLBPc2；
S32、然后分别将第一次卷积获得的所有特征集合Mc1和第二次卷积获得的所有特征集合Mc2做相乘和相加运算，其中，Mc1＝Mmc1+Mwc1+MLBPc1；Mc2＝Mmc2+Mwc2+MLBPc2；最后获得校准后的特征图集合Moc作为解码阶段的输入，具体计算公式如下：
Moc＝Mc1×Mc2+Mc1+Mc2 公式12。
所述步骤S4中，将所述特征集合Mm，所述特征集合Mw，所述特征集合MLBP进行特征融合处理具体是：
S41、将所述特征集合Mm的掩码和所述特征集合Mw的掩码相乘，获得初步融合特征集合Mmw；
S42、将所述初步融合特征集合Mmw的掩码和所述特征集合MLBP的掩码相乘，获得新的特征图集合Msum。
所述步骤5具体为：
S51、将所述融合特征c后接一个K维的全连接层，利用softmax分类器进行分类，由softmax分类器中的softmax函数将一个K维的任意实数向量压缩成另一个K维的实数向量，其中实数向量中的每个元素取值都介于(0，1)之间，计算公式如下：
Figure BDA0003085740230000081
在本实施例中，由于只需将病例诊断为有病与非有病两类，因此K取2。
最后，将所述测试集中的肺部X光图片输入到训练好的模型得出预测值，并将预测值与真实值进行比较算出模型的准确率，模型的输入为图片的像素点(在计算机中读取的是一个矩阵)，输出的是预测结果。预测值指的是通过模型预测出来的值，真实值指的是实际的值。准确率等于预测值与真实值相等的个数/总个数。
基于同一发明构思，本申请还提供了一种儿科肺炎辅助诊断模型，是由上述实施例一中的训练方法训练得到，详见实施例二。
实施例二
在本实施例中提供了一种儿科肺炎辅助诊断模型，如图1所示，是由上述实施例一中的训练方法训练得到。
由于本实施例二中的儿科肺炎辅助诊断模型即为经实施例一中的儿科肺炎辅助诊断模型训练方法训练得到，因此本实施例中未详尽说明之处，请参见实施例一中的说明。
本申请实施例中提供的技术方案，至少具有如下技术效果或优点：
虽然以上描述了本发明的具体实施方式，但是熟悉本技术领域的技术人员应当理解，我们所描述的具体的实施例只是说明性的，而不是用于对本发明的范围的限定，熟悉本领域的技术人员在依照本发明的精神所作的等效的修饰以及变化，都应当涵盖在本发明的权利要求所保护的范围内。

Claims (9)

1.一种儿科肺炎辅助诊断模型的训练方法，其特征在于：包括：
S1、获取儿科肺炎患者的医学图像及对应的医学诊断语句，将获得的医学图像处理为统一规格m×m像素大小的医疗图像集合，m由网络模型大小决定；根据医学诊断语句将医学图像分为有病和没病两类，并混合后形成训练集；
S2、对所述训练集分别提取LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw；
S3、将所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw分别采用注意力模块进行注意力处理；
S4、将注意力处理后的所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw进行特征融合，获得新的特征图集合Msum；
S5、将所获得的特征图集合Msum输入到分类器中进行分类，获得训练好的模型。
2.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤S1中还对所述医学图像进行图像扩充处理，图像扩充处理具体是：
S11、将所述医学图像作为初始图像；
S12、对所述初始图像进行旋转处理，得到旋转医学图像；在所述初始图像上添加噪声处理，得到噪声医学图像；
其中，所述旋转采用下述公式1计算：
Figure FDA0003085740220000011
公式1中，(i2，j2)是原图像F(i2，j2)的像素点的坐标，θ为旋转角度，(i1，j1)是将对应像素点(i2，j2)旋转之后的像素点的坐标；
所述添加噪声处理采用下述公式2计算：
(i1，j1)＝(i2，j2)+XMeans+sigma*G(d) 公式2
公式2中，(i2,j2)是原图像F(i2,j2)的像素点的坐标，(i1,j1)是将对应像素点(i2,j2)加入高斯噪声之后的像素点的坐标，XMeans代表高斯分布的平均值，sigma代表高斯分布的标准方差，d为一个线性的随机数，G(d)是随机数的高斯分布随机值；
S13、将所述初始图像、所述旋转医学图像和所述噪声医学图像进行放缩处理，获得统一规格的医疗诊断图像，并作为所述训练图像集合。
所述缩放处理采用下述公式3计算：
Figure FDA0003085740220000021
公式3中，(i2,j2)是原图像某一像素点的坐标，
Figure FDA0003085740220000022
为水平方向缩放的比率，
Figure FDA0003085740220000023
为垂直方向缩放的比率，(i1,j1)是将对应像素点(i2,j2)缩放之后的像素点的坐标。
3.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：
所述步骤S2中，对所述训练集提取LBP特征图集合MLBP的具体过程是：
(11)将所述训练集中的医学图像划分成像素大小为n×n的子块，m和n为正整数且m能被n整除；
(12)取每一子块作为中心子块，在中心子块中选择一个像素点作为窗口中心像素点(xc,yc)，在邻域半径为R的圆形区域内选取p个相邻子块的像素点作为采样点，所述窗口中心像素点采用
Figure FDA0003085740220000024
的算子分别与p个采样点进行灰度值比较，若所述采样点的灰度值更大，则该窗口中心像素点的位置标记为1，否则标记为0，从而产生p位二进制数，得到该中心子块的窗口中心像素点(xc,yc)的LBP特征值；
(13)计算每个子块的直方图，并对所述直方图进行归一化处理，从而将所述LBP特征值与图形的空间信息结合在一起，得到各子块的统计直方图LBPH；
(14)将各子块的统计直方图LBPH连接成为一张肺部X光图片的LBP特征集合a，所述LBP特征集合a中的LBP特征的维度为2×n2×m/n；
(15)通过PCA对所述LBP特征集合a中的LBP特征进行降维处理，得到h维的LBP特征，h为用户设定值。
4.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤S2中，所述多尺度特征图集合Mm的提取具体是：
(21)采用ResNet50模型通过次残差后提取特征图Mm1，
(21)使用不同的扩张因子对所述特征图Mm1进行空洞卷积处理；同时，采用自适应平均池化操作，获取相同分辨率、相同通道数以及不同尺度的上下文特征信息，将所述上下文特征信息进行级联，得到所述多尺度特征图集合Mm；
其中，所述空洞卷积处理的输出计算公式表示如下：
Figure FDA0003085740220000031
公式10中，Mi为不同的扩张因子卷积后的特征图，Cov(Mm1，r)为对特征图进行扩张率为r的空洞卷积，Mag为平均池化后的特征图，Avgpool(Mm1)为对特征图Mm1进行平均池化运算，C(M1，M3，M5，M7，Meg)为对所有的特征图即M1，M3，M5，M7和Meg进行级联。
5.根据权利要求4所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤S2中，所述多级别特征图集合Mw的提取具体为：
(31)获得所述多尺度特征图集合Mm，并进行全局平均池化操作后得到全局上下文特征信息；
(32)校准所述全局上下文特征信息，生成具有全局信息的掩码；浅层特征经过卷积运算与所述掩码相乘，再经过卷积运算得到特征图的输出；
Figure FDA0003085740220000041
公式11中，Mlow为浅层特征信息，Mhigh为深层特征信息，Mglobal为全局上下文特征信息，Mg为全局上下文信息生成的掩码与浅层特征融合的结果，M1为浅层特征生成的掩码与深层特征融合的结果，Mh为深层特征生成的掩码与浅层特征融合的结果，Coy(·)为卷积运算，Mask(·)为不同级别特征生成对应的掩码，Add(·)为特征图对应元素相加，Ma为Mg、M1和Mh级联后的特征图，Mb为Mgloba1、Mlow和Mhigh经过卷积后联的特征图，Mw为Ma和Mb相加后的特征图。
6.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤S3中，将所述LBP特征图集合MLBP、多尺度特征图集合Mm、以及多级别特征图集合Mw分别采用注意力模块处理具体是：
S31、将上述特征图像集合Mm，Mw，MLBP分别进行两次卷案件操作，分别得到：第一次卷积后获得的多尺度特征集合Mmc1、多层次特征集合Mwc1和LBP特征集合MLBPc1，第二次卷积后获得的多尺度特征集合Mmc2、多层次特征集合Mwc2和LBP特征集合MLBPc2；
S32、然后分别将第一次卷积获得的所有特征集合Mc1和第二次卷积获得的所有特征集合Mc2做相乘和相加运算，其中，Mc1＝Mmc1+Mwc1+MLBPc1；Mc2＝Mmc2+Mwc2+MLBPc2；最后获得校准后的特征图集合Moc作为解码阶段的输入，具体计算公式如下：
Moc＝Mc1×Mc2+Mc1+Mc2 公式12。
7.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤S4中，将所述特征集合Mm，所述特征集合Mw，所述特征集合MLBP进行特征融合处理具体是：
S41、将所述特征集合Mm的掩码和所述特征集合Mw的掩码相乘，获得初步融合特征集合Mmw；
S42、将所述初步融合特征集合Mmw的掩码和所述特征集合MLBP的掩码相乘，获得新的特征图集合Msum。
8.根据权利要求1所述的一种儿科肺炎辅助诊断模型的训练方法，其特征在于：所述步骤5具体为：
S51、将所述融合特征c后接一个K维的全连接层，利用softmax分类器进行分类，由softmax分类器中的softmax函数将一个K维的任意实数向量压缩成另一个K维的实数向量，其中实数向量中的每个元素取值都介于(0，1)之间，计算公式如下：
Figure FDA0003085740220000051
9.一种儿科肺炎辅助诊断模型，其特征在于：由权利要求1至8任一项所述的一种儿科肺炎辅助诊断模型的训练方法训练获得。
CN202110579827.2A 2021-05-26 2021-05-26 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型 Active CN113205153B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110579827.2A CN113205153B (zh) 2021-05-26 2021-05-26 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110579827.2A CN113205153B (zh) 2021-05-26 2021-05-26 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型

Publications (2)

Publication Number Publication Date
CN113205153A true CN113205153A (zh) 2021-08-03
CN113205153B CN113205153B (zh) 2023-05-30

Family

ID=77023338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110579827.2A Active CN113205153B (zh) 2021-05-26 2021-05-26 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型

Country Status (1)

Country Link
CN (1) CN113205153B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114462558A (zh) * 2022-04-13 2022-05-10 南昌工程学院 一种数据增广的监督学习图像缺陷分类方法与系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108647703A (zh) * 2018-04-19 2018-10-12 北京联合大学 一种基于显著性的分类图像库的类型判断方法
CN109815923A (zh) * 2019-01-30 2019-05-28 华侨大学 基于lbp特征与深度学习的金针菇菇头分选识别方法
WO2020215236A1 (zh) * 2019-04-24 2020-10-29 哈尔滨工业大学（深圳） 图像语义分割方法和系统
CN112287940A (zh) * 2020-10-30 2021-01-29 西安工程大学 一种基于深度学习的注意力机制的语义分割的方法
CN112734715A (zh) * 2020-12-31 2021-04-30 同济大学 一种肺部ct图像的肺结节分割方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108647703A (zh) * 2018-04-19 2018-10-12 北京联合大学 一种基于显著性的分类图像库的类型判断方法
CN109815923A (zh) * 2019-01-30 2019-05-28 华侨大学 基于lbp特征与深度学习的金针菇菇头分选识别方法
WO2020215236A1 (zh) * 2019-04-24 2020-10-29 哈尔滨工业大学（深圳） 图像语义分割方法和系统
CN112287940A (zh) * 2020-10-30 2021-01-29 西安工程大学 一种基于深度学习的注意力机制的语义分割的方法
CN112734715A (zh) * 2020-12-31 2021-04-30 同济大学 一种肺部ct图像的肺结节分割方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王德才;: "基于深度学习的肺癌检测方法研究", 数字技术与应用 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114462558A (zh) * 2022-04-13 2022-05-10 南昌工程学院 一种数据增广的监督学习图像缺陷分类方法与系统

Also Published As

Publication number Publication date
CN113205153B (zh) 2023-05-30

Similar Documents

Publication Publication Date Title
Candemir et al. A review on lung boundary detection in chest X-rays
CN113674253B (zh) 基于U-Transformer的直肠癌CT影像自动分割方法
Shao et al. A detection approach for solitary pulmonary nodules based on CT images
CN113283446B (zh) 图像中目标物识别方法、装置、电子设备及存储介质
CN106023151B (zh) 一种开放环境下中医舌象目标检测方法
Gaggion et al. Improving anatomical plausibility in medical image segmentation via hybrid graph neural networks: applications to chest x-ray analysis
CN108010013A (zh) 一种肺ct图像肺结节检测方法
Liu et al. A fully automatic segmentation algorithm for CT lung images based on random forest
CN114582470A (zh) 一种模型的训练方法、训练装置及医学影像报告标注方法
CN106897572A (zh) 基于流形学习的肺结节病例匹配辅助检测系统及其工作方法
Ali et al. Ensemble lung segmentation system using deep neural networks
Gupte et al. Deep learning models for calculation of cardiothoracic ratio from chest radiographs for assisted diagnosis of cardiomegaly
Li et al. Automatic annotation algorithm of medical radiological images using convolutional neural network
CN113205153B (zh) 一种儿科肺炎辅助诊断模型的训练方法及训练所得的模型
Ghali et al. Vision transformers for lung segmentation on CXR images
Malathi et al. Active contour based segmentation and classification for pleura diseases based on Otsu’s thresholding and support vector machine (SVM)
Sharafeldeen et al. Accurate segmentation for pathological lung based on integration of 3d appearance and surface models
Chaisangmongkon et al. External validation of deep learning algorithms for cardiothoracic ratio measurement
CN110837844A (zh) 基于ct图像不相似性特征的胰腺囊性肿瘤良恶性分类方法
Fonseca et al. Automatic orientation identification of pediatric chest x-rays
CN117151162A (zh) 基于自监督与专家门控的跨解剖区域器官增量分割方法
Ibrahim et al. Lung Segmentation Using ResUnet++ Powered by Variational Auto Encoder-Based Enhancement in Chest X-ray Images
Almakady et al. Volumetric texture analysis based on three-dimensional gaussian markov random fields for copd detection
Bannister et al. Detecting 3D syndromic faces as outliers using unsupervised normalizing flow models
CN115546893A (zh) 一种啦啦操视频的评价可视化方法和系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant