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.

Blind beam-hardening correction from Poisson measurements

Polychromatic sparse image reconstruction and mass attenuation spectrum estimation via B-spline basis function expansion

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

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

Penalized-Likelihood Sinogram Restoration for Computed Tomography

Beam Hardening Correction for Sparse-View Ct Reconstruction

Beam hardening correction based on image noise statistics

Multi-Material Blind Beam Hardening Correction Based on Non-Linearity Adjustment of Projections

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

Spectrum Binning Approach for Multi-Material Beam Hardening Correction (MMBHC) in CT

Spectral CT Reconstruction with Image Sparsity and Spectral Mean

Statistical iterative spectral CT imaging method based on blind separation of polychromatic projections

Simple beam hardening correction method (2DCalBH) based on 2D linearization

Statistical Image Reconstruction Using Mixed Poisson-Gaussian Noise Model for X-Ray CT

Optimization-based Reconstruction for Correcting Non-Linear Partial Volume Artifacts in CT

Statistical image reconstruction with beam‐hardening compensation for X‐ray CT by a calibration step (2DIterBH)

Basis-image reconstruction directly from sparse-view data in spectral CT

A Beam Hardening Correction Method Based on HL Consistency

Robust Beam Hardening Artifacts Reduction for Computed Tomography Using Spectrum Modeling.

SBL-LCGL: sparse Bayesian learning based on Laplace distribution for robust cone-beam x-ray luminescence computed tomography

Alternating Minimization Algorithm with Automatic Relevance Determination for Transmission Tomography under Poisson Noise