Photon Counting X-ray Research Articles

Pulse pileup model for spectral resolved X-ray photon-counting detectors with dead time and retrigger capability

We developed an analytical model to evaluate the effect of signal pileup on the recorded energy spectrum in X-ray photon-counting detectors affected by dead time and equipped with retrigger capability. The retrigger function allows the system to work in a specific non-paralyzable counting mode by counting the time-over-threshold of piled-up signals in multiples of a predefined and selectable retrigger time. The model, designed for rectangle-like-shaped signals, allows for arbitrary input energy spectra and can significantly help understand and optimize the behavior of counting detectors with spectral capabilities and retrigger mechanisms in applications involving polychromatic beams, e.g., spectral X-ray imaging and computed tomography (CT), in a time-efficient way. Dedicated numerical simulations were used to validate the model under several conditions of incoming flux and threshold energy, with excellent results.

Empirical optimization of energy bin weights for compressing measurements with realistic photon counting x-ray detectors.

Photon counting detectors (PCDs) provide higher spatial resolution, improved contrast-to-noise ratio (CNR), and energy discriminating capabilities. However, the greatly increased amount of projection data in photon counting computed tomography (PCCT) systems becomes challenging to transmit through the slip ring, process, and store. This study proposes and evaluates an empirical optimization algorithm to obtain optimal energy weights for energy bin data compression. This algorithm is universally applicable to spectral imaging tasks including 2 and 3 material decomposition (MD) tasks and virtual monoenergetic images (VMIs). This method is simple to implement while preserving spectral information for the full range of object thicknesses and is applicable to different PCDs, for example, silicon detectors and CdTe detectors. We used realistic detector energy response models to simulate the spectral response of different PCDs and an empirical calibration method to fit a semi-empirical forward model for each PCD. We numerically optimized the optimal energy weights by minimizing the average relative Cramér-Rao lower bound (CRLB) due to the energy-weighted bin compression, for MD and VMI tasks over a range of material area density (0-40g/cm2 water, 0-2.16g/cm2 calcium). We used Monte Carlo simulation of a step wedge phantom and an anthropomorphic head phantom to evaluate the performance of this energy bin compression method in the projection domain and image domain, respectively. The results show that for 2 MD, the energy bin compression method can reduce PCCT data size by 75% and 60%, with an average variance penalty of less than 17% and 3% for silicon and CdTe detectors, respectively. For 3 MD tasks with a K-edge material (iodine), this method can reduce the data size by 62.5% and 40% with an average variance penalty of less than 12% and 13% for silicon and CdTe detectors, respectively. We proposed an energy bin compression method that is broadly applicable to different PCCT systems and object sizes, with high data compression ratio and little loss of spectral information.

Demonstrating a Novel, Hidden Source of Spectral Distortion in X-ray Photon Counting Detectors and Assessing Novel Trigger Schemes Proposed to Avoid It.

X-ray photon counting spectral imaging (x-CSI) determines a detected photon's energy by comparing the charge it induces with several thresholds, counting how many times each is crossed (the standard method, STD). This paper is the first to demonstrate that this approach can unexpectedly delete counts from the recorded energy spectrum under some clinically relevant conditions: a process we call negative counting. Four alternative counting schemes are proposed and simulated for a wide range of sensor geometries (pixel pitch 100-600 µm, sensor thickness 1-3 mm), number of thresholds (3, 5, 8, 24 and 130) and medically relevant X-ray fluxes (106-109 photons mm-2 s-1). Spectral efficiency and counting efficiency are calculated for each simulation. Performance gains are explained mechanistically and correlated well with the improved suppression of "negative counting". The best performing scheme (Shift Register, SR) entirely eliminates negative counting, remaining close to an ideal scheme at fluxes of up to 108 photons mm-2 s-1. At the highest fluxes considered, the deviation from ideal behaviour is reduced by 2/3 in SR compared with STD. The results have significant implications both for generally improving spectral fidelity and as a possible path toward the 109 photons mm-2 s-1 goal in photon-counting CT.

High-timing-precision detection of single X-ray photons by superconducting nanowires.

