Quantitative material characterization from multi-energy photon counting CT

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.

First experimental evaluation of a high-resolution deep silicon photon-counting sensor

Purpose

Current photon-counting computed tomography detectors are limited to a pixel size of around 0.3 to 0.5 mm due to excessive charge sharing degrading the dose efficiency and energy resolution as the pixels become smaller. In this work, we present measurements of a prototype photon-counting detector that leverages the charge sharing to reach a theoretical sub-pixel resolution in the order of 1    μ m . The goal of the study is to validate our Monte-Carlo simulation using measurements, enabling further development.

Approach

We measure the channel response at the MAX IV Lab, in the DanMAX beamline, with a 35 keV photon beam, and compare the measurements with a 2D Monte Carlo simulation combined with a charge transport model. Only a few channels on the prototype are connected to keep the number of wire bonds low.

Results

The measurements agree generally well with the simulations with the beam close to the electrodes but diverge as the beam is moved further away. The induced charge cloud signals also seem to increase linearly as the beam is moved away from the electrodes.

Conclusions

The agreement between measurements and simulations indicates that the Monte-Carlo simulation can accurately model the channel response of the detector with the photon interactions close to the electrodes, which indicates that the unconnected electrodes introduce unwanted effects that need to be further explored. With the same Monte-Carlo simulation previously indicating a resolution of around 1    μ m with similar geometry, the results are promising that an ultra-high resolution detector is not far in the future.

Leveraging 3D Bioprinting and Photon-Counting Computed Tomography to Enable Noninvasive Quantitative Tracking of Multifunctional Tissue Engineered Constructs

3D bioprinting is revolutionizing the fields of personalized and precision medicine by enabling the manufacturing of bioartificial implants that recapitulate the structural and functional characteristics of native tissues. However, the lack of quantitative and noninvasive techniques to longitudinally track the function of implants has hampered clinical applications of bioprinted scaffolds. In this study, multimaterial 3D bioprinting, engineered nanoparticles (NPs), and spectral photon-counting computed tomography (PCCT) technologies are integrated for the aim of developing a new precision medicine approach to custom-engineer scaffolds with traceability. Multiple CT-visible hydrogel-based bioinks, containing distinct molecular (iodine and gadolinium) and NP (iodine-loaded liposome, gold, methacrylated gold (AuMA), and Gd2 O3 ) contrast agents, are used to bioprint scaffolds with varying geometries at adequate fidelity levels. In vitro release studies, together with printing fidelity, mechanical, and biocompatibility tests identified AuMA and Gd2 O3 NPs as optimal reagents to track bioprinted constructs. Spectral PCCT imaging of scaffolds in vitro and subcutaneous implants in mice enabled noninvasive material discrimination and contrast agent quantification. Together, these results establish a novel theranostic platform with high precision, tunability, throughput, and reproducibility and open new prospects for a broad range of applications in the field of precision and personalized regenerative medicine.

Research on accuracy of material identification based on photon counting spectral CT

BACKGROUND

Photon counting spectral CT is a significant direction in the development of CT technology and material identification is an important application of spectral CT. However, spectrum estimation in photon counting spectral CT is highly complex and may affect quantification accuracy of material identification.

OBJECTIVE

To address the problem of energy spectrum estimation in photon-counting spectral CT, this study investigates empirical material decomposition algorithms to achieve accurate quantitative decomposition of the effective atomic number.

METHODS

The spectrum is first calibrated using the empirical dual-energy calibration (EDEC) method and the effective atomic number is then quantitatively estimated based on the EDEC method. The accuracy of estimating the effective atomic number of materials under different calibration conditions is investigated by designing different calibration phantoms, and accurate quantitation is achieved using suitable calibration settings. Last, the validity of this method is verified through simulations and experimental studies.

RESULTS

The results demonstrate that the error in estimating the effective atomic number is reduced to within 4% for low and medium Z materials, thereby enabling accurate material identification.

CONCLUSION

The empirical dual-energy correction method can solve the problem of energy spectrum estimation in photon counting spectral CT. Accurate effective atomic number estimation can be achieved with suitable calibration.

Photon-counting detector CT improves quality of arterial phase abdominal scans: A head-to-head comparison with energy-integrating CT

PURPOSE

Photon-counting detector (PCD)-CT is expected to have a substantial impact on oncologic abdominal imaging. We compared subjective and objective image quality between PCD-CT and conventional energy-integrating detector (EID-)CT arterial phase abdominal scans.

METHODS

This study included 84 patients undergoing both types of abdominal CT. EID-CT scans were acquired with a tube voltage of 100 kVp. With PCD-CT, acquired with 120-kVp, we reconstructed polychromatic T3D images and virtual monoenergetic images (VMIs) in 10-keV intervals from 40 to 90 keV. Quantitative image analysis included noise and contrast-to-noise ratio (CNR) of hepatic vessels, kidney cortex, and hypervascular liver lesions to liver parenchyma. Three raters used a 5-point Likert scale for qualitative image analysis of image noise and contrast, lesion conspicuity, and overall image quality. Radiation dose exposure (CT dose index) was compared between the two CT types.

