A dynamic attenuator improves spectral imaging with energy-discriminating, photon counting detectors

Spectral Photon Counting CT: Imaging Algorithms and Performance Assessment

Photon counting x-ray detectors (PCDs) with spectral capabilities have the potential to revolutionize computed tomography (CT) for medical imaging. The ideal PCD provides accurate energy information for each incident x-ray, and at high spatial resolution. This information enables material-specific imaging, enhanced radiation dose efficiency, and improved spatial resolution in CT images. In practice, PCDs are affected by non-idealities, including limited energy resolution, pulse pileup, and cross talk due to charge sharing, K-fluorescence, and Compton scattering. In order to maximize their performance, PCDs must be carefully designed to reduce these effects and then later account for them during correction and post-acquisition steps. This review article examines algorithms for using PCDs in spectral CT applications, including how non-idealities impact image quality. Performance assessment metrics that account for spatial resolution and noise such as the detective quantum efficiency (DQE) can be used to compare different PCD designs, as well as compare PCDs with conventional energy integrating detectors (EIDs). These methods play an important role in enhancing spectral CT images and assessing the overall performance of PCDs.

Anomalous edge response of cadmium telluride-based photon counting detectors jointly caused by high-flux radiation and inter-pixel communication

This work reports an edge enhancing effect experimentally observed in cadmium telluride (CdTe)-based photon counting detector (PCD) systems operated under the charge summing (CS) mode and irradiated by high-flux x-rays. Experimental measurements of the edge spread functions (ESFs) of a PCD system (100μm pixel size, 88 ns deadtime) were performed at different input flux levels from 4.5 × 105count per second (cps) mm-2to 1.5 × 109cps mm-2for the single pixel mode (SP) and the CS mode. A theoretical model that incorporates the impacts of inter-pixel communications and the arbitration process involved in the CS mode was developed to help explain the physical origin of the observed edge enhancing effect. Compared with the monotonically increasing ESF of the SP mode, the ESF of the CS mode measured at high-flux levels shows a peak at an intermediate location (50μm from the edge). The peak became more pronounced with increasing flux levels. The theoretically calculated ESFs agreed well with experimental results with relative errors less than 5% at all flux levels and tested. These results indicate that the anomalous edge enhancing effect is jointly caused by the pileup effect and the CS circuit that introduces negative correlations between adjacent pixels. When the input flux is high enough to deliver photons to multiple adjacent pixels within the same deadtime period, the CS mode may treat the coincident x-rays as shared charges and thus introduce count losses in addition to the well-known pileup count loss. When a high contrast object partially blocks certain pixels from x-rays, the adjacent unblocked pixels have an increased probability of registering counts as a result of the negative correlation. This leads to a peak on the ESF at a pixel-to-edge distance half of the pixel pitch.

Assessment of Multienergy Interpixel Coincidence Counters (MEICC) for Charge Sharing Correction or Compensation for Photon Counting Detectors With Boxcar Signals

Recently, multi-energy inter-pixel coincidence counter (MEICC) has been proposed for charge sharing correction and compensation for photon counting detectors (PCDs), which uses energy-dependent coincidence counters to record coincident events between multiple energy windows of a pixel-of-interest and those of neighboring pixels. A Monte Carlo (MC) simulation study was performed to assess the performance of MEICC; however, the performance might have been overestimated in a previous study. The charge sharing increases the number of photons recorded at a PCD pixel at the expense of the spatial resolution, and therefore, when spatially uniform flat-field x-ray signals are used, it gives PCDs with charge sharing more signals than a PCD without charge sharing. In this paper, we propose to use spatially modulated boxcar signals for evaluating the performances for high spatial frequency tasks because they provide consistent signals regardless of the presence of absence of charge sharing. The flat-field signals must be used for low spatial frequency tasks. We assessed the performances of MEICC and other PCDs with both flat-field signals and boxcar signals, with optimal threshold energies, and with two different pixel sizes. As it is expected, normalized Cramér-Rao lower bounds (nCRLBs) measured with the boxcar signals were worse than those with flat-field signals in general. The nCRLBs of MEICC with 225-μm pixel were close to the current 450-μm PCD. We studied a combination of flat-field signals and N×N super-pixels, where the output of N×N pixels were added, using an MC simulation and a simple charge sharing counting model. The study showed that charge sharing had two opposing impacts on the conventional CT imaging-a negative impact with double-counting among N×N pixels and a positive impact with single-counting spill-in and spill-out across the super-pixel boundary-and the positive impact diminished with increasing N. A use of large N×N super-pixels such as N≥25 was suggested to approximate the zero-frequency detection quantum efficiency of PCD with charge sharing.

