CN110840479B - 一种骨成分确定方法、装置、设备及存储介质 - Google Patents
一种骨成分确定方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN110840479B CN110840479B CN201911222174.1A CN201911222174A CN110840479B CN 110840479 B CN110840479 B CN 110840479B CN 201911222174 A CN201911222174 A CN 201911222174A CN 110840479 B CN110840479 B CN 110840479B
- Authority
- CN
- China
- Prior art keywords
- bone
- data
- positioning image
- bone component
- scanning
- 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
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 254
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004088 simulation Methods 0.000 claims abstract description 41
- 239000000203 mixture Substances 0.000 claims abstract description 39
- 238000012937 correction Methods 0.000 claims abstract description 15
- 230000011218 segmentation Effects 0.000 claims description 6
- 238000002591 computed tomography Methods 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 229910000389 calcium phosphate Inorganic materials 0.000 description 7
- 239000001506 calcium phosphate Substances 0.000 description 7
- 235000011010 calcium phosphates Nutrition 0.000 description 7
- QORWJWZARLRLPR-UHFFFAOYSA-H tricalcium bis(phosphate) Chemical compound [Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O QORWJWZARLRLPR-UHFFFAOYSA-H 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 230000004807 localization Effects 0.000 description 3
- 230000037182 bone density Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 208000029725 Metabolic bone disease Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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
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)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明实施例公开了一种骨成分确定方法、装置、设备及存储介质。该方法包括：分别获取受检部位的定位像扫描数据与旋转扫描数据，其中，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的；得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据；对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分以得到校正后骨成分。本发明实施例的技术方案，通过在定位像扫描和旋转扫描中采用不同能量的X射线，实现了骨成分的量化且精准的确定，从而提高了骨硬化校正的准确性，且不会增加硬件成本或改变临床操作工作流程。
Description
技术领域
本发明实施例涉及医学图像处理技术领域，尤其涉及一种骨成分确定方法、装置、设备及存储介质。
背景技术
电子计算机断层扫描(Computed Tomography，CT)中的X射线源通常是具有一定频谱宽度的球管，这使得CT图像只能获得人体组织的等效密度值，而无法区分不同组织的原子序数。这样一来，当受检部位含有大量骨组织时，在CT图像中会出现严重的骨硬化伪影。
现有的骨硬化伪影的校正方案，通常是先定性估计骨成分，其中，骨成分可以是骨组织中的磷酸钙的和/或水的比例；然后，根据已估计的骨成分对CT重建图像进行正投影以实现图像校正。但是，由于骨成分是估计值，而不同个体的骨成分会有差异，这使得依据现有方案校正后的CT重建图像中依然存在残留的骨硬化伪影，这将严重影响临床诊断的精准性。
发明内容
本发明实施例提供了一种骨成分确定方法、装置、设备及存储介质，以实现精准且量化的确定骨成分的效果。
第一方面，本发明实施例提供了一种骨成分确定方法，可以包括：
分别获取受检部位的定位像扫描数据与旋转扫描数据，其中，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的；
得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据；
对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分以得到校正后骨成分。
可选的，根据对比结果调整初始骨成分以得到校正后骨成分，可以包括：
根据对比结果调整初始骨成分，并根据调整结果更新初始骨成分；
重复执行根据初始骨成分和定位像扫描角度的步骤，直至对比结果满足预设阈值条件，将初始骨成分作为校正后骨成分。
可选的，在此基础上，上述方法还可以包括：
从重建图像中分割出骨组织图像，并根据骨组织图像的电子计算机断层扫描值，得到初始骨成分。
可选的，从重建图像中分割出骨组织图像，可以包括：
基于预设分割阈值从重建图像中分割出骨组织图像。
可选的，对比定位像扫描数据和定位像模拟数据，可以包括：
基于预设指标对比定位像扫描数据和定位像模拟数据，其中，预设指标包括均方差、均值和总值中的至少一个。
可选的定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的，可以包括：
定位像扫描数据是基于第一预设管电压扫描得到的，旋转扫描数据是基于第二预设管电压扫描得到的，第一预设管电压和第二预设管电压间的差值大于预设电压阈值。
可选的，上述骨成分确定方法，还可以包括：
基于校正后骨成分对重建图像进行骨硬化校正。
第二方面，本发明实施例还提供了一种骨成分确定装置，该装置可以包括：
获取模块，用于分别获取受检部位的定位像扫描数据与旋转扫描数据，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的；
得到模块，用于得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据；
调整模块，用于对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分，以得到校正后骨成分。
第三方面，本发明实施例还提供了一种设备，该设备可以包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当一个或多个程序被一个或多个处理器执行，使得一个或多个处理器实现本发明任意实施例所提供的骨成分确定方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明任意实施例所提供的骨成分确定方法。
本发明实施例的技术方案，分别获取受检部位的基于不同的预设管电压扫描得到的定位像扫描数据与旋转扫描数据，这一操作无需特殊硬件支持和/或增加扫描剂量；针对旋转扫描数据的重建图像，根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据，由于初始骨成分和骨组织的真实的骨成分间可能存在差异，则定位像模拟数据和定位像扫描数据间亦可能存在差异；根据定位像扫描数据和定位像模拟数据间的对比结果调整初始骨成分以得到校正后骨成分，也就是说，根据两个数据间的差异反向推导两个骨成分间的差异，由此得到的校正后骨成分无限接近于骨组织的真实的骨成分。上述技术方案，通过在定位像扫描和旋转扫描中采用不同能量的X射线，实现了骨成分的量化且精准的确定，从而提高了骨硬化校正的准确性，特别是改善了头部骨硬化伪影的校正结果，且不会增加硬件成本或改变临床操作工作流程。
附图说明
图1是本发明实施例一中的一种骨成分确定方法的流程图；
图2是本发明实施例一中的一种骨成分确定方法的优选实施例图；
图3是本发明实施例二中的一种骨成分确定装置的结构框图；
图4是本发明实施例三中的一种设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
在介绍本发明实施例之前，先对本发明实施例的应用场景进行示例性说明：人体中的骨组织主要是由磷酸钙和水构成的，当能够精准且量化的确定骨成分后，才能更好地校正CT图像中的骨硬化伪影，上述骨成分可以是骨组织中的磷酸钙的比例和/或水的比例。
CT临床扫描过程通常包含定位像扫描与旋转扫描，其中，定位像扫描是扫描机架保持不动、扫描床移动且进行连续放线扫描，获得类似于X射线机的定位像扫描数据；旋转扫描可以分为断层扫描和螺旋扫描，断层扫描的扫描床保持不动，螺旋扫描的扫描床移动，获得旋转扫描数据。因此，定位像扫描数据可以认为是某一个角度下受检部位的投影数据，且该投影数据在旋转扫描数据中存在对应的投影数据。
由上可知，针对定位像扫描数据和在旋转扫描数据中相同的定位像扫描角度和扫描床位置对应的投影数据，当基于相同的预设管电压进行定位像扫描和旋转扫描时，二者是一致的；当基于不同的预设管电压进行定位像扫描和旋转扫描时，二者存在偏差。由于人体主要是由水和骨组织构成，且骨组织的密度和原子序数都比水大，则骨成分的大小对不同能谱下的投影数据的影响更大。因此，可以考虑通过对这一偏差的调整来定量分析骨成分，也就是说，通过两个不同能谱的投影数据间的比较来定量分析骨成分。
实施例一
图1是本发明实施例一中提供的一种骨成分确定方法的流程图。本实施例可适用于定量确定骨成分的情况，尤其适用于基于双能管电压定量确定骨成分的情况。该方法可以由本发明实施例提供的骨成分确定装置来执行，该装置可以由软件和/或硬件的方式实现，该装置可以集成在各种用户设备上。
参见图1，本发明实施例的方法具体包括如下步骤：
S110、分别获取受检部位的定位像扫描数据与旋转扫描数据，其中，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的。
其中，考虑到同一个体的同一区域的骨成分基本一致，不同个体的同一区域的骨成分可能存在差异的情况，可以根据受检部位的基于不同的预设管电压扫描得到的定位像扫描数据与旋转扫描数据，得到受检部位的骨组织的骨成分。而且，两个预设管电压的差值越大，骨成分的定量分析越准确，这是因为磷酸钙和水在不同的预设管电压的X射线下的衰减系数是不同的。比如，若基于第一预设管电压扫描得到定位像扫描数据，且基于第二预设管电压扫描得到旋转扫描数据，则第一预设管电压和第二预设管电压间的差值可以大于预设电压阈值。示例性的，第一预设管电压是80kVp，第二预设管电压是120kVp，或者，第一预设管电压是120kVp，第二预设管电压是80kVp，等等。
需要说明的是，一方面，基于定位像扫描数据与旋转扫描数据确定骨成分，而不是基于两个旋转扫描数据确定骨成分的原因在于，由于同一个体的同一区域的骨成分基本一致，则针对同一预设管电压在不同角度下扫描得到的定位像扫描数据基本一致，换言之，可以基于某一角度下扫描得到的定位像扫描数据，表征各个角度下扫描得到的定位像扫描数据，一个定位像扫描数据就可以表征一个旋转扫描数据。这样一来，基于两个不同的预设管电压扫描得到的定位像扫描数据与旋转扫描数据，已经可以定量分析骨成分。而且，基于定位像扫描数据与旋转扫描数据确定骨成分的实现方案，无需特殊硬件支持，比如双能球管、具有能谱分辨能力的探测器等等，亦无需增加扫描剂量，比如两次不同能量的旋转扫描，这可显著降低需要受检者承受的检查成本和辐射剂量。
另一方面，不同的预设管电压可以通过采用不同的管电压、在定位像扫描和旋转扫描时采用不同的射线滤过装置等方式来实现。而且，定位像扫描并未限定是某一角度下的扫描，比如，其可以采用正位与侧位片，正位片是受检部位的正面的定位像，侧位片是受检部位的侧面的定位像。
S120、得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据。
其中，对旋转扫描数据进行图像重建得到旋转扫描数据的重建图像，该重建图像可以呈现出受检部位的三维密度数据。在此基础上，可以根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据，而且，进行正投影的管电压和进行定位像扫描时的管电压是一致的。即，根据正投影过程中涉及到的物理模型的材料、投影角度和X射线的能谱，对重建图像进行正投影，由此得到唯一的投影数据，即定位像模拟数据。
需要说明的是，上述初始骨成分可以通过多种方式获得，比如，可以根据经验数据获得，以骨成分是骨组织中的磷酸钙的比例为例，儿童的骨成分可以是40％，因为儿童的骨密度比较低，成年人的骨成分可以是80％，因为成年人的骨密度比较高，不同骨密度的骨组织对X射线的响应是不同的。再比如，若对多个标准模体进行预先扫描，当各标准模体的骨成分已知且互不相同时，根据扫描结果可以得到骨成分和CT值的映射关系，这样一来，实际应用中，初始骨成分可以根据重建图像中骨组织图像的CT值和映射关系获得。再比如，若校正后骨成分的确定过程是一个迭代过程，则初始骨成分只是对迭代步数产生影响，而对迭代结果并不会产生影响，因此，可以随机设置一个初始骨成分。
当然，对于上文所述的根据重建图像中骨组织图像的CT值确定初始骨成分的实现方案，可以先从重建图像中分割出骨组织图像，再根据骨组织图像的CT值确定初始骨成分。针对从重建图像中分割出骨组织图像，一种可选的实施方式是，由于重建图像主要由骨组织和水构成，而骨组织的CT值和水的CT值的差异较大，比如，水的CT值是0，骨组织的CT值是几百，骨组织的CT值与骨组织中磷酸钙的比例有关，磷酸钙的比例越大，骨组织的CT值越大，因此，可以基于预设分割阈值从重建图像中分割出骨组织图像。当然，还可以基于其它实施方式从重建图像中分割出骨组织图像，在此未做具体限定。
S130、对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分以得到校正后骨成分。
其中，对比定位像扫描数据和定位像模拟数据的实现方式有多种，比如，可以基于预设指标对比定位像扫描数据和定位像模拟数据，该预设指标可以是均方差、均值、总值和单点对比中的至少一个，也可以是其余的数据指标。预设指标的具体选择可以根据实际的应用场景确定，比如，若需要得到精准度较高的校正后骨成分，可以对定位像扫描数据和定位像模拟数据中的像素点进行一一对比；若需要快速确定校正后骨成分，可以对定位像扫描数据中各像素点的总值和定位像模拟数据中各像素点的总值进行对比。
由于初始骨成分是根据各种已知数据估算得到的，其和骨组织的真实的骨成分间可能存在差异，当定位像模拟数据是在初始骨成分的基础上得到时，这就使得定位像模拟数据和定位像扫描数据(即，骨组织的真实的扫描数据)间可能存在差异，通过两个数据间的差异就可以反向推导出两个骨成分间的差异。因此，可以根据定位像扫描数据和定位像模拟数据间的对比结果，调整初始骨成分以得到校正后骨成分。示例性的，若初始骨成分是40％，根据对比结果可以确定校正后骨成分大于初始骨成分，则可将初始骨成分由40％调整至50％，并将50％作为校正后骨成分。可选的，进一步的，可以基于校正后骨成分对重建图像进行骨硬化校正，精准确定的校正后骨成分可以提高骨硬化校正的准确性。
本发明实施例的技术方案，分别获取受检部位的基于不同的预设管电压扫描得到的定位像扫描数据与旋转扫描数据，这一操作无需特殊硬件支持和/或增加扫描剂量；针对旋转扫描数据的重建图像，根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据，由于初始骨成分和骨组织的真实的骨成分间可能存在差异，则定位像模拟数据和定位像扫描数据间亦可能存在差异；根据定位像扫描数据和定位像模拟数据间的对比结果调整初始骨成分以得到校正后骨成分，也就是说，根据两个数据间的差异反向推导两个骨成分间的差异，由此得到的校正后骨成分无限接近于骨组织的真实的骨成分。上述技术方案，通过在定位像扫描和旋转扫描中采用不同能量的X射线，实现了骨成分的量化且精准的确定，从而提高了骨硬化校正的准确性，特别是改善了头部骨硬化伪影的校正结果，且不会增加硬件成本或改变临床操作工作流程。
当然，若为进一步提高校正后骨成分的精确度，可以迭代执行上述步骤，一种可选的技术方案，根据对比结果调整初始骨成分以得到校正后骨成分，具体可以包括：根据对比结果调整初始骨成分，并根据调整结果更新初始骨成分；重复执行根据初始骨成分和定位像扫描角度的步骤，直至对比结果满足预设阈值条件，将初始骨成分作为校正后骨成分。也就是说，当定位像扫描数据和定位像模拟数据间的差异比较小时，此时得到的初始骨成分就十分接近于骨组织的真实的骨成分，可将初始骨成分作为校正后骨成分。
示例性的，若初始骨成分是40％，根据对比结果可以确定校正后骨成分大于初始骨成分，则可将初始骨成分由40％调整至50％，并将初始骨成分更新为50％；在此基础上，重新计算定位像模拟数据，并再次对比定位像扫描数据和定位像模拟数据，根据对比结果将初始骨成分由50％调整至45％……循环往复，直至根据最新的初始骨成分计算得到的定位像模拟数据和最先得到的定位像扫描数据间的差值满足预设阈值条件，将此时的初始骨成分作为校正后骨成分。
在上述技术方案的基础上，骨成分确定方法的具体实现过程还可以如图2所示：基于80kVp扫描得到受检部位的定位像扫描数据，并基于120kVp扫描得到受检部位的旋转扫描数据；在完成旋转扫描数据的CT图像重建后，提取重建图像中的骨组织，并根据骨组织图像的CT值初始化一个骨成分，得到初始骨成分，其可通过x％-Ca3P2O8表示，即骨组织中磷酸钙的比例是x％；在该初始骨成分下进行80kVp的定位像扫描角度的正投影，得到定位像模拟数据；将定位像模拟数据与对应的定位像扫描数据进行比较，如有偏差则进行迭代调整x，直至两者的偏差小于预设阈值，将此时的初始骨成分作为校正后骨成分。在确定受检部位的校正后骨成分后，可对重建图像进行常规的骨硬化校正。
实施例二
图3为本发明实施例二提供的骨成分确定装置的结构框图，该装置用于执行上述任意实施例所提供的骨成分确定方法。该装置与上述各实施例的骨成分确定方法属于同一个发明构思，在骨成分确定装置的实施例中未详尽描述的细节内容，可以参考上述骨成分确定方法的实施例。参见图3，该装置具体可包括：获取模块310、得到模块320和调整模块330。
其中，获取模块310，用于分别获取受检部位的定位像扫描数据与旋转扫描数据，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的；
得到模块320，用于得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据；
调整模块330，用于对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分，以得到校正后骨成分。
可选的，调整模块330，具体可以包括：
调整更新单元，用于根据对比结果调整初始骨成分，并根据调整结果更新初始骨成分；
重复单元，用于重复执行根据初始骨成分和定位像扫描角度的步骤，直至对比结果满足预设阈值条件，将初始骨成分作为校正后骨成分。
可选的，在上述装置的基础上，该装置还可包括：
分割模块，用于从重建图像中分割出骨组织图像，并根据骨组织图像的电子计算机断层扫描值，得到初始骨成分。
可选的，分割模块，具体可用于：
基于预设分割阈值从重建图像中分割出骨组织图像。
可选的，调整模块330，具体可以包括：
对比单元，用于基于预设指标对比定位像扫描数据和定位像模拟数据，其中，预设指标包括均方差、均值和总值中的至少一个。
可选的，获取模块310，具体可以包括：
获取单元，用于定位像扫描数据是基于第一预设管电压扫描得到的，旋转扫描数据是基于第二预设管电压扫描得到的，第一预设管电压和第二预设管电压间的差值大于预设电压阈值。
可选的，在上述装置的基础上，该装置还可包括：
校正模块，用于基于校正后骨成分对重建图像进行骨硬化校正。
本发明实施例二提供的骨成分确定装置，通过获取模块可以分别获取受检部位的基于不同的预设管电压扫描得到的定位像扫描数据与旋转扫描数据，这一操作无需特殊硬件支持和/或增加扫描剂量；针对旋转扫描数据的重建图像，得到模块可以根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据，由于初始骨成分和骨组织的真实的骨成分间可能存在差异，则定位像模拟数据和定位像扫描数据间亦可能存在差异；调整模块可以根据定位像扫描数据和定位像模拟数据间的对比结果调整初始骨成分以得到校正后骨成分，也就是说，根据两个数据间的差异反向推导两个骨成分间的差异，由此得到的校正后骨成分无限接近于骨组织的真实的骨成分。上述装置，通过在定位像扫描和旋转扫描中采用不同能量的X射线，实现了骨成分的量化且精准的确定，从而提高了骨硬化校正的准确性，特别是改善了头部骨硬化伪影的校正结果，且不会增加硬件成本或改变临床操作工作流程。
本发明实施例所提供的骨成分确定装置可执行本发明任意实施例所提供的骨成分确定方法，具备执行方法相应的功能模块和有益效果。
值得注意的是，上述骨成分确定装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
实施例三
图4为本发明实施例三提供的一种设备的结构示意图，如图4所示，该设备包括存储器410、处理器420、输入装置430和输出装置440。设备中的处理器420的数量可以是一个或多个，图4中以一个处理器420为例；设备中的存储器410、处理器420、输入装置430和输出装置440可以通过总线或其它方式连接，图4中以通过总线450连接为例。
存储器410作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的骨成分确定方法对应的程序指令/模块(例如，骨成分确定装置中的获取模块310、得到模块320和调整模块330)。处理器420通过运行存储在存储器410中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的骨成分确定方法。
存储器410可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据设备的使用所创建的数据等。此外，存储器410可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器410可进一步包括相对于处理器420远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置430可用于接收输入的数字或字符信息，以及产生与装置的用户设置以及功能控制有关的键信号输入。输出装置440可包括显示屏等显示设备。
实施例四
本发明实施例四提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种骨成分确定方法，包括：
分别获取受检部位的定位像扫描数据与旋转扫描数据，其中，定位像扫描数据与旋转扫描数据是基于不同的预设管电压扫描得到的；
得到旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对重建图像进行正投影，得到定位像模拟数据；
对比定位像扫描数据和定位像模拟数据，并根据对比结果调整初始骨成分以得到校正后骨成分。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的骨成分确定方法中的相关操作。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。依据这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种骨成分确定方法，其特征在于，包括：
分别获取受检部位的定位像扫描数据与旋转扫描数据，其中，所述定位像扫描数据与所述旋转扫描数据是基于不同的预设管电压扫描得到的；
得到所述旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对所述重建图像进行正投影，得到定位像模拟数据，其中，进行正投影的管电压和进行定位像扫描时的管电压是一致的；
对比所述定位像扫描数据和所述定位像模拟数据，并根据对比结果调整所述初始骨成分以得到校正后骨成分。
2.根据权利要求1所述的方法，其特征在于，所述根据对比结果调整所述初始骨成分以得到校正后骨成分，包括：
根据对比结果调整所述初始骨成分，并根据调整结果更新所述初始骨成分；
重复执行所述根据初始骨成分和定位像扫描角度的步骤，直至所述对比结果满足预设阈值条件，将所述初始骨成分作为校正后骨成分。
3.根据权利要求1所述的方法，其特征在于，还包括：
从所述重建图像中分割出骨组织图像，并根据所述骨组织图像的电子计算机断层扫描值，得到所述初始骨成分。
4.根据权利要求3所述的方法，其特征在于，所述从所述重建图像中分割出骨组织图像，包括：
基于预设分割阈值从所述重建图像中分割出骨组织图像。
5.根据权利要求1所述的方法，其特征在于，所述对比所述定位像扫描数据和所述定位像模拟数据，包括：
基于预设指标对比所述定位像扫描数据和所述定位像模拟数据，其中，所述预设指标包括均方差、均值和总值中的至少一个。
6.根据权利要求1所述的方法，其特征在于，所述定位像扫描数据与所述旋转扫描数据是基于不同的预设管电压扫描得到的，包括：
所述定位像扫描数据是基于第一预设管电压扫描得到的，所述旋转扫描数据是基于第二预设管电压扫描得到的，所述第一预设管电压和所述第二预设管电压间的差值大于预设电压阈值。
7.根据权利要求1所述的方法，其特征在于，还包括：
基于所述校正后骨成分对所述重建图像进行骨硬化校正。
8.一种骨成分确定装置，其特征在于，包括：
获取模块，用于分别获取受检部位的定位像扫描数据与旋转扫描数据，所述定位像扫描数据与所述旋转扫描数据是基于不同的预设管电压扫描得到的；
得到模块，用于得到所述旋转扫描数据的重建图像，并根据初始骨成分和定位像扫描角度，对所述重建图像进行正投影，得到定位像模拟数据，其中，进行正投影的管电压和进行定位像扫描时的管电压是一致的；
调整模块，用于对比所述定位像扫描数据和所述定位像模拟数据，并根据对比结果调整所述初始骨成分，以得到校正后骨成分。
9.一种设备，其特征在于，所述设备包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-7中任一所述的骨成分确定方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1-7任一所述的骨成分确定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911222174.1A CN110840479B (zh) | 2019-12-03 | 2019-12-03 | 一种骨成分确定方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911222174.1A CN110840479B (zh) | 2019-12-03 | 2019-12-03 | 一种骨成分确定方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110840479A CN110840479A (zh) | 2020-02-28 |
| CN110840479B true CN110840479B (zh) | 2023-05-26 |
Family
ID=69607440
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911222174.1A Active CN110840479B (zh) | 2019-12-03 | 2019-12-03 | 一种骨成分确定方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110840479B (zh) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2000021441A1 (en) * | 1998-09-30 | 2000-04-20 | Demetech Ab | Method and device for measuring the bone mineral content of the skeleton |
| CN1551745A (zh) * | 2001-09-03 | 2004-12-01 | 在放射检查方法中确定密度分布和原子序数分布的方法 | |
| JP2010068832A (ja) * | 2008-09-16 | 2010-04-02 | Hitachi Medical Corp | Ｘ線ｃｔ装置 |
| CN103186883A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 一种ct图像重建中骨硬化伪影的校正方法 |
| CN103559729A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种双能谱ct图像迭代重建方法 |
| CN103690184A (zh) * | 2012-09-27 | 2014-04-02 | 西门子公司 | 用于确定骨骼矿物质密度值的方法和计算机断层造影系统 |
| CN103900931A (zh) * | 2012-12-26 | 2014-07-02 | 首都师范大学 | 一种多能谱ct成像方法及成像系统 |
| CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
| CN105608721A (zh) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
| CN108601571A (zh) * | 2016-12-15 | 2018-09-28 | 皇家飞利浦有限公司 | Ct成像系统和用于ct成像系统的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8139829B2 (en) * | 2005-02-09 | 2012-03-20 | Hitachi Medical Corporation | Diagnostic imaging support system and diagnostic imaging support program |
| US8811709B2 (en) * | 2012-11-14 | 2014-08-19 | General Electric Company | System and method for multi-material correction of image data |
-
2019
- 2019-12-03 CN CN201911222174.1A patent/CN110840479B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2000021441A1 (en) * | 1998-09-30 | 2000-04-20 | Demetech Ab | Method and device for measuring the bone mineral content of the skeleton |
| CN1551745A (zh) * | 2001-09-03 | 2004-12-01 | 在放射检查方法中确定密度分布和原子序数分布的方法 | |
| JP2010068832A (ja) * | 2008-09-16 | 2010-04-02 | Hitachi Medical Corp | Ｘ線ｃｔ装置 |
| CN103186883A (zh) * | 2011-12-29 | 2013-07-03 | 上海联影医疗科技有限公司 | 一种ct图像重建中骨硬化伪影的校正方法 |
| CN103690184A (zh) * | 2012-09-27 | 2014-04-02 | 西门子公司 | 用于确定骨骼矿物质密度值的方法和计算机断层造影系统 |
| CN103900931A (zh) * | 2012-12-26 | 2014-07-02 | 首都师范大学 | 一种多能谱ct成像方法及成像系统 |
| CN103961125A (zh) * | 2013-01-31 | 2014-08-06 | 东北大学 | 一种用于锥束ct的ct值校正方法 |
| CN103559729A (zh) * | 2013-11-18 | 2014-02-05 | 首都师范大学 | 一种双能谱ct图像迭代重建方法 |
| CN105608721A (zh) * | 2016-01-30 | 2016-05-25 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及装置 |
| CN108601571A (zh) * | 2016-12-15 | 2018-09-28 | 皇家飞利浦有限公司 | Ct成像系统和用于ct成像系统的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 黄科峰等.能谱CT有效平均原子序数对离体尿路结石成分分析的前瞻性研究.《临床军医杂志》.2011,第39卷(第4期),第615-618页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110840479A (zh) | 2020-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110570492B (zh) | 一种基于神经网络的ct伪影抑制方法、设备以及介质 | |
| JP6106365B2 (ja) | 乳房の形態学的情報を取得する方法 | |
| US8965092B2 (en) | Method and apparatus to determine acquisition parameters in dual-energy tomosynthesis | |
| JP5715335B2 (ja) | 任意のスペクトルにおける計算機式断層写真法画像の形成のための順投影 | |
| US20060159223A1 (en) | Method and apparatus for correcting for beam hardening in CT images | |
| JP2019516460A (ja) | 空間とスペクトル情報に基づく複数エネルギーのｃｔ画像におけるノイズ制御のためのシステムと方法 | |
| EP2732431A2 (en) | Image processing for spectral ct | |
| CN110599559B (zh) | 多能量金属伪影减少 | |
| US10605933B2 (en) | X-ray spectral calibration technique for cone-beam CT | |
| JP2013013722A (ja) | Ｘ線撮像における散乱補正の方法及びシステム | |
| US11419566B2 (en) | Systems and methods for improving image quality with three-dimensional scout | |
| US20120308101A1 (en) | Method and system for reconstruction of tomographic images | |
| JP7346429B2 (ja) | スペクトルコンピュータ断層撮影（ｃｔ）スキャナによって生成される画質が向上したバーチャル非造影画像 | |
| CN111627083B (zh) | 骨硬化伪影校正方法、装置、计算机设备和可读存储介质 | |
| US20140016847A1 (en) | Multi-phase computed tomography image reconstruction | |
| US20190175133A1 (en) | Variance reduction for monte carlo-based scatter estimation | |
| US11337671B2 (en) | Methods and systems for improved spectral fidelity for material decomposition | |
| CN110458913B (zh) | 一种多阈值分割ct图像校正图像重建中骨硬化伪影的方法 | |
| US20060257011A1 (en) | Method and system for filtering scan data | |
| CN110840479B (zh) | 一种骨成分确定方法、装置、设备及存储介质 | |
| US10552992B2 (en) | Poly-energetic reconstruction method for metal artifacts reduction | |
| JP7346546B2 (ja) | Ｃｔ仮想単色撮像のための自動的適応エネルギ設定を行う方法 | |
| US11954769B2 (en) | System for reconstructing an image of an object | |
| CN107341836B (zh) | 一种ct螺旋扫描图像重建方法及装置 | |
| KR101613391B1 (ko) | 치료 계획용 ct 영상을 이용한 산란선 제거 방법 및 시스템 |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co.,Ltd. Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |