CN107485405B - 一种利用参考模块测量骨密度的装置 - Google Patents

一种利用参考模块测量骨密度的装置 Download PDF

Info

Publication number
CN107485405B
CN107485405B CN201710712476.1A CN201710712476A CN107485405B CN 107485405 B CN107485405 B CN 107485405B CN 201710712476 A CN201710712476 A CN 201710712476A CN 107485405 B CN107485405 B CN 107485405B
Authority
CN
China
Prior art keywords
bone
image
ray
value
bone density
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
Application number
CN201710712476.1A
Other languages
English (en)
Other versions
CN107485405A (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.)
Zhejiang Kangyuan Medical Device Co ltd
Original Assignee
Zhejiang Kangyuan Medical Device Co ltd
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 Zhejiang Kangyuan Medical Device Co ltd filed Critical Zhejiang Kangyuan Medical Device Co ltd
Priority to CN201710712476.1A priority Critical patent/CN107485405B/zh
Publication of CN107485405A publication Critical patent/CN107485405A/zh
Application granted granted Critical
Publication of CN107485405B publication Critical patent/CN107485405B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus 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/505Apparatus 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明公开一种利用参考模块测量骨密度的装置,包括,数字X‑射线影像设备:用于采集被测骨骼投影的图像;图像预处理装置:将包含有被测骨骼的图像进行预处理;图像分割处理装置:将预处理后的被测骨骼的图像进行分割处理;骨骼衰减因子计算装置,依据X‑射线的衰减特性计算得到被照部位的骨骼衰减因子;骨密度等效分析装置:把上述m(x,y)值与一个已知等效骨密度值的参考模块所成的影像中具有相同的m(x,y)值的区域相对应,对应的参考模块区域的等效骨密度值即为被测骨骼部位的骨密度值。本发明与传统的DXA方法相比较,具有测量精度高、便于操作等优点。

Description

