CN107202805B - 基于卷积核的锥束ct散射伪影校正方法 - Google Patents
基于卷积核的锥束ct散射伪影校正方法 Download PDFInfo
- Publication number
- CN107202805B CN107202805B CN201710401166.8A CN201710401166A CN107202805B CN 107202805 B CN107202805 B CN 107202805B CN 201710401166 A CN201710401166 A CN 201710401166A CN 107202805 B CN107202805 B CN 107202805B
- Authority
- CN
- China
- Prior art keywords
- scattering
- photon
- point
- ray
- calculating
- 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
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000012937 correction Methods 0.000 title claims abstract description 23
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000002591 computed tomography Methods 0.000 claims abstract description 9
- 235000013405 beer Nutrition 0.000 claims description 7
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 230000003993 interaction Effects 0.000 abstract description 6
- 238000000333 X-ray scattering Methods 0.000 abstract description 4
- 239000000463 material Substances 0.000 abstract description 4
- 239000000126 substance Substances 0.000 abstract description 3
- 230000009471 action Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000007408 cone-beam computed tomography Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001436672 Bhatia Species 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- 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]
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/05—Investigating materials by wave or particle radiation by diffraction, scatter or reflection
- G01N2223/051—Investigating materials by wave or particle radiation by diffraction, scatter or reflection correcting for scatter
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种锥束CT散射伪影校正方法,特别是一种基于卷积核的锥束CT散射伪影校正方法,包括:对被测物进行CT扫描,得到投影数据;估算X射线的初始的光子数N0;利用初始光子数N0计算探测器上的光子数分布;计算X射线散射分布;从原始投影数据中扣除散射分布。本技术方案提出了一种新的卷积核,其中最关键的部分是散射核函数的求解,针对同材质均匀物体,分析X射线与物质相互作用的每个过程,并用数学公式加以描述,从而求得整个平面的散射分布。
Description
技术领域
本发明涉及一种锥束CT散射伪影校正方法,特别是一种基于卷积核的锥束CT散射伪影校正方法。
背景技术
锥束CT从出现以来就凭借着其精度高、成像速度快、剂量小等优势被广泛的应用于临床医学和工业无损检测中。但是受其大视野和使用大面积探测器的影响,散射问题已经成为了影响CBCT成像质量的主要原因,受到了越来越多的关注。为了解决CBCT系统中的散射引起的图像质量下降的问题,学者们提出了许多方法。根据方法中使用手段的不同,主要可以分为三类方法:硬件校正法、软件校正法、软硬件混合校正法。其中硬件校正法是在CT系统中添加额外的硬件,通过减少X射线来达到减少散射的目的,例如准直器、抗散射光栅等。这些方法对于散射的削弱均有一定的效果,但是存在散射校正不完全、与CBCT空间结构不相适应等问题。
软件法是在获得全部的投影数据后,利用数字图像处理算法估计散射光子的分布,并在投影数据中将其去除。常见的有卷积法、反卷积法和蒙特卡洛模拟法等。软件法不需要添加额外的硬件,实现简单,但是校正精度与卷积核的选取有很大关系。Shaw等人1987年提出卷积法,通过将投影数据与核函数卷积、加权后得到散射估值,再从投影数据中将其减去实现散射校正。现在核函数主流的求解方法是通过蒙特卡洛模拟笔形束穿过不同厚度的体膜,每个笔形束都用点扩展函数来表征,散射分布就是每个点扩展函数的累积响应。谢世鹏、M SUN、Esther Meyer、Nanmina Bhatia等人从2010年到2016年均有相关的文章。
混合法是硬件法与软件法的结合,在向CT系统中添加一些特制的硬件如BSA、调制板等后经过一定的处理达到散射校正的目的。混合法效果一般较好,但是需要特制的硬件和扫描次数。
发明内容
鉴于此,本发明提供一种基于卷积核的锥束CT散射伪影校正方法,提出了一种新的卷积核,其中最关键的部分是散射核函数的求解,针对同材质均匀物体,分析X射线与物质相互作用的每个过程,并用数学公式加以描述,从而求得整个平面的散射分布。
为了达到上述目的,本发明是通过以下技术方案实现的:
本发明提供一种基于卷积核的锥束CT散射伪影校正方法,包括:
对被测物进行CT扫描,得到投影数据;
估算X射线的初始的光子数N0;
利用初始光子数N0计算探测器上的光子数分布;
计算X射线散射分布;
从原始投影数据中扣除散射分布。
进一步地,估算X射线的初始的光子数N0,包括:
利用对被测物体进行CT扫描时的电压、电流和X射线能谱估算X射线的初始的光子数N0。
进一步地,利用初始光子数N0计算探测器上的光子数分布,包括:
计算被测物体的入射点到散射点的光子衰减;
计算经过散射点后光子散射角的概率分布;
计算从散射点到探测器的光子衰减。
进一步地,计算被测物体的入射点到散射点的光子衰减,包括:
根据比尔定律,被测物体的入射点到散射点的光子衰减为:
其中,N1即经过物体对X射线的衰减后光子数,μ为物体对X射线的衰减系数。
进一步地,计算经过散射点后光子散射角的概率分布,包括:
X射线在散射点发生散射,X射线方向发生偏转,偏转方向与原方向的夹角θ根据初始能量的大小呈现一定的概率分布,公式如下:
其中,E0为入射光子能量,E1为入射光子能量,α=E0/m0c2,m0为静态电子质量,c为真空中光速,re为经典电子半径;
以探测器为中心点建立坐标系,散射点在探测器上对应的坐标为(x’,y’),则有
根据上述公式可以得到散射光子偏转θ时的概率
进一步地,计算从散射点到探测器的光子衰减,包括:
根据比尔定律,计算被测物体的入射点到散射点的光子衰减为:
进一步地,计算X射线散射分布,包括:
计算散射点对探测器上任一点造成的散射光子数
Nf=ρ×N3×P×ΔΩ,
其中,ρ为一调节常数,表示电子密度等对光子数的影响;N3表示X射线与物体作用过程中光子数;P表示散射光子偏转θ的概率;ΔΩ表示探测器对a点的张角。
本发明提供一种基于卷积核的锥束CT散射伪影校正方法,具有如下有益效果:
本发明基于锥束CT散射伪影矫正方法的卷积法提出了一种新的卷积核,其中最关键的部分是散射核函数的求解,针对同材质均匀物体,分析X射线与物质相互作用的每个过程,并用数学公式加以描述,从而求得整个平面的散射分布。
附图说明
图1为本发明实施例所提供的基于卷积核的锥束CT散射伪影校正方法的流程示意图;
图2为本发明实施例所提供的基于卷积核的锥束CT散射伪影校正方法另一种实施方式的流程示意图;
图3为X射线与被测物体相互作用示意图;
图4为以探测器中心点为极坐标下散射光子的角分布。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
以下结合具体情况说明本发明的示例性实施例:
请参考图1,图1为本发明实施例所提供的基于卷积核的锥束CT散射伪影校正方法的流程示意图;本实施例提供一种基于卷积核的锥束CT散射伪影校正方法,包括以下步骤:
步骤S101、对被测物进行CT扫描,得到投影数据;
步骤S102、估算X射线的初始的光子数N0;
步骤S103、利用初始光子数N0计算探测器上的光子数分布;
步骤S104、计算X射线散射分布;
步骤S105、从原始投影数据中扣除散射分布。
请参考图2,图2为本发明实施例所提供的基于卷积核的锥束CT散射伪影校正方法另一种实施方式的流程示意图;本实施例提供一种基于卷积核的锥束CT散射伪影校正方法,包括:
步骤S201、对被测物进行CT扫描,得到投影数据。
步骤S202、利用对被测物体进行CT扫描时的电压、电流和X射线能谱估算X射线的初始的光子数N0。
步骤S203、计算被测物体的入射点到散射点的光子衰减;
根据比尔定律,被测物体的入射点到散射点的光子衰减为:
其中,N1即经过物体对X射线的衰减后光子数,μ为物体对X射线的衰减系数。
步骤S204、计算经过散射点后光子散射角的概率分布;
其中,X射线在散射点发生散射,X射线方向发生偏转,偏转方向与原方向的夹角θ根据初始能量的大小呈现一定的概率分布,公式如下:
其中,E0为入射光子能量,E1为入射光子能量,α=E0/m0c2,m0为静态电子质量,c为真空中光速,re为经典电子半径;
以探测器为中心点建立坐标系,散射点在探测器上对应的坐标为(x’,y’),则有
根据上述公式可以得到散射光子偏转θ时的概率
步骤S205、计算从散射点到探测器的光子衰减;
然后,根据比尔定律,计算被测物体的入射点到散射点的光子衰减为:
步骤S206、计算X射线散射分布;
在本实施例中,计算散射点对探测器上任一点造成的散射光子数
Nf=ρ×N3×P×ΔΩ,
其中,ρ为一调节常数,表示电子密度等对光子数的影响;N3表示X射线与物体作用过程中光子数;P表示散射光子偏转θ的概率;ΔΩ表示探测器对a点的张角。
步骤S207、从原始投影数据中扣除散射分布。
该方法的具体操作过程如下:
首先对被测物体进行CT扫描,利用得到的数据粗略的估计出物体的轮廓信息,然后利用扫描物体时的电压和电流以及相关的能谱估算出初始的光子数N0,并用求得物体轮廓信息和光子数依据所给公式计算出探测器上对应的光子分布。结合投影数据,建立线性方程组求解初始光子数,利用初始光子数解得散射分布,从初始的投影数据中口出散射分布。
X射线与物体相互作用的过程可以分为三个部分,具体过程如图3,图中a为散射点,b为所选取的截面。首先是从入射点到发生散射点a之前的部分L1段,此时物体对X射线起衰减作用,该部分遵从比尔定律,假设X射线初始的光子数为N0,那么有:
N1即经过物体对X射线的衰减后光子数,μ为物体对X射线的衰减系数。在同材质的均物体中,μ保持恒定不变。
其次X射线在a点发生散射,此时X射线方向发生偏转,能量也随之变化。X射线偏转方向与原方向的夹角根据初始能量的大小呈现一定的概率分布。Klein-Nishina公式描述的是实际散射物理过程中单个光子发生散射其康普顿散射微分截面,即这一概率分布。其具体公式如下:
其中E0为入射光子能量,α=E0/m0c2,m0为静态电子质量,c为真空中光速,re为经典电子半径,Te=e/m0c2=2.8×10-13cm。图4给出了极坐标表示的单个电子康普顿散射微分截面与入射光子能量与散射角度的关系,从图中可以看出,不同能量下散射发生时,散射光子偏转不同方向的概率,并且随着入射光子能量的增加,散射光子偏转方向趋向于前向分布。
K-N公式描述的就是在入射能量为E0的条件下,散射角度发生变化的概率。
在只考虑一次散射的情况下,在a点发生散射后,角度发生偏转,以探测器中心点建立坐标系,a点在探测器上对应的坐标为(x’,y’)则偏转角:
根据K-N公式,可以计算出散射光子偏转θ时的概率P:
同时发生散射前后时光子数的变化可表示为:
发生散射后,物体对散射光子的衰减同样服从比尔定律:
N=N2e-μL;
散射点a对探测器上任一点造成的散射光子数可以表示为:
Nf=ρ×N3×P×ΔΩ;
其中ρ为一调节常数,表示电子密度等对光子数的影响。N3表示X射线与物体作用过程中光子数。P表示散射光子偏转θ的概率。ΔΩ表示探测器对a点的张角,可以认为是一个常数。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上通过具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (4)
1.一种基于卷积核的锥束CT散射伪影校正方法,其特征在于,包括:
对被测物进行CT扫描,得到投影数据;
估算X射线的初始的光子数N0;
利用初始光子数N0计算探测器上的光子数分布,包括:计算被测物体的入射点到散射点的光子衰减;计算经过散射点后光子散射角的概率分布;计算从散射点到探测器的光子衰减;其中,所述计算经过散射点后光子散射角的概率分布,具体为:以探测器为中心点建立坐标系,散射点在探测器上对应的坐标为(x’,y’),则有θ为X射线在散射点发生散射,X射线方向发生偏转时,偏转方向与原方向的夹角;L1表示入射点到散射点的距离;
根据下式可以得到散射光子偏转θ时的概率:
计算X射线散射分布,包括:计算散射点对探测器上任一点造成的散射光子数:
Nf=ρ×N3×P×ΔΩ,
其中,ρ为一调节常数,表示电子密度等对光子数的影响;N3表示X射线与物体作用过程中光子数;P表示散射光子偏转θ的概率;ΔΩ表示探测器对a点的张角;
从原始投影数据中扣除散射分布。
2.根据权利要求1所述的基于卷积核的锥束CT散射伪影校正方法,其特征在于,估算X射线的初始的光子数N0,包括:
利用对被测物体进行CT扫描时的电压、电流和X射线能谱估算X射线的初始的光子数N0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710401166.8A CN107202805B (zh) | 2017-05-31 | 2017-05-31 | 基于卷积核的锥束ct散射伪影校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710401166.8A CN107202805B (zh) | 2017-05-31 | 2017-05-31 | 基于卷积核的锥束ct散射伪影校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107202805A CN107202805A (zh) | 2017-09-26 |
CN107202805B true CN107202805B (zh) | 2020-05-05 |
Family
ID=59907033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710401166.8A Active CN107202805B (zh) | 2017-05-31 | 2017-05-31 | 基于卷积核的锥束ct散射伪影校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107202805B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3731759A4 (en) | 2017-12-29 | 2021-01-20 | Shanghai United Imaging Healthcare Co., Ltd. | SYSTEMS AND METHODS FOR SCRAP CORRECTION OF AN IMAGE |
CN108606805B (zh) * | 2018-03-15 | 2022-02-11 | 东软医疗系统股份有限公司 | 一种散射校正方法、装置及设备 |
CN109480892B (zh) * | 2018-12-29 | 2022-11-04 | 上海联影医疗科技股份有限公司 | 一种图像的生成方法 |
CN109846501B (zh) * | 2019-04-02 | 2023-02-28 | 深圳市安健科技股份有限公司 | 散射线校正方法及终端 |
CN113804709B (zh) * | 2021-09-22 | 2022-08-12 | 清华大学 | 基于拟蒙特卡罗和强制探测的校正ct散射信号的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102967555A (zh) * | 2012-11-16 | 2013-03-13 | 西安交通大学 | 一种含散射修正的光子成像系统图像重建系统及方法 |
CN103200873A (zh) * | 2011-08-18 | 2013-07-10 | 株式会社东芝 | 光子计数型的x射线计算机断层装置以及散射线校正方法 |
CN105447832A (zh) * | 2015-12-14 | 2016-03-30 | 天津三英精密仪器有限公司 | 一种基于探测器单元标定的ct图像伪影校正方法及应用 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7463712B2 (en) * | 2006-05-18 | 2008-12-09 | The Board Of Trustees Of The Leland Stanford Junior University | Scatter correction for x-ray imaging using modulation of primary x-ray spatial spectrum |
EP2036038A1 (en) * | 2006-06-22 | 2009-03-18 | Philips Intellectual Property & Standards GmbH | Method and system for error compensation |
WO2010133985A2 (en) * | 2009-05-20 | 2010-11-25 | Koninklijke Philips Electronics, N.V. | Continuous time-of-flight scatter simulation method |
US9086366B2 (en) * | 2012-02-15 | 2015-07-21 | L-3 Communications Security And Detection Systems, Inc. | Determining a material property based on scattered radiation |
CN104166962B (zh) * | 2014-07-29 | 2018-06-15 | 南京邮电大学 | 一种使用散射核方法的锥束ct散射校正方法 |
-
2017
- 2017-05-31 CN CN201710401166.8A patent/CN107202805B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103200873A (zh) * | 2011-08-18 | 2013-07-10 | 株式会社东芝 | 光子计数型的x射线计算机断层装置以及散射线校正方法 |
CN102967555A (zh) * | 2012-11-16 | 2013-03-13 | 西安交通大学 | 一种含散射修正的光子成像系统图像重建系统及方法 |
CN105447832A (zh) * | 2015-12-14 | 2016-03-30 | 天津三英精密仪器有限公司 | 一种基于探测器单元标定的ct图像伪影校正方法及应用 |
Non-Patent Citations (3)
Title |
---|
同位素Tl-201的单光子发射CT影像的康普顿散射校正;钱霖;《计算物理》;20020131;第19卷(第1期);第37-42页 * |
工业锥束CT伪影校正方法研究;张峰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20120315(第03期);正文第11-15、25-32页 * |
钱霖.同位素Tl-201的单光子发射CT影像的康普顿散射校正.《计算物理》.2002,第19卷(第1期),第37-42页. * |
Also Published As
Publication number | Publication date |
---|---|
CN107202805A (zh) | 2017-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107202805B (zh) | 基于卷积核的锥束ct散射伪影校正方法 | |
Sun et al. | Improved scatter correction using adaptive scatter kernel superposition | |
US7203276B2 (en) | X-ray scatter image reconstruction by balancing of discrepancies between detector responses, and apparatus therefor | |
US7519143B2 (en) | Method and system for generating a scatter corrected X-ray image | |
JP2007529738A (ja) | 干渉性散乱ctのビーム硬化補正および減衰補正 | |
Dhaene et al. | A realistic projection simulator for laboratory based X-ray micro-CT | |
Bhatia et al. | Scattering correction using continuously thickness-adapted kernels | |
Dixon et al. | An improved analytical model for CT dose simulation with a new look at the theory of CT dose | |
KR101874235B1 (ko) | 방사능의 추정 방법 | |
Sun et al. | Correction for patient table‐induced scattered radiation in cone‐beam computed tomography (CBCT) | |
US9615807B2 (en) | Systems and methods for improving image quality in cone beam computed tomography | |
Shi et al. | A Monte Carlo study of the impact of phosphor optical properties on EPID imaging performance | |
Díaz et al. | A fast scatter field estimator for digital breast tomosynthesis | |
Lindgren et al. | X-ray modeling of realistic synthetic radiographs of thin titanium welds | |
Shi et al. | GPU-accelerated Monte Carlo simulation of MV-CBCT | |
Thanasupsombat et al. | A Simple Scatter Reduction Method in Cone‐Beam Computed Tomography for Dental and Maxillofacial Applications Based on Monte Carlo Simulation | |
CN108065950B (zh) | 一种放射成像方法及其系统 | |
JP6744749B2 (ja) | 放射線撮影装置、放射線撮影システム、放射線撮影方法、及びプログラム | |
US10079078B2 (en) | Method for correcting a spectrum | |
US20120163695A1 (en) | Intra-detector scatter correction | |
Apel et al. | Cosmic ray energy reconstruction from the S (500) observable recorded in the KASCADE-Grande air shower experiment | |
Bhatia et al. | Separable scatter model of the detector and object contributions using continuously thickness-adapted kernels in CBCT | |
US6249568B1 (en) | Process for improving a signal/noise ratio of the image of a moving object | |
Wang et al. | Reconstruction algorithm for point source neutron imaging through finite thickness scintillator | |
CN111685786B (zh) | X射线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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |