CN113222951B - 一种识别髋关节x线的骨质疏松人工智能诊断装置 - Google Patents
一种识别髋关节x线的骨质疏松人工智能诊断装置 Download PDFInfo
- Publication number
- CN113222951B CN113222951B CN202110549854.5A CN202110549854A CN113222951B CN 113222951 B CN113222951 B CN 113222951B CN 202110549854 A CN202110549854 A CN 202110549854A CN 113222951 B CN113222951 B CN 113222951B
- Authority
- CN
- China
- Prior art keywords
- image
- neural network
- network model
- hip joint
- ray
- 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.)
- Active
Links
- 210000004394 hip joint Anatomy 0.000 title claims abstract description 41
- 208000001132 Osteoporosis Diseases 0.000 title claims abstract description 35
- 238000003745 diagnosis Methods 0.000 title claims abstract description 31
- 238000003062 neural network model Methods 0.000 claims abstract description 55
- 238000003709 image segmentation Methods 0.000 claims abstract description 31
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 29
- 238000012549 training Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000010276 construction Methods 0.000 claims abstract description 7
- 210000002436 femur neck Anatomy 0.000 claims description 31
- 238000013527 convolutional neural network Methods 0.000 claims description 26
- 230000037182 bone density Effects 0.000 claims description 25
- 238000002372 labelling Methods 0.000 claims description 19
- 238000013528 artificial neural network Methods 0.000 claims description 11
- 230000007704 transition Effects 0.000 claims description 8
- 238000011176 pooling Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 7
- 238000007689 inspection Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000013473 artificial intelligence Methods 0.000 claims 4
- 238000000034 method Methods 0.000 claims 1
- 238000012216 screening Methods 0.000 abstract description 3
- 210000001624 hip Anatomy 0.000 description 10
- 206010017076 Fracture Diseases 0.000 description 7
- 208000010392 Bone Fractures Diseases 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004195 computer-aided diagnosis Methods 0.000 description 2
- 230000001054 cortical effect Effects 0.000 description 2
- 238000013136 deep learning model Methods 0.000 description 2
- 238000001739 density measurement Methods 0.000 description 2
- 210000000527 greater trochanter Anatomy 0.000 description 2
- 210000000528 lesser trochanter Anatomy 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000037411 Aortic calcification Diseases 0.000 description 1
- 208000020084 Bone disease Diseases 0.000 description 1
- 206010005963 Bone formation increased Diseases 0.000 description 1
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 206010039984 Senile osteoporosis Diseases 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000003187 abdominal effect Effects 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000037180 bone health Effects 0.000 description 1
- 230000002308 calcification Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000009693 chronic damage Effects 0.000 description 1
- 208000037976 chronic inflammation Diseases 0.000 description 1
- 230000006020 chronic inflammation Effects 0.000 description 1
- 230000001447 compensatory effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 230000004821 effect on bone Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 208000020089 femoral neck fracture Diseases 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000652 homosexual effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 210000001519 tissue Anatomy 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/505—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/048—Activation functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- 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/10116—X-ray 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/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- 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)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Software Systems (AREA)
- Veterinary Medicine (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Quality & Reliability (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Physiology (AREA)
- Image Analysis (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种识别髋关节X线的骨质疏松人工智能诊断装置，包括：数据获取模块、模型构建模块、模型训练模块、图像检测模块；其中，数据获取模块用于获取髋关节Ward三角区标注的医学图像，模型构建模块用于构建图像分割神经网络模型和图像分类神经网络模型，模型训练模块用于对两个神经网络模型进行训练，图像检测模块用于对待检测的髋关节X线片通过训练后的图像分割神经网络模型和图像分类神经网络模型的处理，获取骨质状况信息数值。采用本发明装置，可以实现仅识别髋关节X线片即可辅助初筛骨质疏松的目的，简单、有效的提升骨质疏松的诊断率。
Description
技术领域
本发明属于医学图像处理技术领域，特别是涉及一种识别髋关节X线的骨质疏松人工智能诊断装置。
背景技术
骨折风险随年龄呈指数增加，50岁以后，几乎每2名女性或每5名男性中就有1名会遭受脆性骨折影响，其中骨质疏松是易发骨折的一个重要原因。
目前，骨质疏松的早期诊断主要依据双能X线吸收测定法(dualenergy X-rayabsorptiometry，DXA)及进行骨密度测量，主要检测部位包括髋部及腰椎前后位,可以说是评估骨质疏松症和骨折风险的重要标准。定量计算机断层扫描(QCT)是一种比较先进的检测方法，它不仅可以测量皮质骨的骨密度，还可以测量松质骨的骨密度。大量研究表明，用QCT测量BMD得到的结果与DXA高度一致，其优点是没有重叠效应，但患者接受的辐射剂量较大，检查费用相对较高，所以临床应用较少。骨密度通常用T-值(T-Score)表示，T-值＝(实测值－同种族同性别正常青年人峰值骨密度)/同种族同性别正常青年人峰值骨密度的标准差。基于DXA测量的中轴骨(腰椎1-4、股骨颈或全髋)骨密度或桡骨远端1/3骨密度对骨质疏松症的诊断标准是T-值≤－2.5。对于绝经后女性、50岁及以上男性，建议参照WHO推荐的诊断标准，骨密度值低于同性别、同种族健康成人的骨峰值1个标准差及以内属正常；降低1～2.5个标准差为骨量低下(或低骨量)；降低等于和超过2.5个标准差为骨质疏松；骨密度降低程度符合骨质疏松诊断标准，同时伴有一处或多处脆性骨折为严重骨质疏松。对于儿童、绝经前女性和50岁以下男性，其骨密度水平的判断建议用同种族的Z-值表示，Z-值＝(骨密度测定值－同种族同性别同龄人骨密度均值)/同种族同性别同龄人骨密度标准差。将Z-值≤－2.0视为“低于同年龄段预期范围”或低骨量。
临床常规影像学检查对评估骨小梁结构的价值有限，并且缺乏相关的定量参数，对骨量减少的诊断效果并不理想。随着理论的发展和信号处理技术的进步，计算机辅助诊断(CAD)已逐渐应用于医学影像的研究。近年来，越来越多的人开始关注骨骼图像的纹理分析。有研究结果表明，纹理分析对骨病诊断具有实用的潜在价值。最近，以卷积神经网络为代表的深度学习模型在医学图像处理中取得了巨大的成功。回顾大量股骨颈纹理分析中的骨质疏松分类，发现卷积神经网络具有较好的分类效果，融合深层特征和底层特征后，分类效果也会有所提高。
髋部Ward三角是在大、小转子和转子间嵴，也就是在压力曲线的支持束和转子束之间，存在骨小梁的薄弱区，骨密度相对最低，称为Ward三角区，为股骨颈骨折好发区，在髋部X线正位片上亦可识别。尽管Ward三角区松质骨含量高,对骨密度变化敏感性高，但由于其投射面积小,取样和重复性误差导致其精确度差,限制了它的临床应用。但Ward三角区含有80％以上的松质骨，是钙储备区，亦有学者认为其对骨密度变化敏感性高，对骨折的风险性预测较好。
由于DXA是二维重叠图像，其缺点是不能区分松质骨和皮质骨，不能处理脊柱重叠组织引起的伪影，如主动脉钙化、腹腔钙化等，从而影响BMD测量结果的准确性。并且以髋部和腰椎骨密度降低来确诊骨质疏松症是有局限性的。因老年骨质疏松患者病程长，慢性炎症或损伤修复导致不同程度代偿性骨质增生，相应部位的骨密度下降不明显。仅以常规DXA检测和髋部、腰椎骨密度来诊断骨质疏松有很大的漏诊率。
同时现有诊疗体系中，患者在骨质检查时往往需要同时做X线和DXA检查，既增加了患者负担，同时过量影像数据也增加了医生的工作量，容易出现漏诊和误诊。
发明内容
本发明提出一种识别髋关节X线的骨质疏松人工智能诊断装置，通过对髋关节X线图像中股骨颈区及其内部Ward三角区进行标注，经由深度学习模型训练，实现在髋关节X线片上预测髋关节股骨颈区及其内部的Ward三角区的形态、位置和二者的骨密度值，从而实现仅识别髋关节X线片即可初筛骨质疏松的目的，简单、有效的提升骨质疏松的诊断准确度，降低漏诊和误诊，提高临床骨质疏松的诊断效率，降低患者经济负担和医生工作量。
为实现上述目的，本发明提供了如下方案：
一种识别髋关节X线的骨质疏松人工智能诊断装置，包括：数据获取模块、模型构建模块、模型训练模块、图像检测模块；
所述数据获取模块，用于获取股骨颈区及其内部Ward三角区标注的医学图像；
所述模型构建模块，用于构建图像分割神经网络模型和图像分类神经网络模型；所述图像分割神经网络模型用于获取所述医学图像的分割图像；所述分类神经网络模型用于对分割图像进行分类识别，得到骨质状况信息数值；
所述模型训练模块，用于对所述图像分割神经网络模型和所述图像分类神经网络模型进行训练；
所述图像检测模块，用于将待检测的未经股骨颈区及Ward三角区标注的所述医学图像通过训练后的所述图像分割神经网络模型的处理，得到分割图像，利用训练后的所述图像分类神经网络模型对所述分割图像的分类识别，获取骨质状况信息数值。
优选的，所述医学图像为髋关节X线图像。
优选的，所述图像分割神经网络模型通过获取所述医学图像的标注图像,并通过对所述标注图像和所述医学图像的图像处理，得到所述分割图像。
优选的，所述图像分割神经网络模型采用图像分割神经网络对输入的所述医学图像的图像特征进行提取与映射,得到所述标注图像；所述图像分割神经网络包含但不限于U-Net卷积神经网络；
所述图像特征包括股骨颈区及其内部的Ward三角区标注图像。
优选的，所述图像分割神经网络包括相互连接的编码器和解码器；
所述编码器设置有卷积层，用于对所述医学图像的所述图像特征进行提取；
所述解码器设置有反卷积层，用于从所述图像特征中获取股骨颈区及其内部的Ward三角区的所述标注图像；
优选的，所述图像分类神经网络模型采用图像分类神经网络对所述分割图像进行分类判别；所述图像分类神经网络包括但不限于DenseNet密集卷积神经网络，
优选的，所述DenseNet密集卷积神经网络包括若干个稠密块和若干个过渡层；
所述过渡层与所述稠密块相间连接。
优选的，所述DenseNet密集卷积神经网络还包括瓶颈层；
所述瓶颈层用于减少所述分割图像的输入数量；
所述瓶颈层设于所述图像分类神经网络模型的输入端口，所述瓶颈层的输出与所述DenseNet密集卷积神经网络中的首个所述稠密块连接。
优选的，所述DenseNet密集卷积神经网络还包括全局平均池化层和分类器；
所述全局平均池化层和所述分类器与所述DenseNet密集卷积神经网络中的最后一个所述稠密块连接。
本发明的有益效果为：
本发明公开了一种识别髋关节X线的骨质疏松人工智能诊断装置，基于已有的DXA检查数值结果，打造骨质疏松诊断装置，实现了在髋关节X线图像上通过识别股骨颈区及其内部Ward三角区来进行智能的骨质疏松初筛及辅助诊断；通过不断优化装置的识别性能，简单、有效的提升骨质疏松的诊断准确度，降低漏诊和误诊，提高临床骨质疏松的诊断效率。更进一步的，患者只需要普通的髋关节X片检查就可以得到相关的骨质健康信息，减轻患者经济压力，也减轻医生的工作量，提升医院诊疗效率，也减少了二次检查和放射线暴露。
附图说明
为了更清楚地说明本发明的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例的结构示意图；
图2为本发明实施例髋部Ward三角示意图；
图3为本发明实施例中U-Net卷积神经网络模型示意图；
图4为采用本实施例装置诊断骨质疏松的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示，为本发明提供的一种识别髋关节X线的骨质疏松人工智能诊断装置结构示意图，包括：数据获取模块、模型构建模块、模型训练模块、图像检测模块；
数据模块用于获取髋关节X线影像中股骨颈区及其内部Ward三角区标注的医学图像。
图2为髋部Ward三角示意图，位于大、小转子和转子间嵴，松质骨含量高，对骨密度变化敏感性高，对骨折的风险性预测较好。
模型构建模块用于构建图像分割神经网络模型和图像分类神经网络模型；
其中，图像分割神经网络模型用于获取医学图像的分割图像，在本实施例中，采用U-Net卷积神经网络对输入的髋关节X线图像的图像特征进行提取与映射,得到所述标注图像，再通过对该标注图像的图像处理，得到分割图像，作为后续图像分类神经网络模型的输入数据。在本实施例中，通过剪切和缩放方式获取分割图像。
在本实施例中，U-Net卷积神经网络包括相互连接的编码器和解码器；其中，编码器设置有卷积层，用于对所述医学图像的图像特征进行提取；解码器设置有反卷积层，用于从所述图像特征中获取所述标注图像；
在本实施例中，为更好地获取图像细节信息，编码器和解码器的各层之间设有用于获取髋部标注图像信息的跨层连接。
U-Net整个神经网络主要有两部分组成：收缩路径和扩展路径。收缩路径主要是用来捕捉图片中的上下文信息，而与之相对称的扩展路径则是为了对图片中所需要分割出来的部分进行精准定位。在本实施例中，图像分割神经网络模型采用U-Net卷积神经网络，以实现对髋关节X线图像的特征提取与映射，其输入为髋关节X线图像，输出为髋关节X线图像对应的标注图像，可以很好地精确定位图像中的关注区域。
在本实施例中，如图3所示，U-Net卷积神经网络的编码器利用九个卷积层从输入数据中提取图像特征。原始图像大小为i×512×512×1，其中i为训练数据的批量大小。前三个卷积层分别为64、128、256通道，卷积核大小为4×4，重叠步长为2。第四卷积层到第九卷积层共有512个通道，卷积核大小为4×4，重叠步长为2。卷积层的激活函数为LeakyReLU，斜率为0.2。解码器的目的是从所获得的髋关节X线图像特征信息中获取股骨颈区及其内部Ward三角区的标注图像。解码器由相应的八个反卷积层组成。前六个反卷积层有512个通道，卷积核大小为4×4，重叠步长为2。第七个反卷积层至最后的反卷积层分别为256、128、64通道，卷积核大小为4×4，重叠步长为2。在本实施例中，为了帮助解码器更好地获取标注图像的细节，引入跨层连接的操作方式连接相应的编码器和解码器层。卷积神经网络的输出图像大小与输入图像大小相同。
在本实施例中，所述分类神经网络模型采用DenseNet密集卷积神经网络，根据分割图像的骨质状况，实现分割图像的分类识别，得到骨质状况信息数值；
在本实施例中，所述DenseNet密集卷积神经网络包括若干个稠密块和若干个过渡层；稠密块与过渡层相间连接，即稠密块位于两个过渡层之间，因此稠密块数量比过渡层数量多一个。
在本实施例中，为了减少所述分割图像的输入数量，DenseNet密集卷积神经网络还设有瓶颈层，包含卷积层和平均池化层；瓶颈层位于图像分类神经网络模型的输入端口，其输出与首个稠密块连接。
在本实施例中，在所述DenseNet密集卷积神经网络的结构中，最后一个稠密块还连接有全局平均池化层和分类器。
在本实施例中采用上述瓶颈层-DenseNet密集卷积神经网络-全局平均池化层和分类器的结构，能够更加有效减少输入特征图数量，提高计算效率。
所述模型训练模块，利用样本数据对图像分割神经网络模型和图像分类神经网络模型进行训练和优化，提高诊断准确度。
所述图像检测模块，用于将待检测未经股骨颈区及Ward三角区标注的所述医学图像通过训练后的图像分割神经网络模型的输入，获取待检测髋关节X线图像的标注图像，通过对待检测髋关节X线图像和标注图像进行剪切和缩放，得到分割图像，利用训练后的图像分类神经网络模型对得到的分割图像进行分类识别，获取图像对应的骨质状况信息，通过模型输出的数值进行骨质疏松程度评估。
图4为采用本实施例装置诊断骨质疏松的流程示意图，包括如下步骤：
S1.获取髋关节X线图像及骨密度DXA检查图像数据，人工标注髋关节X线图像中股骨颈区及其内部的Ward三角区，利用对应骨密度DXA检查数值结果进行图像赋值；
S2.构建图像分割神经网络模型和图像分类神经网络模型，其中，图像分割神经网络模型用于通过训练学习以获取髋关节X线图像，髋关节X线图像通过剪切和缩放获取分割图像，该分割图像作为图像分类神经网络模型的输入数据；分类神经网络模型用于通过训练学习对输入的分割图像进行骨质状况分类识别，得到骨质状况信息数值；
S3.利用样本数据对图像分割神经网络模型和图像分类神经网络模型进行训练；
S4.将待检测髋关节X线图像作为训练后的图像分割神经网络模型的输入，获取待检测髋关节X线图像的标注图像，通过对待检测髋关节X线图像和标注图像进行剪切和缩放，得到分割图像；
S5.利用训练后的图像分类神经网络模型对得到的分割图像进行分类判别，获取图像对应的骨质状况信息。
以上所述的实施例仅是对本发明优选方式进行的描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。
Claims (5)
1.一种识别髋关节X线的骨质疏松人工智能诊断装置，其特征在于，包括：数据获取模块、模型构建模块、模型训练模块、图像检测模块；
所述数据获取模块，用于获取股骨颈区及其内部Ward三角区标注的医学图像；
所述模型构建模块，用于构建图像分割神经网络模型和图像分类神经网络模型；所述图像分割神经网络模型用于获取所述医学图像的分割图像；所述分类神经网络模型用于对所述分割图像进行分类识别，得到骨质状况信息数值；
所述模型训练模块，用于对所述图像分割神经网络模型和所述图像分类神经网络模型进行训练；
所述图像检测模块，用于将待检测的未经股骨颈区及Ward三角区标注的所述医学图像通过训练后的所述图像分割神经网络模型的处理，得到分割图像，利用训练后的所述图像分类神经网络模型对所述分割图像的分类识别，获取骨质状况信息数值；
所述医学图像为髋关节X线图像；
所述图像分割神经网络模型通过获取所述医学图像的标注图像,并通过对所述标注图像和所述医学图像的图像处理，得到所述分割图像；
所述图像分割神经网络模型采用图像分割神经网络对输入的所述医学图像的图像特征进行提取与映射,得到所述标注图像；所述图像分割神经网络包括但不限于U-Net卷积神经网络；
所述图像特征包括股骨颈区及其内部的Ward三角区标注图像；
所述图像分割神经网络包括相互连接的编码器和解码器；
所述编码器设置有卷积层，用于对所述医学图像的所述图像特征进行提取；
所述解码器设置有反卷积层，用于从所述图像特征中获取股骨颈区及其内部的Ward三角区的所述标注图像；
使用所述诊断装置诊断骨质疏松的方法包括如下步骤：
S1.获取髋关节X线图像及骨密度DXA检查图像数据，人工标注髋关节X线图像中股骨颈区及其内部的Ward三角区，利用对应骨密度DXA检查数值结果进行图像赋值；
S2.构建图像分割神经网络模型和图像分类神经网络模型，其中，图像分割神经网络模型用于通过训练学习以获取髋关节X线图像，髋关节X线图像通过剪切和缩放获取分割图像，该分割图像作为图像分类神经网络模型的输入数据；分类神经网络模型用于通过训练学习对输入的分割图像进行骨质状况分类识别，得到骨质状况信息数值；
S3.利用样本数据对图像分割神经网络模型和图像分类神经网络模型进行训练；
S4.将待检测髋关节X线图像作为训练后的图像分割神经网络模型的输入，获取待检测髋关节X线图像的标注图像，通过对待检测髋关节X线图像和标注图像进行剪切和缩放，得到分割图像；
S5.利用训练后的图像分类神经网络模型对得到的分割图像进行分类判别，获取图像对应的骨质状况信息。
2.根据权利要求1所述的识别髋关节X线的骨质疏松人工智能诊断装置，其特征在于，
所述图像分类神经网络模型采用图像分类神经网络对所述分割图像进行分类判别；所述图像分类神经网络包括但不限于DenseNet密集卷积神经网络。
3.根据权利要求2所述的识别髋关节X线的骨质疏松人工智能诊断装置，其特征在于，
所述DenseNet密集卷积神经网络还包括若干个稠密块和若干个过渡层；
所述过渡层与所述稠密块相间连接。
4.根据权利要求3所述的识别髋关节X线的骨质疏松人工智能诊断装置，其特征在于，
所述DenseNet密集卷积神经网络还包括瓶颈层；
所述瓶颈层用于减少所述分割图像的输入数量；
所述瓶颈层设于所述图像分类神经网络模型的输入端口，所述瓶颈层的输出与所述DenseNet密集卷积神经网络中的首个所述稠密块连接。
5.根据权利要求4所述的识别髋关节X线的骨质疏松人工智能诊断装置，其特征在于，
所述DenseNet密集卷积神经网络还包括全局平均池化层和分类器；
所述全局平均池化层和所述分类器与所述DenseNet密集卷积神经网络中的最后一个所述稠密块连接。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549854.5A CN113222951B (zh) | 2021-05-20 | 2021-05-20 | 一种识别髋关节x线的骨质疏松人工智能诊断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110549854.5A CN113222951B (zh) | 2021-05-20 | 2021-05-20 | 一种识别髋关节x线的骨质疏松人工智能诊断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113222951A CN113222951A (zh) | 2021-08-06 |
| CN113222951B true CN113222951B (zh) | 2023-05-02 |
Family
ID=77093771
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110549854.5A Active CN113222951B (zh) | 2021-05-20 | 2021-05-20 | 一种识别髋关节x线的骨质疏松人工智能诊断装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113222951B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114049315B (zh) * | 2021-10-29 | 2023-04-18 | 北京长木谷医疗科技有限公司 | 关节识别方法、电子设备、存储介质及计算机程序产品 |
| CN113822231A (zh) * | 2021-11-08 | 2021-12-21 | 中国人民解放军陆军特色医学中心 | 一种基于深度学习图像识别的转子间骨折手术辅助*** |
| CN114723763B (zh) * | 2022-05-24 | 2022-09-02 | 博志生物科技(深圳)有限公司 | 一种医学图像分割方法、装置、设备及存储介质 |
| CN116570367A (zh) * | 2023-05-12 | 2023-08-11 | 北京长木谷医疗科技股份有限公司 | 机器人手术操作骨磨削骨质智能感知预测方法装置及设备 |
| CN117635951B (zh) * | 2024-01-24 | 2024-05-03 | 苏州大学附属第二医院 | 基于x线图像自动识别髋部骨质疏松的判定方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1633594A (zh) * | 2002-02-27 | 2005-06-29 | 成像治疗仪股份有限公司 | X射线图像定量分析的方法和装置 |
| CN1682236A (zh) * | 2002-08-20 | 2005-10-12 | 成像治疗仪股份有限公司 | 用于x射线图像分析的方法和设备 |
| CN110033438A (zh) * | 2019-03-14 | 2019-07-19 | 上海市第六人民医院 | 髋关节标记***及其标记方法 |
| CN112396591A (zh) * | 2020-11-25 | 2021-02-23 | 暨南大学附属第一医院（广州华侨医院） | 一种基于腰椎x线图像的骨质疏松智能评估方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070047794A1 (en) * | 2000-10-11 | 2007-03-01 | Philipp Lang | Methods and devices for analysis of x-ray images |
| CN110648337A (zh) * | 2019-09-23 | 2020-01-03 | 武汉联影医疗科技有限公司 | 髋关节分割方法、装置、电子设备和存储介质 |
| CN110796636A (zh) * | 2019-09-25 | 2020-02-14 | 中国人民解放军战略支援部队信息工程大学 | 基于卷积神经网络的ct图像骨质状况检测方法及装置 |
-
2021
- 2021-05-20 CN CN202110549854.5A patent/CN113222951B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1633594A (zh) * | 2002-02-27 | 2005-06-29 | 成像治疗仪股份有限公司 | X射线图像定量分析的方法和装置 |
| CN1682236A (zh) * | 2002-08-20 | 2005-10-12 | 成像治疗仪股份有限公司 | 用于x射线图像分析的方法和设备 |
| CN110033438A (zh) * | 2019-03-14 | 2019-07-19 | 上海市第六人民医院 | 髋关节标记***及其标记方法 |
| CN112396591A (zh) * | 2020-11-25 | 2021-02-23 | 暨南大学附属第一医院（广州华侨医院） | 一种基于腰椎x线图像的骨质疏松智能评估方法 |
Non-Patent Citations (2)
| Title |
|---|
| S. M. Nazia Fathima 等.Estimation of t-score and BMD values from X-ray images for detection of osteoporosis.《ICCSP '19: Proceedings of the 3rd International Conference on Cryptography, Security and Privacy》.2019,220–224. * |
| 高小玲 等.原发性骨淋巴瘤影像学特征及误诊分析.《华南国防医学杂志》.2018,26-30+57. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113222951A (zh) | 2021-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113222951B (zh) | 一种识别髋关节x线的骨质疏松人工智能诊断装置 | |
| US9460506B2 (en) | System and method for predicting future fractures | |
| Pietka et al. | Feature extraction in carpal-bone analysis | |
| US8965087B2 (en) | System and method of predicting future fractures | |
| Areeckal et al. | Current and emerging diagnostic imaging-based techniques for assessment of osteoporosis and fracture risk | |
| Wani et al. | Computer-aided diagnosis systems for osteoporosis detection: a comprehensive survey | |
| JP2007524438A (ja) | 放射線画像処理技術における補償の方法 | |
| Liu et al. | Diagnostic and gradation model of osteoporosis based on improved deep U-Net network | |
| Hussain et al. | Computer-aided osteoporosis detection from DXA imaging | |
| Areeckal et al. | Early diagnosis of osteoporosis using radiogrammetry and texture analysis from hand and wrist radiographs in Indian population | |
| Tang et al. | CNN-based qualitative detection of bone mineral density via diagnostic CT slices for osteoporosis screening | |
| US7539332B1 (en) | Method and system for automatically identifying regions of trabecular bone tissue and cortical bone tissue of a target bone from a digital radiograph image | |
| US7387439B2 (en) | X-ray beam calibration for bone mineral density assessment using mammography system | |
| US20160180520A1 (en) | Quantitative method for 3-d joint characterization | |
| Seo et al. | Deep focus approach for accurate bone age estimation from lateral cephalogram | |
| Nazia Fathima et al. | A Survey on Osteoporosis Detection Methods with a Focus on X-ray and DEXA Images | |
| KR20210054925A (ko) | 골밀도 산출을 위한 관심영역 추출 시스템 및 방법 | |
| Xue et al. | A dual-selective channel attention network for osteoporosis prediction in computed tomography images of lumbar spine | |
| Sha et al. | The improved faster-RCNN for spinal fracture lesions detection | |
| CN116433620A (zh) | 基于ct图像的骨密度预测及骨质疏松智能筛查方法和*** | |
| KR102258070B1 (ko) | 발의 유형 평가 방법 및 이를 이용한 발의 유형 평가용 디바이스 | |
| Patil et al. | Detection of osteoporosis in lumbar spine [L1-L4] trabecular bone: a review article | |
| Taleb-Ahmed et al. | On information fusion to improve segmentation of MRI sequences | |
| Su et al. | A Hierarchical Opportunistic Screening Model for Osteoporosis Using Deep Learning | |
| Shakya et al. | Benchmarking encoder-decoder architectures for biplanar X-ray to 3D shape reconstruction |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Qin Yanguo Inventor after: He Meng Inventor after: Guo Deming Inventor after: Tang Xiongfeng Inventor after: Shen Xianyue Inventor after: Chen Bo Inventor before: Qin Yanguo Inventor before: He Meng Inventor before: Tang Xiongfeng Inventor before: Guo Deming Inventor before: Shen Xianyue Inventor before: Chen Bo |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |