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.

Beaming hardening correction algorithm based on polychromatic system parameters in computed tomography

Beam hardening correction method for Cone-Beam computed tomography based on reprojection of slice contours

Beam hardening artifact correction in CT based on re projection

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

A project of X-ray hardening correction in large ICT

Beam hardening correction based on image noise statistics

The Correction Method Based on Optimal Solution in Digital Imaging system

Pixel-by-pixel correction of beam hardening artifacts by bowtie filter in fan-beam CT

A fast beam hardening correction method incorporated in a filtered back-projection based MAP algorithm.

Eliminating Beam-Hardening Artifacts In High-Energy Industrial Computed Tomography(Ict)

Beam Hardening Correction Based on Monte Carlo Simulation

A Practical Analysis of Scatter-to-primary Ratio after Beam Hardening Correction for X-Ray CT

Beam Hardening Correction for Middle-Energy Industrial Computerized Tomography

Penalized-Likelihood Sinogram Restoration for Computed Tomography

Monte Carlo Simulation for Polychromatic X-ray Fluorescence Computed Tomography with Sheet-Beam Geometry

A patient-specific scatter artifacts correction method

Energy-angle Correlation Correction Algorithm for Monochromatic Computed Tomography Based on Thomson Scattering X-ray Source

Scatter Correction for X-Ray Conebeam CT Using One-Dimensional Primary Modulation

Lighting Deviation Correction for Integrating-Sphere Multispectral Imaging Systems.

Blind beam-hardening correction from Poisson measurements

Correction for photon attenuation in SPECT: analytical framework, average attenuation factors, and a new hybrid approach.