CN106618620B - 骨硬化伪影校正系数计算方法及装置 - Google Patents
骨硬化伪影校正系数计算方法及装置 Download PDFInfo
- Publication number
- CN106618620B CN106618620B CN201611244293.3A CN201611244293A CN106618620B CN 106618620 B CN106618620 B CN 106618620B CN 201611244293 A CN201611244293 A CN 201611244293A CN 106618620 B CN106618620 B CN 106618620B
- Authority
- CN
- China
- Prior art keywords
- substance
- hardening
- projection value
- bone
- different thicknesses
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 88
- 210000000988 bone and bone Anatomy 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 33
- 208000034189 Sclerosis Diseases 0.000 title description 2
- 239000000126 substance Substances 0.000 claims abstract description 78
- 238000002591 computed tomography Methods 0.000 claims abstract description 24
- 238000003325 tomography Methods 0.000 claims abstract description 17
- 239000000463 material Substances 0.000 claims description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 34
- 229910000389 calcium phosphate Inorganic materials 0.000 claims description 14
- 239000001506 calcium phosphate Substances 0.000 claims description 14
- 235000011010 calcium phosphates Nutrition 0.000 claims description 14
- QORWJWZARLRLPR-UHFFFAOYSA-H tricalcium bis(phosphate) Chemical group [Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O QORWJWZARLRLPR-UHFFFAOYSA-H 0.000 claims description 14
- 230000001419 dependent effect Effects 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- VVQNEPGJFQJSBK-UHFFFAOYSA-N Methyl methacrylate Chemical compound COC(=O)C(C)=C VVQNEPGJFQJSBK-UHFFFAOYSA-N 0.000 claims description 2
- 229920005372 Plexiglas® Polymers 0.000 claims description 2
- NTHWMYGWWRZVTN-UHFFFAOYSA-N sodium silicate Chemical compound [Na+].[Na+].[O-][Si]([O-])=O NTHWMYGWWRZVTN-UHFFFAOYSA-N 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 6
- 230000015654 memory Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000001514 detection method Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000013334 tissue model Methods 0.000 description 3
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 229910052791 calcium Inorganic materials 0.000 description 2
- 239000011575 calcium Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 210000004872 soft tissue Anatomy 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000002083 X-ray spectrum Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000004519 manufacturing process 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
- 238000012805 post-processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- 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/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
-
- 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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (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)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Pulmonology (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种骨硬化伪影校正系数计算方法及装置，该方法包括：使用计算机断层设备进行模体扫描；设定骨组织由第一物质及第二物质组成，获取不同厚度的所述第一及第二物质组合在断层扫描系统的理论投影值；对该理论投影值进行第一物质硬化校正，获取第一物质硬化校正后投影值；计算不同厚度的所述第一物质及第二物质组合在该断层扫描系统的理想投影值；根据所述第二物质的厚度、第一及第二物质的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数。利用本发明方案获得的校正系数，能够在具备良好普适性的前提下，高效地去除骨硬化伪影。
Description
本申请是申请日为2016年1月30日、申请号为“201610069409.8”、发明名称为“计算机断层成像伪影校正方法及装置”的发明专利申请的分案申请。
【技术领域】
本发明涉及计算机断层成像技术领域，尤其涉及一种计算机断层成像伪影的校正方法及装置。
【背景技术】
计算机断层成像是用射线对人体的特定部位按一定厚度的层面进行扫描，根据不同的人体组织对射线的吸收能力不同，利用计算机重建出断层面图像的技术。
在利用X射线进行计算机断层扫描及重建过程中，因球管产生的X射线具有一定频谱宽度，物质对X射线的吸收系数随X线能量的增大而减小，连续能谱的X射线穿过如人体等被扫描物体后，低能量射线易被吸收，高能量射线较易穿过，射束平均能量会变高，射线逐渐变硬。该种效应称之为射束硬化效应。射束硬化效应的存在，会使图像重建时出现伪影，影响图像的重建质量。因而现有技术在重建图像之前会对投影数据进行基于水模的射线硬化校正，这种校正会将软组织的X射线硬化现象消除，但无法消除由于人体骨头引起的射线硬化伪影，即骨硬化伪影。
已知存在各种针对骨硬化伪影的校正方法：一种是基于图像后处理技术，通过经验参数来消除骨硬化伪影(如：Jiang Hsieh et al,“An iterative approach to thebeam hardening correction in cone beam CT”,Med.Phys.27 1,January 2000)，其弊端在于矫正系数的来源缺乏理论依据，校正准确性及效率较差；另一种是通过预先扫描特制的骨组织仿体来产生矫正系数(如：专利CN01124649.9-计算机层析X射线摄影设备)，其弊端在于通过扫描骨组织仿体得到的系数往往缺乏普适性(不同年龄人群的骨组织成分差异很大，往往不能用同样的矫正系数)。因而，上述方案均不能很好的解决骨硬化伪影问题。
因此，需要提出一种新的计算机断层成像骨硬化伪影校正方法及实施该方法的装置，能够在具备良好普适性的前提下，高效地去除骨硬化伪影。
【发明内容】
本发明解决的是计算机断层成像图像中出现骨硬化伪影的问题。
为解决上述问题，本发明提出一种计算机断层成像骨硬化伪影校正方法，包括：使用计算机断层设备进行模体扫描；设定骨组织由第一物质及第二物质组成，获取不同厚度的所述第一及第二物质组合在断层扫描系统的理论投影值；对该理论投影值进行第一物质硬化校正，获取第一物质硬化校正后投影值；计算不同厚度的所述第一物质及第二物质组合在该断层扫描系统的理想投影值；根据所述第二物质的厚度、第一及第二物质的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数；使用所述骨硬化校正系数进行伪影校正。
可选地，还包括：获取使所述模体的测量投影值与理论投影值相等时的每个探测单元对应的等效滤过厚度。
可选地，所述模体为厚度及材料已知的均匀模体。
可选地，所述模体的材料为水或有机玻璃。
可选地，所述第一物质硬化校正包括：计算不同厚度第一物质在该断层扫描系统的理论投影值及理想投影值；对该理论投影值及理想投影值进行多项式拟合，获得第一物质硬化校正系数；使用该第一物质硬化校正系数进行校正。
可选地，所述根据所述第二物质厚度、理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数包括：以第二物质厚度为自变量，以理想投影值及第一物质硬化校正后投影值之差为因变量进行多项式拟合，获得骨硬化校正系数。
可选地，所述根据所述第二物质厚度、理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数包括：以第一物质厚度、第二物质厚度为自变量，以理想投影值及第一物质硬化校正后投影值之差为因变量进行曲面拟合，获得骨硬化校正系数。
可选地，所述第一物质为水，所述第二物质为含钙物质。
可选地，所述第二物质为磷酸钙。
本发明还提供一种计算机断层成像骨硬化伪影校正装置，包括：存储单元，存储有使用上述方法所获得的骨硬化伪影校正系数；校正单元，根据扫描数据选择所述校正系数进行骨硬化伪影校正。
本发明对比现有技术有如下的有益效果：
本方案基于基材料理论，将骨组织认为由两种基材料物质组成，通过一次模体扫描，即可产生具有普适性的矫正系数，只需要在重建过程中应用合适的组织模型，即可得到较好的骨硬化伪影校正效果。
【附图说明】
图1是本发明的计算机断层成像系统的结构示意图；
图2是本发明一实施例的骨硬化伪影校正方法流程示意图；
图3是本发明一实施例中利用水硬化校正系数进行校正的流程示意图；
图4是本发明一实施例中校正系数表示例。
【具体实施方式】
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图对本发明的具体实施方式做详细的说明。
图1是一种计算机断层成像系统的结构示意图，如图1所示，计算机断层成像系统100包括机架110，所述机架110具有围绕系统轴线旋转的可旋转的部分130。可旋转的部分130具有相对设置的X射线源131和X射线探测器132的X射线系统。
计算机断层成像系统100还具有检查床120，在进行检查时，患者在该检查床120上可以沿着Z轴方向被推入到扫描腔体中。X射线源131绕S轴旋转，探测器132相对于X射线源131一起运动，以采集投影测量数据，这些数据在之后被用于重建图像。还可以进行螺旋扫描，在螺旋扫描期间，通过患者沿着S轴的连续运动和X射线源131的同时旋转，X射线源131相对于患者产生螺旋轨迹。
所述计算机断层成像系统100还可以包括控制单元和图像重建单元，所述控制单元用于在扫描过程中根据特定的扫描协议控制计算机断层成像系统100的各部件。所述图像重建单元用于根据探测器132采样的待校正数据重建出图像。
以上，仅以示例方式阐释了可使用本发明所提供骨硬化伪影校正方法的计算机断层成像设备，本领域技术人员理解，其它如使用X射线的C型臂系统等设备，或组合式医学成像系统(例如：组合式正电子发射断层成像-计算机断层成像，Positron EmissionTomography-Computed tomography Tomography,PET-CT)，或使用其它类型射线的断层成像设备等，均可适用本发明所述校正方法及装置，本发明对计算机断层成像设备的类型与结构并不做具体限定。
当受检对象在上述任一种类的计算机断层成像设备中进行扫描成像时，因射束硬化效应的存在会导致骨硬化伪影的存在，影响图像成像质量及导致阅图不便，因而需对此类伪影进行校正。以下，是根据本发明公布方案对骨硬化伪影进行校正的具体实施方式举例：
图2是根据本发明一实施例的骨硬化伪影校正方法流程示意图：
执行步骤S1，使用计算机断层设备进行模体扫描。此处的模体可选择使用厚度及材料已知的均匀模体，优选地，模体的材料可选择与人体软组织化学成分相似的材料，例如水或有机玻璃。
执行步骤S2，获取不同厚度的第一物质及第二物质组合在断层扫描系统的理论投影值。在X射线扫描中，根据基材料分解理论，任何组织的线性衰减系数可以表示为两种基材料的质量衰减系数的线性组合。本方案中设定骨组织由两种不同的物质(基材料)组成，例如，将骨组织认为是一定比例水和另外一种成分的物质混合构成。另外一种成分可优选为含钙量较高的物质(本实施例中选择磷酸钙)。
在对物体进行计算机断层扫描过程中，测量投影值表示探测器实际测到的被扫描物体投影值，理论投影值表示考虑X射线光子能量分布而计算得到的该被扫描物体的投影值，而理想投影值则表示光子能量均为E0时(E0为可配置参数，代表光子能量为单一的该特定值)的X射线经过该物体时的投影值。例如，测量投影值projMeas可表示如下：
公式(1)中，I0和I分别表示入射被扫描物体和透射出被扫描物体的X射线强度。
步骤S1中的模体的理论投影值ProjCal计算公式可表示如下：
公式(2)中，E代表X射线光子能量，S(E)为球管发出的X射线谱，D(E)为断层扫描系统的探测器响应，μfilter(E)和μphan(E)分别为等效滤过材料和模体材料的线性衰减系数，Lfilter为每一探测单元对应的等效滤过厚度，Lphan为步骤S1中扫描的模体厚度。
进一步地，该实施例中，Lfilter为使模体的测量投影值与理论投影值相等时(可在允许误差范围内)的每个探测单元所对应的等效滤过厚度。其可通过迭代的方法求取，例如，迭代修改该厚度值并根据公式(2)计算模体的理论投影值，直到该理论投影值与测量投影值相等时(允许误差范围内)为止。
在步骤S2中，定义ProjCali,j为不同厚度的第一物质及第二物质组合在断层扫描系统的理论投影值，则对于本实施例中的水与磷酸钙组合，其不同厚度组合的理论投影值可通过下述公式计算获得：
其中，μH2O(E)表示水的线性衰减系数，LH2O,i(i＝0,1,2,…)表示不同的水厚度，μphospca(E)表示磷酸钙的线性衰减系数，Lphospca,j(j＝0,1,2,…)表示不同的磷酸钙厚度(本说明书公式中，如未做特殊说明，相同符号变量含义均与其它公式相同)。
执行步骤S3，对该理论投影值进行第一物质硬化校正，获取第一物质硬化校正后投影值。对于本实施例中的水与磷酸钙组合，则进行水硬化校正，其可通过产生水硬化校正系数的方法进行。
图3示出了本实施例中利用水硬化校正系数进行水硬化校正的流程：
首先，执行步骤S301，确定不同厚度水的理论投影值。该值ProjCalH2O,i可通过以下公式获得：
其中,μH2O(E)表示水的线性衰减系数，LH2O,i(i＝0,1,2,…)表示不同的水厚度。
执行步骤S302，确定不同厚度水的理想投影值。该理想投影值ProjIdealH2O,i可通过以下公式获得：
ProjIdealH2O,i＝μH2O(E0)LH2O,i (5)
式(5)中，μH2O(E0)表示水对于能量为E0的X光子的线性衰减系数，LH2O,i(i＝0,1,2,…)表示不同的水厚度。
执行步骤S303，对不同厚度水的理想投影值及理论投影值进行多项式拟合，获得水硬化校正系数。该水硬化校正系数αk可通过以下公式获得：
式(6)中，N1表示多项式阶数，αk(k＝0,1,…)。
最后，执行步骤S304，使用该水硬化校正系数进行校正。针对步骤S2所获得的不同组合的水及磷酸钙组合的理论投影值ProjCali,j，其校正后投影值ProjCorrectedi,j可通过以下公式获得：
其中，αk为步骤S303中得到的水硬化矫正系数，k＝0,1,…。
完成第一物质硬化校正后，继续如图2所示，执行步骤S4，计算不同厚度的第一物质及第二物质组合在该断层扫描系统的理想投影值。本实施例中，该理想投影值ProjIdeali,j可通过以下公式获得：
ProjIdeali,j＝μH2O(E0)LH2O,i+μphospca(E0)Lphospca,j (8)
其中，μH2O(E0)表示水对于能量为E0的X光子的线性衰减系数，μphospca(E0)表示磷酸钙对于能量为E0的X光子的线性衰减系数，LH2O,i(i＝0,1,2,…)表示不同的水厚度，Lphospca,j(j＝0,1,2,…)表示不同的磷酸钙厚度。
执行步骤S5，根据所述第二物质的厚度、第一及第二物质的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数。本实施例中，针对水与磷酸钙的组合，分别以水的厚度LH2O,i(i＝0,1,2,…)，磷酸钙的厚度Lphospca,j(j＝0,1,2,…)为自变量，以步骤S4中获取的不同厚度的水和磷酸钙组合的理想投影值和步骤S3中获得的水硬化矫正后投影值之差为因变量，进行曲面拟合，获得拟合参数作为骨硬化矫正系数，该过程可依据以下公式实现：
ProjErrori,j＝ProjIdeali,j-ProjCorrectedi,j＝f(LH2O,i,Lphospca,j) (9)
式(9)中，f(LH2O,i,Lphospca,j)为以LH2O,i(i＝0,1,2,…)和Lphospca,j(j＝0,1,2,…)为自变量的曲面函数。
根据本实施例的一个变化例，也可以磷酸钙厚度Lphospca,j(j＝0,1,2,…)为自变量，以步骤S4中获取的不同厚度的水和磷酸钙组合的理想投影值和步骤S3中获得的水硬化矫正后投影值之差为因变量进行多项式拟合，产生骨硬化矫正系数，该变化例过程可依据以下公式实现：
式(10)中，下标i0表示水的厚度为固定值LH2O,i0，N2为多项式阶数，βk(k＝0,1,…,N2)为多项式系数。
执行步骤S6，使用该骨硬化伪影校正系数进行伪影校正。根据步骤S5获取的校正系数可以图4所示形式存储于计算机断层扫描设备重建单元中，扫描重建时，根据扫描对象选择相应的骨组织模型(不同成分比例组合的基材料模型)，再根据骨组织模型提取原始图像中X射线经过的基材料的等效厚度，来选择对应的骨伪影校正系数进行校正，即可得到较好的骨硬化伪影去除效果。
例如，使用该伪影校正系数进行校正的一种实施方式可以为：接收扫描数据，并基于该扫描数据重建待校正图像及该待校正图像的参照图像；对该参照图像的像素点赋予第一物质比例，并基于该第一物质比例获取该参照图像的第一物质基图；对该第一物质基图及该参照图像进行投影，获取该投影操作每一投影射线对应的第一物质等效长度；根据所述每一投影射线对应的第一物质等效长度，确定伪影校正系数表中的伪影校正系数；使用该伪影校正系数对所述待校正图像进行伪影校正。
相应地，本发明也提出了一种骨硬化伪影校正装置，包括：存储单元，存储有根据上述方法所获得的骨硬化伪影校正系数；校正单元，根据扫描数据选择所述校正系数进行骨硬化伪影校正。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于计算机可读存储介质中，存储介质可以包括但不限于：软盘、光盘、CD-ROM、磁光盘、ROM(只读存储器)、RAM(随机存取存储器)、EPROM(可擦除可编程只读存储器)、EEPROM(电可擦除可编程只读存储器)、磁卡或光卡、闪存、或适于存储机器可执行指令的其他类型的介质/机器可读介质。
本发明中，各实施例采用递进式写法，重点描述与前述实施例的不同之处，各实施例中的相同方法或结构参照前述实施例的相同部分。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。
Claims (8)
1.一种骨硬化伪影校正系数计算方法，其特征在于，包括：
使用计算机断层设备进行模体扫描；
设定骨组织由第一物质及第二物质组成，获取不同厚度的所述第一及第二物质组合在断层扫描系统的理论投影值；
对该理论投影值进行第一物质硬化校正，获取第一物质硬化校正后投影值；
计算不同厚度的所述第一物质及第二物质组合在该断层扫描系统的理想投影值；
根据所述第二物质的厚度、不同厚度的所述第一物质及第二物质组合的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数。
2.根据权利要求1所述的计算方法，其特征在于，所述模体为厚度及材料已知的均匀模体。
3.根据权利要求2所述的计算方法，其特征在于，所述模体的材料为水或有机玻璃。
4.根据权利要求1所述的计算方法，其特征在于，所述第一物质硬化校正包括：
计算不同厚度第一物质在该断层扫描系统的理论投影值及不同厚度第一物质在该断层扫描系统的理想投影值；
对所述不同厚度第一物质在该断层扫描系统的理论投影值及不同厚度第一物质在该断层扫描系统的理想投影值进行多项式拟合，获得第一物质硬化校正系数；
使用该第一物质硬化校正系数进行校正。
5.根据权利要求1所述的计算方法，其特征在于，所述根据所述第二物质厚度、不同厚度的所述第一物质及第二物质组合的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数包括：
以第二物质厚度为自变量，以不同厚度的所述第一物质及第二物质组合的理想投影值及第一物质硬化校正后投影值之差为因变量进行多项式拟合，获得骨硬化校正系数。
6.根据权利要求1所述的计算方法，其特征在于，所述根据所述第二物质厚度、不同厚度的所述第一物质及第二物质组合的理想投影值及第一物质硬化校正后投影值，获得骨硬化校正系数包括：
以第一物质厚度、第二物质厚度为自变量，以不同厚度的所述第一物质及第二物质组合的理想投影值及第一物质硬化校正后投影值之差为因变量进行曲面拟合，获得骨硬化校正系数。
7.根据权利要求1所述的计算方法，其特征在于，所述第一物质为水，所述第二物质为含钙物质。
8.根据权利要求7所述的计算方法，其特征在于，所述第二物质为磷酸钙。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611244293.3A CN106618620B (zh) | 2016-01-30 | 2016-01-30 | 骨硬化伪影校正系数计算方法及装置 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611244293.3A CN106618620B (zh) | 2016-01-30 | 2016-01-30 | 骨硬化伪影校正系数计算方法及装置 |
| CN201610069409.8A CN105608721B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像伪影校正方法及装置 |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610069409.8A Division CN105608721B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像伪影校正方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106618620A CN106618620A (zh) | 2017-05-10 |
| CN106618620B true CN106618620B (zh) | 2020-08-04 |
Family
ID=55988636
Family Applications (4)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610069409.8A Active CN105608721B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像伪影校正方法及装置 |
| CN201611247279.9A Active CN106725569B (zh) | 2016-01-30 | 2016-01-30 | 骨硬化伪影校正方法及装置 |
| CN201611244293.3A Active CN106618620B (zh) | 2016-01-30 | 2016-01-30 | 骨硬化伪影校正系数计算方法及装置 |
| CN201611244252.4A Active CN106618619B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像设备 |
Family Applications Before (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610069409.8A Active CN105608721B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像伪影校正方法及装置 |
| CN201611247279.9A Active CN106725569B (zh) | 2016-01-30 | 2016-01-30 | 骨硬化伪影校正方法及装置 |
Family Applications After (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611244252.4A Active CN106618619B (zh) | 2016-01-30 | 2016-01-30 | 计算机断层成像设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (4) | CN105608721B (zh) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017128892A1 (zh) * | 2016-01-30 | 2017-08-03 | 上海联影医疗科技有限公司 | 计算机断层成像伪影校正方法及系统 |
| CN108109183B (zh) * | 2016-11-25 | 2021-06-11 | 上海东软医疗科技有限公司 | 射束硬化校正方法及装置 |
| CN106683144B (zh) * | 2016-12-30 | 2020-07-14 | 上海联影医疗科技有限公司 | 一种图像迭代重建方法及装置 |
| US10977843B2 (en) | 2017-06-28 | 2021-04-13 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining parameters for medical image processing |
| CN107330949B (zh) * | 2017-06-28 | 2020-11-03 | 上海联影医疗科技有限公司 | 一种伪影校正方法及系统 |
| CN107730569B (zh) * | 2017-11-30 | 2021-01-01 | 上海联影医疗科技股份有限公司 | 一种医学图像伪影校正方法及装置 |
| CN108109185B (zh) * | 2017-12-18 | 2021-07-20 | 上海联影医疗科技股份有限公司 | 一种生成用于消除ct伪影的校正系数的方法，以及一种基于校正系数消除ct伪影的方法 |
| US10922855B2 (en) | 2017-11-30 | 2021-02-16 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for determining at least one artifact calibration coefficient |
| CN108670282B (zh) * | 2018-03-28 | 2021-12-17 | 上海联影医疗科技股份有限公司 | 一种骨硬化伪影校正方法 |
| CN108918559B (zh) * | 2018-07-28 | 2021-08-17 | 北京纳米维景科技有限公司 | 一种实现图像自校正的x射线图像探测器及其方法 |
| WO2020077592A1 (zh) * | 2018-10-18 | 2020-04-23 | 清华大学 | Ct系统能谱不一致性的校正方法 |
| CN110097517B (zh) * | 2019-04-28 | 2022-12-27 | 东软医疗系统股份有限公司 | 去除图像伪影的方法及装置 |
| CN111053566B (zh) * | 2019-12-30 | 2021-09-07 | 上海联影医疗科技股份有限公司 | X射线图像重建方法、装置、x射线成像系统及介质 |
| CN110811660B (zh) * | 2019-10-25 | 2023-05-12 | 赛诺威盛科技(北京)股份有限公司 | 一种校正ct射线束硬化伪影的方法 |
| CN110840479B (zh) * | 2019-12-03 | 2023-05-26 | 上海联影医疗科技股份有限公司 | 一种骨成分确定方法、装置、设备及存储介质 |
| CN111080740B (zh) * | 2019-12-27 | 2023-06-16 | 上海联影医疗科技股份有限公司 | 一种图像校正方法、装置、设备及介质 |
| CN111627083B (zh) * | 2020-05-26 | 2023-11-21 | 上海联影医疗科技股份有限公司 | 骨硬化伪影校正方法、装置、计算机设备和可读存储介质 |
| CN112053329B (zh) * | 2020-08-21 | 2023-04-14 | 东软医疗系统股份有限公司 | 骨硬化系数的获取方法、装置、存储介质和电子设备 |
| CN113989400B (zh) * | 2021-09-26 | 2024-08-30 | 清华大学 | 一种ct图像生成方法、装置、电子设备及计算机存储介质 |
| CN113961124B (zh) * | 2021-09-27 | 2024-02-27 | 上海联影医疗科技股份有限公司 | 医学图像显示方法、装置、计算机设备和存储介质 |
| CN114041811A (zh) * | 2021-11-03 | 2022-02-15 | 海南大学 | 模体和模体在扫描设备中的应用方法 |
| CN117796827A (zh) * | 2022-09-26 | 2024-04-02 | 同方威视技术股份有限公司 | 用于成像设备的标定方法、装置、成像设备 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1296011C (zh) * | 2003-07-16 | 2007-01-24 | 东软飞利浦医疗设备系统有限责任公司 | 一种ct机射束硬化的校正方法 |
| US7391844B2 (en) * | 2005-01-14 | 2008-06-24 | General Electric Company | Method and apparatus for correcting for beam hardening in CT images |
| JP2009531109A (ja) * | 2006-03-29 | 2009-09-03 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 温度アーチファクト訂正 |
| DE102007028660B3 (de) * | 2007-06-21 | 2009-01-29 | Siemens Ag | Verfahren zur Korrektur von Bewegungsartefakten bei der Aufnahme von MR-Bildern |
| CN101226642A (zh) * | 2008-01-25 | 2008-07-23 | 西安交通大学 | 基于ct数据一致性的投影射束硬化校正方法 |
| CN101510298B (zh) * | 2009-03-17 | 2010-12-29 | 西北工业大学 | 一种ct伪影的综合校正方法 |
| US8611627B2 (en) * | 2009-12-23 | 2013-12-17 | General Electric Company | CT spectral calibration |
| CN103186883B (zh) * | 2011-12-29 | 2015-09-30 | 上海联影医疗科技有限公司 | 一种ct图像重建中骨硬化伪影的校正方法 |
| US8811709B2 (en) * | 2012-11-14 | 2014-08-19 | General Electric Company | System and method for multi-material correction of image data |
-
2016
- 2016-01-30 CN CN201610069409.8A patent/CN105608721B/zh active Active
- 2016-01-30 CN CN201611247279.9A patent/CN106725569B/zh active Active
- 2016-01-30 CN CN201611244293.3A patent/CN106618620B/zh active Active
- 2016-01-30 CN CN201611244252.4A patent/CN106618619B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106618619B (zh) | 2021-01-08 |
| CN106725569B (zh) | 2020-08-04 |
| CN106618619A (zh) | 2017-05-10 |
| CN106618620A (zh) | 2017-05-10 |
| CN105608721A (zh) | 2016-05-25 |
| CN105608721B (zh) | 2017-04-05 |
| CN106725569A (zh) | 2017-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106618620B (zh) | 骨硬化伪影校正系数计算方法及装置 | |
| US10540766B2 (en) | Systems and methods for artifact correction of computerized tomography | |
| US8483471B2 (en) | Method and system for scatter correction in X-ray imaging | |
| JP6492005B2 (ja) | Ｘ線ｃｔ装置、再構成演算装置、及び再構成演算方法 | |
| JP5220309B2 (ja) | Ｘ線断層撮影装置 | |
| JP2009112627A (ja) | Ｘ線ｃｔ装置 | |
| CN104408758A (zh) | 一种低剂量能谱ct图像处理方法 | |
| CN103961125A (zh) | 一种用于锥束ct的ct值校正方法 | |
| CN110458762B (zh) | 一种基于可调双因子的ct图像射束硬化伪影校正系统 | |
| US20120308101A1 (en) | Method and system for reconstruction of tomographic images | |
| CN105608720B (zh) | 计算机断层成像伪影校正方法及装置 | |
| CN108670282B (zh) | 一种骨硬化伪影校正方法 | |
| JP2004357969A (ja) | Ｘ線計測装置 | |
| CN100365664C (zh) | 基于hl一致性条件的ct投影数据射束硬化效应校正方法 | |
| JP4887132B2 (ja) | Ｘ線ｃｔ装置 | |
| Xiu et al. | An Innovative Beam Hardening Correction Method for Computed Tomography Systems. | |
| CN107730569A (zh) | 一种医学图像伪影校正方法及装置 | |
| Kim et al. | Accuracy improvement of interior region of interest reconstruction using field modulation CT acquisition | |
| JP2015019863A (ja) | 放射線断層撮影装置及びプログラム | |
| US7155045B2 (en) | Method and apparatus for correcting the contrast density of a radiography image | |
| CN107341836B (zh) | 一种ct螺旋扫描图像重建方法及装置 | |
| Jacobson et al. | A kV–MV approach to CBCT metal artifact reduction using multi-layer MV-CBCT | |
| Tang et al. | X-ray projection simulation based on physical imaging model | |
| CN110840479B (zh) | 一种骨成分确定方法、装置、设备及存储介质 | |
| Usui et al. | Cycle-Generative Adversarial Network-Based Image Correction of 4-Dimensional Cone-Beam Computed Tomography for Lung Cancer Adaptive Radiation Therapy |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai city Jiading District Industrial Zone Jiading Road No. 2258 Patentee before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |