CN104166962B - 一种使用散射核方法的锥束ct散射校正方法 - Google Patents
一种使用散射核方法的锥束ct散射校正方法 Download PDFInfo
- Publication number
- CN104166962B CN104166962B CN201410365880.2A CN201410365880A CN104166962B CN 104166962 B CN104166962 B CN 104166962B CN 201410365880 A CN201410365880 A CN 201410365880A CN 104166962 B CN104166962 B CN 104166962B
- Authority
- CN
- China
- Prior art keywords
- image
- scattering
- correction
- scattered
- shutter
- 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
Abstract
本发明公开了一种使用散射核方法的锥束CT散射校正方法，其步骤为：步骤一：测量遮挡板下的散射信号；步骤二：自适应散射核参数的计算；步骤三：估计扩散函数的参数；步骤四：通过扩散函数对图像进行散射校正；步骤五：图像遮挡区域的修复；步骤六：锥束CT的重建。本发明不需要获取模型参数的先验数据，不增加图像获取及重建的复杂度，添加的设备复杂度不高，且不需要重复扫描被测物体，能高效的进行锥束CT散射校正。
Description
技术领域
本发明涉及一种医学图像处理技术领域，尤其涉及锥束CT散射校正领域。
背景技术
锥束CT具有更高的扫描速度和辐射利用率，能有效的减少X射线管的负载输出，降低扫描成本，能快速获得高分辨率三维断层图像数据。影响锥束CT重建图像质量的因素有很多，如X线散射、噪声、几何误差、能谱、探测单元响应不一致等，但由于锥束平板CT使用大范围的X射线平板探测器，这使得成像质量与传统CT相比较更易受到X射线散射的影响，因散射而形成的伪影、CT数的不准确等严重影响对重建图像的分析与判断，因此，针对X射线散射校正的研究显得尤为重要。
国内外有许多方法提出用于解决锥束CT的散射校正，按照类型我们可分为基于散射测量的方法和基于非散射测量的方法。非散射测量的方法包括空气间隙滤air-gap，反散射网格anti-scatter grids，分析方法和蒙特卡洛方法。这些方法通常会各种限制，如：设备结构、被测物密度分布和计算量巨大，只能降低一部分因散射造成的伪影。基于散射测量方法是在X射线管与被扫描物体之间加一个射线遮挡设备，来测量部分的散射射线，通过插值能方法来估计整个区域的散射分布，基于散射测量方法中，有些方法需要多次扫描来获取散射信号和主射线信号,这样增加了扫描时间和病人的X射线照射剂量；有些方法只通过一次扫描就可以了，但这些方法中有些增加了图像获取和重建的复杂度；有些是增加了设备的复杂度和扫描模式的难度。
在散射校正的方法中，蒙特卡洛模拟方法非常耗费时间；模型分析计算方法需要知道X射线谱的信息、被照物体的物质组成和几何结构，实际这些参数不太容易获得；初级射线调制方法中校正结果受限于调制板自身的结构；基于部分散射射线测量方法，有些需要增加照射剂量，有些方法对散射分布的估计准确度不高；Star-Lack等使用散射核的研究方法进行锥束CT散射校正，在不增加照射剂量的情况下，可以方便的解决不同能级的X射线散射对锥束CT重建图像的影响，但其需要使用先验知识估计模型参数，且其没有较好的构建散射模型，一定程度上阻碍了方法使用的广度和精度。
发明内容
针对现有技术中存在的问题，本发明提供了一种使用散射核方法的锥束CT散射校正方法，是依靠加入遮挡板来获取部分区域散射分布，自适应的调整散射分布模型参数，进而求解散射核的参数。
本发明的技术解决方案是：一种使用散射核方法的锥束CT散射校正方法，其步骤为：
步骤一：测量遮挡板下的散射信号；
步骤二：自适应散射核参数的计算；
步骤三：估计扩散函数的参数；
步骤四：通过扩散函数对图像进行散射校正；
步骤五：图像遮挡区域的修复；
步骤六：锥束CT的重建。
其中，遮挡板形状为镂空“十”字型，由两个铅条十字交叉构成，铅条宽2mm、高3mm、长120mm，铅条放置在圆形塑料板上；所述遮挡板放置在X射线源与被照物体之间，且围绕圆心按一定方向旋转；所述遮挡板离射线源230mm，离被照物770mm，离平板探测器1170mm；遮挡板放置在X射线源与被照物体之间是防止重建时因相同区域连续遮挡造成重建误差过大；遮挡板的形状可以多样化，主要目的是测出部分区域的散射分布情况。
其中估计扩散函数的参数和通过扩散函数对图像进行散射校正是通过迭代方法进行的。
在通过扩散函数对图像进行散射校正中：散射分布函数PSF公式为：
估计(X,Y)点处的散射公式为：
这里xm,yn是大小为(M*N)图像坐标系中的坐标，为了增加计算效率，其图像可以是降采样后的图像；I'P(xm,yn)是初始射线强度估计。
在获取遮挡板下方的散射分布后，进行计算α,β,A,B,σ1,σ2，计算公式为：
Af(xm,yn)＝A.(Ip(xm,yn)/I0(xm,yn))α.(ln(I0(xm,yn)/Ip(xm,yn)))β，
对应的不同区域可以求出不同的α,β,A,B,σ1,σ2，在初始散射校正后图像通过与散射核做卷积计算得到进一步的散射校正。
同时，图像遮挡区域的修复是使用基于样例纹理的遮挡区域修复算法来进一步抑制因为遮挡部分区域对断层重建造成的影响，其步骤为：首先，标记要修复的区域；其次，使用被修复区域临近的像素的纹理计算纹理和结构；最后，填充。
本发明于现有技术相比具有积极的有益效果：本发明的使用散射核方法的锥束CT散射校正方法采用了自适应散射核的方式进行自适应的获取构建的模型参数，依靠加入遮挡板来获取部分区域散射分布，自适应的调整散射分布模型参数，进而求解散射核的参数，不需要获取模型参数的先验数据，不增加图像获取及重建的复杂度，添加的设备复杂度不高，也不需要重复扫描被测物体，能高效的进行锥束CT散射校正。
附图说明
图1为本发明中遮挡板结构状态示意图；
图2为本发明流程图；
图3为蒙特卡洛模拟投影图；
图4为使用本发明方法计算的散射分布图；
图5为去除散射的投影图；
图6为修复后的投影图一；
图7为模体Ⅱ的几何结构示意图；
图8为蒙特卡洛模拟结果；
图9为修复后的投影图二；
图10为应用Star-Lack方法校正结果示意图；
图11为模体CatPhan 500及锥束CT实验室系统；
图12为本发明散射校正前后对比图一；
图13为本发明散射校正前后对比图二。
具体实施方式
以下结合附图，通过具体实施例对本发明技术方案做进一步的说明。
本实施例的使用散射核方法的锥束CT散射校正方法是使用蒙特卡洛粒子进行模拟。图1为本发明中遮挡板结构状态示意图；图2为本发明流程图；
如图3、4、5、6，图3为蒙特卡洛模拟投影图；图4为使用本发明方法计算的散射分布图；图5为去除散射的投影图；图6为修复后的投影图一。模体Ⅰ是由PMMA构成的长方体40cm×40cm×8cm，使用5×107个光子模拟，结果如图3所示，图4是使用本发明方法计算的散射分布图；图5为去除散射的投影图；图6是对遮挡板下方的图像使用最邻近插值的方法进行插补，进行散射校正后，并使用基于样本纹理的图像修复技术后的投影图。
对散射校正效果分析见表1：
表1 杯状伪影(τcup)与投影噪声的分析
本发明是杯状伪影τcup从26.09％下降到1.26％，但投影噪声有一定的上升，从11上升到16，未进行散射校正的投影图表示为PI_NONE,散射校正后的投影图表示为PI_SC。
如图7、8、9、10，图7为模体的几何结构示意图；图8为蒙特卡洛模拟结果；图9为修复后的投影图二；图10为应用Star-Lack方法校正结果示意图。摸体Ⅱ的几何结构如图7所示，材料由PPMA树脂组成，使用5×108个光子模拟，120KV，蒙特卡洛模拟结果如图8所示，从图中能较明显看出，本发明能降低散射带来的伪影。
表2 为模体Ⅱ的定量分析表：
表2 模体Ⅱ投影图的散射校正定理分析
如图11、12、13所示，图11为模体CatPhan 500及锥束CT实验室系统；图12为本发明散射校正前重建断层图；图13为本发明散射校正后重建断层图。从图13中能看出本发明方法能使杯状伪影和条纹状伪影明显减少，其中杯状伪影从36％减少到10％。图13中方框用来强调与图12中相应位置条纹状伪影的减少。
虽然本发明已以较佳实施例公开如上，但它们并不是用来限定本发明，任何熟习此技艺者，在不脱离本发明之精神和范围内，自当可作各种变化或润饰，但同样在本发明的保护范围之内。
Claims (1)
1.一种使用散射核方法的锥束CT散射校正方法，其特征在于，其步骤为：
步骤一：测量遮挡板下的散射信号；所述步骤一中遮挡板形状为遮挡中心是镂空的“十”字型，由两个铅条十字交叉构成，铅条宽2mm、高3mm、长120mm，铅条放置在圆形塑料板上；所述遮挡板放置在X射线源与被照物体之间，且围绕圆心按一定方向旋转；所述遮挡板离射线源230mm，离被照物770mm，离平板探测器1170mm；遮挡板放置在X射线源与被照物体之间是防止重建时因相同区域连续遮挡造成重建误差过大；遮挡板是测出部分区域的散射分布情况；
步骤二：自适应散射核参数的计算；
步骤三：估计扩散函数的参数；
步骤四：通过扩散函数对图像进行散射校正；所述估计扩散函数的参数和所述通过扩散函数对图像进行散射校正是通过迭代方法进行；
所述通过扩散函数对图像进行散射校正，散射分布函数PSF公式为：
PSF(Ip，Io，r)＝A·(Ip/Io)α·(1n(Io/Ip))β·[exp(-r2/2σ1 2)+Bexp(-r2/2σ2 2))]；
估计(x，y)点处的散射公式为：
获取遮挡板下方的散射分布后，进行计算α，β，A，B，σ1，σ2，计算公式为：
Af(xm，yn)＝A·(Ip(xm，yn)/Io(xm，yn))α·(ln(Io(xm，yn)/Ip(xm，yn)))β，
对应的不同区域可以求出不同的α，β，A，B，σ1，σ2，在初始散射校正后图像通过与散射核做卷积计算得到进一步的散射校正；
式中，xm，yn是大小为M*N图像坐标系中的坐标，其图像是降采样后的图像；I′P(xm，yn)是初始射线强度估计；
步骤五：图像遮挡区域的修复；所述图像遮挡区域的修复是使用基于样例纹理的遮挡区域修复算法来进一步抑制因为遮挡部分区域对断层重建造成的影响，其步骤为：首先，标记要修复的区域；其次，使用被修复区域临近的像素的纹理计算纹理和结构；最后，填充；
步骤六：锥束CT的重建。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410365880.2A CN104166962B (zh) | 2014-07-29 | 2014-07-29 | 一种使用散射核方法的锥束ct散射校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410365880.2A CN104166962B (zh) | 2014-07-29 | 2014-07-29 | 一种使用散射核方法的锥束ct散射校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104166962A CN104166962A (zh) | 2014-11-26 |
| CN104166962B true CN104166962B (zh) | 2018-06-15 |
Family
ID=51910755
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410365880.2A Active CN104166962B (zh) | 2014-07-29 | 2014-07-29 | 一种使用散射核方法的锥束ct散射校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104166962B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504656B (zh) * | 2014-12-10 | 2017-06-27 | 浙江大学 | 锥束ct图像域快速散射修正方法 |
| DE102016217984A1 (de) * | 2016-09-20 | 2018-04-05 | Siemens Healthcare Gmbh | Sinogrammbasierte Streustrahlenkorrektur in der Computertomographie |
| CN106361367B (zh) * | 2016-12-01 | 2019-10-08 | 上海联影医疗科技有限公司 | 一种检测器的校正方法和使用该校正方法的装置及设备 |
| GB201703291D0 (en) * | 2017-03-01 | 2017-04-12 | Ibex Innovations Ltd | Apparatus and method for the correction of scatter in a radiographic system |
| CN107202805B (zh) * | 2017-05-31 | 2020-05-05 | 中国人民解放军信息工程大学 | 基于卷积核的锥束ct散射伪影校正方法 |
| CN107802280B (zh) * | 2017-11-08 | 2021-02-19 | 上海联影医疗科技股份有限公司 | 校正曲线生成方法、投影图像的校正方法、系统及存储介质 |
| CN108804834A (zh) * | 2018-06-14 | 2018-11-13 | 中科超精（安徽）科技有限公司 | 复杂非均匀结构散射核获取方法 |
| CN109846501B (zh) * | 2019-04-02 | 2023-02-28 | 深圳市安健科技股份有限公司 | 散射线校正方法及终端 |
| CN113409414B (zh) * | 2021-06-08 | 2024-03-26 | 江苏一影医疗设备有限公司 | X线图像的散射修正方法、装置、电子设备、存储介质 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102068270A (zh) * | 2010-12-16 | 2011-05-25 | 西安交通大学 | 一种使用静止环状射线阻挡阵列的锥束ct扫描和散射校正方法 |
| CN103578082A (zh) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | 一种锥束ct散射校正方法及系统 |
-
2014
- 2014-07-29 CN CN201410365880.2A patent/CN104166962B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102068270A (zh) * | 2010-12-16 | 2011-05-25 | 西安交通大学 | 一种使用静止环状射线阻挡阵列的锥束ct扫描和散射校正方法 |
| CN103578082A (zh) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | 一种锥束ct散射校正方法及系统 |
Non-Patent Citations (3)
| Title |
|---|
| Region filling and object removal by exemplar-based image inpainting;A. Criminisi, P. Perez, K. Toyama;《IEEE Transactions on Image Processing》;20040930;第13卷(第9期);摘要、第III部分、图4-5 * |
| Scatter correction for cone-beam computed tomography using self-adaptive scatter kernel superposition;谢世朋, 罗立名;《中国物理C》;20120630;第36卷(第6期);摘要、第567页第2段,第2.1-2.3节,第3.2节、图1,3,5 * |
| 基于衰减板的锥束CT散射校正;谢世朋，罗立民;《电子学报》;20110715;第39卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104166962A (zh) | 2014-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104166962B (zh) | 一种使用散射核方法的锥束ct散射校正方法 | |
| Maldera et al. | Digital breast tomosynthesis: dose and image quality assessment | |
| US10175181B2 (en) | X-ray imaging system allowing the correction of the scatter radiation and precise detection of the distance between the source and the detector | |
| Telsemeyer et al. | Quantitative carbon ion beam radiography and tomography with a flat-panel detector | |
| Zhang et al. | Comparison of patient specific dose metrics between chest radiography, tomosynthesis, and CT for adult patients of wide ranging body habitus | |
| Loy Rodas et al. | Seeing is believing: increasing intraoperative awareness to scattered radiation in interventional procedures by combining augmented reality, Monte Carlo simulations and wireless dosimeters | |
| Dauvergne et al. | Monte Carlo comparison of x-ray and proton CT for range calculations of proton therapy beams | |
| Yang et al. | A breast-specific, negligible-dose scatter correction technique for dedicated cone-beam breast CT: a physics-based approach to improve Hounsfield Unit accuracy | |
| CN108292428A (zh) | 图像重建的系统和方法 | |
| CN103987320A (zh) | X线ct装置以及散射x线校正方法 | |
| Rinaldi et al. | Experimental investigations on carbon ion scanning radiography using a range telescope | |
| Wang et al. | Measurement of the spatial resolution and the relative density resolution in an industrial cone-beam micro computed tomography system | |
| Saunders Jr et al. | The effect of breast compression on mass conspicuity in digital mammography | |
| Mettivier et al. | Cone‐beam breast computed tomography with a displaced flat panel detector array | |
| Altunbas et al. | A unified scatter rejection and correction method for cone beam computed tomography | |
| Quiñones et al. | Filtered back-projection reconstruction for attenuation proton CT along most likely paths | |
| US9615807B2 (en) | Systems and methods for improving image quality in cone beam computed tomography | |
| Díaz et al. | A fast scatter field estimator for digital breast tomosynthesis | |
| Sølie et al. | Image quality of list-mode proton imaging without front trackers | |
| Darne et al. | A proton imaging system using a volumetric liquid scintillator: a preliminary study | |
| Tendler et al. | Image quality evaluation of projection-and depth dose-based approaches to integrating proton radiography using a monolithic scintillator detector | |
| Hayati et al. | Monte Carlo modeling of a conventional X-ray computed tomography scanner for gel dosimetry purposes | |
| CN105528771B (zh) | 一种使用能量函数方法的锥束ct中杯状伪影的校正方法 | |
| Szafraniec et al. | Synchrotron based planar imaging and digital tomosynthesis of breast and biopsy phantoms using a CMOS active pixel sensor | |
| Shah et al. | Characterization of X-ray scattering for various phantoms and clinical breast geometries using breast CT on a dedicated hybrid system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20141126 Assignee: Nanjing causal Artificial Intelligence Research Institute Co., Ltd Assignor: Nanjing Post & Telecommunication Univ. Contract record no.: X2019320000168 Denomination of invention: Cone beam CT scattering correction method by use of scattering nucleus method Granted publication date: 20180615 License type: Common License Record date: 20191028 |