Abstract:We formulate computed tomography (CT) sinogram preprocessing as a statistical restoration problem in which the goal is to obtain the best estimate of the line integrals needed for reconstruction from the set of noisy, degraded measurements. CT measurement data are degraded by a number of factors-including beam hardening and off-focal radiation-that produce artifacts in reconstructed images unless properly corrected. Currently, such effects are addressed by a sequence of sinogram-preprocessing steps, including deconvolution corrections for off-focal radiation, that have the potential to amplify noise. Noise itself is generally mitigated through apodization of the reconstruction kernel, which effectively ignores the measurement statistics, Although in high-noise situations adaptive filtering methods that loosely model data statistics are sometimes applied. As an alternative, we present a general imaging model relating the degraded measurements to the sinogram of ideal line integrals and propose to estimate these line integrals by iteratively optimizing a statistically based objective function. We consider three different strategies for estimating the set of ideal line integrals, one based on direct estimation of ideal "monochromatic" line integrals that have been corrected for single-material beam hardening, one based on estimation of ideal "polychromatic" line integrals that can be readily mapped to monochromatic line integrals, and one based on estimation of ideal transmitted intensities, from which ideal, monochromatic line integrals can be readily estimated. The first two approaches involve maximization of a penalized Poisson-likelihood objective function while the third involves minimization of a quadratic penalized weighted least squares (PWLS) objective applied in the transmitted intensity domain. We find that at low exposure levels typical of those being considered for screening CT, the Poisson-likelihood based approaches outperform the PWLS objective as well as a standard approach based on adaptive filtering followed by deconvolution. At higher exposure levels, the approaches all perform similarly.

Piecewise Tikhonov Regularization Scheme To Reconstruct Discontinuous Density In Computerized Tomography

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Convergent Data-Driven Regularizations for CT Reconstruction

A Numerical Differentiation Method and Its Application to Reconstruction of Discontinuity

X-Ray CT Image Reconstruction Via Wavelet Frame Based Regularization and Radon Domain Inpainting

Numerical Differentiation by a Tikhonov Regularization Method Based on the Discrete Cosine Transform

Tikhonov Regularization Based on Near-Optimal Regularization Parameter with Application to Capacitance Tomography Image Reconstruction

Window Function-Based Regularization for Electrical Capacitance Tomography Image Reconstruction

Penalized-Likelihood Sinogram Restoration for Computed Tomography

A cascading l0 regularization reconstruction method in nonsubsampled contourlet domain for limited-angle CT

CT image reconstruction algorithms based on the Hanke Raus parameter choice rule

Low-dose CT reconstruction via edge-preserving total variation regularization

A hybrid image reconstruction algorithm for electrical capacitance tomography

GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

A relaxed iterated Tikhonov regularization for linear ill-posed inverse problems

Better than the Total Variation Regularization

Tikhonov regularization with conjugate gradient least squares method for large-scale discrete ill-posed problem in image restoration

Well-posedness and Tikhonov regularization of an inverse source problem for a parabolic equation with an integral constraint

A Wavelet-Based Multiresolution Regularized Least Squares Reconstruction Approach for Optical Tomography

A General-Thresholding Solution for ℓp (0 < p <1) Regularized CT Reconstruction