RESULTS

Mean CT dose index and effective dose were respectively 18 % and 26 % lower with PCD-CT versus EID-CT. Compared with EID-CT, CNRs of kidney cortex and vessel to liver parenchyma were significantly higher in PCD-CT VMIs at energies ≤ 60 keV and in polychromatic T3D images (p < 0.004). Overall image quality of PCD-CT VMIs at 50 and 60 keV was rated as significantly better (p < 0.01) than the EID-CT images (inter-reader agreement alpha = 0.80). Lesion conspicuity was significantly better in low-keV VMIs (p < 0.03) and worse in > 70-keV VMIs.

CONCLUSIONS

With low-keV VMI, PCD-CT yields significantly improved objective and subjective quality of arterial phase oncological imaging compared with EID-CT. This advantage may translate into higher diagnostic confidence and lower radiation dose protocols.

Synergistic multi-spectral CT reconstruction with directional total variation

This work considers synergistic multi-spectral CT reconstruction where information from all available energy channels is combined to improve the reconstruction of each individual channel. We propose to fuse these available data (represented by a single sinogram) to obtain a polyenergetic image which keeps structural information shared by the energy channels with increased signal-to-noise ratio. This new image is used as prior information during a channel-by-channel minimization process through the directional total variation. We analyse the use of directional total variation within variational regularization and iterative regularization. Our numerical results on simulated and experimental data show improvements in terms of image quality and in computational speed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 2'.

Regularization by denoising sub-sampled Newton method for spectral CT multi-material decomposition

Spectral Computed Tomography (CT) is an emerging technology that enables us to estimate the concentration of basis materials within a scanned object by exploiting different photon energy spectra. In this work, we aim at efficiently solving a model-based maximum-a-posterior problem to reconstruct multi-materials images with application to spectral CT. In particular, we propose to solve a regularized optimization problem based on a plug-in image-denoising function using a randomized second order method. By approximating the Newton step using a sketching of the Hessian of the likelihood function, it is possible to reduce the complexity while retaining the complex prior structure given by the data-driven regularizer. We exploit a non-uniform block sub-sampling of the Hessian with inexact but efficient conjugate gradient updates that require only Jacobian-vector products for denoising term. Finally, we show numerical and experimental results for spectral CT materials decomposition. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Photon-counting normalized metal artifact reduction (NMAR) in diagnostic CT

PURPOSE

Metal artifacts can drastically reduce the diagnostic value of computed tomography (CT) images. Even the state-of-the-art algorithms cannot remove them completely. Photon-counting CT inherently provides spectral information, similar to dual-energy CT. Many applications, such as material decomposition, are not possible when metal artifacts are present. Our aim is to develop a prior-based metal artifact reduction specifically for photon-counting CT that can correct each bin image individually or in their combinations.

METHODS

Photon-counting CT sorts incoming photons into several energy bins, producing bin and threshold images containing spectral information. We use this spectral information to obtain a better prior image for the state-of-the-art metal artifact reduction algorithm FSNMAR. First, we apply a non-linear transformation to the bin images to obtain bone-emphasized images. Subsequently, we forward-project the bin images and bone-emphasized images and multiply the resulting sinograms with each other element-wise to mimic beam hardening effects. These sinograms are reconstructed and linearly combined to produce an artifact-reduced image. The coefficients of this linear combination are automatically determined by minimizing a threshold-based cost function in the image domain. After thresholding, we obtain the prior image for FSNMAR, which is applied to the individual bin images and the lowest threshold image. We test our photon-counting normalized metal artifact reduction (PCNMAR) on forensic CT data and compare it to conventional FSNMAR, where the prior is generated via linear sinogram inpainting. For numerical analysis, we compute both the standard deviation in an ROI with metal artifacts and the CNR of soft tissue and fat.

RESULTS

PCNMAR can effectively reduce metal artifacts without sacrificing the overall image quality. Compared to FSNMAR, our method produces fewer secondary artifacts and is more consistent with the measurements. Areas that contain metal, air, and soft tissue are more accurate in PCNMAR. In some cases, the standard deviation in the artifact ROI is reduced by more than 50% relative to FSNMAR, while the CNR values are similar. If extreme artifacts are present, PCNMAR is unable to outperform FSNMAR. Using either two, four, or only the highest energy bin to produce the prior image yielded comparable results.

CONCLUSIONS

PCNMAR is an effective method of reducing metal artifacts in photon-counting CT. The spectral information available in photon-counting CT is highly beneficial for metal artifact reduction, especially the high-energy bin, which inherently contains fewer artifacts. While scanning with four instead of two bins does not provide a better artifact reduction, it allows for more freedom in the selection of energy thresholds.

Photon-counting x-ray detectors for CT

The introduction of photon-counting detectors is expected to be the next major breakthrough in clinical x-ray computed tomography (CT). During the last decade, there has been considerable research activity in the field of photon-counting CT, in terms of both hardware development and theoretical understanding of the factors affecting image quality. In this article, we review the recent progress in this field with the intent of highlighting the relationship between detector design considerations and the resulting image quality. We discuss detector design choices such as converter material, pixel size, and readout electronics design, and then elucidate their impact on detector performance in terms of dose efficiency, spatial resolution, and energy resolution. Furthermore, we give an overview of data processing, reconstruction methods and metrics of imaging performance; outline clinical applications; and discuss potential future developments.

A Clinically Driven Task-Based Comparison of Photon Counting and Conventional Energy Integrating CT for Soft Tissue, Vascular, and High-Resolution Tasks

Photon-counting CT detectors are the next step in advancing CT system development and will replace the current energy integrating detectors (EID) in CT systems in the near future. In this context, the performance of PCCT was compared to EID CT for three clinically relevant tasks: abdominal soft tissue imaging, where differentiating low contrast features is important; vascular imaging, where iodine detectability is critical; and, high-resolution skeletal and lung imaging. A multi-tiered phantom was imaged on an investigational clinical PCCT system (Siemens Healthineers) across different doses using three imaging modes: macro and ultra-high resolution (UHR) PCCT modes and EID CT. Images were reconstructed using filtered backprojection and soft tissue (B30f), vascular (B46f), or high-resolution (B70f; U70f for UHR) kernels. Noise power spectra, task transfer functions, and detectability index were evaluated. For a soft tissue task, PCCT modes showed comparable noise and resolution with improved contrast-to-noise ratio. For a vascular task, PCCT modes showed lower noise and improved iodine detectability. For a high resolution task, macro mode showed lower noise and comparable resolution while UHR mode showed higher noise but improved spatial resolution for both air and bone. PCCT offers competitive advantages to EID CT for clinical tasks.

Dual-energy CT and coronary imaging

Dual-energy computed tomography has been proposed for enhancing the evaluation of coronary artery disease in many fronts. However, the clinical translation of such applications has followed a slower pace of clinical translation. This paper will review the evidence supporting the use of dual-energy computed tomography in coronary artery disease (CAD) and provide some practical illustrations, while underscoring the challenges and gaps in knowledge that have contributed to this phenomenon.

Feasibility of unconstrained three-material decomposition: imaging an excised human heart using a prototype silicon photon-counting CT detector

Rationale and objectives

The purpose of this study was to evaluate the feasibility of unconstrained three-material decomposition in a human tissue specimen containing iodinated contrast agent, using an experimental multi-bin photon-counting silicon detector. It was further to evaluate potential added clinical value compared to a 1st-generation state-of-the-art dual-energy computed tomography system.

Materials and methods

A prototype photon-counting silicon detector in a bench-top setup for x-ray tomographic imaging was calibrated using a multi-material calibration phantom. A heart with calcified plaque was obtained from a deceased patient, and the coronary arteries were injected with an iodinated contrast agent mixed with gelatin. The heart was imaged in the experimental setup and on a 1st-generation state-of-the-art dual-energy computed tomography system. Projection-based three-material decomposition without any constraints was performed with the photon-counting detector data, and the resulting images were compared with those obtained from the dual-energy system.

Results

The photon-counting detector images show better separation of iodine and calcium compared to the dual-energy images. Additional experiments confirmed that unbiased estimates of soft tissue, calcium, and iodine could be achieved without any constraints.

Conclusion

The proposed experimental system could provide added clinical value compared to current dual-energy systems for imaging tasks where mix-up of iodine and calcium is an issue, and the anatomy is sufficiently small to allow iodine to be differentiated from calcium. Considering its previously shown count rate capability, these results show promise for future integration of this detector in a clinical CT scanner.

Key Points

• Spectral photon-counting detectors can solve some of the fundamental problems with conventional single-energy CT.

• Dual-energy methods can be used to differentiate iodine and calcium, but to do so must rely on constraints, since solving for three unknowns with only two measurements is not possible. Photon-counting detectors can improve upon these methods by allowing unconstrained three-material decomposition.

• A prototype photon-counting silicon detector with high count rate capability allows performing unconstrained three-material decomposition and qualitatively shows better differentiation of iodine and calcium than dual-energy CT.

Novel multi-energy X-ray imaging methods: Experimental results of new image processing techniques to improve material separation in computed tomography and direct radiography

We present novel multi-energy X-ray imaging methods for direct radiography and computed tomography. The goal is to determine the contribution of thickness, mass density and atomic composition to the measured X-ray absorption in the sample. Algorithms have been developed by our own to calculate new X-ray images using data from an unlimited amount of scans/images of different tube voltages by pixelwise fitting of the detected gray levels. The resulting images then show a contrast that is influenced either by the atomic number of the elements in the sample (photoelectric interactions) or by the mass density (Compton scattering). For better visualization, those images can be combined to a color image where different materials can easily be distinguished. In the case of computed tomography no established true multi-energy methodology that does not require an energy sensitive detector is known to the authors. The existing dual-energy methods often yield noisy results that need spatial averaging for clear interpretation. The goal of the method presented here is to qualitatively calculate atomic number and mass density images without loosing resolution while reducing the noise by the use of more than two X-ray energies. The resulting images are generated without the need of calibration stan-dards in an automatic and fast data processing routine. They provide additional information that might be of special interest in cases like archaeology where the destruction of a sample to determine its composition is no option, but a increase in measurement time is of little concern.

