CN108542414A - 计算机断层图像的散射校正方法以及装置 - Google Patents
计算机断层图像的散射校正方法以及装置 Download PDFInfo
- Publication number
- CN108542414A CN108542414A CN201810456437.4A CN201810456437A CN108542414A CN 108542414 A CN108542414 A CN 108542414A CN 201810456437 A CN201810456437 A CN 201810456437A CN 108542414 A CN108542414 A CN 108542414A
- Authority
- CN
- China
- Prior art keywords
- image
- projection
- data
- projected image
- target object
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 56
- 238000002591 computed tomography Methods 0.000 title claims abstract description 29
- 239000002184 metal Substances 0.000 claims description 84
- 229910052751 metal Inorganic materials 0.000 claims description 84
- 238000004590 computer program Methods 0.000 claims description 22
- 238000011157 data evaluation Methods 0.000 claims description 5
- 238000007408 cone-beam computed tomography Methods 0.000 abstract description 21
- 230000005855 radiation Effects 0.000 abstract description 11
- 238000010586 diagram Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 6
- 238000003709 image segmentation Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 239000007769 metal material Substances 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000003702 image correction Methods 0.000 description 1
- WABPQHHGFIMREM-UHFFFAOYSA-N lead(0) Chemical compound [Pb] WABPQHHGFIMREM-UHFFFAOYSA-N 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000005295 random walk Methods 0.000 description 1
- 238000000790 scattering method Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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/40—Arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4064—Arrangements for generating radiation specially adapted for radiation diagnosis specially adapted for producing a particular type of beam
- A61B6/4085—Cone-beams
-
- 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
- A61B6/5282—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to scatter
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (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)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
本申请涉及一种计算机断层图像的散射校正方法、系统、计算机设备和存储介质。采用计算机断层图像的散射校正方法,当前标准投影图像利用前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,CBCT系统中在扫描一周均匀采集多张的标准投影图像后,对于每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,能够有效避免重建的CT图像中形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量,在保证CBCT图像质量的同时最大限度的降低目标物体的辐射剂量。
Description
技术领域
本申请涉及图像处理技术领域,特别是涉及一种计算机断层图像的散射校正方法、装置、计算机设备和存储介质。
背景技术
锥形束计算机断层成像(Cone-beam Computed Tomography,CBCT)在医学成像领域起着非常重要的作用,但是,由于CBCT采用的探测器为平板探测器,无法采用准直器遮挡散射线,所以重建的CBCT图像存在严重的散射伪影。CBCT散射严重制约了CBCT图像在临床的深入应用和发展。
传统的基于硬件的散射校正方法通常是在X射线影像系统中添加若干硬件工具,以减少达到平板探测器的散射射线实现散射校正的目的,或者是为了计算出散射射线的量从而从总的射线量中减去散射射线的量以实现散射校正的目的。基于硬件的散射方法主要包括反散射栅、空气隙和衰减板校正方法,但是,反散射栅、空气隙方法往往都存在原发信号损失的情况,往往需要增加病人的受照剂量补偿被衰减的源射线强度,而衰减板校正方法需要对病人进行二次扫描导致病人受照剂量增大;基于硬件的散射校正方法普遍存在X射线照射剂量增大的缺点,难以应用于临床应用中。
发明内容
基于此,有必要针对上述散射校正方法存在X射线照射剂量增大的技术问题,提供一种计算机断层图像的散射校正方法、装置、计算机设备和存储介质。
一种计算机断层图像的散射校正方法,包括以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各所述标准投影图像中的第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像第一侧的散射数据以及所述前一张标准投影图像和所述后一张标准投影图像中第二侧的散射数据;
根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据所述第一目标投影图像以及所述第二目标投影图像获取所述当前标准投影图像散射校正后的目标投影图像。
一种计算机断层图像的散射校正装置,包括:
标准投影图像获取模块,用于获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
第一散射数据获取模块,用于分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
第二散射数据获取模块,用于分别根据各所述第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像第一侧的散射数据和所述前一张标准投影图像以及所述后一张标准投影图像中第二侧的散射数据;
当前投影数据估算模块,用于根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
目标投影图像获取模块,用于将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减,获得所述当前标准投影图像散射校正后的目标投影图像。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各所述标准投影图像中的第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像第一侧的散射数据以及所述前一张标准投影图像和所述后一张标准投影图像中第二侧的散射数据;
根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据所述第一目标投影图像以及所述第二目标投影图像获取所述当前标准投影图像散射校正后的目标投影图像。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各所述标准投影图像中的第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像第一侧的散射数据以及所述前一张标准投影图像和所述后一张标准投影图像中第二侧的散射数据;
根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据所述第一目标投影图像以及所述第二目标投影图像获取所述当前标准投影图像散射校正后的目标投影图像。
上述计算机断层图像的散射校正方法、装置、计算机设备和存储介质,当前标准投影图像利用前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,CBCT系统中在扫描一周均匀采集多张的标准投影图像后,对于每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,能够有效避免重建的CBCT图像中形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量,在保证CBCT图像质量的同时最大限度的降低目标物体的辐射剂量。
附图说明
图1为本发明一个实施例中计算机断层图像的散射校正方法的应用环境图;
图2为本发明一个实施例中散射校正板的结构示意图;
图3为本发明一个实施例中计算机断层图像的散射校正方法的流程图;
图4为本发明一个实施例中根据第一散射数据估算其他位置坐标上的第二散射数据的示意图;
图5为本发明另一个实施例中计算机断层图像的散射校正方法的流程图;
图6为本发明一个实施例中计算机断层图像的散射校正装置的结构示意图;
图7为本发明另一个实施例中计算机断层图像的散射校正装置的结构示意图;
图8为本发明另一个实施例中计算机断层图像的散射校正装置的结构示意图;
图9为本发明一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
需要说明的是,本发明实施例所涉及的术语“第一\第二\第三\第四”仅仅是是区别类似的对象,不代表针对对象的特定排序,可以理解地,“第一\第二\第三”在允许的情况下可以互换特定的顺序或先后次序,以使这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。例如,本发明实施例所涉及的术语“第一侧”以及“第二侧”仅仅是是区别投影图像中的左半侧的图像以及右半侧的图像,或者是区别投影图像中的上半侧的图像以及下半侧的图像;例如,若“第一侧”表示为投影图像中左半侧的图像,则“第二侧”表示为投影图像中右半侧的图像,若“第一侧”表示投影图像中为右半侧的图像,则“第二侧”表示为投影图像中左半侧图像。
本申请提供的计算机断层图像的散射校正方法,可以应用于如图1所示的应用环境中,散射校正板120放置于射线源110与探测器130之间,探测器130与计算机设备140连接;其中,散射校正板120在探测器130的投影占探测器130面积的一半,且散射校正板120可以通过机械传动的方式进行水平偏移,计算机设备140可以是服务器,也可以是个人计算机、笔记本电脑、智能手机和平板电脑。
当需要对被扫描对象进行计算机断层成像扫描时,被扫描对象可以放置于散射校正板120与射线源110之间,也可以放置于散射校正板120与探测器130之间;射线源110发出锥束X线束对被扫描对象以及散射校正板120进行扫描,由探测器130接收透过被扫描对象以及散射校正板120的X线获得投影图像,探测器130将投影图像转变转为数字信号输入计算机设备140中,计算机设备140对获得的投影图像进行散射校正。
可选的,散射校正板120可以是具有一定厚度的矩形薄板,其材料可以选择对X射线吸收量较小的材料;矩形薄板上镶嵌着一系列的金属球,金属球呈交叉点状分布,如图2所示;金属球的材料选择对X射线吸收量较大的金属材料,使得X射线无法传过金属材料。散射校正板120可以设置为其所在平面与探测器130所在平面平行。
在其中一个实施例中,如图3所示,图3为本发明一个实施例中计算机断层图像的散射校正方法的流程图,提供了一种计算机断层图像的散射校正方法,以该方法应用于图1中的应用环境中,目标物体放置于散射校正板120与探测器130之间为例进行说明,包括以下步骤:
步骤S310:获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,当前标准投影图像、前一张标准投影图像以及后一张标准投影图像间的扫描角度间隔相等;当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,当前标准投影图像中第二侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第一侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据。
本步骤中,当前标准投影图像是散射校正板120在探测器130上的投影占据探测器130面积的第一侧时,锥束X线束对目标物体以及散射校正板120进行扫描获得的标准投影图像,所以,当前标准投影图像中第一侧的散射校正板加目标物体投影数据为X射线透过目标对象以及散射校正板120获得的投影图像,第二侧的目标物体投影数据为X射线透过目标对象获得的投影图像;
前一张标准投影图像以及后一张标准投影图像是散射校正板120在探测器130上的投影占据探测器130面积的第二侧时,锥束X线束对目标物体以及散射校正板120进行扫描获得的标准投影图像,所以,前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据为X射线透过目标对象获得的投影图像,第二侧的散射校正板加目标物体投影数据为X射线透过目标对象以及散射校正板120获得的投影图像。
步骤S320:分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据。
本步骤中,各标准投影图像包括当前标准投影图像、前一张标准投影图像以及后一张标准投影图像。由于散射校正板加目标物体投影数据是X射线透过目标对象以及散射校正板120获得的投影图像,X射线先穿过散射校正板120,穿过金属球部分的射线被吸收,剩余的X射线再穿过目标物体产生散射,则金属球在探测器130上相应的位置得到的是散射数据。同样的,若目标物体放置于散射校正板120与射线源110之间,X射线穿过目标物体以后,散射分布已经形成,再穿过散射校正板120,金属球吸收初始X射线,则金属球在探测器130上相应的位置得到的是散射X射线强度,即散射数据。
根据散射校正板中的金属球在各个标准投影图像中的第一中心位置,获取各个标准投影图像中该第一中心位置的第一散射数据;例如,对于当前标准投影图像,根据散射校正板中的金属球在当前标准投影图像的中心位置,获取当前标准投影图像第一侧的散射校正板加目标物体投影数据在该中心位置上的散射数据;对于前一张标准投影图像,根据散射校正板中的金属球在前一张标准投影图像的中心位置,获取前一张标准投影图像第二侧的散射校正板加目标物体投影数据在该中心位置上的散射数据;对于后一张标准投影图像,根据散射校正板中的金属球在后一张标准投影图像的中心位置,获取后一张标准投影图像中散射校正板加目标物体投影数据在该中心位置上的散射数据。
步骤S330:分别根据各标准投影图像中的第一散射数据估算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各第一散射数据以及各第二散射数据获取当前标准投影图像中第一侧的散射数据以及前一张标准投影图像和后一张标准投影图像中第二侧的散射数据。
本步骤中,各标准投影图像包括当前标准投影图像、前一张标准投影图像以及后一张标准投影图像。对于每张标准投影图像,根据散射校正板加目标物体投影数据中已知的第一中心位置的第一散射数据,估算散射校正板加目标物体投影数据所在一侧的其他位置的散射数据,从而获得标准投影图像上散射校正板加目标物体投影数据所在一侧的所有位置的散射数据;例如,对于当前标准投影图像,根据当前标准投影图像中第一侧上已知的第一中心位置上的散射数据,估算当前标准投影图像中第一侧其他位置上的散射数据,从而获得当前标准投影图像中第一侧上所有位置的散射数据;对于前一张标准投影图像,根据前一张标准投影图像中第二侧上已知的第一中心位置的散射数据,估算前一张标准投影图像中第二侧上其他位置的散射数据,从而获得前一张标准投影图像中第二侧上所有位置的散射数据,对于后一张标准投影图像执行同样的操作,不再次赘述。
步骤S340:根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据;根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据。
本步骤中,由于散射信号强度变化缓慢,频率较低,可以根据当前标准投影图像前后相邻的两张标准投影图像,即根据前一张标准投影图像中第一侧的目标物体投影数据以及后一张标准投影图像中第一侧的目标物体投影数据,估算当前标准投影图像中第一侧的目标物体投影数据,根据前一张标准投影图像中第二侧的散射数据以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据。
步骤S350:将当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据第一目标投影图像以及第二目标投影图像获取当前标准投影图像散射校正后的目标投影图像。
本步骤中,将当前标准投影图像中第一侧的目标物体投影数据与当前标准投影图像中第一侧的散射数据相减获得当前标准投影图像中第一侧散射校正后的目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与当前标准投影图像中第二侧的散射数据相减获得当前标准投影图像中第二侧散射校正后的目标投影图像,从而获得当前标准投影图像完整的散射校正后的目标投影图像。
上述计算机断层图像的散射校正方法中,当前标准投影图像可以利用前一张标准投影图像以及后一张标准投影图像的标准投影图像数据进行散射校正,在CBCT系统中,在扫描一周均匀采集多张的标准投影图像后,对与每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,能够有效避免重建的CT图像形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量。
在其中一个实施例中,获取散射校正板以及目标物体的当前标准投影图像的步骤之前,还包括以下步骤:获取散射校正板的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像;其中,当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像分别与当前标准投影图像、前一张标准投影图像以及后一张标准投影图像的扫描角度相同;当前散射校正板投影图像中第一侧包括金属球投影图像,前一张散射校正板投影图像以及后一张散射校正板投影图像中第二侧包括金属球投影图像;根据各散射校正板投影图像中的金属球投影图像获取金属球在各散射校正板投影图像中的第二中心位置;分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据的步骤之前,还包括以下步骤:根据各第二中心位置确定金属球在各标准投影图像中的第一中心位置。
本实施例是第一中心位置S的获取方法,散射校正板投影图像是指在不放置目标对象的情况下,锥束X线束仅对散射校正板120进行扫描获得的标准投影图像;分别获取与当前标准投影图像、前一张标准投影图像以及后一张标准投影图像的扫描角度以及散射校正板放置位置均相同的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像,当前散射校正板投影图像中的金属球的投影位置与当前标准投影图像中的金属球的投影位置是一样的,同样的,前一张散射校正板投影图像与前一张标准投影图像的中金属球的投影位置是一样的,后一张散射校正板投影图像与后一张标准投影图像的中金属球的投影位置是一样的,通过计算散射校正板图像中散射校正板中金属球的第二中心位置,可估计金属球在各标准投影图像中的第一中心位置。
进一步的,在其中一个实施例中,根据各散射校正板投影图像中的金属球投影数据获取金属球在各散射校正板投影图像中的第二中心位置的步骤,包括以下步骤:从散射校正板投影图像中分割出金属球投影图像;检测金属球投影图像的边缘值,根据边缘值获取第二中心位置。
本实施例中,对散射校正板投影图像进行图像分割处理,分割散射校正板投影图像中的金属球投影图像,对金属球投影图像进行边缘检测,根据边缘值确定第二中心位置。
可选的,图像分割处理可以利用阈值分割法或者随机游走分割法进行图像分割处理。
以阈值分割法为例,对本实施例中的第二中心位置获取方法进行说明。检测每一张散射校正板投影图像中像素点的像素值,若像素值小于等于预设阈值,则该像素点判定为金属信息,并将该像素点的像素值重新设置像素值为1;像素值大于预设阈值,则该像素点判定为非金属信息,并将该像素点的像素值重新设置像素值为0,这样一张散射校正板投影图像中每一个像素点的像素值只含有0和1值,1值代表金属区域,0代表非金属区域。然后,对进行图像分割处理后的散射校正板投影图像进行逐行检测,若检测到相邻像素点的像素值分别为0值和1值,记录1值所对应的像素点的列位置c1,再继续逐行检测,直到检测到相邻的像素点的像素值分别为1值和0值,记录1值所对应的像素点的列位置c2,将c1和c2相加并除以2,得到该金属球的中心点位于散射校正板投影图像的横坐标。同理,对进行图像分割处理后的散射校正板投影图像进行逐列检测,直到检测到相邻像素点的像素值分别为0值和1值,记录1值所对应的像素点的行位置r1,再继续逐列检测,直到检测到相邻像素点的像素值为1值和0值,记录1值所对应的像素点的行位置r2,将r1和r2相加并除以2,得到金属球的中心点位于散射校正板投影图像像素矩阵的纵坐标,这样就得到了一个金属球的在散射校正板投影图像中的中心点位置坐标,其它金属球的中心点位置坐标同理计算,此处不一一赘述。
在其中一个实施例中,分别根据各第一散射数据估算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据的步骤,包括以下步骤:分别计算其他位置坐标到与其他位置坐标相邻的第一中心位置的距离值;根据各距离值计算各相邻的第一中心位置对应的散射数据的权重因子;根据各相邻的第一中心位置对应的散射数据以及对应的权重因子计算其他位置坐标的第二散射数据。
本实施例中,对于任一待估算的其他位置,利用与其相邻的已知散射数据计算待估算的其他位置坐标上的散射数据;具体的,获取与待估算的其他位置相邻的第一中心位置,计算待估算的其他位置到与其相邻的第一中心位置的距离值,根据距离值获取相邻的第一中心位置上散射数据的权重因子,然后将不同权重的相邻的第一中心位置上散射数据相加,获得待估算的其他位置的第二散射数据,剩余的其他位置坐标的第二散射数据计算同理,从而可以获得标准投影图像某一侧全部位置上的散射数据。
参见图4,其中点410为一个待估算的其他位置,获取与点410相邻的4个第一中心位置(点421、点422、点423以及点424),分别计算点410到点421、点422、点423以及点424的距离值,根据距离值获取这4个第一中心位置上散射数据的权重因子,其中,权重因子与距离值成反比关系,最后将不同权重的4个第一中心位置上散射数据相加,即为该待估算的其他位置的散射数据的值。
在其中一个实施例中,根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据的步骤,包括以下步骤:获取前一张标准投影图像中第一侧的第一待估算位置的第一目标物体投影数据,以及后一张标准投影图像中第一侧的第一待估算位置的第二目标物体投影数据;根据第一目标物体投影数据以及第二目标物体投影数据获取当前标准投影图像中第一侧的第一待估算位置的目标物体投影数据。
本实施例中,对于当前标准投影图像中第一侧的待估算位置的目标物体投影数据,可以对前一张标准投影图像中第一侧的待估算位置的目标物体投影数据以及后一张标准投影图像中第一侧的待估算位置的目标物体投影数据利用插值法估算而得,其中,插值法可以为线性插值法、多项式插值法或者B样条插值法。
以线性插值法为例,获取在前一张标准投影图像中第一侧的待估算位置的目标物体投影数据,以及后一张标准投影图像中第一侧的待估算位置的目标物体投影数据,赋予这两个目标投影数据数值是0.5的权重因子,估算当前标准投影图像中第一侧的待估算位置的目标物体投影数据。
在其中一个实施例中,根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据的步骤,包括以下步骤:获取前一张标准投影图像中第二侧的第二待估算位置的第三散射数据,以及后一张标准投影图像中第二侧的第二待估算位置的第四散射数据;根据第三散射数据以及第四散射数据获取当前标准投影图像中第二侧的第二待估算位置的散射数据。
本实施例中,对于当前标准投影图像中第二侧的待估算位置的散射数据,可以对前一张标准投影图像中第二侧的待估算位置的散射数据以及后一张标准投影图像中第二侧的待估算位置的散射数据利用插值法估算而得,其中,插值法可以为线性插值法、多项式插值法或者B样条插值法。
以线性插值法为例,获取在前一张标准投影图像中第二侧的待估算位置的散射数据,以及后一张标准投影图像中第二侧的待估算位置的散射数据,赋予这两个目标投影数据数值是0.5的权重因子,估算当前标准投影图像中第二侧的待估算位置的散射数据。
在其中一个实施例中,获得当前标准投影图像散射校正后的目标投影图像的步骤之后,还包括以下步骤:对多张目标投影图像进行CT图像重建。
本实施例中,CBCT系统中在扫描一周均匀采集多张的标准投影图像后,对与每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正后,对多张散射校正后的目标投影图像进行CT重建,获得无散射伪影的CT图像。
参见图5,图5为本发明另一个实施例中计算机断层图像的散射校正方法的流程图,计算机断层图像的散射校正方法,包括以下步骤:
步骤S510:获取散射校正板的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像;其中,当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像分别与当前标准投影图像、前一张标准投影图像以及后一张标准投影图像的扫描角度相同;当前散射校正板投影图像中第一侧包括金属球投影图像,前一张散射校正板投影图像以及后一张散射校正板投影图像中第二侧包括金属球投影图像。
步骤S520:从散射校正板投影图像中的金属球投影图像中分割出金属球投影图像;检测金属球投影图像的边缘值,根据边缘值获取金属球在各散射校正板投影图像中的第二中心位置。
步骤S530:获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,当前标准投影图像、前一张标准投影图像以及后一张标准投影图像间的扫描角度间隔相等;当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,当前标准投影图像中第二侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第一侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据。
步骤S540:根据各第二中心位置确定金属球在各标准投影图像中的第一中心位置。
步骤S550:分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据。
步骤S560:分别计算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标到与其他位置坐标相邻的第一中心位置的距离值;根据各距离值计算各相邻的第一中心位置对应的散射数据的权重因子;根据各相邻的第一中心位置对应的散射数据以及对应的权重因子计算其他位置坐标的第二散射数据。
步骤S570:分别根据各第一散射数据以及各第二散射数据获取当前标准投影图像第一侧的散射数据以及前一张标准投影图像和后一张标准投影图像中第二侧的散射数据。
步骤S580:根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据;根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据。
步骤S590:将当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据第一目标投影图像以及第二目标投影图像获取当前标准投影图像散射校正后的目标投影图像。
上述计算机断层图像的散射校正方法中,当前标准投影图像可以利用前一张标准投影图像以及后一张标准投影图像的投影图像数据进行散射校正,在CBCT系统中,在扫描一周均匀采集多张的标准投影图像后,对与每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,能够有效避免重建的CT图像形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量,在保证CBCT图像质量的同时最大限度的降低目标物体的辐射剂量。相较于传统的基于软件的校正方法,例如反卷积法以及蒙特卡罗模拟算法,反卷积法计算速度虽快,但散射核设计难度大导致散射估计不准确,而蒙特卡罗模拟算法对计算量要求高,散射校正效率低,难以在临床应用;而本实施例中计算机断层图像的散射校正方法,散射校正过程简单,计算量低,散射校正效率。
为了更加清晰本发明技术方案的效果,下面结合一个在锥形束计算机断层成像系统中使用本发明的应用实例来进一步阐述。
(1)设置曝光条件,其中包括设置X射线球管电压、X射线球管电流、锥形束计算机断层成像系统扫描一周需要均匀采集投影数据的张数,每张投影数据对应一个扫描角度。
具体的,X射线球管电压可以设置为90kV,X线球管电流可以设置为5mA,锥形束计算机断层成像系统扫描一周均匀采集投影数据的张数可以设置为360张,即每间隔1°采集1张投影数据。
(2)制作散射校正板并将散射校正板置于扫描物与探测器之间。
具体的,散射校正板平面与探测器平面平行,散射校正板的设计如图2所示,其中,散射校正板可以设计为厚度2mm的矩形薄板,水平长为6.5cm和竖直宽为13cm,材料为有机玻璃,散射校正板中嵌入一系列直径为2mm的金属铅球,在薄板中呈网格交叉点状分布,相邻两个金属铅球心之间的距离为5mm。
(3)在不放置目标物体的情况下,对散射校正板进行扫描获取散射校正板投影图像;其中,在扫描奇数张散射校正板投影图像时,散射校正板沿着其平面水平方向偏置到第一侧,使得散射校正板在探测器上的投影占据探测器面积的一半;在扫描偶数张散射校正板投影图像时,散射校正板沿着其平面水平方向偏置到第二侧,使得散射校正板在探测器上的投影占据探测器面积的另一半。
具体的,扫描奇数张散射校正板投影图像时,比如第1,3,5,…,359张,散射校正板沿着其平面水平方向左侧偏置,使得散射校正板在探测器上的投影占据探测器面积的左一半;扫描偶数张散射校正板投影图像时,比如第2,4,6,…,360张,散射校正板沿着其平面水平方向右侧偏置,使得散射校正板在探测器上的投影占据探测器面积的右一半。
可选的,扫描奇数张散射校正板投影图像和扫描偶数张散射校正板投影图像时可以交替进行,这时需要散射校正板交替水平偏移;也可以扫描完所有奇数张的散射校正板投影图像后,再进行散射校正板水平偏移,再扫描完所有偶数张的散射校正板投影图像。
(4)对所有散射校正板投影图像中的金属球投影图像进行图像分割处理,分割出所有金属球投影图像,并计算每张散射校正板投影图像中每个金属球投影图像的第一中心点位置;
(5)放置目标物体,对散射校正板以及目标物体进行扫描获取标准投影图像;扫描奇数张标准投影图像时,散射校正板沿着其平面水平方向偏置到第一侧,使得散射校正板在探测器上的投影占据探测器面积的一半;扫描偶数张标准投影图像时,散射校正板沿着其平面水平方向偏置到第一侧,使得散射校正板在探测器上的投影占据探测器面积的另一半。
具体的,本步骤和步骤(3)相似,只是扫描对象为散射校正板以及目标物体。扫描奇数张的标准投影图像时,比如第1,3,5,…,359张,散射校正板沿着其平面水平方向左侧偏置,使得散射校正板在探测器上的投影占据探测器面积的左一半,获得的标准投影图像中左侧为散射校正板加目标物体投影数据,右侧为目标物体投影数据。扫描偶数张的标准投影图像时,比如第2,4,6,…,360张,散射校正板沿着其平面水平方向右侧偏置,使得散射校正板在探测器上的投影占据探测器面积的右一半,获得的标准投影图像中右侧为散射校正板加目标物体投影数据,左侧为目标物体投影数据。
可选的,扫描奇数张的标准投影图像和扫描偶数张的标准投影图像时可以交替进行,这时需要散射校正板交替水平偏移;也可以扫描完所有奇数张的标准投影图像后,再进行散射校正板水平偏移,再扫描完所有偶数张的标准投影图像。
(6)根据散射校正板投影图像中每个金属球投影图像的第一中心点位置,确定在对应角度的标准投影图像中金属球在该标准投影图像的第一中心位置,并获取该第一中心位置上的散射数据,例如,对于第3张标准投影图像,它的左半侧为散射校正板加目标物体投影数据,右半侧为目标物体投影数据;对于第3张散射校正板投影图像,它的左半侧为金属球投影图像,且金属球投影图像的中心点位置S已知,根据中心点位置S可以确定在第3张标准投影图像左半侧上金属球投影的第一中心位置;对于其他标准投影图像同理,此处不一一赘述。
(7)对每张标准投影图像中散射校正板加目标物体投影数据的一侧,根据该标准投影图像的第一中心位置上的第一散射数据,确定该标准投影图像中散射校正板加目标物体投影数据的一侧其他位置上的第二散射数据,从而能够获得标准投影图像中某一侧的全部散射数据,形成标准投影图像中一侧为散射数据另一侧为目标物体投影数据。例如,对于第3张标准投影图像左半侧上金属球投影的第一中心位置以及散射数据已知,根据已知的散射数据估算第3张标准投影图像左半侧其他位置上的散射数据,从而获得第3张标准投影图像左半侧全部的散射数据,此时,形成了第3张标准投影图像左半侧为散射数据,右半侧为目标物体投影数据;对于其他标准投影图像同理,此处不一一赘述。
(8)对于任意一张标准投影图像,该标准投影图像上只含有散射数据的一侧的目标物体投影数据,由该标准投影图像的前一张标准投影图像以及后一张标准投影图像的同侧的目标物体投影数据估算获得。例如,对于第3张标准投影图像,它的左半侧为散射数据,右半侧为目标物体投影数据,此时,第3张标准投影图像的左半侧的目标物体投影数据可以根据第2张标准投影图像左半侧的目标物体投影数据和第4张标准投影图像左半侧的目标物体投影数据估算获得;同理的,对于第2张标准投影图像,它的左半侧为目标物体投影数据,右半侧为散射数据,此时,第2张标准投影图像的右半侧的目标物体投影数据可以根据第1张标准投影图像右半侧的目标物体投影数据和第3张标准投影图像右半侧的目标物体投影数据估算获得;对于其他标准投影图像同理,此处不一一赘述,应该说明的是,第1张标准投影图像的前一张标准投影图像为第360张标准投影图像,第360张标准投影图像的后一张标准投影图像为第1张标准投影图像。
(9)对于任意一张标准投影图像,该标准投影图像上只含有目标物体投影数据的一侧的散射数据,由该标准投影图像的前一张标准投影图像以及后一张标准投影图像同侧的散射数据估算获得。本步骤中与步骤(8)类似,例如,对于第3张标准投影图像,它的左半侧为散射数据,右半侧为目标物体投影数据,此时,第3张标准投影图像的右半侧的散射数据可以根据第2张标准投影图像右半侧的散射数据和第4张标准投影图像右半侧的散射数据估算获得,对于其他标准投影图像同理,此处不一一赘述。
(10)对每张标准投影图像中的目标物体投影数据与对应的散射数据进行相减,得到校正后的目标投影图像。
(11)对校正后的目标投影图像进行CT图像重建,得到无散射伪影的CT图像。
本实施例中,CBCT系统中在扫描一周均匀采集多张的标准投影图像后,对于每一张标准投影图像利用其前一张标准投影图像以及后一张标准投影图像的图像数据进行散射校正,例如,对于第2张标准投影图像可以利用第1张标准投影图像以及第3张标准投影图像的图像数据进行散射校正,利用本方法能够有效避免重建的CT图像中形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量,在保证CBCT图像质量的同时最大限度的降低目标物体的辐射剂量。
应该理解的是,虽然图3或图5的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图3或图5中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在其中一个实施例中,如图6所示,图6为本发明一个实施例中计算机断层图像的散射校正装置的结构示意图,提供了一种计算机断层图像的散射校正装置,包括:
标准投影图像获取模块610,用于获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,当前标准投影图像、前一张标准投影图像以及后一张标准投影图像间的扫描角度间隔相等;当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,当前标准投影图像中第二侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第一侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
第一散射数据获取模块620,用于分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
第二散射数据获取模块630,用于分别根据各标准投影图像中的第一散射数据估算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各第一散射数据以及各第二散射数据获取当前标准投影图像中第一侧的散射数据以及前一张标准投影图像和后一张标准投影图像中第二侧的散射数据;
当前投影数据估算模块640,用于根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据;根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据;
目标投影图像获取模块650,用于将当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据第一目标投影图像以及第二目标投影图像获取当前标准投影图像散射校正后的目标投影图像。
上述计算机断层图像的散射校正装置,能够有效避免重建的CBCT图像中形成散射伪影,且只需要对目标物体扫描一次,不增加目标物体的辐射剂量,在保证CBCT图像质量的同时最大限度的降低目标物体的辐射剂量。
如图7所示,图7为本发明另一个实施例中计算机断层图像的散射校正装置的结构示意图,本实施例中,计算机断层图像的散射校正装置还包括第一中心位置获取模块660,用于获取散射校正板的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像;其中,当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像分别与当前标准投影图像、前一张标准投影图像以及后一张标准投影图像的扫描角度相同;当前散射校正板投影图像中第一侧包括金属球投影图像,前一张散射校正板投影图像以及后一张散射校正板投影图像中第二侧包括金属球投影图像;根据各散射校正板投影图像中的金属球投影图像获取金属球在各散射校正板投影图像中的第二中心位置;根据各第二中心位置确定金属球在各标准投影图像中的第一中心位置。
在其中一个实施例中,第一中心位置获取模块660用于从散射校正板投影图像中分割出金属球投影图像;检测金属球投影图像的边缘值,根据边缘值获取第二中心位置。
在其中一个实施例中,第二散射数据获取模块630用于分别计算其他位置坐标到与其他位置坐标相邻的第一中心位置的距离值;根据各距离值计算各相邻的第一中心位置对应的散射数据的权重因子;根据各相邻的第一中心位置对应的散射数据以及对应的权重因子计算其他位置坐标的第二散射数据。
在其中一个实施例中,当前投影数据估算模块640用于获取前一张标准投影图像中第一侧的第一待估算位置的第一目标物体投影数据,以及后一张标准投影图像中第一侧的第一待估算位置的第二目标物体投影数据;根据第一目标物体投影数据以及第二目标物体投影数据获取当前标准投影图像中第一侧的第一待估算位置的目标物体投影数据。
在其中一个实施例中,当前投影数据估算模块640用于获取前一张标准投影图像中第二侧的第二待估算位置的第三散射数据,以及后一张标准投影图像中第二侧的第二待估算位置的第四散射数据;根据第三散射数据以及第四散射数据获取当前标准投影图像中第二侧的第二待估算位置的散射数据。
如图8所示,图8为本发明又一个实施例中计算机断层图像的散射校正装置的结构示意图,本实施例中,计算机断层图像的散射校正装置还包括CT图像重建模块670,用于对多张目标投影图像进行CT图像重建。
关于计算机断层图像的散射校正装置的具体限定可以参见上文中对于计算机断层图像的散射校正方法的限定,在此不再赘述。上述计算机断层图像的散射校正装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是终端,其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种计算机断层图像的散射校正方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏,该计算机设备的输入装置可以是显示屏上覆盖的触摸层,也可以是计算机设备外壳上设置的按键、轨迹球或触控板,还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解,图9中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在其中一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,当前标准投影图像、前一张标准投影图像以及后一张标准投影图像间的扫描角度间隔相等;当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,当前标准投影图像中第二侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第一侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各标准投影图像中的第一散射数据估算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各第一散射数据以及各第二散射数据获取当前标准投影图像中第一侧的散射数据以及前一张标准投影图像和后一张标准投影图像中第二侧的散射数据;
根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据;根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据;
将当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据第一目标投影图像以及第二目标投影图像获取当前标准投影图像散射校正后的目标投影图像。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,当前标准投影图像、前一张标准投影图像以及后一张标准投影图像间的扫描角度间隔相等;当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,当前标准投影图像中第二侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第一侧包括目标物体投影数据,前一张标准投影图像以及后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各标准投影图像中的第一散射数据估算各标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各第一散射数据以及各第二散射数据获取当前标准投影图像中第一侧的散射数据以及前一张标准投影图像和后一张标准投影图像中第二侧的散射数据;
根据前一张标准投影图像以及后一张标准投影图像中第一侧的目标物体投影数据估算当前标准投影图像中第一侧的目标物体投影数据;根据前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算当前标准投影图像中第二侧的散射数据;
将当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据第一目标投影图像以及第二目标投影图像获取当前标准投影图像散射校正后的目标投影图像。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:获取散射校正板的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像;其中,当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像分别与当前标准投影图像、前一张标准投影图像以及后一张标准投影图像的扫描角度相同;当前散射校正板投影图像中第一侧包括金属球投影图像,前一张散射校正板投影图像以及后一张散射校正板投影图像中第二侧包括金属球投影图像;根据各散射校正板投影图像中的金属球投影图像获取金属球在各散射校正板投影图像中的第二中心位置;分别根据散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各标准投影图像中散射校正板加目标物体投影数据的第一散射数据的步骤之前,还包括以下步骤:根据各第二中心位置确定金属球在各标准投影图像中的第一中心位置。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:从散射校正板投影图像中分割出金属球投影图像;检测金属球投影图像的边缘值,根据边缘值获取第二中心位置。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:分别计算其他位置坐标到与其他位置坐标相邻的第一中心位置的距离值;根据各距离值计算各相邻的第一中心位置对应的散射数据的权重因子;根据各相邻的第一中心位置对应的散射数据以及对应的权重因子计算其他位置坐标的第二散射数据。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:获取前一张标准投影图像中第一侧的第一待估算位置的第一目标物体投影数据,以及后一张标准投影图像中第一侧的第一待估算位置的第二目标物体投影数据;根据第一目标物体投影数据以及第二目标物体投影数据获取当前标准投影图像中第一侧的第一待估算位置的目标物体投影数据。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:获取前一张标准投影图像中第二侧的第二待估算位置的第三散射数据,以及后一张标准投影图像中第二侧的第二待估算位置的第四散射数据;根据第三散射数据以及第四散射数据获取当前标准投影图像中第二侧的第二待估算位置的散射数据。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:对多张目标投影图像进行CT图像重建。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种计算机断层图像的散射校正方法,其特征在于,包括以下步骤:
获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
分别根据各所述标准投影图像中的第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像中第一侧的散射数据以及所述前一张标准投影图像和所述后一张标准投影图像中第二侧的散射数据;
根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减获得第一目标投影图像,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减获得第二目标投影图像,根据所述第一目标投影图像以及所述第二目标投影图像获取所述当前标准投影图像散射校正后的目标投影图像。
2.根据权利要求1所述的计算机断层图像的散射校正方法,其特征在于,所述获取散射校正板以及目标物体的当前标准投影图像的步骤之前,还包括以下步骤:
获取所述散射校正板的当前散射校正板投影图像、前一张散射校正板投影图像以及后一张散射校正板投影图像;其中,所述当前散射校正板投影图像、所述前一张散射校正板投影图像以及所述后一张散射校正板投影图像分别与所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像的扫描角度相同;所述当前散射校正板投影图像中第一侧包括金属球投影图像,所述前一张散射校正板投影图像以及所述后一张散射校正板投影图像中第二侧包括金属球投影图像;
根据各所述散射校正板投影图像中的金属球投影图像获取所述金属球在各所述散射校正板投影图像中的第二中心位置;
所述分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据的步骤之前,还包括以下步骤:
根据各所述第二中心位置确定所述金属球在各标准投影图像中的第一中心位置。
3.根据权利要求2所述的计算机断层图像的散射校正方法,其特征在于,所述根据各所述散射校正板投影图像中的金属球投影数据获取所述金属球在各所述散射校正板投影图像中的第二中心位置的步骤,包括以下步骤:
从所述散射校正板投影图像中分割出金属球投影图像;
检测所述金属球投影图像的边缘值,根据所述边缘值获取所述第二中心位置。
4.根据权利要求1所述的计算机断层图像的散射校正方法,其特征在于,所述分别根据各所述标准投影图像中的第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据的步骤,包括以下步骤:
分别计算所述其他位置坐标到与所述其他位置坐标相邻的第一中心位置的距离值;
根据各所述距离值计算各所述相邻的第一中心位置对应的散射数据的权重因子;
根据各所述相邻的第一中心位置对应的散射数据以及对应的权重因子计算所述其他位置坐标的第二散射数据。
5.根据权利要求1所述的计算机断层图像的散射校正方法,其特征在于,所述根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据的步骤,包括以下步骤:
获取所述前一张标准投影图像中第一侧的第一待估算位置的第一目标物体投影数据,以及所述后一张标准投影图像中第一侧的所述第一待估算位置的第二目标物体投影数据;
根据所述第一目标物体投影数据以及第二目标物体投影数据获取所述当前标准投影图像中第一侧的所述第一待估算位置的目标物体投影数据。
6.根据权利要求1所述的计算机断层图像的散射校正方法,其特征在于,所述根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据的步骤,包括以下步骤:
获取所述前一张标准投影图像中第二侧的第二待估算位置的第三散射数据,以及所述后一张标准投影图像中第二侧的所述第二待估算位置的第四散射数据;
根据所述第三散射数据以及所述第四散射数据获取所述当前标准投影图像中第二侧的所述第二待估算位置的散射数据。
7.根据权利要求1所述的计算机断层图像的散射校正方法,其特征在于,所述获得所述当前标准投影图像散射校正后的目标投影图像的步骤之后,还包括以下步骤:
对多张目标投影图像进行CT图像重建。
8.一种计算机断层图像的散射校正装置,其特征在于,包括:
标准投影图像获取模块,用于获取散射校正板以及目标物体的当前标准投影图像、前一张标准投影图像以及后一张标准投影图像;其中,所述当前标准投影图像、所述前一张标准投影图像以及所述后一张标准投影图像间的扫描角度间隔相等;所述当前标准投影图像中第一侧包括散射校正板加目标物体投影数据,所述当前标准投影图像中第二侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第一侧包括目标物体投影数据,所述前一张标准投影图像以及所述后一张标准投影图像中第二侧包括散射校正板加目标物体投影数据;
第一散射数据获取模块,用于分别根据所述散射校正板中的金属球在各标准投影图像中的第一中心位置,获取各所述标准投影图像中散射校正板加目标物体投影数据的第一散射数据;
第二散射数据获取模块,用于分别根据各所述第一散射数据估算各所述标准投影图像中散射校正板加目标物体投影数据中其他位置坐标上的第二散射数据,根据各所述第一散射数据以及各所述第二散射数据获取所述当前标准投影图像第一侧的散射数据和所述前一张标准投影图像以及所述后一张标准投影图像中第二侧的散射数据;
当前投影数据估算模块,用于根据所述前一张标准投影图像以及所述后一张标准投影图像中第一侧的目标物体投影数据估算所述当前标准投影图像中第一侧的目标物体投影数据;根据所述前一张标准投影图像以及后一张标准投影图像中第二侧的散射数据估算所述当前标准投影图像中第二侧的散射数据;
目标投影图像获取模块,用于将所述当前标准投影图像中第一侧的目标物体投影数据与第一侧的散射数据相减,将所述当前标准投影图像中第二侧的目标物体投影数据与第二侧的散射数据相减,获得所述当前标准投影图像散射校正后的目标投影图像。
9.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述的计算机断层图像的散射校正方法的步骤。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的计算机断层图像的散射校正方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810456437.4A CN108542414B (zh) | 2018-05-14 | 2018-05-14 | 计算机断层图像的散射校正方法以及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810456437.4A CN108542414B (zh) | 2018-05-14 | 2018-05-14 | 计算机断层图像的散射校正方法以及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108542414A true CN108542414A (zh) | 2018-09-18 |
CN108542414B CN108542414B (zh) | 2019-08-23 |
Family
ID=63494817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810456437.4A Expired - Fee Related CN108542414B (zh) | 2018-05-14 | 2018-05-14 | 计算机断层图像的散射校正方法以及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108542414B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111001097A (zh) * | 2019-12-28 | 2020-04-14 | 上海联影医疗科技有限公司 | 一种放疗剂量评估系统、装置及存储介质 |
WO2020082280A1 (en) * | 2018-10-25 | 2020-04-30 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for scatter correction |
CN112184569A (zh) * | 2020-09-08 | 2021-01-05 | 杭州三坛医疗科技有限公司 | 图像修复方法和图像修复装置 |
CN113100802A (zh) * | 2021-04-16 | 2021-07-13 | 上海联影医疗科技股份有限公司 | 一种校正机械偏差的方法和系统 |
CN116188617A (zh) * | 2023-04-21 | 2023-05-30 | 有方(合肥)医疗科技有限公司 | Ct图像数据处理方法、装置及ct系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1635850A (zh) * | 2002-02-21 | 2005-07-06 | 罗切斯特大学 | X射线散射校正 |
US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
CN103578082A (zh) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | 一种锥束ct散射校正方法及系统 |
CN104161536A (zh) * | 2014-07-30 | 2014-11-26 | 西北工业大学 | 一种基于互补光栅的锥束ct散射校正方法及其装置 |
CN105615911A (zh) * | 2014-10-20 | 2016-06-01 | 北卡罗来纳大学教堂山分校 | 用于固定的数字化胸部断层融合成像的系统及相关方法 |
-
2018
- 2018-05-14 CN CN201810456437.4A patent/CN108542414B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1635850A (zh) * | 2002-02-21 | 2005-07-06 | 罗切斯特大学 | X射线散射校正 |
US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
CN103578082A (zh) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | 一种锥束ct散射校正方法及系统 |
CN104161536A (zh) * | 2014-07-30 | 2014-11-26 | 西北工业大学 | 一种基于互补光栅的锥束ct散射校正方法及其装置 |
CN105615911A (zh) * | 2014-10-20 | 2016-06-01 | 北卡罗来纳大学教堂山分校 | 用于固定的数字化胸部断层融合成像的系统及相关方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020082280A1 (en) * | 2018-10-25 | 2020-04-30 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for scatter correction |
US11087507B2 (en) | 2018-10-25 | 2021-08-10 | Shanghai United Imaging Healthcare. Co., Ltd. | System and method for scatter correction |
US11875433B2 (en) | 2018-10-25 | 2024-01-16 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for scatter correction |
CN111001097A (zh) * | 2019-12-28 | 2020-04-14 | 上海联影医疗科技有限公司 | 一种放疗剂量评估系统、装置及存储介质 |
CN111001097B (zh) * | 2019-12-28 | 2022-09-16 | 上海联影医疗科技股份有限公司 | 一种放疗剂量评估系统、装置及存储介质 |
CN112184569A (zh) * | 2020-09-08 | 2021-01-05 | 杭州三坛医疗科技有限公司 | 图像修复方法和图像修复装置 |
CN113100802A (zh) * | 2021-04-16 | 2021-07-13 | 上海联影医疗科技股份有限公司 | 一种校正机械偏差的方法和系统 |
CN113100802B (zh) * | 2021-04-16 | 2023-07-28 | 上海联影医疗科技股份有限公司 | 一种校正机械偏差的方法和系统 |
CN116188617A (zh) * | 2023-04-21 | 2023-05-30 | 有方(合肥)医疗科技有限公司 | Ct图像数据处理方法、装置及ct系统 |
CN116188617B (zh) * | 2023-04-21 | 2023-08-08 | 有方(合肥)医疗科技有限公司 | Ct图像数据处理方法、装置及ct系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108542414B (zh) | 2019-08-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108542414B (zh) | 计算机断层图像的散射校正方法以及装置 | |
Wang et al. | Acuros CTS: A fast, linear Boltzmann transport equation solver for computed tomography scatter–Part II: System modeling, scatter correction, and optimization | |
US9848844B2 (en) | Iterative reconstruction process | |
JP6438752B2 (ja) | 二重エネルギーctスキャン画像における物体識別方法 | |
US9589373B2 (en) | Monte carlo modeling of field angle-dependent spectra for radiographic imaging systems | |
JP6028804B2 (ja) | デジタル画像処理方法および撮影装置 | |
Abdoli et al. | Assessment of metal artifact reduction methods in pelvic CT | |
CN108606805B (zh) | 一种散射校正方法、装置及设备 | |
JP7254322B2 (ja) | 医用画像処理システムとその画像再構成方法のためのコリメータ | |
Altunbas et al. | A unified scatter rejection and correction method for cone beam computed tomography | |
Sun et al. | Correction for patient table‐induced scattered radiation in cone‐beam computed tomography (CBCT) | |
JP2015112482A (ja) | 二重エネルギーct造影剤強調されたスキャン画像内の石灰化部分を識別する方法 | |
JP2013003145A (ja) | 核画像におけるアーティファクト除去 | |
Xu et al. | Statistical iterative reconstruction to improve image quality for digital breast tomosynthesis | |
Michielsen et al. | Patchwork reconstruction with resolution modeling for digital breast tomosynthesis | |
CN108903964B (zh) | 计算机断层图像的散射校正方法以及装置 | |
JP6526428B2 (ja) | 医用画像処理装置、医用画像処理方法および医用画像診断装置 | |
Thanasupsombat et al. | A simple scatter reduction method in cone-beam computed tomography for dental and maxillofacial applications based on Monte Carlo simulation | |
US11538160B2 (en) | Systems and methods for scatter correction of image | |
Kim et al. | An additional tilted‐scan‐based CT metal‐artifact‐reduction method for radiation therapy planning | |
US7184062B2 (en) | Statistically-based image blending methods and systems for pasting multiple digital sub-images together | |
CN108154541A (zh) | 医学成像设备及其操作方法 | |
US10552992B2 (en) | Poly-energetic reconstruction method for metal artifacts reduction | |
CN113759412B (zh) | 获取束流形状和能量探测单元响应特征的方法、装置 | |
JP6147512B2 (ja) | 核医学診断装置、画像処理装置および画像再構成プログラム |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190823 |