Image quality improvement of a one-step spectral CT reconstruction on a prototype photon-counting scanner

Objective. X-ray spectral computed tomography (CT) allows for material decomposition (MD). This study compared a one-step material decomposition MD algorithm with a two-step reconstruction MD algorithm using acquisitions of a prototype CT scanner with a photon-counting detector (PCD).Approach. MD and CT reconstruction may be done in two successive steps, i.e. decompose the data in material sinograms which are then reconstructed in material CT images, or jointly in a one-step algorithm. The one-step algorithm reconstructed material CT images by maximizing their Poisson log-likelihood in the projection domain with a spatial regularization in the image domain. The two-step algorithm maximized first the Poisson log-likelihood without regularization to decompose the data in material sinograms. These sinograms were then reconstructed into material CT images by least squares minimization, with the same spatial regularization as the one step algorithm. A phantom simulating the CT angiography clinical task was scanned and the data used to measure noise and spatial resolution properties. Low dose carotid CT angiographies of 4 patients were also reconstructed with both algorithms and analyzed by a radiologist. The image quality and diagnostic clinical task were evaluated with a clinical score.Main results. The phantom data processing demonstrated that the one-step algorithm had a better spatial resolution at the same noise level or a decreased noise value at matching spatial resolution. Regularization parameters leading to a fair comparison were selected for the patient data reconstruction. On the patient images, the one-step images received higher scores compared to the two-step algorithm for image quality and diagnostic.Significance. Both phantom and patient data demonstrated how a one-step algorithm improves spectral CT image quality over the implemented two-step algorithm but requires a longer computation time. At a low radiation dose, the one-step algorithm presented good to excellent clinical scores for all the spectral CT images.

Constrained one-step material decomposition reconstruction of head CT data from a silicon photon-counting prototype

BACKGROUND

Spectral CT material decomposition provides quantitative information but is challenged by the instability of the inversion into basis materials. We have previously proposed the constrained One-Step Spectral CT Image Reconstruction (cOSSCIR) algorithm to stabilize the material decomposition inversion by directly estimating basis material images from spectral CT data. cOSSCIR was previously investigated on phantom data.

PURPOSE

This study investigates the performance of cOSSCIR using head CT datasets acquired on a clinical photon-counting CT (PCCT) prototype. This is the first investigation of cOSSCIR for large-scale, anatomically complex, clinical PCCT data. The cOSSCIR decomposition is preceded by a spectrum estimation and nonlinear counts correction calibration step to address nonideal detector effects.

METHODS

Head CT data were acquired on an early prototype clinical PCCT system using an edge-on silicon detector with eight energy bins. Calibration data of a step wedge phantom were also acquired and used to train a spectral model to account for the source spectrum and detector spectral response, and also to train a nonlinear counts correction model to account for pulse pileup effects. The cOSSCIR algorithm optimized the bone and adipose basis images directly from the photon counts data, while placing a grouped total variation (TV) constraint on the basis images. For comparison, basis images were also reconstructed by a two-step projection-domain approach of Maximum Likelihood Estimation (MLE) for decomposing basis sinograms, followed by filtered backprojection (MLE + FBP) or a TV minimization algorithm (MLE + TVmin ) to reconstruct basis images. We hypothesize that the cOSSCIR approach will provide a more stable inversion into basis images compared to two-step approaches. To investigate this hypothesis, the noise standard deviation in bone and soft-tissue regions of interest (ROIs) in the reconstructed images were compared between cOSSCIR and the two-step methods for a range of regularization constraint settings.

RESULTS

cOSSCIR reduced the noise standard deviation in the basis images by a factor of two to six compared to that of MLE + TVmin , when both algorithms were constrained to produce images with the same TV. The cOSSCIR images demonstrated qualitatively improved spatial resolution and depiction of fine anatomical detail. The MLE + TVmin algorithm resulted in lower noise standard deviation than cOSSCIR for the virtual monoenergetic images (VMIs) at higher energy levels and constraint settings, while the cOSSCIR VMIs resulted in lower noise standard deviation at lower energy levels and overall higher qualitative spatial resolution. There were no statistically significant differences in the mean values within the bone region of images reconstructed by the studied algorithms. There were statistically significant differences in the mean values within the soft-tissue region of the reconstructed images, with cOSSCIR producing mean values closer to the expected values.

CONCLUSIONS

The cOSSCIR algorithm, combined with our previously proposed spectral model estimation and nonlinear counts correction method, successfully estimated bone and adipose basis images from high resolution, large-scale patient data from a clinical PCCT prototype. The cOSSCIR basis images were able to depict fine anatomical details with a factor of two to six reduction in noise standard deviation compared to that of the MLE + TVmin two-step approach.

Technical Note: The nearest neighborhood-based approach for estimating basis line-integrals using photon-counting detector

PURPOSE

This study aims to develop a calibration-based estimator for the photon-counting detector (PCD)-based x-ray computed tomography.

METHODS

We propose the nearest neighborhood (NN)-based estimator, which searches for the nearest calibration data for a given PCD output and sets the associated basis line-integrals as the estimate. Searching for the nearest neighbors can be accelerated using the pre-calculated k-d tree for the data.

RESULTS

The proposed method is compared to the model-based maximum likelihood (ML) estimator. For slab phantom study, both ML and NN-based methods achieve the Cramér-Rao lower bound and are unbiased for various combinations of three basis materials (water, bone, and gold). The proposed method is also validated for K-edge imaging and presents almost unbiased Au concentrations in the region of interest.

CONCLUSIONS

The proposed NN-based method is demonstrated to be as accurate as the model-based ML estimator, but it is computationally efficient and requires only calibration measurements.