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

Approaches, advantages, and challenges to photon counting detector and multi-energy CT

Photon counting detector CT (PCD-CT) is the newest major development in CT technology and has been commercially available since 2021. It offers major technological advantages over current standard-of-care energy integrating detector CT (EID-CT) including improved spatial resolution, improved iodine contrast to noise ratio, multi-energy imaging, and reduced noise. This article serves as a foundational basis to the technical approaches and concepts of PCD-CT technology with primary emphasis on detector technology in direct comparison to EID-CT. The article also addresses current technological challenges to PCD-CT with particular attention to cross talk and its causes (e.g., Compton scattering, fluorescence, charge sharing, K-escape) as well as pile-up.

Influence of helical pitch and gantry rotation time on image quality and file size in ultrahigh-resolution photon-counting detector CT

The goal of this experimental study was to quantify the influence of helical pitch and gantry rotation time on image quality and file size in ultrahigh-resolution photon-counting CT (UHR-PCCT). Cervical and lumbar spine, pelvis, and upper legs of two fresh-frozen cadaveric specimens were subjected to nine dose-matched UHR-PCCT scan protocols employing a collimation of 120 × 0.2 mm with varying pitch (0.3/1.0/1.2) and rotation time (0.25/0.5/1.0 s). Image quality was analyzed independently by five radiologists and further substantiated by placing normed regions of interest to record mean signal attenuation and noise. Effective mAs, CT dose index (CTDIvol), size-specific dose estimate (SSDE), scan duration, and raw data file size were compared. Regardless of anatomical region, no significant difference was ascertained for CTDIvol (p ≥ 0.204) and SSDE (p ≥ 0.240) among protocols. While exam duration differed substantially (all p ≤ 0.016), the lowest scan time was recorded for high-pitch protocols (4.3 ± 1.0 s) and the highest for low-pitch protocols (43.6 ± 15.4 s). The combination of high helical pitch and short gantry rotation times produced the lowest perceived image quality (intraclass correlation coefficient 0.866; 95% confidence interval 0.807-0.910; p < 0.001) and highest noise. Raw data size increased with acquisition time (15.4 ± 5.0 to 235.0 ± 83.5 GByte; p ≤ 0.013). Rotation time and pitch factor have considerable influence on image quality in UHR-PCCT and must therefore be chosen deliberately for different musculoskeletal imaging tasks. In examinations with long acquisition times, raw data size increases considerably, consequently limiting clinical applicability for larger scan volumes.

Three-Dimensional Visualization of Shunt Valves with Photon Counting CT and Comparison to Traditional X-ray in a Simple Phantom Model

This study introduces an application of innovative medical technology, Photon Counting Computer Tomography (PC CT) with novel detectors, for the assessment of shunt valves. PC CT technology offers enhanced visualization capabilities, especially for small structures, and opens up new possibilities for detailed three-dimensional imaging. Shunt valves are implanted under the skin and redirect excess cerebrospinal fluid, for example, to the abdominal cavity through a catheter. They play a vital role in regulating cerebrospinal fluid drainage in various pathologies, which can lead to hydrocephalus. Accurate imaging of shunt valves is essential to assess the rate of drainage, as their precise adjustment is a requirement for optimal patient care. This study focused on two adjustable shunt valves, the proGAV 2.0® and M. blue® (manufactured by Miethke, Potsdam, Germany). A comprehensive comparative analysis of PC CT and traditional X-ray techniques was conducted to explore this cutting-edge technology and it demonstrated that routine PC CT can efficiently assess shunt valves' adjustments. This technology shows promise in enhancing the accurate management of shunt valves used in settings where head scans are already frequently required, such as in the treatment of hydrocephalus.

Dual energy CT reconstruction using the constrained one step spectral image reconstruction algorithm

BACKGROUND

The constrained one-step spectral CT Image Reconstruction method (cOSSCIR) has been developed to estimate basis material maps directly from spectral CT data using a model of the polyenergetic x-ray transmissions and incorporating convex constraints into the inversion problem. This 'one-step' approach has been shown to stabilize the inversion in the case of photon-counting CT, and may provide similar benefits to dual-kV systems that utilize integrating detectors. Since the approach does not require the same rays be acquired for every spectral measurement, cOSSCIR can apply to dual energy protocols and systems used clinically, such as fast and slow kV switching systems and dual source scanning.

PURPOSE

The purpose of this study is to investigate the use of cOSSCIR applied to dual-kV data, using both registered and unregistered spectral acquisitions, specifically slow and fast kV switching imaging protocols. For this application, cOSSCIR is investigated using inverse crime simulations and dual-kV experiments. This study is the first demonstration of cOSSCIR on the dual-kV reconstruction problem.

METHODS