一种利用参考模块测量骨密度的装置
技术领域
本发明涉及医疗器械技术领域,具体涉及应用X-射线测量技术定量测量骨密度的装置。
背景技术
目前市场上已有多种多样的骨密度检测方法及产品【1】,如双能量X-射线吸收法(dual X-ray absorptiometry,DXA),定量CT法(quantitative computed tomography,QCT)、定量磁共振法(quantitative magnetic resonance imaging,QMRI)和定量超声法(quantitative ultrasound,QUS),其中双能量X-射线吸收法(DXA)是WHO建议的用于临床诊断骨质疏松的金标准方法,但是其存在着设备相对昂贵、维修成本高、需要经过专门培训的专业人士操作的等诸多弱点。在这种情况下,人工智能与传统检测技术的结合,使得其他新的检测技术被不断地被人们所关注、开发。
发明内容
本发明所要解决的技术问题是,提供一种利用参考模块测量骨密度的装置,降低成本,便于操作。
为解决上述技术问题,本发明采用如下技术方案:一种利用参考模块测量骨密度的装置,包括,
数字X-射线影像设备:用于采集被测骨骼投影的图像;
图像预处理装置:将包含有被测骨骼的图像进行预处理;
图像分割处理装置:将预处理后的被测骨骼的图像进行分割处理,分割出骨骼和骨骼周围的软组织区域,并分别计算出这些区域的几何参数,其中的几何参数包括大小、相对位置,以及这些区域中所包含的图像灰度值I(x,y);
骨骼衰减因子计算装置,将X-射线曝光参数以及上述描述的几何参数作为输入变量,依据X-射线的衰减特性,透过被照骨骼区域后,X-射线探测器所接收到的图像灰度信号I(x,y)表示为,
Figure BDA0001383015550000021
m(x,y)=ln[I0(x,y)/I(x,y)]+ms-ma
其中,m(x,y)是被照部位的骨骼衰减因子,I0(x,y)是X-射线探测器在没有经过任何物体衰减所接受到的图像信号强度,S(x,y)是射线锥形束散射修正因子,ms=lnS(x,y)是射线散射修正值,m0则是软组织吸收的修正值;
骨密度等效分析装置:把上述m(x,y)值与一个已知等效骨密度值的参考模块所成的影像中具有相同的m(x,y)值的区域相对应,对应的参考模块区域的等效骨密度值即为被测骨骼部位的骨密度值。
优选的,图像预处理包括对X-射线探测器的每个像素单元的的信号降噪处理、X-射线探测器灵敏度以及X-射线分布的不均匀性矫正,以及对照射投影几何失真进行有效的修正。
优选的,所述骨骼衰减因子计算装置具有深度学习模块,先计算出软组织的吸收衰减以及射线锥形束射线散射对图像骨骼区域的灰度信号的修正量,最后根据修正量得出有效骨骼衰减因子。
优选的,所述参考模块由一组不同厚度的铝片和丙烯酸树脂、以及聚乙烯组合成梯形状。
优选的,还包括骨质疏松症诊断参数计算装置,依据人群特定部位的骨密度值数据库,得出高出或低于健康年轻人的标准方差数作为骨质疏松症诊断的参数。
本发明通过普通的数字X-射线影像设备(DR)来采集人体骨骼投影的图像,采用智能计算与X-射线放射吸收原理相结合的方法来对包含该骨骼的图像进行处理,计算出该骨骼中感兴趣的骨骼区域的骨密度值。由于与传统的DXA方法相比较,具有测量精度高、便于操作等优点,因此有望成为最具有发展前景的诊断骨质疏松症的技术之一。
附图说明
图1为本发明中参考模块的设计示意图。
具体实施方式
本发明通过普通的数字X-射线影像设备(DR)来采集人体骨骼投影的图像,采用智能计算与X-射线放射吸收原理相结合的方法来对包含该骨骼的图像进行处理,计算出该骨骼中感兴趣的骨骼区域的骨密度值。以下对本发明的具体实施方式做出具体说明。
一种利用参考模块测量骨密度的装置,包括,
数字X-射线影像设备:用于采集被测骨骼投影的图像;
图像预处理装置:将包含有被测骨骼的图像进行预处理;
图像分割处理装置:将预处理后的被测骨骼的图像进行分割处理,分割出骨骼和骨骼周围的软组织区域,并分别计算出这些区域的几何参数,其中的几何参数包括大小、相对位置,以及这些区域中所包含的图像灰度值I(x,y);
骨骼衰减因子计算装置,将X-射线曝光参数以及上述描述的几何参数作为输入变量,依据X-射线的衰减特性,透过被照骨骼区域后,X-射线探测器所接收到的图像灰度信号I(x,y)表示为,
Figure BDA0001383015550000041
m(x,y)=ln[I0(x,y)/I(x,y)]+ms-m0
其中,m(x,y)是被照部位的骨骼衰减因子,I0(x,y)是X-射线探测器在没有经过任何物体衰减所接受到的图像信号强度,S(x,y)是射线锥形束散射修正因子,ms=lnS(x,y)是射线散射修正值,m0则是软组织吸收的修正值;
骨密度等效分析装置:把上述m(x,y)值与一个已知等效骨密度值的参考模块所成的影像中具有相同的m(x,y)值的区域相对应,对应的参考模块区域的等效骨密度值即为被测骨骼部位的骨密度值。
利用上述装置对骨密度进行测量的测量方法说明如下,智能放射吸收测量法是以普通数字X射线影像设备为图像采集平台,首先采集被测量的骨骼投影的图像,接下来将该图像输入到图像处理工作站进行预处理,接下来将该图像输入到图像处理工作站进行预处理。这预处理可包括对平板探测器的每个像素单元的灵敏度,以及平板探测器表面X-射线分布的不均匀性进行矫正,同时也对射线投影造成的几何失真进行有效的修正。然后采用常用的图像分割算法,比如GVF Snake自动分割方法,将预处理后的骨骼的图像进行分割处理,有效地再投影图像中分割出感兴趣的骨骼和其周围的软组织区域,以及这些区域中所包含的图像灰度值I0,及其分布。然后,依据X-射线的衰减特性计算得出骨骼衰减因子。最后,我们采用“参考模块法”,把上述m(x,y)值与一个已知等效骨密度值的参考模块所成的影像中具有相同的m(x,y)值的区域相对应,对应的参考模块区域的等效骨密度值为我们被测骨骼部位的骨密度值。
一种更全面的方法可以通过深度学习的方法,骨骼衰减因子计算装置具有深度学习模块,先计算出软组织的吸收衰减以及射线锥形束射线散射对图像骨骼区域的灰度信号的修正量,最后得出有效骨骼衰减因子m(x,y)。
在测量骨密度值时,是对某一个特定的部位区域来进行测量的。而这些区域则是需要依据骨骼特定的解剖学部位来定义的,比方说在测量手臂时,采用尺骨茎突尖端、桡骨远端关节软骨下终板、尺骨桡侧与桡骨尺侧交界处作为参考点指标来定位被测区域的骨密度值。
骨密度测量对测量系统精度要求特别高,因此,必须对测量系统采用特殊的矫正方法。在实际测量过程中,操作人员需要对所使用的数字X-射线影像设备进行必要的矫正。这些矫正包含多方面的内容,包括X-射线矫正、系统矫正、图像矫正。采用矫正确认后处理测量系统的所测得的结果的稳定性和一致性,如果发现测量结果有偏差,算法就会对测量系统进行相应的自动矫正。
在本实施例中,图像分割处理装置可以采用GVF Snake图像分割算法,将预处理后的被测骨骼的图像进行分割处理。当然,本领域技术人员可以理解是,成熟的图像分割算法目前有很多种,具体使用包括但并不限于GVF Snake图像分割算法。
参考模块包含有一系列不同等效骨密度值的规则区域。采用参考模块法来进行定量分析计算被测骨骼的骨密度值。这里特殊加工的“参考模块”,如图1所示。图中A:模拟前臂的部分,黑色部分是铝材料用来模拟骨骼,灰色部分是丙烯酸树脂用来模拟皮肤软组织。图中B:不同厚度铝片所组成的铝梯,用于模拟不同面密度的骨骼。图中C:不同厚度的聚乙烯(Polyethyleneh或Acrylic)所组成的阶梯,用于模拟不同厚度的脂肪组织,肌肉组织。它包含可用于跟踪原始入射射线特性的由一组不同厚度的铝片和丙烯酸树脂、以及聚乙烯所组合成的梯形状。这里可分别采用三种材料、纯铝、丙烯酸树脂、聚乙烯来模拟标准骨骼、脂肪、肌肉的射线吸收衰减、射线的散射和软组织的吸收,当然其他类似等效的材料也可以用于分析模块中,这些材料的射线的等效衰减特性,包括等效密度值都已经被确定。
进一步地,该装置还包括骨质疏松症诊断参数计算装置,可以依据人群特定部位的骨密度值数据库,得出高出(+)或低于(-)健康年轻人的标准方差数作为骨质疏松症诊断的参数,即T-值。
除上述优选实施例外,本发明还有其他的实施方式,本领域技术人员可以根据本发明作出各种改变和变形,只要不脱离本发明的精神,均应属于本发明权利要求书中所定义的范围。