Spectral Differential Phase Contrast X-Ray Radiography

We investigate the combination of two emerging X-ray imaging technologies, namely spectral imaging and differential phase contrast imaging. By acquiring spatially and temporally registered images with several different X-ray spectra, spectral imaging can exploit differences in the energy-dependent attenuation to generate material selective images. Differential phase contrast imaging uses an entirely different contrast generation mechanism: The phase shift that an X-ray wave exhibits when traversing an object. As both methods can determine the (projected) electron density, we propose a novel material decomposition algorithm that uses the spectral and the phase contrast information simultaneously. Numerical experiments show that the combination of these two imaging techniques benefits from the strengths of the individual methods while the weaknesses are mitigated: Quantitatively accurate basis material images are obtained and the noise level is strongly reduced, compared to conventional spectral X-ray imaging.

Robust Energy Calibration Technique for Photon Counting Spectral Detectors

This paper describes the implementation of a novel and robust threshold energy calibration method for photon counting detectors using polychromatic X-ray tubes. Methods often used for such energy calibration may require re-orientation of the detector or introduce calibration errors that are flux and acquisition time-dependent. Our newly proposed "differential intensity ratios" (DIR) method offers a practical and robust alternative to existing methods. We demonstrate this robustness against photon flux used in calibration, spectral errors such as pulse pile-up as well as the detector's inherent spectral resolution limits. The demonstrated significant insensitivity of the proposed DIR signature to detector spectral distortions and energy resolution is a key finding. The proposed DIR calibration method is demonstrated using Medipix3RX detectors with a CdTe sensor under varying flux conditions. A per pixel calibration using the DIR method has also been implemented to demonstrate an improvement over the global energy resolution of the PCD.

Quantification of multiple mixed contrast and tissue compositions using photon-counting spectral computed tomography

Quantitative material decomposition of multiple mixed, or spatially coincident, contrast agent (gadolinium and iodine) and tissue (calcium and water) compositions is demonstrated using photon-counting spectral computed tomography (CT). Material decomposition is performed using constrained maximum likelihood estimation (MLE) in the image domain. MLE is calibrated by multiple linear regression of all pure material compositions, which exhibits a strong correlation ( R 2 > 0.91 ) between the measured x-ray attenuation in each photon energy bin and known concentrations in the calibration phantom. Material decomposition of mixed compositions in the sample phantom provides color material concentration maps that clearly identify and differentiate each material. The measured area under the receiver operating characteristic curve is > 0.95 , indicating highly accurate material identification. Material decomposition also provides accurate quantitative estimates of material concentrations in mixed compositions with a root-mean-squared error < 12 % of the maximum concentration for each material. Thus, photon-counting spectral CT enables quantitative molecular imaging of multiple spatially coincident contrast agent (gadolinium and iodine) and tissue (calcium and water) compositions, which is not possible with current clinical molecular imaging modalities, such as nuclear imaging and magnetic resonance imaging.

Joint Statistical Iterative Material Image Reconstruction for Spectral Computed Tomography Using a Semi-Empirical Forward Model

By acquiring tomographic measurements with several distinct photon energy spectra, spectral computed tomography (spectral CT) is able to provide additional material-specific information compared with conventional CT. This information enables the generation of material selective images, which have found various applications in medical imaging. However, material decomposition typically leads to noise amplification and a degradation of the signal-to-noise ratio. This is still a fundamental problem of spectral CT, especially for low-dose medical applications. Inspired by the success for low-dose conventional CT, several statistical iterative reconstruction algorithms for spectral CT have been developed. These algorithms typically rely on detailed knowledge about the spectrum and the detector response. Obtaining this knowledge is often difficult in practice, especially if photon counting detectors are used to acquire the energy specific information. In this paper, a new algorithm for joint statistical iterative material image reconstruction is presented. It relies on a semi-empirical forward model which is tuned by calibration measurements. This strategy allows to model spatially varying properties of the imaging system without requiring detailed prior knowledge of the system parameters. We employ an efficient optimization algorithm based on separable surrogate functions to accelerate convergence and reduce the reconstruction time. Numerical as well as real experiments show that our new algorithm leads to reduced statistical bias and improved image quality compared with projection-based material decomposition followed by analytical or iterative image reconstruction.

Experimental study to optimize configurations of PCD Spectral CT

BACKGROUND

High dose efficiency of photon counting detector based spectral CT (PCD-SCT) and its value in some clinical diagnosis have been well acknowledged. However, it has not been widely adopted in practical use for medical diagnosis and security inspection.

OBJECTIVE

