Abstract: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.

An algorithm for constrained one-step inversion of spectral CT data

Non-convex primal-dual algorithm for image reconstruction in spectral CT

Direct inversion of spectral CT data into a materials decomposition and the effect of multiple soft tissues

One-step inversion algorithms for spectral CT, with application to dynamic Cone Beam CT

An inversion algorithm for P-functions with applications to Multi-energy CT

An extended primal-dual algorithm framework for nonconvex problems: Application to image reconstruction in spectral CT

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

Three material decomposition for spectral computed tomography enabled by block-diagonal step-preconditioning

Optimization-based Algorithm for Solving the Discrete X-Ray Transform with Nonlinear Partial Volume Effect

Derivative-Free iterative One-Step Reconstruction for Multispectral CT

A Fast Regularized Iterative Algorithm for Fan-Beam CT Reconstruction.

High resolution image reconstruction with constrained, total-variation minimization

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

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

Preliminary Study of Image Reconstruction from Sparse-View Data in Phase-Contrast CT

Alternating Minimization Based Framework For Simultaneous Spectral Calibration And Image Reconstruction In Spectral Ct

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

Tv-Constrained Incremental Algorithms For Low-Intensity Ct Image Reconstruction

One-step inverse generation network for sparse-view dual-energy CT reconstruction and material imaging

An Investigation of Direct Image Reconstruction in DECT with Physical Data

Sparse2Inverse: Self-supervised inversion of sparse-view CT data