Claims (4)

1.一种利用参考模块测量骨密度的装置,其特征在于包括,
数字X-射线影像设备:用于采集被测骨骼投影的图像;
图像预处理装置:将包含有被测骨骼的图像进行预处理;
图像分割处理装置:将预处理后的被测骨骼的图像进行分割处理,分割出骨骼区域和骨骼周围的软组织区域,并分别计算出这些区域的几何参数,其中的几何参数包括大小、相对位置,以及这些区域中所包含的图像灰度信号强度I(x,y);
骨骼衰减因子计算装置,将X-射线曝光参数以及上述描述的几何参数作为输入变量,依据X-射线的衰减特性,透过被照射部位后,X-射线探测器所接收到的图像灰度信号强度I(x,y)表示为,
Figure FDA0002830796410000011
m(x,y)=ln[I0(x,y)/I(x,y)]+ms-m0
其中,m(x,y)是被照部位的骨骼衰减因子,I0(x,y)是X-射线探测器在没有经过任何物体衰减所接受到的图像灰度信号强度,S(x,y)是射线锥形束散射修正因子,ms=lnS(x,y)是射线散射修正值,m0则是软组织吸收的修正值,所述骨骼衰减因子计算装置具有深度学习模块,先计算出软组织的吸收衰减以及射线锥形束射线散射对图像骨骼区域的灰度信号的修正量,最后根据修正量得出有效骨骼衰减因子;
骨密度等效分析装置:把上述m(x,y)值与一个已知等效骨密度值的参考模块所成的影像中具有相同的m(x,y)值的区域相对应,对应的参考模块区域的等效骨密度值即为被测骨骼部位的骨密度值。
2.根据权利要求1所述的一种利用参考模块测量骨密度的装置,其特征在于,图像预处理包括对X-射线探测器的每个像素单元的的信号降噪处理、X-射线探测器灵敏度以及X-射线分布的不均匀性矫正,以及对照射投影几何失真进行有效的修正。
3.根据权利要求1所述的一种利用参考模块测量骨密度的装置,其特征在于,所述参考模块由一组不同厚度的铝片和丙烯酸树脂、以及聚乙烯组合成梯形状。
4.根据权利要求1所述的一种利用参考模块测量骨密度的装置,其特征在于,还包括骨质疏松症诊断参数计算装置,依据人群特定部位的骨密度值数据库,得出高出或低于健康年轻人的标准方差数作为骨质疏松症诊断的参数。
CN201710712476.1A 2017-08-18 2017-08-18 一种利用参考模块测量骨密度的装置 Active CN107485405B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710712476.1A CN107485405B (zh) 2017-08-18 2017-08-18 一种利用参考模块测量骨密度的装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710712476.1A CN107485405B (zh) 2017-08-18 2017-08-18 一种利用参考模块测量骨密度的装置