To evaluate the influence on PCD-SCT from multiple aspects including the number of energy channels, k-edge materials, energy thresholding, basis functions in spectral information decomposition, and the combined optimal setting for these parameters and configurations.

METHODS

Basis material decomposition after spatial reconstruction is applied for PCD-SCT. A "one-step" synthesis method, merging decomposition with synthesis, is proposed to obtain virtual monochromatic images. An I-RMSE is computed using the bias part of I-RMSE to describe the difference of a synthesized signal from ground truth and the standard deviation part of I-RMSE to express the noise level. In addition, virtual monochromatic images commonly used in the medical area are also synthesized. Both numerical simulations and practical experiments are conducted for validation.

RESULTS

Results indicated that the I-RMSE for matters significantly reduced with an increased number of energy channels compared with dual-energy channel. The maximum reduction is 6% for triple-, 18% for quadruple-and 24% for quintuple-energy, respectively. However, the improvement is not linear, and also slows down after the number of energy channels reaches a certain number. Contrast agents of high concentration can introduce up to 50% error to surrounding matters. Moreover, different energy partitions influence the total error, which demonstrates the necessity of energy threshold optimization. Last, the optimal basis-material combination varies according to targeted imaging matters and the interested monochromatic energies.

CONCLUSIONS

Gain from more energy channels could be significant with the increase of energy channel number. Introduction of contrast agents in scanned objects will increase overall error in spectral CT imaging. Energy thresholding optimization is beneficial for information recovery. Moreover, the choice of basis materials could also be important to obtain low noise results. With these studies of the effect from various configurations for PCD-SCT, one may optimize the configuration of PCD-SCT accordingly.

Tutorial on X-ray photon counting detector characterization

BACKGROUND

Recent advances in photon counting detection technology have led to significant research interest in X-ray imaging.

OBJECTIVE

As a tutorial level review, this paper covers a wide range of aspects related to X-ray photon counting detector characterization.

METHODS

The tutorial begins with a detailed description of the working principle and operating modes of a pixelated X-ray photon counting detector with basic architecture and detection mechanism. Currently available methods and techniques for charactering major aspects including energy response, noise floor, energy resolution, count rate performance (detector efficiency), and charge sharing effect of photon counting detectors are comprehensively reviewed. Other characterization aspects such as point spread function (PSF), line spread function (LSF), contrast transfer function (CTF), modulation transfer function (MTF), noise power spectrum (NPS), detective quantum efficiency (DQE), bias voltage, radiation damage, and polarization effect are also remarked.

RESULTS

A cadmium telluride (CdTe) pixelated photon counting detector is employed for part of the characterization demonstration and the results are presented.

CONCLUSIONS

This review can serve as a tutorial for X-ray imaging researchers and investigators to understand, operate, characterize, and optimize photon counting detectors for a variety of applications.

Dual Energy CT Kidney Stone Differentiation in Photon Counting Computed Tomography

This study evaluates the capabilities of a whole-body photon counting CT system to differentiate between four common kidney stone materials, namely uric acid (UA), calcium oxalate monohydrate (COM), cystine (CYS),and apatite (APA) ex vivo. Two different x-ray spectra (120 kV and 140 kV) were applied and two acquisition modes were investigated; The macro-mode generates two energy threshold based image-volumes and two energy bin based image-volumes. In the chesspattern-mode, however, four energy thresholds are applied. A virtual low energy image, as well as a virtual high energy image are derived from initial threshold-based images, while considering their statistically correlated nature. The energy bin based images of the macro-mode, as well as the virtual low and high energy image of the chesspattern-mode serve as input for our dual energy evaluation. The dual energy ratio of the individually segmented kidney stones were utilized to quantify the discriminability of the different materials. The dual energy ratios of the two acquisition modes showed high correlation for both applied spectra. Wilcoxon-rank sum tests and the evaluation of the area under the receiver operating characteristics curves suggest that the UA kidney stones are best differentiable from all other materials (AUC = 1.0), followed by CYS (AUC ≈ 0.9 compared against COM and APA). COM and APA, however, are hardly distinguishable (AUC between 0.63 and 0.76). The results hold true for the measurements of both spectra and both acquisition modes.

Numerical observer for atherosclerotic plaque classification in spectral computed tomography

