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.

Resolution Properties Of Non-Parametric Regression Sinogram Smoothing Using An Explicit Poisson Model

Resolution Properties of Non-Parametric Regression Sinogram Smoothing Using an Explicit Poisson Model

Nonparametric Regression Sinogram Smoothing Using a Roughness-Penalized Poisson Likelihood Objective Function

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

Potential Equivalence of Sinogram and Image-Domain Penalized Likelihood Methods

A Noise Reduction Method for Non-Stationary Noise Model of SPECT Sinogram Based on Kalman Filter

Few-view Tomography Using Roughness-Penalized Nonparametric Regression and Periodic Spline Interpolation

Robust Bayesian Super-Resolution Approach Via Sparsity Enforcing a Priori for Near-Field Aeroacoustic Source Imaging

3D Distance-Weighted Wiener Filter for Poisson Noise Reduction in Sinogram Space for SPECT Imaging

Analytical Noise Treatment for Low-Dose Ct Projection Data by Penalized Weighted Least-Square Smoothing in the K-L Domain

Penalized-Likelihood Sinogram Restoration for Computed Tomography

Periodic Splines and Gaussian Processes for the Resolution of Linear Inverse Problems

Combined Transformation of Ordering SPECT Sinograms for Signal Extraction from Measurements with Poisson Noise.

PET image reconstruction: a robust state space approach.

TU‐G‐110‐09: Spectral CT (SCT) Reconstruction Using Tight Wavelet Frames

Favorable Noise Uniformity Properties of Fourier-Based Interpolation and Reconstruction Approaches in Single-Slice Helical Computed Tomography

Graph-based joint denoising and super-resolution of generalized piecewise smooth images

Spline-based Inverse Radon Transform in Two and Three Dimensions

Nonlinear tomographic reconstruction via nonsmooth optimization

Development of a high-performance noise-reduction filter for tomographic reconstruction

A Detailed Investigation of Non-Iterative Methods for 3d Spect Image Reconstruction