Precisely acquiring the timing information of individual X-ray photons is important in both fundamental research and practical applications. The timing precision of commonly used X-ray single-photon detectors remains in the range of one hundred picoseconds to microseconds. In this work, we report on high-timing-precision detection of single X-ray photons through the fast transition to the normal state from the superconductive state of superconducting nanowires. We successfully demonstrate a free-running X-ray single-photon detector with a timing resolution of 20.1ps made of 100-nm-thick niobium nitride film with an active area of 50μm by 50μm. By using a repeated differential timing measurement on two adjacent X-ray single-photon detectors, we demonstrate a precision of 0.87ps in the arrival-time difference of X-ray photon measurements. Therefore, our work significantly enhances the timing precision in X-ray photon counting, opening a new niche for ultrafast X-ray photonics and many associated applications.

Background rates in several hard X-ray photon counting pixel array detectors

Background rates in several pixel array detectors are investigated with an eye toward using them with hard, >6 keV, X-rays in very low-rate experiments - e.g., at signal rates <0.1/s/cm2. Measured background event rates for a detector with an unshielded 0.75 mm thick CdTe sensor on the experimental floor at SPring-8 varied from 0.4/s/cm2 with a 6 keV lower-level discriminator (LLD) threshold to 0.2/s/cm2 with a 75 keV LLD threshold. The background for a detector with a 1 mm thick silicon sensor was smaller, ∼0.08/s/cm2 for a 3keV threshold dropping to ∼0.07/s/cm2 at a 17 keV threshold. These rates are dominated by terrestrial sources, such as gamma rays emitted from trace impurities in concrete, with only a small contribution, <0.01/s/cm2, from direct detection of cosmic ray muons (CRMs). 15 mm of Pb shielding reduces the measured rates to < 0.05/s/cm2 in CdTe and to < 0.02/s/cm2 in 1 mm silicon. Additional processing, “time slicing”, as may be used in low-rate experiments where backgrounds may be most problematic, is shown to reduce the background rates further, to < 0.004/s/cm2 with the silicon sensor, and to between 0.002 and 0.02/s/cm2 for CdTe, where the exact value for CdTe depends sensitively on the detector threshold, and the use of the retriggering and/or the use of dual discriminator thresholds. Background rates are presented as functions of discriminator threshold, shielding, and processing. We also discuss the magnitude of, and the correction for, the limitation of the detector dynamic range that can be introduced by the time slicing. Finally, we present one example where time slicing was used.

Degradation of signal-to-noise ratio in counting detectors due to pile-up effects

This work proposes a methodology to quantify how the non-linearity of counting detectors due to pulse pile-up effects impacts detrimentally the signal-to-noise ratio of measurements and the resulting DQE. The method is presented from a conceptual point of view, justified and validated by both analytical and Monte Carlo methods applied to the ideal theoretical pile-up models usually discussed in the literature. Although the motivation of this work is the application to photon counting X-ray detectors, the study may be relevant for other type of particle or radiation detection systems operating at high count rate regimes, when pile-up is not negligible. The high count rates can be managed in actual experiments by implementing pile-up compensation methods in the detector readout chain, by applying correction algorithms to the measured data or by a combination of both. What has been less discussed in literature is the impact of those pile-up compensation techniques on the final statistical properties of the resulting data, and the proposed method can be potentially used for that purpose. The application of the method to more realistic systems is illustrated with an example of a simulated 2D X-ray photon counting detector that includes all the primary physical effects and the full readout chain.

Characterization of 3.2 Gbps readout in 65 nm CMOS technology

A new class of photon-counting pixel detectors allows for capturing of an image in several photon energy bins in one shot. A decreased pixel pitch and an increased number of energy bins are needed to enhance the spatial and spectral resolution of the detector. This led to new requirements for the readout systems and their bandwidths, as more data is generated for the same detection area. Fast differential serial communication enables high-speed data rates, thus providing an ideal solution to transfer large amounts of data generated by the detector’s front-end electronics. However, its implementation provides extra challenges. This work introduces a novel high-speed serial readout designed in a 65 nm CMOS technology that will be used in the future photon-counting X-ray imaging detectors. The design of the serial transmitter is presented together with the characterization of jitter and channel performance.