Multi-energy inter-pixel coincidence counters for charge sharing correction and compensation in photon counting detectors

PURPOSE

Smaller pixel sizes of x-ray photon counting detectors (PCDs) are advantageous for count rate capabilities but disadvantageous for charge sharing. With charge sharing, the energy of an x-ray photon may be split and one photon may produce two or more counts at adjacent pixels, both at lower energies than the incident energy. This "double-counting" increases noise variance and degrades the spectral response. Overall, it has a significantly negative impact on the performance of PCD-based computed tomography (CT). Charge sharing is induced by the detection physics and occurs regardless of count rates; thus, it is impossible to avoid. We propose in this paper a method that has a potential to address both noise and bias added by charge sharing.

METHODS

We propose applying a multi-energy inter-pixel coincidence counter (MEICC) technique, which uses energy-dependent coincidence counters, keeps the book of charge sharing events during data acquisition, and provides the exact number of charge sharing occurrences, which can be used to either correct or compensate for them after the acquisition is completed. MEICC does not interfere with the primary counting process; therefore, PCDs with MEICC will remain as fast as those without MEICC. MEICC can be implemented using current electronics technology because its inter-pixel coincidence counters used to handle digital data are rather simple. We evaluated Cramér-Rao lower bound (CRLB) of PCDs with and without MEICC using a Monte Carlo simulation.

RESULTS

When the number of energy windows was four or larger and eight neighboring pixels were used, the CRLBs of 225-µm PCD with MEICC normalized by those of the current PCD with the same number of windows were 0.361-0.383 for water density images of two basis functions, which was only 5.7-16.4% worse than those of a PCD without charge sharing (which were at 0.329-0.358). In contrast, the normalized CRLBs of the PCD with one coincidence counter were 0.466-0.499, which were 37.3-45.6% worse than the PCD without charge sharing. The use of eight neighboring pixels provided ~10% better CRLB values than four neighboring pixels for MEICC. With four energy windows, decreasing the number of coincidence counters from 16 to 9 only slightly increased the CRLB from 0.255 to 0.269 (which corresponded to as little as a 5.5% change). The normalized CRLBs of MEICC for K-edge imaging (gold) were 0.295-0.426, while those of the one coincidence counter were 0.926-0.959 and the ideal PCDs were 0.126-0.146.

CONCLUSIONS

The proposed MEICC provides spectral information that can be used to address charge sharing problems in PCDs and is expected to satisfy the requirements for clinical x-ray CT. MEICC is very effective, especially for K-edge imaging, which requires accurate spectral information.

Fluid-filled dynamic bowtie filter: Description and comparison with other modulators

PURPOSE

A dynamic bowtie filter can modulate flux along both fan and view angles for reduced patient dose, scatter, and required photon flux, which is especially important for photon counting detectors (PCDs). Among the proposed dynamic bowtie designs, the piecewise-linear attenuator (Hsieh and Pelc, Med Phys. 2013;40:031910) offers more flexibility than conventional filters, but relies on analog positioning of a limited number of wedges. In this work, we study our previously proposed dynamic attenuator design, the fluid-filled dynamic bowtie filter (FDBF) that has digital control. Specifically, we use computer simulations to study fluence modulation, reconstructed image noise, and radiation dose and to compare it to other attenuators. FDBF is an array of small channels each of which, if it can be filled with dense fluid or emptied quickly, has a binary effect on the flux. The cumulative attenuation from each channel along the x-ray path contributes to the FDBF total attenuation.

METHODS

An algorithm is proposed for selecting which FDBF channels should be filled. Two optimization metrics are considered: minimizing the maximum-count-rate for PCDs and minimizing peak-variance for energy-integrating detectors (EIDs) at fixed radiation dose (for optimizing dose efficiency). Using simulated chest, abdomen, and shoulder data, the performance is compared with a conventional bowtie and a piecewise-linear attenuator. For minimizing peak-variance, a perfect-attenuator (hypothetical filter capable of adjusting the fluence of each ray individually) and flat-variance attenuator are also included in the comparison. Two possible fluids, solutions of zinc bromide and gadolinium chloride, were tested.

