Department of Medical Physics, University of Wisconsin-Madison School of Medicine and Public Health, 1111 Highland Avenue, Madison, WI, 53705, USA.
Department of Radiology, University of Wisconsin-Madison School of Medicine and Public Health, 600 Highland Avenue, Madison, WI, 53792, USA.
Med Phys. 2018 May;45(5):1942-1956. doi: 10.1002/mp.12856. Epub 2018 Apr 6.
Low-signal correction (LSC) in the raw counts domain has been shown to effectively reduce noise streaks in CT because the data inconsistency associated with photon-starved regions may be mitigated prior to the log transformation step. However, a systematic study of the performance of these raw data correction methods is still missing in literature. The purpose of this work was to provide such a systematic study for two well-known low-signal correction schemes using either the adaptive trimmed mean (ATM) filter or the anisotropic diffusion (AD) filter in the raw counts domain.
Image data were acquired experimentally using an anthropomorphic chest phantom and a benchtop cone-beam CT (CBCT) imaging system. Phantom scans were repeated 50 times at a reduced dose level of 0.5 mGy and a reference level of 1.9 mGy. The measured raw counts at 0.5 mGy underwent LSC using the ATM and AD filters. Two relevant parameters were identified for each filter and approximately one hundred operating points in each parameter space were analyzed. Following LSC and log transformation, FDK reconstruction was performed for each case. Noise and spatial resolution properties were assessed across the parameter spaces that define each LSC filter; the results were summarized through 2D contour maps to better understand the trade-offs between these competing image quality features. 2D noise power spectrum (NPS) and modulation transfer function (MTF) were measured locally at two spatial locations in the field-of-view (FOV): a posterior region contaminated by noise streaks and an anterior region away from noise streaks. An isotropy score metric was introduced to characterize the directional dependence of the NPS and MTF (viz., ϵ and ϵ , respectively), with a range from 0 for highly anisotropic to 1 for perfectly isotropic. The noise magnitude and coarseness were also measured.
(a) Both the ATM and AD LSC methods were successful in reducing noise streaks, but their noise and spatial resolution properties were found to be highly anisotropic and shift-variant. (b) NPS isotropy scores in the posterior region were generally improved from ϵ = 0.09 for the images without LSC to the range ϵ = (0.11, 0.67) for ATM and ϵ = (0.06, 0.67) for AD, depending on the filter parameters used. (c) The noise magnitude was reduced across the parameter space of either LSC filter whenever a change along the axis of the controlling parameter led to stronger raw data filtration. Changes in noise magnitude were inversely related to changes in spatial resolution along the direction perpendicular to the streaks. No correlation was found, however, between the contour maps of noise magnitude and the NPS isotropy. (d) Both filters influenced the noise coarseness anisotropically, with coarser noise occurring along directions perpendicular to the noise streaks. The anisotropic noise coarseness was intrinsically and directly related to resolution losses in a given direction: coarseness plots mimic the topography of the 2D MTF, i.e., the coarser the noise, the lower the resolution.
Both AD and ATM LSC schemes enable low-dose CBCT imaging. However, it was found that noise magnitude and overall spatial resolution vary considerably across the parameter space for each filter, and more importantly these image quality features are highly anisotropic and shift-variant.
在原始计数域中进行低信号校正(LSC)已被证明可以有效地减少 CT 中的噪声条纹,因为与光子饥饿区域相关的数据不一致性可以在对数变换步骤之前得到缓解。然而,文献中仍然缺少对这些原始数据校正方法性能的系统研究。本研究的目的是提供一种系统的研究,使用自适应修剪均值(ATM)滤波器或各向异性扩散(AD)滤波器在原始计数域中对两种著名的低信号校正方案进行研究。
使用人体胸部模拟体模和台式锥形束 CT(CBCT)成像系统进行实验性图像数据采集。在降低剂量水平 0.5 mGy 和参考剂量水平 1.9 mGy 下,重复扫描体模 50 次。0.5 mGy 下的测量原始计数经过 ATM 和 AD 滤波器的 LSC。为每个滤波器确定了两个相关参数,并在每个参数空间中分析了大约一百个工作点。在 LSC 和对数变换后,对每个病例进行 FDK 重建。在定义每个 LSC 滤波器的参数空间中评估噪声和空间分辨率特性;通过二维轮廓图总结结果,以便更好地理解这些竞争的图像质量特性之间的权衡。在视场(FOV)中的两个空间位置局部测量二维噪声功率谱(NPS)和调制传递函数(MTF):受噪声条纹污染的后区和远离噪声条纹的前区。引入各向异性分数度量来描述 NPS 和 MTF 的方向依赖性(即,ϵ和ϵ,分别),范围从 0 到完全各向同性的 1。还测量了噪声幅度和粗糙度。
(a)ATM 和 AD LSC 方法都成功地减少了噪声条纹,但发现它们的噪声和空间分辨率特性具有高度各向异性和位移变化。(b)在 LSC 之后,后区的 NPS 各向异性分数通常从没有 LSC 的图像的 ϵ=0.09 提高到 ATM 的 ϵ=(0.11,0.67)和 AD 的 ϵ=(0.06,0.67),这取决于使用的滤波器参数。(c)在任意 LSC 滤波器的参数空间中,只要沿着控制参数的轴发生变化导致更强的原始数据滤波,噪声幅度就会降低。噪声幅度的变化与条纹垂直方向的空间分辨率的变化成反比。然而,在噪声幅度轮廓图和 NPS 各向异性之间没有发现相关性。(d)两种滤波器都以各向异性的方式影响噪声粗糙度,噪声粗糙度沿垂直于噪声条纹的方向更粗糙。各向异性噪声粗糙度与给定方向的分辨率损失直接相关:粗糙度图模仿二维 MTF 的地形,即噪声越粗糙,分辨率越低。
AD 和 ATM LSC 方案都可用于低剂量 CBCT 成像。然而,研究发现,对于每个滤波器,噪声幅度和整体空间分辨率在参数空间中变化很大,更重要的是,这些图像质量特征具有高度的各向异性和位移变化。