Effect of superhydrophobic surfaces on bubble column flow dynamics

Bubbly flow in bubble column reactors promotes mixing necessary for many chemical processes. We show that if superhydrophobic-coated material is introduced into a bubble column, there can be a substantial difference in gas holdup and earlier initiation of churn-turbulent flow which can alter larger-scale mixing without a need to change the superficial gas velocity. Addition of superhydrophobic surface can also cause bubbles to (directly or indirectly assisted by the surface) escape faster to the free surface resulting in a reduced void fraction (i.e., reduced gas holdup). As the flow becomes optically opaque at few percent gas phase volume fraction, we utilize two dual plane wire mesh sensors to obtain velocity profiles and bubble size distributions, in addition to the traditional pressure and level based gas holdup measurements to calculate average phase fraction. Additionally, a custom build photon-counting dual energy threshold X-ray computed tomography system is employed to get a higher resolution measurement of the time average phase fraction non-intrusively. We report satisfactory agreement between these techniques with differences arising for understood reasons, and use the insight thus yielded to discuss the effect of superhydrophobic surfaces on bubble column flow dynamics.

Multi-voltage threshold digitizer using a time-varied threshold for photon-counting X-ray security inspection imaging

In this paper, we present a multi-voltage threshold (MVT) digitizer using a time-varied threshold and investigate the feasibility of using this digitizer to develop the data acquisition system for photon-counting X-ray security inspection imaging. This digitizer compares the scintillation pulse with a preset time-varied voltage signal and records the time information of the intersections. As the amplitude of the threshold varies as a function of time, the voltage values of the intersections could be calculated using the recorded time, and the energy information could be extracted via the acquired time–voltage samples. Unlike the constant voltage thresholds MVT digitizer requiring multiple voltage comparators, this method needs only one voltage comparator, which notably reduces the complexity in configuration. We demonstrate the development of such an MVT digitizer using a sine signal threshold and evaluate its energy performance with a yttrium orthosilicate (YSO) scintillator and Silicon Photomultiplier (SiPM) detector (designed for X-ray detection). The digitizer obtains energy resolutions of 41.4%, 26.6%, and 25.5% at 30, 59.5, and 80.9 keV, respectively. For comparison, the YSO/SiPM detector is also evaluated with a high-speed oscilloscope, and measured energy resolutions are 40.6%, 26.5%, and 25.1% at 30, 59.5, and 80.9 keV, respectively. These results indicate that the proposed MVT digitizer using a time-varied threshold is a promising method for developing a YSO/SiPM-based photon-counting X-ray detector.

Characterization of charge sharing and fluorescence effects by multiple counts analysis in a Pixie-II based detection system

In CdTe X-ray photon counting detectors (PCD) featuring small pixel sizes (<100µm), charge-sharing and fluorescences emitted by Cd and Te are responsible for multiple counts from a single interacting photon. These effects can impair the imaging and spectroscopic performance of PCDs. Multiple counts can be partially or totally discriminated by properly setting the energy threshold implemented by the PCD system. Using monochromatic radiation and the Pixirad-1/Pixie-II CdTe PCD, this work characterizes and quantifies clusters of multiple counts as a function of energy and threshold.

Experimental optimization of single-exposure dual-energy angiography with photon-counting x-ray detectors.

Photon-counting x-ray detectors may enable single-exposure dual-energy (DE) x-rayangiography. The purpose of this paper is to experimentally optimize the energy thresholds and tube voltage for single-exposure DE x-rayangiography. We optimized single-exposure DE x-ray angiography using the iodine signal-difference-to-noise ratio (SDNR) per root patient air kerma (κ) as a figureof merit. We measured the iodine SDNR by imaging an iodine stepwedge immersed in a water tank with a depth of 30cm in the direction of x-ray propagation. The stepwedge was imaged using tube voltages ranging from 90 to 150kV and a cadmium telluride (CdTe) x-ray detector with two energy bins and analog charge summing for charge sharing suppression. The energy threshold that separates the two energy bins was varied from approximately 35keV to approximately 75% of the maximum energy of the x-ray beam. Curve fitting was used to determine the threshold that maximized . The effect of scatter was determined from measurements of the scatter-to-primary ratios (SPRs) of the low-energy and high-energy images and a semi-empirical model of the relationship between SDNR and SPR. Using the optimal parameters, we imaged a phantom with vessel-simulating structures and backgroundclutter. The optimal energy thresholds increased monotonically from ∼50 to ∼85keV over the range of tube voltages considered. For tube voltages greater than 90kV, the optimal energy thresholds consistently allocated approximately two thirds of all detected primary photons to the low energy bin; this ratio was preserved without scatter. Consistent with prior modeling studies, increased monotonically with tube voltage from 90 to 150kV; at 150kV was approximately 38% higher than that at 90kV for an iodine area density of ∼50mg/cm2 . Scatter reduced SDNR by approximately 25% for SPRs of ∼1 and 0.4 in low-energy and high-energy images, respectively. Achieving optimal image quality in single-exposure DE angiography with photon-counting x-ray detectors will require high tube voltages (i.e.,>130kV) and, for thick patients, energy thresholds that allocate approximately two thirds of all primary photons to the low-energy image. Future work will compare the image quality of singe-exposure photon-counting and kV-switching approaches to DE x-rayangiography.

Development of 64-channel LSI with ultrafast analog and digital signal processing dedicated for photon-counting computed tomography with multi-pixel photon counter

X-ray computed tomography (CT) is non destructive visualization inside the human body. Recently, X-ray photon-counting CT (PC-CT) has drawn attention for reducing the high doses used for patients and acquiring spectral information to identify materials. To enable photon-counting imaging with a wide region (∼60 mm long) of a small animal, we developed a photon-counting system using a 64-channel multi-pixel photon counter (MPPC) array with a fast temporal response (a few nanoseconds) and a very large signal gain (∼106), combined with a 64-channel yttrium–gadolinium–aluminum–gallium garnet scintillator array. In particular, to realize ultra-fast analog and digital signal processing (>10 MHz/channel), we developed a 64-channel large-scale integrated circuit (LSI) named MPPC-CT64. We had previously developed a 16-channel PC-CT system with a 16-channel LSI (MPPC-CT16). Although the MPPC-CT16 realized photon-counting imaging for a ∼16-mm-long phantom, there were some energy uncertainties in the LSI, which degraded the obtained CT image quality. The MPPC-CT64 implements a function for correcting the threshold energies and also increases the number of energy thresholds from four to six, which provides more precise measurements of CT values dependent on X-ray energy. In this paper, we briefly present the electric architecture and performance evaluation of the LSI combined with MPPCs.

Development and validation of a simulation tool for K-edge Subtraction imaging with polychromatic spectra and X-ray photon counting detectors

In K-Edge Subtraction (KES) imaging the K-edge of a high-Z contrast medium is used to generate a differential image in which, in principle, only the contrast agent is visible, suppressing any background. A possible approach to KES imaging is by using polychromatic spectra and spectral X-ray photon counting detectors (XPCD). In this paper we present the development and validation of a simulation tool for KES imaging with spectral XPCD. In addition an analytical formula, based on theoretical considerations, is derived, that further supports the simulation results. The validation is made on experimental data obtained with an X-ray tube, an iodine-based contrast agent, and Pixirad-1/Pixie-III detection system.

Precise phase retrieval for propagation-based images using discrete mathematics

The ill-posed problem of phase retrieval in optics, using one or more intensity measurements, has a multitude of applications using electromagnetic or matter waves. Many phase retrieval algorithms are computed on pixel arrays using discrete Fourier transforms due to their high computational efficiency. However, the mathematics underpinning these algorithms is typically formulated using continuous mathematics, which can result in a loss of spatial resolution in the reconstructed images. Herein we investigate how phase retrieval algorithms for propagation-based phase-contrast X-ray imaging can be rederived using discrete mathematics and result in more precise retrieval for single- and multi-material objects and for spectral image decomposition. We validate this theory through experimental measurements of spatial resolution using computed tomography (CT) reconstructions of plastic phantoms and biological tissues, using detectors with a range of imaging system point spread functions (PSFs). We demonstrate that if the PSF substantially suppresses high spatial frequencies, the potential improvement from utilising the discrete derivation is limited. However, with detectors characterised by a single pixel PSF (e.g. direct, photon-counting X-ray detectors), a significant improvement in spatial resolution can be obtained, demonstrated here at up to 17%.

Perovskite CsPbBr3 Single Crystal Detector for High Flux X-Ray Photon Counting.

X-ray photon-counting detectors capable of resolving the energies of single X-ray photons are critical in medical imaging, and a high count rate is essential for photon-counting detectors. Here, we report the performance of the perovskite CsPbBr3 single-crystal semiconductor detector for X-ray photon counting. The CsPbBr3 detector noise floor, energy response linearity, energy resolution, count rate, and polarization were evaluated. By fine-tuning the detector working conditions, our CsPbBr3 detector with a planar electrode can count an incident X-ray photon rate of ~0.099 and ~0.336 Mcps/pixel at 10% and 30% deadtime loss, respectively, with corresponding energy resolutions of ~18% at 59.5 keV and ~12% at 122 keV at same pulse processing conditions. We also demonstrated that our CsPbBr3 detectors show negligible polarization under an X-ray flux of ~0.45 M photons/s/mm2 for the typical timescale of multiple clinical X-ray scans, such as 1 s- 100 s. Our evaluation demonstrates the high potential of CsPbBr3 detectors for X-ray photon-counting applications in medical and industrial diagnostics with a lower cost than the current state-of-the-art.

Experimental characterization of an X-ray photon counting detector

This paper presents the experimental test results of the X-ray photon counting detectors based on our pixel readout ASIC. The chip integrates a 64 × 128 array of pixels in the size of 150 µm × 150 µm and each pixel consists of four 12-bit energy windows. Two types of CdTe detectors have been bump bonded to the ASICs and have been tested.The electronic characteristics of the readout chip were evaluated first by injecting charge through the calibration capacitors. The non-uniformity of the energy thresholds among pixels was measured and then were compensated by tuning the local DACs.The detector performance was then characterized using X-ray tubes. The energy responses of the detectors were measured using the characteristic X-ray of the target materials. The global energy thresholds were then calibrated with these specific energies of photons. Polarization effects under different photon fluxes were also studied. Preliminary imaging was conducted. The detector response under uniform irradiation was investigated. The spatial uniformity of the detectors was analyzed and the flat field correction was then conducted to improve image quality. CT detailed test results will be discussed in this paper.

Estimating optical parameters of biological tissues with photon-counting micro-CT.

Wavelength-dependent absorption and scattering properties determine the fluorescence photon transport in biological tissues and image resolution of optical molecular tomography. Currently, these parameters are computed from optically measured data. For small animal imaging, estimation of optical parameters is a large-scale optimization problem, which is highly ill-posed. In this paper, we propose a new, to the best of our knowledge, approach to estimate optical parameters of biological tissues with photon-counting micro-computed tomography (micro-CT). From photon-counting x-ray data, multi-energy micro-CT images can be reconstructed to perform multi-organ segmentation and material decomposition in terms of tissue constituents. The concentration and characteristics of major tissue constituents can be utilized to calculate the optical absorption and scattering coefficients of the involved tissues. In our study, we perform numerical simulation, phantom experiments, and in vivo animal studies to calculate the optical parameters using our proposed approach. The results show that our approach can estimate optical parameters of tissues with a relative error of <10%, accurately mapping the optical parameter distributions in a small animal.

Synchrotron X-ray Radiation (SXR) in Medical Imaging: Current Status and Future Prospects

Synchrotron X-ray radiation (SXR) has been widely studied to explore the structure of matter. Recently, there has been an intense focus on the medical application of SXR in imaging. This review is intended to explore the latest applications of SXR in medical imaging and to shed light on the advantages and drawbacks of this modality. The article highlights the latest developments in other fields that can greatly enhance the capability and applicability of SXR. The potentials of using machine and deep learning (DL)-based methods to generate synthetic images to use in regular clinics along with the use of photon counting X-ray detectors for spectral medical imaging with SXR are also discussed.