An integrating detector model was developed for the purpose of reconstructing dual-kV data, and an inverse crime study was used to validate the detector model within the cOSSCIR framework using a simulated pelvic phantom. Experiments were also used to evaluate cOSSCIR on the dual energy problem. Dual-kV data was obtained from a physical phantom containing analogs of adipose, bone, and liver tissues, with the aim of recovering the material coefficients in the bone and adipose basis material maps. cOSSCIR was applied to acquisitions where all rays performed both spectral measurements (registered) and fast and slow kV switching acquisitions (unregistered). cOSSCIR was also compared to two image-domain decomposition approaches, where image-domain methods are the conventional approach for decomposing unregistered spectral data.

RESULTS

Simulations demonstrate the application of cOSSCIR to the dual-kV inversion problem by successfully recovering the material basis maps on ideal data, while further showing that unregistered data presents a more challenging inversion problem. In our experimental reconstructions, the recovered basis material coefficient errors were found to be less than 6.5% in the bone, adipose, and liver regions for both registered and unregistered protocols. Similarly, the errors were less than 4% in the 50 keV virtual mono-energetic images, and the recovered material decomposition vectors nearly overlap their corresponding ground-truth vectors. Additionally, a preliminary two material decomposition study of iodine quantification recovered an average concentration of 9.2 mg/mL from a 10 mg/mL experimental iodine analog.

CONCLUSIONS

Using our integrating detector and spectral models, cOSCCIR is capable of accurately recovering material basis maps from dual-kV data for both registered and unregistered data. The material decomposition quantification compare favorably to the image domain approaches, and our results were not affected by the imaging protocol. Our results also suggest the extension of cOSSCIR to iodine quantification using two material decomposition.

Simultaneous activity and attenuation estimation in TOF-PET with TV-constrained nonconvex optimization

An alternating direction method of multipliers (ADMM) framework is developed for nonsmooth biconvex optimization for inverse problems in imaging. In particular, the simultaneous estimation of activity and attenuation (SAA) problem in time-of-flight positron emission tomography (TOF-PET) has such a structure when maximum likelihood estimation (MLE) is employed. The ADMM framework is applied to MLE for SAA in TOF-PET, resulting in the ADMM-SAA algorithm. This algorithm is extended by imposing total variation (TV) constraints on both the activity and attenuation map, resulting in the ADMM-TVSAA algorithm. The performance of this algorithm is illustrated using the penalized maximum likelihood activity and attenuation estimation (P-MLAA) algorithm as a reference. Additional results on step-size tuning and on the use of unconstrained ADMM-SAA are presented in the previous arXiv submission: arXiv:2303.17042v1.

Advances in metal artifact reduction in CT images: A review of traditional and novel metal artifact reduction techniques

Metal artifacts degrade CT image quality, hampering clinical assessment. Numerous metal artifact reduction methods are available to improve the image quality of CT images with metal implants. In this review, an overview of traditional methods is provided including the modification of acquisition and reconstruction parameters, projection-based metal artifact reduction techniques (MAR), dual energy CT (DECT) and the combination of these techniques. Furthermore, the additional value and challenges of novel metal artifact reduction techniques that have been introduced over the past years are discussed such as photon counting CT (PCCT) and deep learning based metal artifact reduction techniques.

Convergence for nonconvex ADMM, with applications to CT imaging

The alternating direction method of multipliers (ADMM) algorithm is a powerful and flexible tool for complex optimization problems of the form m i n { f ( x ) + g ( y ) : A x + B y = c } . ADMM exhibits robust empirical performance across a range of challenging settings including nonsmoothness and nonconvexity of the objective functions f and g , and provides a simple and natural approach to the inverse problem of image reconstruction for computed tomography (CT) imaging. From the theoretical point of view, existing results for convergence in the nonconvex setting generally assume smoothness in at least one of the component functions in the objective. In this work, our new theoretical results provide convergence guarantees under a restricted strong convexity assumption without requiring smoothness or differentiability, while still allowing differentiable terms to be treated approximately if needed. We validate these theoretical results empirically, with a simulated example where both f and g are nondifferentiable-and thus outside the scope of existing theory-as well as a simulated CT image reconstruction problem.

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.

Direct Multi-Material Reconstruction via Iterative Proximal Adaptive Descent for Spectral CT Imaging

