Stability and Robustness Analysis of Regularized Reconstruction Methods for Low-Dose Computed Tomography in Parallel-Beam Geometry
Authors
Abstract
Low-dose computed tomography (LDCT) reduces radiation exposure but increases the ill-posedness of the reconstruction problem due to noise and sparse data. While regularized methods like Tikhonov and Total Variation (TV) improve image quality, their performance depends heavily on noise characteristics, sampling conditions, and parameter selection. This study presents a systematic stability and robustness analysis of Filtered Back Projection (FBP), Tikhonov regularization, and TV minimization within a 2D parallel-beam CT framework. A unified simulation pipeline based on the Radon transform is developed and evaluated using both the modified Shepp-Logan phantom and a clinical thorax image. Reconstruction behavior is investigated under multiple degradation scenarios involving Gaussian, Poisson, and mixed noise models, across baseline (180 projections) and sparse-view (60 projections) acquisition geometries. To ensure a fair comparison, regularization parameters are optimized for each scenario through an exhaustive SSIM-based grid-search. Quality is assessed via RMSE, PSNR, and SSIM, while robustness is quantified through an empirical Stability Factor S measuring perturbation amplification from measurement to image space. The results show that FBP is highly sensitive to noise and undersampling. Tikhonov regularization improves structural fidelity compared with FBP but remains more sensitive to perturbation than TV. Conversely, TV provides the best compromise between noise suppression, edge preservation, accuracy, and numerical stability. These findings highlight the stability-resolution trade-off in LDCT and demonstrate that the proposed Stability Factor S offers valuable complementary information to conventional metrics.