RESULTS

To obtain the same SNR as routine clinical protocols, the proposed FDBF reduces the maximum-count-rate (across projection data, averaged over the test objects) of PCDs to 1.2 Mcps/mm² , which is 55.8 and 3.3 times lower than the max-count-rate of the conventional bowtie and the piecewise-linear bowtie, respectively. (Averaged across objects for FDBF, the max-count-rate without object and FDBF is 2063.5 Mcps/mm² , and the max-count-rate with object without FDBF is 749.8 Mcps/mm² .) Moreover, for the peak-variance analysis, the FDBF can reduce entrance-energy-fluence (sum of energy incident on objects, used as a surrogate for dose) to 34% of the entrance-energy-fluence from the conventional filter on average while achieving the same peak noise level. Its entrance-energy-fluence reduction performance is only 7% worse than the perfect-attenuator on average and is 13% better than the piecewise-linear filter for chest and shoulder. Furthermore, the noise-map in reconstructed image domain from the FDBF is more uniform than the piecewise-linear filter, with 3 times less variation across the object. For the dose reduction task, the zinc bromide solution performed slightly poorer than stainless steel but was better than the gadolinium chloride solution.

CONCLUSIONS

The FDBF allows finer control over flux distribution compared to piecewise-linear and conventional bowtie filters. It can reduce the required maximum-count-rate for PCDs to a level achievable by current detector designs and offers a high dose reduction factor.

Digital count summing vs analog charge summing for photon counting detectors: A performance simulation study

PURPOSE

Charge sharing is a significant problem for CdTe-based photon counting detectors (PCDs) and can cause high-energy photons to be misclassified as one or more low-energy events. Charge sharing is especially problematic in PCDs for CT because the high flux necessitates small pixels, which increase the magnitude of charge sharing. Analog charge summing (ACS) is a powerful solution to reduce spectral distortion arising from charge sharing but may be difficult to implement. We investigate correction of the signal after digitization by the comparator ("digital count summing"), which is only able to correct a subset of charge sharing events but may have implementation advantages. We compare and quantify the relative performance of digital and ACS in simulations.

METHODS

Transport of photons in CdTe was modeled using Monte Carlo simulations. Energy deposited in the CdTe substrate was converted to electrical charges of a predetermined shape, and all charges within the detector pixel are assumed to be perfectly collected. In ACS, the maximum charge received over any 2 × 2 block of pixels was grouped together prior to digitization. In digital count summing (DCS), the charge was digitized in each pixel, and subsequently, adjacent pixels that detected events grouped their charge to record a single, higher energy event. All simulations were performed at the limit of low flux (no pileup). The default tube voltage was 120 kVp, object thickness was 20 cm of water, pixel pitch was 250 μm, and charge cloud modeled as a Gaussian with σ = 40 μm. Variation of these parameters was examined in a sensitivity analysis.

RESULTS

Detectors that used no correction, DCS, and ACS misclassified 51%, 39%, and 15% of incident photons, respectively. For iodine basis material imaging, DCS exhibited 100% greater dose efficiency compared to uncorrected, and ACS exhibited an additional 111% greater dose efficiency compared to digital charge summing. For a nonspectral task, the dose efficiency improvement as estimated by improvement of zero-frequency detective quantum efficiency, DQE(0) was 10% for DCS compared to uncorrected and 10% for ACS compared to DCS. A sensitivity analysis showed that DCS generally achieved half the benefit of ACS over a range of conditions, although the benefit was markedly less if the charge cloud was instead modeled as a small sphere.

CONCLUSIONS

Summing of counts after digitization may be a simpler alternative to summing of charge prior to digitization due to the relative complexity of analog circuit design. Over most conditions studied, it provides roughly half the benefit of ACS and may offer certain implementation advantages.

Spectral resolution and high-flux capability tradeoffs in CdTe detectors for clinical CT

PURPOSE

Photon-counting detectors using CdTe or CZT substrates are promising candidates for future CT systems but suffer from a number of nonidealities, including charge sharing and pulse pileup. By increasing the pixel size of the detector, the system can improve charge sharing characteristics at the expense of increasing pileup. The purpose of this work is to describe these considerations in the optimization of the detector pixel pitch.

METHODS