Spectral computed tomography (SCT) generates better image quality than conventional computed tomography (CT). It has overcome several limitations for imaging atherosclerotic plaque. However, the literature evaluating the performance of SCT based on objective image assessment is very limited for the task of discriminating plaques. We developed a numerical-observer method and used it to assess performance on discrimination vulnerable-plaque features and compared the performance among multienergy CT (MECT), dual-energy CT (DECT), and conventional CT methods. Our numerical observer was designed to incorporate all spectral information and comprised two-processing stages. First, each energy-window domain was preprocessed by a set of localized channelized Hotelling observers (CHO). In this step, the spectral image in each energy bin was decorrelated using localized prewhitening and matched filtering with a set of Laguerre-Gaussian channel functions. Second, the series of the intermediate scores computed from all the CHOs were integrated by a Hotelling observer with an additional prewhitening and matched filter. The overall signal-to-noise ratio (SNR) and the area under the receiver operating characteristic curve (AUC) were obtained, yielding an overall discrimination performance metric. The performance of our new observer was evaluated for the particular binary classification task of differentiating between alternative plaque characterizations in carotid arteries. A clinically realistic model of signal variability was also included in our simulation of the discrimination tasks. The inclusion of signal variation is a key to applying the proposed observer method to spectral CT data. Hence, the task-based approaches based on the signal-known-exactly/background-known-exactly (SKE/BKE) framework and the clinical-relevant signal-known-statistically/background-known-exactly (SKS/BKE) framework were applied for analytical computation of figures of merit (FOM). Simulated data of a carotid-atherosclerosis patient were used to validate our methods. We used an extended cardiac-torso anthropomorphic digital phantom and three simulated plaque types (i.e., calcified plaque, fatty-mixed plaque, and iodine-mixed blood). The images were reconstructed using a standard filtered backprojection (FBP) algorithm for all the acquisition methods and were applied to perform two different discrimination tasks of: (1) calcified plaque versus fatty-mixed plaque and (2) calcified plaque versus iodine-mixed blood. MECT outperformed DECT and conventional CT systems for all cases of the SKE/BKE and SKS/BKE tasks (all [Formula: see text]). On average of signal variability, MECT yielded the SNR improvements over other acquisition methods in the range of 46.8% to 65.3% (all [Formula: see text]) for FBP-Ramp images and 53.2% to 67.7% (all [Formula: see text]) for FBP-Hanning images for both identification tasks. This proposed numerical observer combined with our signal variability framework is promising for assessing material characterization obtained through the additional energy-dependent attenuation information of SCT. These methods can be further extended to other clinical tasks such as kidney or urinary stone identification applications.

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

We develop a primal-dual algorithm that allows for one-step inversion of spectral CT transmission photon counts data to a basis map decomposition. The algorithm allows for image constraints to be enforced on the basis maps during the inversion. The derivation of the algorithm makes use of a local upper bounding quadratic approximation to generate descent steps for non-convex spectral CT data discrepancy terms, combined with a new convex-concave optimization algorithm. Convergence of the algorithm is demonstrated on simulated spectral CT data. Simulations with noise and anthropomorphic phantoms show examples of how to employ the constrained one-step algorithm for spectral CT data.

Volumetric CT with sparse detector arrays (and application to Si-strip photon counters)

Novel x-ray medical imaging sensors, such as photon counting detectors (PCDs) and large area CCD and CMOS cameras can involve irregular and/or sparse sampling of the detector plane. Application of such detectors to CT involves undersampling that is markedly different from the commonly considered case of sparse angular sampling. This work investigates volumetric sampling in CT systems incorporating sparsely sampled detectors with axial and helical scan orbits and evaluates performance of model-based image reconstruction (MBIR) with spatially varying regularization in mitigating artifacts due to sparse detector sampling. Volumetric metrics of sampling density and uniformity were introduced. Penalized-likelihood MBIR with a spatially varying penalty that homogenized resolution by accounting for variations in local sampling density (i.e. detector gaps) was evaluated. The proposed methodology was tested in simulations and on an imaging bench based on a Si-strip PCD (total area 5 cm  ×  25 cm) consisting of an arrangement of line sensors separated by gaps of up to 2.5 mm. The bench was equipped with translation/rotation stages allowing a variety of scanning trajectories, ranging from a simple axial acquisition to helical scans with variable pitch. Statistical (spherical clutter) and anthropomorphic (hand) phantoms were considered. Image quality was compared to that obtained with a conventional uniform penalty in terms of structural similarity index (SSIM), image uniformity, spatial resolution, contrast, and noise. Scan trajectories with intermediate helical width (~10 mm longitudinal distance per 360° rotation) demonstrated optimal tradeoff between the average sampling density and the homogeneity of sampling throughout the volume. For a scan trajectory with 10.8 mm helical width, the spatially varying penalty resulted in significant visual reduction of sampling artifacts, confirmed by a 10% reduction in minimum SSIM (from 0.88 to 0.8) and a 40% reduction in the dispersion of SSIM in the volume compared to the constant penalty (both penalties applied at optimal regularization strength). Images of the spherical clutter and wrist phantoms confirmed the advantages of the spatially varying penalty, showing a 25% improvement in image uniformity and 1.8  ×  higher CNR (at matched spatial resolution) compared to the constant penalty. The studies elucidate the relationship between sampling in the detector plane, acquisition orbit, sampling of the reconstructed volume, and the resulting image quality. They also demonstrate the benefit of spatially varying regularization in MBIR for scenarios with irregular sampling patterns. Such findings are important and integral to the incorporation of a sparsely sampled Si-strip PCD in CT imaging.

Theoretical comparison of a dual energy system and photon counting silicon detector used for material quantification in spectral CT