Publications (2)

Publication Number Publication Date
CN107485405A CN107485405A (zh) 2017-12-19
CN107485405B true CN107485405B (zh) 2021-02-19

Family

ID=60646723

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710712476.1A Active CN107485405B (zh) 2017-08-18 2017-08-18 一种利用参考模块测量骨密度的装置

Country Status (1)

Country Link
CN (1) CN107485405B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108682001A (zh) * 2018-03-30 2018-10-19 深圳市深图医学影像设备有限公司 骨密度的测量方法及装置
CN109060597A (zh) * 2018-08-28 2018-12-21 中国工程物理研究院激光聚变研究中心 一种金属冲击破碎过程密度诊断方法及系统
JP6744614B1 (ja) * 2018-09-10 2020-08-19 京セラ株式会社 推定装置、推定システム及び推定プログラム
CN113491526B (zh) * 2020-04-07 2023-12-05 辽宁开普医疗系统有限公司 一种基于dr系统的骨密度校正和测量的方法
CN112089429B (zh) * 2020-09-18 2023-09-26 重庆邮电大学 一种基于深度学习算法骨密度仪诊断系统

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06261894A (ja) * 1993-01-18 1994-09-20 Hitachi Ltd 骨塩定量方法
CN102217944A (zh) * 2010-04-14 2011-10-19 深圳市艾克瑞电气有限公司 应用数字x线图像进行骨密度测量的方法
JP2012165983A (ja) * 2011-02-17 2012-09-06 Shimadzu Corp X線撮影装置
CN102805637A (zh) * 2011-05-31 2012-12-05 赵建国 一种三维骨密度的计算方法和测量仪
CN103442638A (zh) * 2011-03-18 2013-12-11 日立阿洛卡医疗株式会社 骨密度测量装置
CN103892856A (zh) * 2012-12-31 2014-07-02 深圳先进技术研究院 一种获取人体骨密度值的方法及系统
CN104873213A (zh) * 2015-04-13 2015-09-02 浙江康源医疗器械有限公司 基于x-射线数字影像的骨骼病变评估方法和装置
CN105388170A (zh) * 2015-12-17 2016-03-09 中国科学技术大学 一种骨密度的测量方法及装置
CN105575455A (zh) * 2015-12-14 2016-05-11 天津三英精密仪器有限公司 一种x射线衰减器设计方法和应用以及利用该方法设计的带有衰减器的ct装置
CN107007294A (zh) * 2016-01-28 2017-08-04 株式会社日立制作所 X射线图像诊断装置以及骨密度测量方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8000512B2 (en) * 2007-10-03 2011-08-16 General Electric Company Slit collimator scatter correction
US7724875B2 (en) * 2007-10-19 2010-05-25 General Electric Company Image guided acquisition of quantitative dual energy data
CN103690187A (zh) * 2013-08-13 2014-04-02 清华大学 基于x射线测量口腔骨密度的设备及方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06261894A (ja) * 1993-01-18 1994-09-20 Hitachi Ltd 骨塩定量方法
CN102217944A (zh) * 2010-04-14 2011-10-19 深圳市艾克瑞电气有限公司 应用数字x线图像进行骨密度测量的方法
JP2012165983A (ja) * 2011-02-17 2012-09-06 Shimadzu Corp X線撮影装置
CN103442638A (zh) * 2011-03-18 2013-12-11 日立阿洛卡医疗株式会社 骨密度测量装置
CN102805637A (zh) * 2011-05-31 2012-12-05 赵建国 一种三维骨密度的计算方法和测量仪
CN103892856A (zh) * 2012-12-31 2014-07-02 深圳先进技术研究院 一种获取人体骨密度值的方法及系统
CN104873213A (zh) * 2015-04-13 2015-09-02 浙江康源医疗器械有限公司 基于x-射线数字影像的骨骼病变评估方法和装置
CN105575455A (zh) * 2015-12-14 2016-05-11 天津三英精密仪器有限公司 一种x射线衰减器设计方法和应用以及利用该方法设计的带有衰减器的ct装置
CN105388170A (zh) * 2015-12-17 2016-03-09 中国科学技术大学 一种骨密度的测量方法及装置
CN107007294A (zh) * 2016-01-28 2017-08-04 株式会社日立制作所 X射线图像诊断装置以及骨密度测量方法

