CN114782375A - 骨质密度测量方法、装置和设备 - Google Patents
骨质密度测量方法、装置和设备 Download PDFInfo
- Publication number
- CN114782375A CN114782375A CN202210448972.1A CN202210448972A CN114782375A CN 114782375 A CN114782375 A CN 114782375A CN 202210448972 A CN202210448972 A CN 202210448972A CN 114782375 A CN114782375 A CN 114782375A
- Authority
- CN
- China
- Prior art keywords
- value
- bone
- value information
- region
- 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.)
- Granted
Links
- 230000037182 bone density Effects 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 46
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 69
- 238000002591 computed tomography Methods 0.000 claims abstract description 63
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 22
- 239000011707 mineral Substances 0.000 claims abstract description 22
- 238000001739 density measurement Methods 0.000 claims abstract description 19
- 210000003205 muscle Anatomy 0.000 claims description 36
- 238000005192 partition Methods 0.000 claims 1
- 238000000638 solvent extraction Methods 0.000 claims 1
- 230000035945 sensitivity Effects 0.000 abstract description 8
- 238000003745 diagnosis Methods 0.000 abstract description 7
- 238000003908 quality control method Methods 0.000 abstract description 5
- 230000005855 radiation Effects 0.000 abstract description 4
- 210000001519 tissue Anatomy 0.000 abstract description 4
- 230000000875 corresponding effect Effects 0.000 description 35
- 208000001132 Osteoporosis Diseases 0.000 description 9
- 238000009547 dual-energy X-ray absorptiometry Methods 0.000 description 9
- 206010017076 Fracture Diseases 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 208000010392 Bone Fractures Diseases 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 208000001164 Osteoporotic Fractures Diseases 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000013399 early diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000011272 standard treatment Methods 0.000 description 2
- 208000008558 Osteophyte Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000326 densiometry Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 238000012795 verification 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种骨质密度测量方法、装置和设备，属于骨密度测量技术领域，该方法、装置和设备通过采用目标患者的常规CT图像，对其进行分析划分，得到不同区域，通过不同区域的CT值和预先构建的骨密度与CT值对应规则，计算得到骨密度值，方便、准确。本申请无需特定体模质控，采用目标患者自身不同组织的CT值和预设计算公式进行骨密度的计算，可以最大程度将用于其他诊断目的的CT扫描数据利用起来，提高效率，节约成本的，减少被试者的辐射量；同时，本发明的结果为三维体积骨密度，敏感性及特异性高。
Description
技术领域
本发明属于骨密度测量技术领域，具体涉及一种骨质密度测量方法、装置和设备。
背景技术
骨质疏松骨折是骨质疏松症的重要并发症，除引起持续疼痛、导致活动能力下降并严重影响患者生活质量外，还具有较高的致残率及病死率。早发现骨质疏松症并进行早期干预，是避免骨质疏松骨折，提高中老年人生活质量的重要途径。但是，目前仅有约9％-14.6％的骨质疏松患者可以被诊断识别，更遑论早期干预。因此，准确、便利的骨密度定量测定极大的有利于骨质疏松诊断或骨密度减少的诊断，做到早期诊断，及时规范治疗，减少脆性骨折的发生，不仅能够极大减轻患者痛苦，还能极大地降低卫生经济负担，具有重要的科学意义。
相关技术中，通常通过测量骨密度(Bone Mass Density，BMD)来分辨骨质疏松症。临床的骨密度测量工具主要包括双能骨密度测量仪(Dual energy X-rayabsorptiometry，DXA)和定量CT(Quantitative computed tomography，QCT)。其中，DXA是利用骨组织对放射物质的吸收测量相对的骨密度减低的程度，该方法基于二维图像完成，容易受到重叠影的影响，且得到的并非绝对骨密度，易受几何形状及扫描野内骨赘、大血管钙化斑块和骨折等影响，不同型号的DXA骨密度仪之间的检测结果会有较大差异，互认度差，测量的BMD与骨折风险也并不具有良好的相关性，临床诊断敏感性和特异性差。而目前国内应用较广泛的Mindways-QCT检测需要每一台机器安装指定的测量软件，同时需要特定定体模进行质控，多数情况下患者要专门序列进行扫描，受到辐射；有成本高、效率较低的缺陷，不利于广泛推广。
因此，如何便捷、准确测量得到骨密度，成为现有技术中亟待解决的技术问题。
发明内容
有鉴于此，本发明的目的在于提供一种骨质密度测量方法、装置和设备方法、装置和设备，以克服目前诊断敏感性和特异性差、准确性低的技术问题。
为实现以上目的，本发明采用如下技术方案：
一方面，一种骨质密度测量方法，包括：
获取目标患者的CT扫描图像；
按照预设规则将所述CT扫描图像划分为不同的区域；
获取每个区域对应的灰度值信息；
根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息；
根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度。
可选的，所述按照预设规则将所述CT扫描图像划分为不同的区域，包括：
识别所述CT扫描图像的不同边界，将所述CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
可选的，还包括：
通过DICOM标签，获取所述CT扫描图像的参数；所述CT扫描图像的参数，包括：斜率和截距值；
所述根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息，包括：
基于所述预设CT值计算规则，根据每个区域中的所述斜率、所述截距值和所述灰度值信息，计算每个所述区域的CT值信息。
可选的，所述根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度，包括：
获取所述松质骨区域的CT值信息中的松质骨CT值、所述脂肪区域的CT值信息中的脂肪CT值和所述肌肉区域的CT值信息中的肌肉CT值；
分别将所述松质骨CT值、所述脂肪CT值和所述肌肉CT值输入至所述预设骨密度与CT值对应规则中，得到所述目标患者的骨密度。
可选的，所述预设骨密度与CT值对应规则，包括：
骨密度＝33.61328+0.81786*VerPeak-0.16364*MusPeak+0.17357*FatPeak；
其中，VerPeak为所述松质骨CT值，MusPeak为所述肌肉CT值，FatPeak为所述脂肪CT值。
可选的，所述获取目标患者的CT扫描图像，包括：
打开目标患者的以DICOM格式存储的所述CT扫描图像。
又一方面，一种骨密度测量装置，包括：
第一获取模块，用于获取目标患者的CT扫描图像；
划分模块，用于按照预设规则将所述CT扫描图像划分为不同的区域；
第二获取模块，用于获取每个区域对应的灰度值信息；
第一计算模块，用于根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息；
第二计算模块，用于根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度。
可选的，所述划分模块，具体用于识别所述CT扫描图像的不同边界，将所述CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
可选的，所述第二计算模块，具体用于获取所述松质骨区域的CT值信息中的松质骨CT值、所述脂肪区域的CT值信息中的脂肪CT值和所述肌肉区域的CT值信息中的肌肉CT值；分别将所述松质骨CT值、所述脂肪CT值和所述肌肉CT值输入至所述预设骨密度与CT值对应规则中，得到所述目标患者的骨密度。
又一方面，一种骨质密度测量设备，包括处理器和存储器，所述处理器与存储器相连：
其中，所述处理器，用于调用并执行所述存储器中存储的程序；
所述存储器，用于存储所述程序，所述程序至少用于执行上述任一所述的骨质密度测量方法。
本发明实施例提供的骨质密度测量方法、装置和设备，通过采用目标患者的常规CT图像，对其进行分析划分，得到不同区域，通过不同区域的CT值和预先构建的骨密度与CT值对应规则，计算得到骨密度值，方便、准确。本申请无需特定体模质控，采用目标患者自身不同组织的CT值和预设计算公式进行骨密度的计算，可以最大程度将用于其他诊断目的的CT扫描数据利用起来，提高效率，节约成本的，减少被试者的辐射量；同时，本发明的结果为三维体积骨密度，敏感性及特异性高。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种骨密度测量方法的流程示意图；
图2是本发明实施例提供的一种骨质密度测量装置的结构示意图；
图3是本发明实施例提供的一种骨质密度测量设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将对本发明的技术方案进行详细的描述。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所得到的所有其它实施方式，都属于本发明所保护的范围。
据估计我国50岁以上中老年人即将达到4.69亿，到2050年时将高达5.71亿。对于中老年人来说，骨质疏松症是众多身体隐患之一。骨质疏松骨折是骨质疏松症的重要并发症，除引起持续疼痛、导致活动能力下降并严重影响患者生活质量外，还具有较高的致残率及病死率。准确、便利的骨密度定量测定极大的有利于骨质疏松诊断或骨密度减少的诊断，做到早期诊断，及时规范治疗，减少脆性骨折的发生，不仅能够极大减轻患者痛苦，还能极大地降低卫生经济负担，具有重要的科学意义。
目前的常规应用于临床的骨密度测量工具主要包括双能骨密度测量仪(Dualenergy X-ray absorptiometry，DXA)和定量CT(Quantitative computed tomography，QCT)。DXA基于二维图像完成骨密度测量，容易受到重叠影的影响，且得到的并非绝对骨密度；另外，DXA基本特定人群数据库的比值，临床诊断敏感性和特异性差。
目前用于骨密度的测量的QCT分为两类，一类是采用专用体模及专用的软件对扫描得到的骨的CT值进行标化得到体积骨密度,需要安装有专门软件的CT机器进行专门序列的扫描，该方法收到体模及是否绑定软件的限制。一类是不需要体模的定量CT，目前是通过计算二维公式将二维的CT值转化而来，尚处于科研阶段，且准确度差，无国内数据库。
因此，如何便捷、准确测量得到骨密度，成为现有技术中亟待解决的技术问题。基于此，本发明实施例提供一种骨质密度测量方法、装置和设备，以实现无体模、方便、准确地骨密度测量。
实施例：
图1是本发明实施例提供的一种骨密度测量方法的流程示意图，请参阅图1，本实施例可以包括以下步骤：
S1、获取目标患者的CT扫描图像。
在具体的实现过程中，可以定义任意一个患者为目标患者，对目标患者进行骨密度的测量。目标患者可以在普通CT上进行检查，得到CT扫描图像。在骨密度测量的过程中，可以直接导入目标患者在普通CT上检查得到的CT扫描图像，并对CT扫描图像进行分析。
在一些实施例中，获取目标患者的CT扫描图像，包括：打开目标患者的以DICOM(Digital Imaging and Communications in Medicine，医学数字成像和通信,是医学图像和相关信息的国际标准)格式存储的CT扫描图像。
例如，骨密度测量时，可以配置DICOM格式打开浏览功能，获取DICOM格式的CT扫描图像，从而进行浏览和读取。
S2、按照预设规则将CT扫描图像划分为不同的区域。
在一些实施例中，按照预设规则将CT扫描图像划分为不同的区域，包括：识别CT扫描图像的不同边界，将CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
在获取到目标患者的CT扫描图像后，可以手动或者自动，或部分手动部分自动，对CT扫描图像进行划分。例如，可以根据选择的不同算法，对CT扫描图像进行自动勾画，如，根据不同区域的特征设定阈值范围，从而划分CT空间分布，可以根据函数求导法基于变化率识别不同区域的边界等。也可以由用户根据工作经验，识别到脂肪区域、肌肉区域和椎体的松质骨区域的边界，进行手动勾画，从而区分出患者的CT扫描图像中的脂肪区域、肌肉区域和椎体的松质骨区域。
S3、获取每个区域对应的灰度值信息。
在划分到不同区域后，获取每个区域的灰度值信息。其中，每个区域的灰度值信息包括该区域内每个像素点的灰度值。每个像素点的灰度值获取为现有技术，可以根据成熟的计算方法进行获取，本申请中不做赘述，请参阅现有技术。
S4、根据灰度值信息、CT扫描图像的参数以及预设CT值计算规则，计算与每个区域的灰度值信息对应的CT值信息。
在一些实施例中，还包括：通过DICOM标签，获取CT扫描图像的参数；CT扫描图像的参数，包括：斜率和截距值；根据灰度值信息、CT扫描图像的参数以及预设CT值计算规则，计算与每个区域的灰度值信息对应的CT值信息，包括：基于预设CT值计算规则，根据每个区域中的斜率、截距值和灰度值信息，计算每个区域的CT值信息。
其中，DICOM标签自动记录与采集决定的CT扫描图像的相关参数，例如，斜率(slope)和截距(intercept)值。可以通过DICOM标签读取到CT扫描图像的slope和intercept值。结合灰度值，利用公式可以计算出与之对应的CT值。
每个区域的CT值信息包括多个CT值，每个像素点可对应一个CT值。
在一些实施例中，标准化反应线性系数CT值的计算公式可以为：
CT＝slope*greyscale+intercept，
即：CT值＝斜率*灰度值+截距值公式(1)
根据上述公式，将获取到的灰度与CT值关系坐标里每个点两值互相转换的斜率、截距值，以及每个点的灰度值输入到公式(1)中，可以计算得到每个区域的CT值信息。
S5、根据预设骨密度与CT值对应规则和CT值信息，确定目标患者的骨密度。
在得到每个区域的CT值信息(即，得到每个区域内每个像素点对应的CT值)后，可以根据目标患者身体不同区域的CT值信息，以及，预设骨密度与CT值对应规则，来计算与CT值对应的骨密度，方便、快捷。
在一些实施例中，根据预设骨密度与CT值对应规则和CT值信息，确定目标患者的骨密度，包括：获取松质骨区域的CT值信息中的松质骨CT值、脂肪区域的CT值信息中的脂肪CT值和肌肉区域的CT值信息中的肌肉CT值；分别将松质骨CT值、脂肪CT值和肌肉CT值输入至预设骨密度与CT值对应规则中，得到目标患者的骨密度。
在一些实施例中，预设骨密度与CT值对应规则，包括：
骨密度＝33.61328+0.81786*VerPeak-0.16364*MusPeak+0.17357*FatPeak
公式(2)；
其中，VerPeak为松质骨CT值，MusPeak为肌肉CT值，FatPeak为脂肪CT值。
为了提高骨密度的计算准确性，在确定松质骨CT值时，可以将松质骨区域内最多像素点对应的CT值作为松质骨CT值，或者，采取松质骨区域内像素点对应的CT值的峰值作为松质骨CT值；可以将肌肉区域内最多像素点对应的CT值作为肌肉CT值，或者，采取肌肉区域内像素点对应的CT值的峰值作为肌肉CT值；可以将脂肪区域内最多像素点对应的CT值作为脂肪CT值，或者，采取脂肪区域内像素点对应的CT值的峰值作为脂肪CT值。
在本申请中，由于从脂肪区域、肌肉区域和椎体的松质骨区域三大不同维度进行CT值的获取，并通过三维度CT值进行骨密度的计算，使得最终计算骨密度为三维体积骨密度，敏感性及特异性高。
采用脂肪区域、肌肉区域和椎体的松质骨区域的数据，可以得到目标患者的全椎体骨密度测量结果。
本申请实施例中，在制定预设骨密度与CT值对应规则，即制定公式(2)时，通过联合694例病人的DXA骨密度数据与对应的CT值参数进行计算获得。例如，通过识别不同区域提取的CT值，和DXA骨密度数据，进行骨密度计算公式的拟合。采用机器学习的方法，在线性转换的基础上，抓取CT值和DXA骨密度的关系，在计算所得骨密度与DXA骨密度差异最小的条件下，自动获取最佳的特征组合和CT值计算骨密度公式。
在获得预设骨密度的计算公式后，对计算公式进行验证和修正，在验证时，采用173例病人进行验证，最终基于该计算公式得到的三维骨密度值与MINDWAN-QCT参照骨密度的标准差为12.143，计算准确度较高。
本发明实施例提供的无体模全椎体骨质密度测量方法，通过采用目标患者的常规CT图像，对其进行分析划分，得到不同区域，通过不同区域的CT值和预先构建的骨密度与CT值对应规则，计算得到骨密度值，方便、准确。本申请无需特定体模质控，采用目标患者自身不同组织的CT值和预设计算公式进行骨密度的计算，可以最大程度将用于其他诊断目的的CT扫描数据利用起来，提高效率，节约成本的，减少被试者的辐射量；同时，本发明的结果为三维体积骨密度，敏感性及特异性高。
基于一个总的发明构思，本发明实施例还提供一种骨质密度测量装置。
本发明还提供了一种骨质密度测量装置，用于实现上述方法实施例。图2是本发明实施例提供的一种骨质密度测量装置的结构示意图。如图2所示，本发明提供的装置可以包括以下结构：
第一获取模块21，用于获取目标患者的CT扫描图像；
划分模块22，用于按照预设规则将CT扫描图像划分为不同的区域；
第二获取模块23，用于获取每个区域对应的灰度值信息；
第一计算模块24，用于根据灰度值信息、CT扫描图像的参数以及预设CT值计算规则，计算与每个区域的灰度值信息对应的CT值信息；
第二计算模块25，用于根据预设骨密度与CT值对应规则和CT值信息，确定目标患者的骨密度。
可选的，划分模块，具体用于识别CT扫描图像的不同边界，将CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
可选的，第一获取模块，还用于通过DICOM标签，获取CT扫描图像的参数；CT扫描图像的参数，包括：斜率和截距值；第一计算模块，具体用于基于预设CT值计算规则，根据每个区域中的斜率、截距值和灰度值信息，计算每个区域的CT值信息。
可选的，第二计算模块，具体用于获取松质骨区域的CT值信息中的松质骨CT值、脂肪区域的CT值信息中的脂肪CT值和肌肉区域的CT值信息中的肌肉CT值；分别将松质骨CT值、脂肪CT值和肌肉CT值输入至预设骨密度与CT值对应规则中，得到目标患者的骨密度。
可选的，预设骨密度与CT值对应规则，包括：
骨密度＝33.61328+0.81786*VerPeak-0.16364*MusPeak+0.17357*FatPeak；
其中，VerPeak为松质骨CT值，MusPeak为肌肉CT值，FatPeak为脂肪CT值。
关于上述实施例中的装置，其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述，此处将不做详细阐述说明。
本发明实施例提供的无体模全椎体骨质密度测量装置，通过采用目标患者的常规CT图像，对其进行分析划分，得到不同区域，通过不同区域的CT值和预先构建的骨密度与CT值对应规则，计算得到骨密度值，方便、准确。本申请无需特定体模质控，采用目标患者自身不同组织的CT值和预设计算公式进行骨密度的计算，可以最大程度将用于其他诊断目的的CT扫描数据利用起来，提高效率，节约成本的，减少被试者的辐射量；同时，本发明的结果为三维体积骨密度，敏感性及特异性高。
基于一个总的发明构思，本发明实施例还提供一种骨质密度测量设备。
本发明还提供了一种骨质密度测量设备，用于实现上述方法实施例。图3是本发明实施例提供的一种骨质密度测量设备的结构示意图。如图3所示，本实施例的骨质密度测量设备包括处理器31和存储器32，处理器31与存储器32相连。其中，处理器31用于调用并执行存储器32中存储的程序；存储器32用于存储程序，程序至少用于执行以上实施例中的骨质密度测量方法。
本申请实施例提供的骨质密度测量设备的具体实施方案可以参考以上任意实施例的骨质密度测量方法的实施方式，此处不再赘述。
可以理解的是，上述各实施例中相同或相似部分可以相互参考，在一些实施例中未详细说明的内容可以参见其他实施例中相同或相似的内容。
需要说明的是，在本发明的描述中，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性。此外，在本发明的描述中，除非另有说明，“多个”的含义是指至少两个。
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本发明的实施例所属技术领域的技术人员所理解。
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(PGA)，现场可编程门阵列(FPGA)等。
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
此外，在本发明各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
上述提到的存储介质可以是只读存储器，磁盘或光盘等。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种骨质密度测量方法，其特征在于，包括：
获取目标患者的CT扫描图像；
按照预设规则将所述CT扫描图像划分为不同的区域；
获取每个区域对应的灰度值信息；
根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息；
根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度。
2.根据权利要求1所述的方法，其特征在于，所述按照预设规则将所述CT扫描图像划分为不同的区域，包括：
识别所述CT扫描图像的不同边界，将所述CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
3.根据权利要求1所述的方法，其特征在于，还包括：
通过DICOM标签，获取所述CT扫描图像的参数；所述CT扫描图像的参数，包括：斜率和截距值；
所述根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息，包括：
基于所述预设CT值计算规则，根据每个区域中的所述斜率、所述截距值和所述灰度值信息，计算每个所述区域的CT值信息。
4.根据权利要求2所述的方法，其特征在于，所述根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度，包括：
获取所述松质骨区域的CT值信息中的松质骨CT值、所述脂肪区域的CT值信息中的脂肪CT值和所述肌肉区域的CT值信息中的肌肉CT值；
分别将所述松质骨CT值、所述脂肪CT值和所述肌肉CT值输入至所述预设骨密度与CT值对应规则中，得到所述目标患者的骨密度。
5.根据权利要求4所述的方法，其特征在于，所述预设骨密度与CT值对应规则，包括：
骨密度＝33.61328+0.81786*VerPeak-0.16364*MusPeak+0.17357*FatPeak；
其中，VerPeak为所述松质骨CT值，MusPeak为所述肌肉CT值，FatPeak为所述脂肪CT值。
6.根据权利要求1所述的方法，其特征在于，所述获取目标患者的CT扫描图像，包括：
打开目标患者的以DICOM格式存储的所述CT扫描图像。
7.一种骨密度测量装置，其特征在于，包括：
第一获取模块，用于获取目标患者的CT扫描图像；
划分模块，用于按照预设规则将所述CT扫描图像划分为不同的区域；
第二获取模块，用于获取每个区域对应的灰度值信息；
第一计算模块，用于根据所述灰度值信息、所述CT扫描图像的参数以及预设CT值计算规则，计算与每个所述区域的所述灰度值信息对应的CT值信息；
第二计算模块，用于根据预设骨密度与CT值对应规则和所述CT值信息，确定所述目标患者的骨密度。
8.根据权利要求7所述的装置，其特征在于，所述划分模块，具体用于识别所述CT扫描图像的不同边界，将所述CT扫描图像划分为脂肪区域、肌肉区域和椎体的松质骨区域。
9.根据权利要求8所述的装置，其特征在于，所述第二计算模块，具体用于获取所述松质骨区域的CT值信息中的松质骨CT值、所述脂肪区域的CT值信息中的脂肪CT值和所述肌肉区域的CT值信息中的肌肉CT值；分别将所述松质骨CT值、所述脂肪CT值和所述肌肉CT值输入至所述预设骨密度与CT值对应规则中，得到所述目标患者的骨密度。
10.一种骨质密度测量设备，其特征在于，包括处理器和存储器，所述处理器与存储器相连：
其中，所述处理器，用于调用并执行所述存储器中存储的程序；
所述存储器，用于存储所述程序，所述程序至少用于执行权利要求1-6任一项所述的骨质密度测量方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210448972.1A CN114782375B (zh) | 2022-04-27 | 2022-04-27 | 骨质密度测量方法、装置和设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210448972.1A CN114782375B (zh) | 2022-04-27 | 2022-04-27 | 骨质密度测量方法、装置和设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114782375A true CN114782375A (zh) | 2022-07-22 |
| CN114782375B CN114782375B (zh) | 2023-06-23 |
Family
ID=82433849
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210448972.1A Active CN114782375B (zh) | 2022-04-27 | 2022-04-27 | 骨质密度测量方法、装置和设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114782375B (zh) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080013813A1 (en) * | 2006-07-12 | 2008-01-17 | General Electric Company | Methods and apparatus for BMD measuring |
| US20150164454A1 (en) * | 2013-12-18 | 2015-06-18 | Siemens Medical Solutions Usa, Inc. | Method and computing unit for measuring and displaying the bone density of a patient |
| CN106447787A (zh) * | 2016-09-18 | 2017-02-22 | 北京理工大学 | 骨骼ct值与弹性模量关系确定方法及其试验加载装置 |
| US20170098315A1 (en) * | 2015-10-06 | 2017-04-06 | Rigaku Corporation | Analyzer, analysis method and analysis program of bone mineral density |
| CN108805950A (zh) * | 2018-05-28 | 2018-11-13 | 牙博士医疗控股集团有限公司 | 口腔牙片信息转换方法及装置 |
| US20190192099A1 (en) * | 2017-12-21 | 2019-06-27 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for medical imaging of intervertebral discs |
| CN110363765A (zh) * | 2019-07-24 | 2019-10-22 | 苏州科技城医院 | 基于ct影像的骨质疏松参数自动测量方法 |
| CN110660044A (zh) * | 2019-08-30 | 2020-01-07 | 博志生物科技(深圳)有限公司 | 骨组织结构形态异常的快速检测方法及电子装置 |
| CN110840473A (zh) * | 2019-11-22 | 2020-02-28 | 江研伟 | 一种基于CT薄层扫描Hu值的骨密度测量系统 |
| CN111212680A (zh) * | 2017-08-23 | 2020-05-29 | 西门子医疗有限公司 | 用于提供适合于在对患者的辐照的规划中使用的结果数据的方法 |
| CN111881910A (zh) * | 2020-07-31 | 2020-11-03 | 杭州依图医疗技术有限公司 | 基于椎骨影像的信息处理方法及计算机可读存储介质 |
| CN113679404A (zh) * | 2021-07-06 | 2021-11-23 | 湖南浚朗科技有限公司 | 一种骨密度测量系统及其测量方法 |
-
2022
- 2022-04-27 CN CN202210448972.1A patent/CN114782375B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080013813A1 (en) * | 2006-07-12 | 2008-01-17 | General Electric Company | Methods and apparatus for BMD measuring |
| US20150164454A1 (en) * | 2013-12-18 | 2015-06-18 | Siemens Medical Solutions Usa, Inc. | Method and computing unit for measuring and displaying the bone density of a patient |
| US20170098315A1 (en) * | 2015-10-06 | 2017-04-06 | Rigaku Corporation | Analyzer, analysis method and analysis program of bone mineral density |
| CN106447787A (zh) * | 2016-09-18 | 2017-02-22 | 北京理工大学 | 骨骼ct值与弹性模量关系确定方法及其试验加载装置 |
| CN111212680A (zh) * | 2017-08-23 | 2020-05-29 | 西门子医疗有限公司 | 用于提供适合于在对患者的辐照的规划中使用的结果数据的方法 |
| US20190192099A1 (en) * | 2017-12-21 | 2019-06-27 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for medical imaging of intervertebral discs |
| CN108805950A (zh) * | 2018-05-28 | 2018-11-13 | 牙博士医疗控股集团有限公司 | 口腔牙片信息转换方法及装置 |
| CN110363765A (zh) * | 2019-07-24 | 2019-10-22 | 苏州科技城医院 | 基于ct影像的骨质疏松参数自动测量方法 |
| CN110660044A (zh) * | 2019-08-30 | 2020-01-07 | 博志生物科技(深圳)有限公司 | 骨组织结构形态异常的快速检测方法及电子装置 |
| CN110840473A (zh) * | 2019-11-22 | 2020-02-28 | 江研伟 | 一种基于CT薄层扫描Hu值的骨密度测量系统 |
| CN111881910A (zh) * | 2020-07-31 | 2020-11-03 | 杭州依图医疗技术有限公司 | 基于椎骨影像的信息处理方法及计算机可读存储介质 |
| CN113679404A (zh) * | 2021-07-06 | 2021-11-23 | 湖南浚朗科技有限公司 | 一种骨密度测量系统及其测量方法 |
Non-Patent Citations (6)
| Title |
|---|
| PENG ZHANG等: "Visceral adiposity is negatively associated with bone density and muscle attenuation", vol. 101, no. 2, pages 337 * |
| THOMAS M. LINK等: "In Vitro and in Vivo Spiral CT to Determine Bone Mineral Density: Initial Experience in Patients at Risk for Osteoporosis", vol. 231, no. 3, pages 805 - 811 * |
| 吴亚楠等: "老年2型糖尿病患者骨量与骨代谢生化指标的变化", no. 07, pages 847 - 850 * |
| 泮宸帅等: "肱骨近端骨质疏松对肩袖损伤修复影响的研究进展", 《中华肩肘外科电子杂志》, vol. 6, no. 1, pages 73 - 76 * |
| 许岩: "无体膜法定量CT骨密度测量在临床上的应用价值", 《中国优秀硕士学位论文全文数据库医药 卫生科技辑》, no. 02, pages 065 - 249 * |
| 郭永杰等: "河南地区健康成人定量CT测量腰椎骨密度的分析", vol. 22, no. 5, pages 442 - 444 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114782375B (zh) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3021753B1 (en) | Systems and methods for determining hepatic function from liver scans | |
| CN100463655C (zh) | 一种腹部ct图像自动测量装置 | |
| WO2006002312A2 (en) | Systems and methods for quantifying symmetry to evaluate medical images | |
| JP2004522465A (ja) | Ｘ線画像解析のための方法および装置 | |
| JP2013542046A (ja) | 超音波画像処理のシステムおよび方法 | |
| CN113222951B (zh) | 一种识别髋关节x线的骨质疏松人工智能诊断装置 | |
| EP4069354A1 (en) | System and method for determining radiation parameters | |
| CN212853503U (zh) | 肝肿瘤智能分析装置 | |
| US20160180520A1 (en) | Quantitative method for 3-d joint characterization | |
| EP1629434A2 (en) | Computation of wall thickness | |
| Xiao et al. | Using 3D medical modeling to evaluate the accuracy of single-photon emission computed tomography (SPECT) bone scintigraphy in diagnosing condylar hyperplasia | |
| JP5579535B2 (ja) | 胎児の肋骨数を測定する超音波システムおよび方法 | |
| CN117455931A (zh) | 基于深度学习mr体部脂肪组织的分割和定量测量方法 | |
| KR20150072910A (ko) | 대상체의 측정에 사용될 캘리퍼 타입에 따라 위치가 교정된 포인트를 의료 영상에서 표시하기 위한 방법 및 장치 | |
| CN114782375A (zh) | 骨质密度测量方法、装置和设备 | |
| CN111583219B (zh) | 颅颌面软硬组织的分析方法及装置、电子设备 | |
| CN108720870A (zh) | 一种基于超声衰减系数的脂肪肝检测系统 | |
| JP2001029343A (ja) | 体脂肪計測装置 | |
| US7848551B2 (en) | Method and system for analysis of bone density | |
| EP3995081A1 (en) | Diagnosis assisting program | |
| CN109993754B (zh) | 从图像上进行颅骨分割的方法和系统 | |
| US9271690B2 (en) | Method and system to estimate visceral adipose tissue by restricting subtraction of subcutaneous adipose tissue to coelom projection region | |
| CN111557681B (zh) | 一种螺旋ct图像分析方法 | |
| CN116843648B (zh) | 基于锥形束ct影像的髁突骨改建三维自动定量测量系统 | |
| US20240225588A1 (en) | Method and Apparatus of Intelligent Analysis for Liver Tumor |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20230321 Address after: No. 524, 515, 516, 517, 518, Floor 5, Unit 2, Building 5, No. 42, Shuxi Road, Jinniu Hi-tech Industrial Park, Chengdu, Sichuan, 610000 Applicant after: Yingtu (Chengdu) Medical Technology Co.,Ltd. Address before: No. 37, Wuhou District National School Lane, Chengdu, Sichuan Province Applicant before: WEST CHINA HOSPITAL OF SICHUAN University |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |