A signal detection model for quantifying overregularization in nonlinear image reconstruction

Visualization of the distortion induced by nonlinear noise reduction in computed tomography

Purpose

We developed a method to visualize the image distortion induced by nonlinear noise reduction algorithms in computed tomography (CT) systems.

Approach

Nonlinear distortion was defined as the induced residual when testing a reconstruction algorithm by the criteria for a linear system. Two types of images were developed: a nonlinear distortion of an object (NLD object ) image and a nonlinear distortion of noise (NLD noise ) image to visualize the nonlinear distortion induced by an algorithm. Calculation of the images requires access to the sinogram data, which is seldomly fully provided. Hence, an approximation of the NLD object image was estimated. Using simulated CT acquisitions, four noise levels were added onto forward projected sinograms of a typical CT image; these were noise reduced using a median filter with the simultaneous iterative reconstruction technique or a total variation filter with the conjugate gradient least-squares algorithm. The linear reconstruction technique filtered back-projection was also analyzed for comparison.

Results

Structures in the NLD object image indicated contrast and resolution reduction of the nonlinear denoising. Although the approximated NLD object image represented the original NLD object image well, it had a higher random uncertainty. The NLD noise image for the median filter indicated both stochastic variations and structures reminding of the object while for the total variation filter only stochastic variations were indicated.

Conclusions

The developed images visualize nonlinear distortions of denoising algorithms. The object may be distorted by the noise and vice versa. Analyzing the distortion correlated to the object is more critical than analyzing a distortion of stochastic variations. The absence of nonlinear distortion may measure the robustness of the denoising algorithm.

Addressing CT metal artifacts using photon-counting detectors and one-step spectral CT image reconstruction

PURPOSE

The constrained one-step spectral CT image reconstruction (cOSSCIR) algorithm with a nonconvex alternating direction method of multipliers optimizer is proposed for addressing computed tomography (CT) metal artifacts caused by beam hardening, noise, and photon starvation. The quantitative performance of cOSSCIR is investigated through a series of photon-counting CT simulations.

METHODS

cOSSCIR directly estimates basis material maps from photon-counting data using a physics-based forward model that accounts for beam hardening. The cOSSCIR optimization framework places constraints on the basis maps, which we hypothesize will stabilize the decomposition and reduce streaks caused by noise and photon starvation. Another advantage of cOSSCIR is that the spectral data need not be registered, so that a ray can be used even if some energy window measurements are unavailable. Photon-counting CT acquisitions of a virtual pelvic phantom with low-contrast soft tissue texture and bilateral hip prostheses were simulated. Bone and water basis maps were estimated using the cOSSCIR algorithm and combined to form a virtual monoenergetic image for the evaluation of metal artifacts. The cOSSCIR images were compared to a "two-step" decomposition approach that first estimated basis sinograms using a maximum likelihood algorithm and then reconstructed basis maps using an iterative total variation constrained least-squares optimization (MLE+TVmin $_{\text{min}}$ ). Images were also compared to a nonspectral TVmin $_{\text{min}}$ reconstruction of the total number of counts detected for each ray with and without normalized metal artifact reduction (NMAR) applied. The simulated metal density was increased to investigate the effects of increasing photon starvation. The quantitative error and standard deviation in regions of the phantom were compared across the investigated algorithms. The ability of cOSSCIR to reproduce the soft-tissue texture, while reducing metal artifacts, was quantitatively evaluated.

RESULTS

Noiseless simulations demonstrated the convergence of the cOSSCIR and MLE+TVmin $_{\text{min}}$ algorithms to the correct basis maps in the presence of beam-hardening effects. When noise was simulated, cOSSCIR demonstrated a quantitative error of -1 HU, compared to 2 HU error for the MLE+TVmin $_{\text{min}}$ algorithm and -154 HU error for the nonspectral TVmin $_{\text{min}}$ +NMAR algorithm. For the cOSSCIR algorithm, the standard deviation in the central iodine region of interest was 20 HU, compared to 299 HU for the MLE+TVmin $_{\text{min}}$ algorithm, 41 HU for the MLE+TVmin $_{\text{min}}$ +Mask algorithm that excluded rays through metal, and 55 HU for the nonspectral TVmin $_{\text{min}}$ +NMAR algorithm. Increasing levels of photon starvation did not impact the bias or standard deviation of the cOSSCIR images. cOSSCIR was able to reproduce the soft-tissue texture when an appropriate regularization constraint value was selected.

CONCLUSIONS

By directly inverting photon-counting CT data into basis maps using an accurate physics-based forward model and a constrained optimization algorithm, cOSSCIR avoids metal artifacts due to beam hardening, noise, and photon starvation. The cOSSCIR algorithm demonstrated improved stability and accuracy compared to a two-step method of decomposition followed by reconstruction.