Also Published As

Publication number Publication date
CN107485405A (zh) 2017-12-19

Similar Documents

Publication Publication Date Title
CN107485405B (zh) 一种利用参考模块测量骨密度的装置
EP1324695B1 (en) Methods and devices for analysis of x-ray images
US6173038B1 (en) Plain x-ray bone densitometry apparatus and method
US20070086560A1 (en) Scatter correction
Lefebvre et al. Automatic detection of the boundary of the calcaneus from ultrasound parametric images using an active contour model; clinical assessment
US20070047794A1 (en) Methods and devices for analysis of x-ray images
CN107115119B (zh) Pet图像衰减系数的获取方法、衰减校正的方法及系统
CA2580726A1 (en) System and method of predicting future fractures
EP1530775A2 (en) Methods and devices for analysis of x-ray images
US6570955B1 (en) Digital x-ray material testing and bone densitometry apparatus and method
EP1605824A2 (en) Methods for the compensation of imaging technique in the processing of radiographic images
US11288848B2 (en) Three-dimensional ultrasound image display method
GB2512876A (en) Methods and apparatus for quantifying inflammation
CN104873213B (zh) 基于x‑射线数字影像的骨骼病变评估方法和装置
Tsujimoto et al. Automated segmentation and detection of increased uptake regions in bone scintigraphy using SPECT/CT images
Korstjens et al. Reliability of an image analysis system for quantifying the radiographic trabecular pattern
Pavan et al. Quality and dose optimization in hand computed radiography
Wilkinson et al. Digital image analysis protocol for determining the radiocarpal joint space in the rheumatoid arthritic wrist
Berger et al. A software program to measure the three-dimensional length of the spine from radiographic images: validation and reliability assessment for adolescent idiopathic scoliosis
US20220296192A1 (en) Bone trabeculae index for x-ray dark-field radiography
KR102578943B1 (ko) 엑스레이(X-ray) 촬영된 고관절 방사선 이미지를 이용한 기계학습 기반 골밀도 도출방법
CN114782375B (zh) 骨质密度测量方法、装置和设备
RU2684567C2 (ru) Способ реконструктивного дозиметрического контроля в протонной терапии сканирующим пучком
Díez et al. Analysis and evaluation of periodic physiological organ motion in radiotherapy treatments
Kuczynski Characterizing the structure-function relationship of hand osteoarthritis using dynamic and high resolution CT imaging

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