The transport of x rays through the CdTe substrate was simulated in a Monte Carlo fashion using GEANT4. Deposited energy was converted into charges distributed as a Gaussian function with size dependent on interaction depth to capture spreading from diffusion and Coulomb repulsion. The charges were then collected in a pixelated fashion. Pulse pileup was incorporated separately with Monte Carlo simulation. The Cramér-Rao lower bound (CRLB) of the measurement variance was numerically estimated for the basis material projections. Noise in these estimates was propagated into CT images. We simulated pixel pitches of 250, 350, and 450 microns and compared the results to a photon counting detector with pileup but otherwise ideal energy response and an ideal dual-energy system (80/140 kVp with tin filtration). The modeled CdTe thickness was 2 mm, the incident spectrum was 140 kVp and 500 mA, and the effective dead time was 67 ns. Charge summing circuitry was not modeled. We restricted our simulations to objects of uniform thickness and did not consider the potential advantage of smaller pixels at high spatial frequencies.

RESULTS

At very high x-ray flux, pulse pileup dominates and small pixel sizes perform best. At low flux or for thick objects, charge sharing dominates and large pixel sizes perform best. At low flux and depending on the beam hardness, the CRLB of variance in basis material projections tasks can be 32%-55% higher with a 250 micron pixel pitch compared to a 450 micron pixel pitch. However, both are about four times worse in variance than the ideal photon counting detector. The optimal pixel size depends on a number of factors such as x-ray technique and object size. At high technique (140 kVp/500 mA), the ratio of variance for a 450 micron pixel compared to a 250 micron pixel size is 2126%, 200%, 97%, and 78% when imaging 10, 15, 20, and 25 cm of water, respectively. If 300 mg/cm² of iodine is also added to the object, the variance ratio is 117%, 91%, 74%, and 72%, respectively. Nonspectral tasks, such as equivalent monoenergetic imaging, are less sensitive to spectral distortion.

CONCLUSIONS

The detector pixel size is an important design consideration in CdTe detectors. Smaller pixels allow for improved capabilities at high flux but increase charge sharing, which in turn compromises spectral performance. The optimal pixel size will depend on the specific task and on the charge shaping time.

A prototype piecewise-linear dynamic attenuator

The piecewise-linear dynamic attenuator has been proposed as a mechanism in CT scanning for personalizing the x-ray illumination on a patient- and application-specific basis. Previous simulations have shown benefits in image quality, scatter, and dose objectives. We report on the first prototype implementation. This prototype is reduced in scale and speed and is integrated into a tabletop CT system with a smaller field of view (25 cm) and longer scan time (42 s) compared to a clinical system. Stainless steel wedges were machined and affixed to linear actuators, which were in turn held secure by a frame built using rapid prototyping technologies. The actuators were computer-controlled, with characteristic noise of about 100 microns. Simulations suggest that in a clinical setting, the impact of actuator noise could lead to artifacts of only 1 HU. Ring artifacts were minimized by careful design of the wedges. A water beam hardening correction was applied and the scan was collimated to reduce scatter. We scanned a 16 cm water cylinder phantom as well as an anthropomorphic pediatric phantom. The artifacts present in reconstructed images are comparable to artifacts normally seen with this tabletop system. Compared to a flat-field reference scan, increased detectability at reduced dose is shown and streaking is reduced. Artifacts are modest in our images and further refinement is possible. Issues of mechanical speed and stability in the challenging clinical CT environment will be addressed in a future design.

Improving pulse detection in multibin photon-counting detectors

Energy-discriminating, photon-counting (EDPC) detectors are attractive for their potential for improved detective quantum efficiency and for their spectral imaging capabilities. However, at high count rates, counts are lost, the detected spectrum is distorted, and the advantages of EDPC detectors disappear. Existing EDPC detectors identify counts by analyzing the signal with a bank of comparators. We explored alternative methods for pulse detection for multibin EDPC detectors that could improve performance at high count rates. The detector signal was simulated in a Monte Carlo fashion assuming a bipolar shape and analyzed using several methods, including the conventional bank of comparators. For example, one method recorded the peak energy of the pulse along with the width (temporal extent) of the pulse. The Cramer-Rao lower bound of the variance of basis material estimates was numerically found for each method. At high count rates, the variance in water material (bone canceled) measurements could be reduced by as much as an order of magnitude. Improvements in virtual monoenergetic images were modest. We conclude that stochastic noise in spectral imaging tasks could be reduced if alternative methods for pulse detection were utilized.