Any method using dual energy computed tomography (CT) has to make prior assumptions in order to quantify k-edge contrast agents. This work estimates the mean square error (MSE) in contrast agent quantification employing a method based on assigning each reconstructed voxel a ratio of soft tissue and fat using dual energy CT. The results are compared to the MSE using a photon counting silicon detector with multiple bins. The square root of the MSEs of the quantifications of iodine and gadolinium for an object consisting of soft tissue and fat using the silicon detector and dual energy CT range from below 2% and 1% of the contrast agent content for 100 mg/cm(3) of iodine and gadolinium, up to approximately 10% and 13%, and 6% and 4%, for 5 mg/cm(3) of iodine and gadolinium, respectively. When adding bone with a voxel volume fraction of 2.2%, the square root of the MSEs of the quantifications of iodine and gadolinium using dual energy CT increases to 25% and 6%, respectively, for 5 mg/cm(3) of contrast agent. In conclusion, results indicate that the noise levels of the material quantification using the silicon detector are higher than the noise levels using a dual energy CT when the composition of the object is known. However, using a dual energy CT increases the risk of model specification error and subsequently a large bias in contrast agent quantification, a problem which does not exist when using a multi-bin CT where the number of energy bins is larger than two.

Efficacy of fixed filtration for rapid kVp-switching dual energy x-ray systems

PURPOSE

Dose efficiency of dual kVp imaging can be improved if the two beams are filtered to remove photons in the common part of their spectra, thereby increasing spectral separation. While there are a number of advantages to rapid kVp-switching for dual energy, it may not be feasible to have two different filters for the two spectra. Therefore, the authors are interested in whether a fixed added filter can improve the dose efficiency of kVp-switching dual energy x-ray systems.

METHODS

The authors hypothesized that a K-edge filter would provide the energy selectivity needed to remove overlap of the spectra and hence increase the precision of material separation at constant dose. Preliminary simulations were done using calcium and water basis materials and 80 and 140 kVp x-ray spectra. Precision of the decomposition was evaluated based on the propagation of the Poisson noise through the decomposition function. Considering availability and cost, the authors chose a commercial Gd2O2S screen as the filter for their experimental validation. Experiments were conducted on a table-top system using a phantom with various thicknesses of acrylic and copper and 70 and 125 kVp x-ray spectra. The authors kept the phantom exposure roughly constant with and without filtration by adjusting the tube current. The filtered and unfiltered raw data of both low and high energy were decomposed into basis material and the variance of the decomposition for each thickness pair was calculated. To evaluate the filtration performance, the authors measured the ratio of material decomposition variance with and without filtration.

RESULTS

Simulation results show that the ideal filter material depends on the object composition and thickness, and ranges across the lanthanide series, with higher atomic number filters being preferred for more attenuating objects. Variance reduction increases with filter thickness, and substantial reductions (40%) can be achieved with a 2× loss in intensity. The authors' experimental results validate the simulations, yet were overall slightly worse than expectation. For large objects, conventional (non-K-edge) beam hardening filters perform well.

CONCLUSIONS

This study demonstrates the potential of fixed K-edge filtration to improve the dose efficiency and material decomposition precision for rapid kVp-switching dual energy systems.

Cascaded systems analysis of photon counting detectors

PURPOSE

Photon counting detectors (PCDs) are an emerging technology with applications in spectral and low-dose radiographic and tomographic imaging. This paper develops an analytical model of PCD imaging performance, including the system gain, modulation transfer function (MTF), noise-power spectrum (NPS), and detective quantum efficiency (DQE).

METHODS

A cascaded systems analysis model describing the propagation of quanta through the imaging chain was developed. The model was validated in comparison to the physical performance of a silicon-strip PCD implemented on an experimental imaging bench. The signal response, MTF, and NPS were measured and compared to theory as a function of exposure conditions (70 kVp, 1-7 mA), detector threshold, and readout mode (i.e., the option for coincidence detection). The model sheds new light on the dependence of spatial resolution, charge sharing, and additive noise effects on threshold selection and was used to investigate the factors governing PCD performance, including the fundamental advantages and limitations of PCDs in comparison to energy-integrating detectors (EIDs) in the linear regime for which pulse pileup can be ignored.

RESULTS

The detector exhibited highly linear mean signal response across the system operating range and agreed well with theoretical prediction, as did the system MTF and NPS. The DQE analyzed as a function of kilovolt (peak), exposure, detector threshold, and readout mode revealed important considerations for system optimization. The model also demonstrated the important implications of false counts from both additive electronic noise and charge sharing and highlighted the system design and operational parameters that most affect detector performance in the presence of such factors: for example, increasing the detector threshold from 0 to 100 (arbitrary units of pulse height threshold roughly equivalent to 0.5 and 6 keV energy threshold, respectively), increased the f50 (spatial-frequency at which the MTF falls to a value of 0.50) by ∼30% with corresponding improvement in DQE. The range in exposure and additive noise for which PCDs yield intrinsically higher DQE was quantified, showing performance advantages under conditions of very low-dose, high additive noise, and high fidelity rejection of coincident photons.

CONCLUSIONS

