Liu Ruirui, Zhang Shuangyue, Zhao Tianyu, O'Sullivan Joseph A, Williamson Jeffrey F, Webb Tyler, Porras-Chaverri Mariela, Whiting Bruce
Department of Radiation Oncology, Washington University School of Medicine, St. Louis, MO, 63110, USA.
Department of Electrical and Systems Engineering, Washington University, St. Louis, MO, USA.
Med Phys. 2021 Feb;48(2):852-870. doi: 10.1002/mp.14652. Epub 2020 Dec 23.
To investigate via Monte Carlo simulations, the impact of scan subject size, antiscatter grid (ASG), collimator size, and bowtie filter on the distribution of scatter radiation in a typical realistically modeled third generation 16 slice diagnostic computed tomography (CT) scanner.
Full radiation transport was simulated with Geant4 in a realistic CT scanner geometric model, including the imaging phantom, bowtie filter (BTF), collimators and detector assembly, except for the ASGs. An analytical method was employed to quantify the probable transmission through the ASG of each photon intersecting the detector array. Normalized scatter profiles (NSP) and scatter-to-primary-ratio (SPR) profiles were simulated for 90 and 140 kVp beams for different size phantoms and slice thicknesses. The impact of CT scatter on the reconstructed attenuation coefficient factor was also studied as were the modulating effects of phantom- and patient-tissue heterogeneities on scatter profiles. A method to characterize the relative spatial frequency content of sinogram signals was developed to assess the latter.
For the 21.4-cm diameter phantom, NSP and SPR increase linearly with collimator opening for both tube potentials, with the 90 kVp scan exhibiting slightly larger NSP and SPR. The BTF modestly modulates scatter under the phantom center, reducing the prominent off-axis lobes by factors of 1.1-1.3. The ASG reduces scatter on the central axis NSP threefold, and reduces scatter at the detectors outside the phantom shadow by factors of 25 to 500. For the phantoms with diameters of 27 and 32 cm, the scatter increases roughly three- and fourfold, respectively, demonstrating that scatter monotonically increases with phantom size, despite deployment of the ASG and BTF. In the absence of a scan subject, the ASG reduces the signal profile arising photons scattered by the BTF. Without ASG, the in-air scatter profile is relatively flat compared to the scatter profile when the ASG is present. For both 90 and 140 kVp photon spectra, the calculated attenuation coefficient decreases linearly with increasing collimation size. For both homogeneous and heterogeneous objects, NSPs are dominated by low spatial frequency content compared to the primary signal. However, the SPR, which quantifies the local magnitude of nonlinear detector response and is dominated by the high frequency content of the primary profile, can contribute strongly to high-spatial frequency streaking artifacts near high-density structures in reconstructed image artifacts.
Public-domain Monte Carlo codes, Geant-4 in particular, is a feasible method for characterizing CT detector response to scattered- and off-focal radiation. Our study demonstrates that the ASG substantially reduces the scatter radiation and reshapes scatter-radiation profiles and affects the accuracy with which the detector array can measure narrow-beam attenuation due its inability to distinguish between true uncollided primary and narrow-angle coherently scattered photons. Hence, incorporating the impact of detector array collimation into the forward-projection signal formation models used by iterative reconstruction algorithms is necessary to use CT for accurately characterizing material properties. While tissue heterogeneities exercise a modest influence on local NPS shape and magnitude, they do not add significant high spatial frequency content.
通过蒙特卡罗模拟,研究扫描对象尺寸、防散射格栅(ASG)、准直器尺寸和蝴蝶结滤波器对典型的真实建模第三代16层诊断计算机断层扫描(CT)扫描仪中散射辐射分布的影响。
在一个真实的CT扫描仪几何模型中,使用Geant4模拟全辐射传输,该模型包括成像体模、蝴蝶结滤波器(BTF)、准直器和探测器组件,但不包括ASG。采用一种分析方法来量化每个与探测器阵列相交的光子通过ASG的可能传输。针对不同尺寸的体模和切片厚度,模拟了90和140 kVp束的归一化散射剖面(NSP)和散射与原发射线比(SPR)剖面。还研究了CT散射对重建衰减系数因子的影响,以及体模和患者组织不均匀性对散射剖面的调制作用。开发了一种表征正弦图信号相对空间频率内容的方法来评估后者。
对于直径为21.4 cm的体模,两种管电压下NSP和SPR均随准直器开口线性增加,90 kVp扫描的NSP和SPR略大。BTF适度调制体模中心下方的散射,将突出的离轴叶减少1.1 - 1.3倍。ASG将中心轴NSP上的散射减少三倍,并将体模阴影外探测器处的散射减少25至500倍。对于直径为27和32 cm的体模,散射分别大致增加三倍和四倍,表明尽管部署了ASG和BTF,散射仍随体模尺寸单调增加。在没有扫描对象的情况下,ASG会降低由BTF散射产生的光子的信号剖面。没有ASG时,与存在ASG时的散射剖面相比,空气中的散射剖面相对平坦。对于90和140 kVp光子光谱,计算得到的衰减系数随准直尺寸增加而线性降低。对于均匀和不均匀物体,与原发射线相比,NSP主要由低空间频率内容主导。然而,SPR量化了非线性探测器响应的局部幅度,并且由原发射线剖面的高频内容主导,它会对重建图像伪影中高密度结构附近的高空间频率条纹伪影产生强烈影响。
公共领域的蒙特卡罗代码,特别是Geant - 4,是表征CT探测器对散射和离焦辐射响应的可行方法。我们的研究表明,ASG显著降低了散射辐射,重塑了散射辐射剖面,并影响探测器阵列测量窄束衰减的准确性,因为它无法区分真正未碰撞的原发射线光子和窄角相干散射光子。因此,将探测器阵列准直的影响纳入迭代重建算法使用的前向投影信号形成模型中,对于使用CT准确表征材料特性是必要的。虽然组织不均匀性对局部NPS形状和幅度有适度影响,但它们不会增加显著的高空间频率内容。