IDeF-X HDBD: Low-Noise ASIC for Imaging Spectroscopy With Semiconductor Detectors in Space Science Applications

We present the IDeF-X HDBD integrated circuit, our latest 32-channel application specified integrated circuit devoted to semiconductor photon counting X-ray detectors, designed to read charges ranging from −40 to 40 fC. The chip reaches an equivalent noise charge (ENC) floor of 17 electrons rms and is optimized for low input capacitance (< 10 pF). Each channel is based on an optimized charge sensitive amplifier (CSA) followed by a CR-RC <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> filter and a peak detector with a power consumption of <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"> <tex-math notation="LaTeX">$850 ~\mu \text{W}$ </tex-math></inline-formula> /channel. Gain and shaping times are tunable. This circuit has been designed with radiation mitigation techniques to meet space application requirements. This ASIC can read either cadmium telluride or silicon detectors for imaging-spectroscopy applications. An energy resolution of 230-eV FHWM at 6 keV and a low-level threshold of 1.2 keV were demonstrated with a single pixel silicon drift detector connected to one channel.

Photon-counting CT via interleaved/gapped spectral channels: Feasibility and imaging performance.

Compared to energy integration, photon-counting x-ray detection facilitates the spectral channelization (energy binning) in spectral CT and thus offers the opportunity to implement data acquisition via sophisticated schemes, for example, gapping or interleaving in spectral channels. In this article, we report our investigation of the performance of material decomposition based spectral imaging in photon-counting CT implemented in such data acquisition schemes, and their comparison with the benchmark scheme and other schemes without spectral gapping or interleaving. Using a deliberately designed anthropomorphic head phantom that mimics the intracranial soft tissues and bony structures, a simulation study is carried out with the focus on two-material decomposition based spectral imaging in photon-counting computed tomography (CT), under both ideal and realistic detector spectral responses. The projection data are acquired in four spectral channels, and then are sorted to implement the schemes of gapping ((ch1 , ch3 ); (ch2 , ch4 ); (ch1 , ch4 )) and interleaving ((ch1 , ch3 ) + (ch2 , ch4 ); (ch1 , ch4 ) + (ch2 , ch3 ); ((ch1 + ch3 ), (ch2 + ch4 )); ((ch1 + ch4 ), (ch2 + ch3 ))) in spectral channels, in addition to the benchmark scheme ((ch1 + ch2 ), (ch3 + ch4 )) and other conventional schemes (ch1 , ch2 ), (ch2 , ch3 ) and (ch3 , ch4 ), where ''ch'' denotes channel, ''+'' denote addition, and (·,·) the operation of material decomposition and image reconstruction. Using the contrast-to-noise ratio between targeted regions of interest as the figure of merit, we study the performance of spectral imaging (material specific and virtual monochromatic) associated with these spectral channelization schemes. Under ideal detector spectral response, the scheme (ch1 , ch4 ) outperforms the benchmark scheme ((ch1 + ch2 ), (ch3 + ch4 )) and others in gapped and/or interleaved spectral channelization in material specific imaging, while the interleaved scheme (ch1 , ch4 ) + (ch2 , ch3 ) performs the best in virtual monochromatic imaging. Notably, only about half x-ray dose is utilized in the scheme (ch1 , ch4 ) for image formation. Under realistic detector spectral response, the difference in imaging performance over all spectral channelization schemes diminishes, along with degradation in each scheme's individual performance. The results suggest that (i) different strategy in spectral channelization may have to be exercised in material specific imaging and virtual monochromatic imaging, respectively, and (ii) the spectral distortion in realistic detector's response due to charge-sharing, Compton scattering, and fluorescent escaping should be mitigated as much as possible. The spectral channelization schemes and associated imaging performance reported herein are novel and thus informative to the community, which may further the understanding of physical fundamentals and design principles for material decomposition based spectral imaging in photon-counting CT and other x-ray related imaging modalities.