Khobragade P, Fan Jiahua, Rupcich Franco, Crotty Dominic J, Schmidt Taly Gilat
Department of Biomedical Engineering, Marquette University and Medical College of Wisconsin, Milwaukee, WI, 53233, USA.
GE Healthcare, Waukesha, WI, 53188, USA.
Med Phys. 2018 Jun 9. doi: 10.1002/mp.13040.
Evaluation of noise texture information in CT images is important for assessing image quality. Noise texture is often quantified by the noise power spectrum (NPS), which requires numerous image realizations to estimate. This study evaluated fractal dimension for quantifying noise texture as a scalar metric that can potentially be estimated using one image realization.
The American College of Radiology CT accreditation phantom (ACR) was scanned on a clinical scanner (Discovery CT750, GE Healthcare) at 120 kV and 25 and 90 mAs. Images were reconstructed using filtered back projection (FBP/ASIR 0%) with varying reconstruction kernels: Soft, Standard, Detail, Chest, Lung, Bone, and Edge. For each kernel, images were also reconstructed using ASIR 50% and ASIR 100% iterative reconstruction (IR) methods. Fractal dimension was estimated using the differential box-counting algorithm applied to images of the uniform section of ACR phantom. The two-dimensional Noise Power Spectrum (NPS) and one-dimensional-radially averaged NPS were estimated using established techniques. By changing the radiation dose, the effect of noise magnitude on fractal dimension was evaluated. The Spearman correlation between the fractal dimension and the frequency of the NPS peak was calculated. The number of images required to reliably estimate fractal dimension was determined and compared to the number of images required to estimate the NPS-peak frequency. The effect of Region of Interest (ROI) size on fractal dimension estimation was evaluated. Feasibility of estimating fractal dimension in an anthropomorphic phantom and clinical image was also investigated, with the resulting fractal dimension compared to that estimated within the uniform section of the ACR phantom.
Fractal dimension was strongly correlated with the frequency of the peak of the radially averaged NPS curve, having a Spearman rank-order coefficient of 0.98 (P-value < 0.01) for ASIR 0%. The mean fractal dimension at ASIR 0% was 2.49 (Soft), 2.51 (Standard), 2.52 (Detail), 2.57 (Chest), 2.61 (Lung), 2.66 (Bone), and 2.7 (Edge). A reduction in fractal dimension was observed with increasing ASIR levels for all investigated reconstruction kernels. Fractal dimension was found to be independent of noise magnitude. Fractal dimension was successfully estimated from four ROIs of size 64 × 64 pixels or one ROI of 128 × 128 pixels. Fractal dimension was found to be sensitive to non-noise structures in the image, such as ring artifacts and anatomical structure. Fractal dimension estimated within a uniform region of an anthropomorphic phantom and clinical head image matched that estimated within the ACR phantom for filtered back projection reconstruction.
Fractal dimension correlated with the NPS-peak frequency and was independent of noise magnitude, suggesting that the scalar metric of fractal dimension can be used to quantify the change in noise texture across reconstruction approaches. Results demonstrated that fractal dimension can be estimated from four, 64 × 64-pixel ROIs or one 128 × 128 ROI within a head CT image, which may make it amenable for quantifying noise texture within clinical images.
评估CT图像中的噪声纹理信息对于评估图像质量很重要。噪声纹理通常通过噪声功率谱(NPS)进行量化,这需要大量图像实例来估计。本研究评估分形维数作为一种标量指标来量化噪声纹理,它有可能通过一个图像实例来估计。
使用美国放射学会CT认证体模(ACR)在临床扫描仪(Discovery CT750,GE医疗)上进行扫描,管电压为120 kV,管电流分别为25和90 mAs。使用滤波反投影(FBP/ASIR 0%)并结合不同的重建核进行图像重建:软组织、标准、细节、胸部、肺部、骨和边缘核。对于每个核,还使用ASIR 50%和ASIR 100%迭代重建(IR)方法进行图像重建。使用差分盒计数算法对ACR体模均匀部分的图像估计分形维数。使用既定技术估计二维噪声功率谱(NPS)和一维径向平均NPS。通过改变辐射剂量,评估噪声幅度对分形维数的影响。计算分形维数与NPS峰值频率之间的斯皮尔曼相关性。确定可靠估计分形维数所需的图像数量,并与估计NPS峰值频率所需的图像数量进行比较。评估感兴趣区域(ROI)大小对分形维数估计的影响。还研究了在仿真人体体模和临床图像中估计分形维数的可行性,并将得到的分形维数与在ACR体模均匀部分内估计的分形维数进行比较。
分形维数与径向平均NPS曲线峰值频率密切相关,对于ASIR 0%,斯皮尔曼等级系数为0.98(P值<0.01)。ASIR 0%时的平均分形维数分别为:软组织核2.49、标准核2.51、细节核2.52、胸部核2.57、肺部核2.61、骨核2.66和边缘核2.7。对于所有研究的重建核,随着ASIR水平的增加,分形维数降低。发现分形维数与噪声幅度无关。从四个64×64像素的ROI或一个128×128像素的ROI成功估计了分形维数。发现分形维数对图像中的非噪声结构敏感,如环形伪影和解剖结构。在仿真人体体模和临床头部图像的均匀区域内估计的分形维数与滤波反投影重建时在ACR体模内估计的分形维数相匹配。
分形维数与NPS峰值频率相关且与噪声幅度无关,这表明分形维数这一标量指标可用于量化不同重建方法下噪声纹理的变化。结果表明,在头部CT图像中,可从四个64×64像素的ROI或一个128×128的ROI估计分形维数,这可能使其适用于量化临床图像中的噪声纹理。