The model for PCD signal and noise performance agreed with measurements of detector signal, MTF, and NPS and provided a useful basis for understanding complex dependencies in PCD imaging performance and the potential advantages (and disadvantages) in comparison to EIDs as well as an important guide to task-based optimization in developing new PCD imaging systems.

Dual-energy cone-beam CT with a flat-panel detector: effect of reconstruction algorithm on material classification

PURPOSE

Cone-beam CT (CBCT) with a flat-panel detector (FPD) is finding application in areas such as breast and musculoskeletal imaging, where dual-energy (DE) capabilities offer potential benefit. The authors investigate the accuracy of material classification in DE CBCT using filtered backprojection (FBP) and penalized likelihood (PL) reconstruction and optimize contrast-enhanced DE CBCT of the joints as a function of dose, material concentration, and detail size.

METHODS

Phantoms consisting of a 15 cm diameter water cylinder with solid calcium inserts (50-200 mg/ml, 3-28.4 mm diameter) and solid iodine inserts (2-10 mg/ml, 3-28.4 mm diameter), as well as a cadaveric knee with intra-articular injection of iodine were imaged on a CBCT bench with a Varian 4343 FPD. The low energy (LE) beam was 70 kVp (+0.2 mm Cu), and the high energy (HE) beam was 120 kVp (+0.2 mm Cu, +0.5 mm Ag). Total dose (LE+HE) was varied from 3.1 to 15.6 mGy with equal dose allocation. Image-based DE classification involved a nearest distance classifier in the space of LE versus HE attenuation values. Recognizing the differences in noise between LE and HE beams, the LE and HE data were differentially filtered (in FBP) or regularized (in PL). Both a quadratic (PLQ) and a total-variation penalty (PLTV) were investigated for PL. The performance of DE CBCT material discrimination was quantified in terms of voxelwise specificity, sensitivity, and accuracy.

RESULTS

Noise in the HE image was primarily responsible for classification errors within the contrast inserts, whereas noise in the LE image mainly influenced classification in the surrounding water. For inserts of diameter 28.4 mm, DE CBCT reconstructions were optimized to maximize the total combined accuracy across the range of calcium and iodine concentrations, yielding values of ∼ 88% for FBP and PLQ, and ∼ 95% for PLTV at 3.1 mGy total dose, increasing to ∼ 95% for FBP and PLQ, and ∼ 98% for PLTV at 15.6 mGy total dose. For a fixed iodine concentration of 5 mg/ml and reconstructions maximizing overall accuracy across the range of insert diameters, the minimum diameter classified with accuracy >80% was ∼ 15 mm for FBP and PLQ and ∼ 10 mm for PLTV, improving to ∼ 7 mm for FBP and PLQ and ∼ 3 mm for PLTV at 15.6 mGy. The results indicate similar performance for FBP and PLQ and showed improved classification accuracy with edge-preserving PLTV. A slight preference for increased smoothing of the HE data was found. DE CBCT discrimination of iodine and bone in the knee was demonstrated with FBP and PLTV at 6.2 mGy total dose.

CONCLUSIONS

For iodine concentrations >5 mg/ml and detail size ∼ 20 mm, material classification accuracy of >90% was achieved in DE CBCT with both FBP and PL at total doses <10 mGy. Optimal performance was attained by selection of reconstruction parameters based on the differences in noise between HE and LE data, typically favoring stronger smoothing of the HE data, and by using penalties matched to the imaging task (e.g., edge-preserving PLTV in areas of uniform enhancement).

Dual-energy CT spectral and energy weighted data sets: carotid stenosis and plaque detection

RATIONALE AND OBJECTIVES

To evaluate dual-energy computed tomography spectral and energy weighted (W) image data sets concerning carotid stenosis and calcified plaque detection.

MATERIALS AND METHODS

Ninety patients were evaluated using contrast media injection. Energy used for dual-energy computed tomography was tin filter with 140/80 kiloVoltage and effective milli Ampere second was 147.4/270.6. Image reconstruction was performed using D30f kernel and 0.0, 0.3, 0.6, 0.8, 1.0 weightings. Data sets were analyzed using both qualitative and quantitative methods.

RESULTS

The signal-to-noise ratio, contrast-to-noise ratio, and figure-of-merit were significantly higher in 0.6-W compared to 140-kV or 80-kV data (all P < .05). Plaque thickness, span, and longitudinal diameters were different for 140-kV, 0.6-W, and 80-kV data (all P < .05). Stenotic intra-luminal diameter was significantly different among 140 kV, 0.6 W, and 80 kV data (all P < .05). A comparison between 0.6 W and digital subtraction angiography was nonsignificant (P > .05) in normal lumen measurement.

CONCLUSIONS

The dimension of calcified plaque and carotid artery with contrast media decreased with increased energy. The percentage of carotid artery stenosis does not vary with different energy. Care must be taken for procedural planning like sizing of stents. Measured diameters of the 0.6 W were close to digital subtraction angiography; we suggest that planning should be based on the images acquired using 0.6 weighting.