Spectral computed tomography (spectral CT) is a promising medical imaging technology because of its ability to provide information on material characterization and quantification. However, with an increasing number of basis materials, the nonlinearity of measurements causes difficulty in decomposition. In addition, noise amplification and beam hardening further reduce image quality. Thus, improving the accuracy of material decomposition while suppressing noise is pivotal for spectral CT imaging. This paper proposes a one-step multi-material reconstruction model as well as an iterative proximal adaptive decent method. In this approach, a proximal step and a descent step with adaptive step size are designed under the forward-backward splitting framework. The convergence analysis of the algorithm is further discussed according to the convexity of the optimization objective function. For simulation experiments with different noise levels, the peak signal-to-noise ratio (PSNR) obtained by the proposed method increases approximately 23 dB, 14 dB, and 4 dB compared to those of other algorithms. Magnified areas of thorax data further demonstrated that the proposed method has a better ability to preserve details in tissues, bones, and lungs. Numerical experiments verify that the proposed method efficiently reconstructed the material maps, and reduced noise and beam hardening artifacts compared with the state-of-the-art methods.

Low-Dose CT Imaging of the Pelvis in Follow-up Examinations-Significant Dose Reduction and Impact of Tin Filtration: Evaluation by Phantom Studies and First Systematic Retrospective Patient Analyses

OBJECTIVES

Low-dose (LD) computed tomography (CT) is still rarely used in musculoskeletal (MSK) radiology. This study evaluates the potentials of LD CT for follow-up pelvic imaging with special focus on tin filtration (Sn) technology for normal and obese patients with and without metal implants.

MATERIALS AND METHODS

In a phantom study, 5 different LD and normal-dose (ND) CT protocols with and without tin filtration were tested using a normal and an obese phantom. Iterative reconstruction (IR) and filtered back projection (FBP) were used for CT image reconstruction. In a subsequent retrospective patient study, ND CT images of 45 patients were compared with follow-up tin-filtered LD CT images with a 90% dose reduction. Sixty-four percent of patients contained metal implants at the follow-up examination. Computed tomography images were objectively (image noise, contrast-to-noise ratio [CNR], dose-normalized contrast-to-noise ratio [CNRD]) and subjectively, using a 6-point Likert score, evaluated. In addition, the figure of merit was calculated. For group comparisons, paired t tests, Wilcoxon signed rank test, analysis of variance, or Kruskal-Wallis tests were used, where applicable.

RESULTS

The LD Sn protocol with 67% dose reduction resulted in equal values in qualitative (Likert score) and quantitative image analysis (image noise) compared with the ND protocol in the phantom study. For follow-up examinations, dose could be reduced up to 90% by using Sn LD CT scans without impairment in the clinical study. However, metal implants resulted in a mild impairment of Sn LD as well as ND CT images. Cancellous bone ( P < 0.001) was assessed worse and cortical bone ( P = 0.063) equally in Sn LD CT images compared with ND CT images. Figure of merit values were significant ( P ≤ 0.02) lower and hence better in Sn LD as in ND protocols. Obese patients benefited in particular from tin filtration in LD MSK imaging in terms of image noise and CNR ( P ≤ 0.05).

CONCLUSIONS

Low-dose CT scans with tin filtration allow maximum dose reduction while maintaining high image quality for certain clinical purposes, for example, follow-up examinations, especially metal implant position, material loosening, and consolidation controls. Overweight patients benefit particularly from tin filter technology. Although metal implants decrease image quality in ND as well as in Sn LD CT images, this is not a relevant limitation for assessability.

Spectral calibration of photon-counting detectors at high photon flux

BACKGROUND

Calibration of photon-counting detectors (PCDs) is necessary for quantitatively accurate spectral computed tomography (CT), but the calibration process can be complicated by nonlinear flux-dependent physical factors such as pulse pile-up.

PURPOSE

This work develops a method for spectral sensitivity calibration of a PCD-based spectral CT system that incorporates nonlinear flux dependence and can thus be employed at high photon flux.

METHODS

A calibration model for the spectral response and polynomial flux dependence is proposed, which incorporates prior x-ray source spectrum and PCD models and that has a small set of parameters for adjusting to the spectral CT system of interest. The model parameters are determined by fitting transmission data from a known object of known composition: a step-wedge phantom composed of different thicknesses of aluminum, a bone equivalent, and polymethyl methacrylate (PMMA), a soft-tissue equivalent. This fitting employs Tikhonov regularization, and the regularization strength and the polynomial order for the intensity modeling are determined by bias and variance analysis. The spectral calibration and nonlinear intensity correction is validated on transmission measurements through a third material, Teflon, at different x-ray photon flux levels.

RESULTS

The nonlinear intensity dependence is determined to be accurately accounted for with a third-order polynomial. The calibrated spectral CT model accurately predicts Teflon transmission to within 1% for flux levels up to 50% of the detector maximum.

CONCLUSIONS

The proposed PCD calibration method enables accurate physical modeling necessary for quantitative imaging in spectral CT. Furthermore, the model applies to high flux settings so that acquisition times will not be limited by restricting the spectral CT system to low flux levels.