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.

Variable Sinograms And Redundant Information In Single-Photon Emission Computed Tomography With Non-Uniform Attenuation

Mathematical Formulation of the Potato Peeler Perspective

[Analytical Algorithms Improving the Reconstruction of SPECT Image].

Variable sinograms and redundant information in tomographic imaging

Image Reconstruction for Discontinuous Objects with Half-Data Knowledge in SPECT with Nonuniform Attenuation

Short-Scan Spect Imaging With Non-Uniform Attenuation And 3d Distance-Dependent Spatial Resolution

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

Preliminary investigation of a novel reconstruction algorithm based upon the potato peeler perspective

/Spl Pi/-Scheme Short-Scan SPECT and Image Reconstruction with Nonuniform Attenuation

Nonlinear Dual Reconstruction of SPECT Activity and Attenuation Images

Penalized-Likelihood Sinogram Restoration for Computed Tomography

Π-Scheme Short-Scan SPECT and Image Reconstruction

Simultaneous reconstruction of tissue attenuation and radioactivity maps in SPECT.

Non-iterative Methods and Their Noise Characteristics in 2D SPECT Image Reconstruction

Potential Equivalence of Sinogram and Image-Domain Penalized Likelihood Methods

Image Reconstruction in 3D Short-Scan SPECT

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

Evaluation of imaging protocol for ECT based on CS image reconstruction algorithm

Interior SPECT- Exact and Stable ROI Reconstruction from Uniformly Attenuated Local Projections.

Image Reconstruction in SPECT with a Half Detector

Efficient Calculation of Resolution and Variance in 2-D Circular-Orbit SPECT