CN115239720A - 基于经典Graf的DDH超声影像人工智能诊断***及方法 - Google Patents
基于经典Graf的DDH超声影像人工智能诊断***及方法 Download PDFInfo
- Publication number
- CN115239720A CN115239720A CN202211155349.3A CN202211155349A CN115239720A CN 115239720 A CN115239720 A CN 115239720A CN 202211155349 A CN202211155349 A CN 202211155349A CN 115239720 A CN115239720 A CN 115239720A
- Authority
- CN
- China
- Prior art keywords
- feature
- image
- ddh
- graf
- standard section
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013473 artificial intelligence Methods 0.000 title claims abstract description 18
- 238000003745 diagnosis Methods 0.000 title abstract description 24
- 239000013598 vector Substances 0.000 claims abstract description 23
- 210000000845 cartilage Anatomy 0.000 claims abstract description 22
- 238000011176 pooling Methods 0.000 claims abstract description 21
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 19
- 238000000605 extraction Methods 0.000 claims abstract description 16
- 230000006870 function Effects 0.000 claims description 34
- 238000010586 diagram Methods 0.000 claims description 33
- 230000004913 activation Effects 0.000 claims description 12
- 238000002604 ultrasonography Methods 0.000 claims description 9
- 238000011478 gradient descent method Methods 0.000 claims description 8
- 238000013528 artificial neural network Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000003672 processing method Methods 0.000 claims description 4
- 210000000588 acetabulum Anatomy 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000002405 diagnostic procedure Methods 0.000 claims 1
- 239000000284 extract Substances 0.000 claims 1
- 210000004394 hip joint Anatomy 0.000 abstract description 6
- 239000012634 fragment Substances 0.000 abstract description 2
- 238000002372 labelling Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000013527 convolutional neural network Methods 0.000 description 2
- 210000001624 hip Anatomy 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 206010073767 Developmental hip dysplasia Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 208000014906 developmental dysplasia of the hip Diseases 0.000 description 1
- 208000035790 developmental dysplasia of the hip 1 Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
- G06V10/443—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by matching or filtering
- G06V10/449—Biologically inspired filters, e.g. difference of Gaussians [DoG] or Gabor filters
- G06V10/451—Biologically inspired filters, e.g. difference of Gaussians [DoG] or Gabor filters with interaction between the filter responses, e.g. cortical complex cells
- G06V10/454—Integrating the filters into a hierarchical structure, e.g. convolutional neural networks [CNN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/54—Extraction of image or video features relating to texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/764—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/774—Generating sets of training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/82—Arrangements for image or video recognition or understanding using pattern recognition or machine learning using neural networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了基于经典Graf的DDH超声影像人工智能诊断***及方法,该诊断***包括分类器,用于从超声影像中识别标准截面;回归器,用于定位标准截面中关键点的位置信息,并利用关键点的位置计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级;所述分类器识别标准截面的方法为:对从超声影像中截取的特征图像进行特征提取与编码,编码后的特征图进行池化操作转化为特征向量,该特征向量经过全连接层后判断该特征图是否为标准截面。本发明的诊断***能够对髋关节未成熟儿童标准髋关节超声影像片段进行筛查,且诊断更快速、客观、准确。
Description
技术领域
本发明涉及超声影像技术领域,尤其涉及基于经典Graf的DDH超声影像人工智能诊断***及方法。
背景技术
超声图像对DDH(儿童发育性髋关节发育不良)严重程度进行全谱分类的Graf 法被儿童骨科医生广为接受认可。而目前大多数医院的髋部超声筛查主要是由经验丰富的临床医生进行操作和评估的,在DDH诊断中常常存在着较低准确性及可重复性等问题。人工智能通过超声图像诊断DDH是一个新兴的研究问题。目前,针对DDH超声图像诊断,Hareendranathan等人设计了一种新颖的舍入索引,在临床医生手动标记的帮助下分析DDH。Quader等人进而提出了一种用于DDH分析的自动骨边界检测***。Sezer等人提出了一种‘分割’的两步法来分析髋关节超声。利用卷积神经网络(Convolutional neuralnetwork,CNN),可构建一个更为快速准确的DDH自动分析诊断***。但这些技术均需要在临床医生手动标记的帮助下分析DDH,仍重度依赖于医生手工标注,并未完全实现自动化。
发明内容
基于背景技术存在的技术问题,本发明提出了基于经典Graf的DDH超声影像人工智能诊断***及方法,诊断***能够对髋关节未成熟儿童标准髋关节超声影像片段进行筛查,且诊断更快速、客观、准确。
本发明提出的基于经典Graf的DDH超声影像人工智能诊断***,包括:
分类器,用于从超声影像中识别标准截面;
回归器,用于定位标准截面中关键点的位置信息,并利用关键点的位置计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
优选地,所述分类器识别标准截面的方法为:对从超声影像中截取的特征图像进行特征提取与编码,编码后的特征图像进行池化操作转化为特征向量,该特征向量经过全连接层后判断该特征图像是否为标准截面。
优选地,所述特征图像的特征提取的处理方法为:
F_{n+1}=Pool(BN(Conv(F_{n})))
式中:F_{n}表示第n阶段的特征;Pool表示2*2的最大值池化;BN表示批正则化,Conv表示卷积操作;
编码后的特征图进行池化操作转化为特征向量的方法为:
f=AvgPool(F_m)
式中:AvgPool表示全局平均池化;m为输入图像进行阶段提取的总数;F_m表示输入图像进行m个阶段的特征提取后最终获得的特向特征。
优选地,所述全连接层通过权重和特征向量完成点积,并得到代表置信系数的数值:
p’=Sigmoid(w*f)
式中:Sigmoid表示sigmoid激活函数,w是全连接权重,p’是最终预测的置信度。
优选地,所述sigmoid激活函数中的参数采用梯度下降法确定取值,所述梯度下降法依据的损失函数为:
Loss=pLog(p’)+(1-p)Log(1-p’)
式中:p是代表该特征图像是否为标准截面,为标准截面时取值1,为非标准截面时取值0。
优选地,所述回归器定位标准截面中关键点的位置信息的方法为:
S11:对从超声影像中截取的特征图像进行特征提取后得到的特征F进行3层反卷积:
F^{i+1}=ReLU(BN(Deconv(F^{i})))
式中:F^i表示经过i次反卷积层获得的特征图,ReLU是激活函数,BN表示批正则化,Deconv表示反卷积层;
S12:对经3次反卷积后获得的特征图F^3进行1*1的卷积得到每个关键点的热力图H’;
S13:对于标准截面的关键点的位置,利用高斯函数生成对应的目标热力图:
H(i,j)=exp{-((i-i’)^2+(j-j’)^2)/(2*δ^2)}
式中:i,j表示目标热力图上的坐标,i’,j’表示标准截面的关键点位置,δ表示高斯函数的衰减程度,H表示目标的热力图;
S14:使用目标的热力图H和回归器输出的热力图H’可以计算损失函数来优化神经网络:
Loss=MSE(H,H’)
式中:MSE表示均方误差函数。
优选地,所述软骨顶线夹角α为基线和软骨顶线之间的夹角,所述骨顶线夹角β为基线和骨顶线之间的夹角,所述基线为软骨膜顶点和髂骨缘点所在的直线,所述骨顶线为髂骨最低点和髋臼骨性顶点所在的直线,所述软骨顶线为软骨顶由凹变凸的转折点和盂唇中心所在的直线。
本发明提出的上述基于经典Graf的DDH超声影像人工智能诊断***的诊断方法,其特征在于,方法步骤如如下:
S21:对从超声影像中截取的特征图进行特征提取与编码;
S22:将编码后的特征图进行池化操作转化为特征向量;
S23:将特征向量通过经过全连接层后判断该特征图是否为标准截面;
S24:对S23中判断为标准截面的特征图进行上采样;
S25:对上采样的结果用1*1卷积得到热力图,并将热力图上的响应值进行加权投片,得到最终关键点的位置坐标;
S26:利用关键点的位置坐标计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
本发明提出的一种处理设备,包括存储器及处理器,所述存储器中存储有至少一种指令,所述至少一种指令由所述处理器加载并执行,以实现上述方法。
本发明提出的一种计算机可读存储介质,其特征在于,所述存储介质中存储有至少一条指令,所述指令由处理器加载并执行以实现上述方法。
本发明的有益技术效果:
本发明通过提取特征的方法是采用由卷积层,激活层,池化层,共同组合成的神经网络,处理输入的图像,本方法采用ResNet架构,提取到表征图像的特征图,提取的特征包含髋关节超声影像的边缘,纹理,相对位置等信息。可以比人工标注、诊断更为敏感地识别出未成熟的髋关节,α角和β角平均误差分别小于2.22°和2.28°。
附图说明
图1为本发明提出的基于经典Graf的DDH超声影像人工智能诊断***的流程图;
图2为本发明提出的标准截面中标注关键点后的示意图;其中1为软骨膜顶点,2为髂骨缘点,3为髂骨最低点,4为髋臼骨性顶点,5为软骨顶由凹变凸的转折点,6为盂唇中心;
图3为本发明提出的按照经典Graf法划线的示意图;
图4为本发明提出的诊断***检测的结果和人工诊断的结果对比图;其中○为人工标注的关键点,×为诊断***标注的关键点。
具体实施方式
实施例1
参照图1,本发明提出的基于经典Graf的DDH超声影像人工智能诊断***,包括:
1)分类器,用于从超声影像中识别标准截面;
作为本发明的一个实施例,该分类器识别标准截面的方法为:对从超声影像中截取的特征图像进行特征提取与编码,编码后的特征图进行池化操作转化为特征向量,该特征向量经过全连接层后判断该特征图是否为标准截面。
具体地,其中特征图像的特征提取的处理方法为:
F_{n+1}=Pool(BN(Conv(F_{n})))
式中:F_{n}表示第n阶段的特征;Pool表示2*2的最大值池化;BN表示批正则化,Conv表示卷积操作;
编码后的特征图进行池化操作转化为特征向量的方法为:
f=AvgPool(F_m)
式中:AvgPool表示全局平均池化,进行全局池化操作后得到的特征向量可以表达图像的分类信息;m为输入图像进行阶段提取的总数;F_m表示输入图像进行m个阶段的特征提取后最终获得的特向特征,对于本实施例来说,对输入图像进行5个阶段的特征提取,其中F_0表示输入图像,F_5表示最终获得的图像特征。
全连接层通过权重和特征向量完成点积,并得到代表置信系数的数值:
p’=Sigmoid(w*f)
式中:Sigmoid表示sigmoid激活函数,w是全连接权重,p’是最终预测的置信度。
sigmoid激活函数中的参数如卷积核、全连接层权重等随机初始化后,采用梯度下降法确定最终的取值,梯度下降法依据是一个由交叉熵表达的的损失函数:
Loss=pLog(p’)+(1-p)Log(1-p’)
式中:p是人工标注的该帧超声是否属于标准帧,标准帧为1,不标准帧为0;通过梯度下降法来优化该损失函数,可以让模型学习超声影像中区分标准帧的特征。
2)回归器,用于定位标准截面中关键点的位置信息,并利用关键点的位置计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
对于回归器来说,在训练阶段按,利用神经网络处理后的图像特征与医生在超声图像上标注的关键点位置构造损失函数,并对该网络进行优化与约束;作为本发明的一个优选的方案,该回归器定位标准截面中关键点的位置信息的方法为:
S11:对从超声影像中截取的特征图像进行特征提取后得到的特征F进行3层反卷积:
F^{i+1}=ReLU(BN(Deconv(F^{i})))
式中:F^i表示经过i次反卷积层获得的特征图,ReLU是激活函数,BN表示批正则化,Deconv表示反卷积层;
S12:对经3次反卷积后获得的特征图F^3进行1*1的卷积得到每个关键点的热力图H’;
S13:对于标准截面的关键点的位置,利用高斯函数生成对应的目标热力图:
H(i,j)=exp{-((i-i’)^2+(j-j’)^2)/(2*δ^2)}
式中:i,j表示目标热力图上的坐标,i’,j’表示标准截面的关键点位置,δ表示高斯函数的衰减程度,一般取值为3个像素,H表示目标的热力图;
S14:使用目标的热力图H和回归器输出的热力图H’可以计算损失函数来优化神经网络:
Loss=MSE(H,H’)
式中:MSE表示均方误差函数。对于本申请的诊断***来说,在训练整体网络的过程中,输入非关键帧是不更新回归网络的反卷积部分参数的。
此外,对于本发明中的软骨顶线夹角α和骨顶线夹角β来说,参照图2和图3软骨顶线夹角α为基线和软骨顶线之间的夹角,骨顶线夹角β为基线和骨顶线之间的夹角,基线为软骨膜顶点和髂骨缘点所在的直线,骨顶线为髂骨最低点和髋臼骨性顶点所在的直线,软骨顶线为软骨顶由凹变凸的转折点和盂唇中心所在的直线。
此外,根据图4的结果可以看出,本发明提出的诊断***检测的结果和人工诊断的结果基本相同,软骨顶线夹角α和骨顶线夹角β平均误差分别小于2.22°和2.28°,进而说明了本发明提出的人工智能诊断***的诊断结果的准确度高。
实施例2
本发明提出了基于的基于经典Graf的DDH超声影像人工智能诊断***的诊断方法,方法步骤如下:
S21:对从超声影像中截取的特征图进行特征提取与编码;
S22:将编码后的特征图进行池化操作转化为特征向量;
S23:将特征向量通过经过全连接层后判断该特征图是否为标准截面;
S24:对S23中判断为标准截面的特征图进行上采样;
S25:对上采样的结果用1*1卷积得到热力图,并将热力图上的响应值进行加权投片,得到最终关键点的位置坐标;
S26:利用关键点的位置坐标计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
其中:特征图像的特征提取的处理方法为:
F_{n+1}=Pool(BN(Conv(F_{n})))
式中:F_{n}表示第n阶段的特征;Pool表示2*2的最大值池化;BN表示批正则化,Conv表示卷积操作。
编码后的特征图进行池化操作转化为特征向量的方法为:
f=AvgPool(F_m)
式中:AvgPool表示全局平均池化;m为输入图像进行阶段提取的总数;F_m表示输入图像进行m个阶段的特征提取后最终获得的特向特征。
全连接层通过权重和特征向量完成点积,并得到代表置信系数的数值:
p’=Sigmoid(w*f)
式中:Sigmoid表示sigmoid激活函数,w是全连接权重,p’是最终预测的置信度。
sigmoid激活函数中的参数采用梯度下降法确定取值,梯度下降法依据的损失函数为:
Loss=pLog(p’)+(1-p)Log(1-p’)
式中:p是代表该特征图像是否为标准截面,为标准截面时取值1,为非标准截面时取值0。
最终关键点的位置坐标定位的方法为:
S11:对从超声影像中截取的特征图像进行特征提取后得到的特征F进行3层反卷积:
F^{i+1}=ReLU(BN(Deconv(F^{i})))
式中:F^i表示经过i次反卷积层获得的特征图,ReLU是激活函数,BN表示批正则化,Deconv表示反卷积层;
S12:对经3次反卷积后获得的特征图F^3进行1*1的卷积得到每个关键点的热力图H’;
S13:对于标准截面的关键点的位置,利用高斯函数生成对应的目标热力图:
H(i,j)=exp{-((i-i’)^2+(j-j’)^2)/(2*δ^2)}
式中:i,j表示目标热力图上的坐标,i’,j’表示标准截面的关键点位置,δ表示高斯函数的衰减程度,H表示目标的热力图;
S14:使用目标的热力图H和回归器输出的热力图H’可以计算损失函数来优化神经网络:
Loss=MSE(H,H’)
式中:MSE表示均方误差函数。
实施例3
本发明提出的一种处理设备,包括存储器及处理器,所述存储器中存储有至少一种指令,所述至少一种指令由所述处理器加载并执行,以实现实施例2中的诊断方法。
实施例4
本发明提出的一种计算机可读存储介质,其特征在于,所述存储介质中存储有至少一条指令,所述指令由处理器加载并执行以实现实施例2中的诊断方法。
以上仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,包括:
分类器,用于从超声影像中识别标准截面;
回归器,用于定位标准截面中关键点的位置信息,并利用关键点的位置信息计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
2.根据权利要求1所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述分类器识别标准截面的方法为:对从超声影像中截取的特征图像进行特征提取与编码,编码后的特征图像进行池化操作转化为特征向量,该特征向量经过全连接层后判断该特征图像是否为标准截面。
3.根据权利要求2所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述特征图像的特征提取的处理方法为:
F_{n+1}=Pool(BN(Conv(F_{n})))
式中:F_{n}表示第n阶段的特征;Pool表示2*2的最大值池化;BN表示批正则化,Conv表示卷积操作;
编码后的特征图进行池化操作转化为特征向量的方法为:
f=AvgPool(F_m)
式中:AvgPool表示全局平均池化;m为输入图像进行阶段提取的总数;F_m表示输入图像进行m个阶段的特征提取后最终获得的特向特征。
4.根据权利要求2所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述全连接层通过权重和特征向量完成点积,并得到代表置信系数的数值:
p’=Sigmoid(w*f)
式中:Sigmoid表示sigmoid激活函数,w是全连接权重,p’是最终预测的置信度。
5.根据权利要求4所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述sigmoid激活函数中的参数采用梯度下降法确定取值,所述梯度下降法依据的损失函数为:
Loss=pLog(p’)+(1-p)Log(1-p’)
式中:p是代表该特征图像是否为标准截面,为标准截面时取值1,为非标准截面时取值0。
6.根据权利要求1所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述回归器定位标准截面中关键点的位置信息的方法为:
S11:对从超声影像中截取的特征图像进行特征提取后得到的特征F进行3层反卷积:
F^{i+1}=ReLU(BN(Deconv(F^{i})))
式中:F^i表示经过i次反卷积层获得的特征图,ReLU是激活函数,Deconv表示反卷积层;
S12:对经3层反卷积后获得的特征图F^3进行1*1的卷积得到每个关键点的热力图H’;
S13:对于标准截面的关键点的位置,利用高斯函数生成对应的目标热力图:
H(i,j)=exp{-((i-i’)^2+(j-j’)^2)/(2*δ^2)}
式中:i,j表示目标热力图上的坐标,i’,j’表示标准截面的关键点位置,δ表示高斯函数的衰减程度,H表示目标的热力图;
S14:使用目标的热力图H和回归器输出的热力图H’来计算损失函数来优化神经网络:
Loss=MSE(H,H’)
式中:MSE表示均方误差函数。
7.根据权利要求1所述的基于经典Graf的DDH超声影像人工智能诊断***,其特征在于,所述软骨顶线夹角α为基线和软骨顶线之间的夹角,所述骨顶线夹角β为基线和骨顶线之间的夹角,所述基线为软骨膜顶点和髂骨缘点所在的直线,所述骨顶线为髂骨最低点和髋臼骨性顶点所在的直线,所述软骨顶线为软骨顶由凹变凸的转折点和盂唇中心所在的直线。
8.如权利要求1-7任一项所述的基于经典Graf的DDH超声影像人工智能诊断***的诊断方法,其特征在于,方法步骤如下:
S21:分类器对从超声影像中截取的特征图像进行特征提取与编码;
S22:分类器将编码后的特征图像进行池化操作转化为特征向量;
S23:分类器将特征向量通过经过全连接层后判断该特征图像是否为标准截面;
S24:对S23中判断为标准截面的特征图像进行上采样;
S25:对上采样的结果用1*1卷积得到热力图,并将热力图上的响应值进行加权投片,得到最终关键点的位置信息;
S26:利用关键点的位置信息计算Graf方法需要的软骨顶线夹角α、骨顶线夹角β,判断DDH的等级。
9.一种处理设备,其特征在于,包括存储器及处理器,所述存储器中存储有至少一种指令,所述至少一种指令由所述处理器加载并执行,以实现权利要求8中所述的方法。
10.一种计算机可读存储介质,其特征在于,所述存储介质中存储有至少一条指令,所述指令由处理器加载并执行以实现权利要求8中所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211155349.3A CN115239720A (zh) | 2022-09-22 | 2022-09-22 | 基于经典Graf的DDH超声影像人工智能诊断***及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211155349.3A CN115239720A (zh) | 2022-09-22 | 2022-09-22 | 基于经典Graf的DDH超声影像人工智能诊断***及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115239720A true CN115239720A (zh) | 2022-10-25 |
Family
ID=83667214
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211155349.3A Pending CN115239720A (zh) | 2022-09-22 | 2022-09-22 | 基于经典Graf的DDH超声影像人工智能诊断***及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115239720A (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20170006946A (ko) * | 2015-07-10 | 2017-01-18 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 초음파 진단 장치의 동작 방법 |
US20180259608A1 (en) * | 2015-11-29 | 2018-09-13 | Arterys Inc. | Automated cardiac volume segmentation |
CN109919943A (zh) * | 2019-04-16 | 2019-06-21 | 广东省妇幼保健院 | 婴儿髋关节角度自动检测方法、***和计算设备 |
CN110895809A (zh) * | 2019-10-18 | 2020-03-20 | 中国科学技术大学 | 准确提取髋关节影像中关键点的方法 |
CN111882531A (zh) * | 2020-07-15 | 2020-11-03 | 中国科学技术大学 | 髋关节超声图像自动分析方法 |
CN112215829A (zh) * | 2020-10-21 | 2021-01-12 | 深圳度影医疗科技有限公司 | 一种髋关节标准切面的定位方法和计算机设备 |
US20210056293A1 (en) * | 2019-08-19 | 2021-02-25 | Zhuhai Eeasy Technology Co., Ltd. | Face detection method |
CN112907507A (zh) * | 2021-01-14 | 2021-06-04 | 杭州米迪智能科技有限公司 | Graf法髋关节超声图像测量方法、装置、设备及存储介质 |
CN113724328A (zh) * | 2021-08-31 | 2021-11-30 | 瓴域影诺(北京)科技有限公司 | 一种髋关节关键点检测方法和*** |
CN113951925A (zh) * | 2021-10-28 | 2022-01-21 | 扬州市妇幼保健院 | 基于Graf超声技术的髋关节标准超声图像获取方法及智能*** |
CN114795258A (zh) * | 2022-04-18 | 2022-07-29 | 浙江大学 | 一种儿童髋关节发育不良诊断*** |
CN114862789A (zh) * | 2022-05-04 | 2022-08-05 | 上海交通大学医学院附属仁济医院 | 一种基于深度学习的小儿髋关节超声图像质量检测方法 |
-
2022
- 2022-09-22 CN CN202211155349.3A patent/CN115239720A/zh active Pending
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20170006946A (ko) * | 2015-07-10 | 2017-01-18 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 초음파 진단 장치의 동작 방법 |
US20180259608A1 (en) * | 2015-11-29 | 2018-09-13 | Arterys Inc. | Automated cardiac volume segmentation |
CN109919943A (zh) * | 2019-04-16 | 2019-06-21 | 广东省妇幼保健院 | 婴儿髋关节角度自动检测方法、***和计算设备 |
US20210056293A1 (en) * | 2019-08-19 | 2021-02-25 | Zhuhai Eeasy Technology Co., Ltd. | Face detection method |
CN110895809A (zh) * | 2019-10-18 | 2020-03-20 | 中国科学技术大学 | 准确提取髋关节影像中关键点的方法 |
CN111882531A (zh) * | 2020-07-15 | 2020-11-03 | 中国科学技术大学 | 髋关节超声图像自动分析方法 |
CN112215829A (zh) * | 2020-10-21 | 2021-01-12 | 深圳度影医疗科技有限公司 | 一种髋关节标准切面的定位方法和计算机设备 |
CN112907507A (zh) * | 2021-01-14 | 2021-06-04 | 杭州米迪智能科技有限公司 | Graf法髋关节超声图像测量方法、装置、设备及存储介质 |
CN113724328A (zh) * | 2021-08-31 | 2021-11-30 | 瓴域影诺(北京)科技有限公司 | 一种髋关节关键点检测方法和*** |
CN113951925A (zh) * | 2021-10-28 | 2022-01-21 | 扬州市妇幼保健院 | 基于Graf超声技术的髋关节标准超声图像获取方法及智能*** |
CN114795258A (zh) * | 2022-04-18 | 2022-07-29 | 浙江大学 | 一种儿童髋关节发育不良诊断*** |
CN114862789A (zh) * | 2022-05-04 | 2022-08-05 | 上海交通大学医学院附属仁济医院 | 一种基于深度学习的小儿髋关节超声图像质量检测方法 |
Non-Patent Citations (6)
Title |
---|
JINGYUAN XU 等: "Hip Landmark Detection With Dependency Mining in Ultrasound Image", 《IEEE TRANSACTIONS ON MEDICAL IMAGING》 * |
JINGYUAN XU 等: "Multi-task hourglass network for online automatic diagnosis of developmental dysplasia of the hip", 《WORLD WIDE WEB》 * |
SICHENG ZHANG 等: "Collagen I in the Hip Capsule Plays a Role in Postoperative Clinical Function in Patients With Developmental Dysplasia of the Hip", 《FRONTIERS IN PEDIATRICS》 * |
YAOSHENG LU 等: "Multitask Deep Neural Network for the Fully Automatic Measurement of the Angle of Progression", 《COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE》 * |
ZHOUBING XU 等: "Less is More: Simultaneous View Classification", 《ARXIV:1805.10376V2 [CS.CV]》 * |
孙锡玮 等: "超声人工智能辅助诊断发育性髋关节发育不良", 《中华骨科杂志》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609525B (zh) | 基于剪枝策略构建卷积神经网络的遥感图像目标检测方法 | |
WO2019200747A1 (zh) | 分割股骨近端的方法、装置、计算机设备和存储介质 | |
CN111259930A (zh) | 自适应注意力指导机制的一般性目标检测方法 | |
CN112967243A (zh) | 一种基于yolo的深度学习芯片封装裂纹缺陷检测方法 | |
CN112085024A (zh) | 一种罐表面字符识别方法 | |
CN114972213A (zh) | 一种基于机器视觉的两阶段主板图像缺陷检测及定位方法 | |
CN113378676A (zh) | 基于多特征融合的图像中人物交互检测方法 | |
CN109801305B (zh) | 基于深度胶囊网络的sar图像变化检测方法 | |
CN109949280B (zh) | 图像处理方法、装置、设备存储介质及生长发育评估*** | |
CN112926652B (zh) | 一种基于深度学习的鱼类细粒度图像识别方法 | |
CN117854072B (zh) | 一种工业视觉缺陷自动标注方法 | |
CN112598031A (zh) | 一种蔬菜病害检测方法及*** | |
CN116665011A (zh) | 一种基于机器视觉的煤矿带式输送机煤流异物识别方法 | |
CN114548253A (zh) | 一种基于图像识别及动态匹配的数字孪生模型构建*** | |
CN115294033A (zh) | 一种基于语义分割网络的轮胎带束层差级和错边缺陷检测方法 | |
CN112991280A (zh) | 视觉检测方法、***及电子设备 | |
CN116452899A (zh) | 一种基于深度学习的超声心动图标准切面识别及评分方法 | |
CN115995040A (zh) | 一种基于多尺度网络的sar图像小样本目标识别方法 | |
CN115861226A (zh) | 基于特征值梯度变化的深度神经网络智能识别表面缺陷方法 | |
Wu et al. | Fast particle picking for cryo-electron tomography using one-stage detection | |
CN112396580B (zh) | 一种圆形零件缺陷检测方法 | |
CN117745708A (zh) | 基于深度学习算法的木地板表面瑕疵检测方法 | |
CN111223113B (zh) | 基于双重密集上下文感知网络的核磁共振海马体分割算法 | |
CN116245855B (zh) | 作物品种鉴定方法、装置、设备及存储介质 | |
CN115239720A (zh) | 基于经典Graf的DDH超声影像人工智能诊断***及方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20221025 |
|
RJ01 | Rejection of invention patent application after publication |