1
|
Hsieh SS, Taguchi K. Spectral information content of Compton scattering events in silicon photon counting detectors. Med Phys 2024; 51:2386-2397. [PMID: 38353409 DOI: 10.1002/mp.16990] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/18/2023] [Revised: 09/14/2023] [Accepted: 01/22/2024] [Indexed: 04/05/2024] Open
Abstract
BACKGROUND Silicon (Si) is a possible sensor material for photon counting detectors (PCDs). A major drawback of Si is that roughly two-thirds of x-ray interactions in the diagnostic energy range are Compton scattering. Because Compton scattering is an energy-insensitive process, it is commonly assumed that Compton events retain little spectral information. PURPOSE To quantify how much information can be recovered from Compton scattering events in models of Si PCDs. METHODS We built a simplified model of Si interactions including two interaction mechanisms: photoelectric effect and Compton scattering. We considered three different binning options that represent strategies for handling Compton events: in Compton censoring, all events under 38 keV (the maximum energy possible from Compton scattering for a 120 keV incident photon) were discarded; in Compton counting, all events between 1 and 38 keV were placed into a single bin; in Compton binning, all events were placed into energy bins of uniform width. These were compared to the ideal detector, which always recorded the correct energy (i.e., 100% photoelectric effect). Every photon was assumed to interact once and only once with Si, and the energy bin width was 5 keV. In the primary analysis, the Si detector was irradiated with a 120 kV spectrum filtered by 30 cm of water, with 99.5% of the arriving spectrum above 38 keV so that there was good separation between photoelectric effect and Compton scattering, and the figures of merit were the Cramér-Rao lower bound (CRLB) of the variance of iodine and water basis material decomposition images, as well as the CRLB of virtual monoenergetic images (i.e., linear combinations of material images) that maximize iodine CNR or water CNR. We also constructed a local linear estimator that attains the CRLB. In secondary analyses, we applied other sources of spectral distortion: (1) a nonzero minimum energy threshold; (2) coarser, 10 keV energy bins; and (3) a model of charge sharing. RESULTS With our chosen spectrum, 67% of the interactions were Compton scattering. Consistent with this, the material decomposition variance for the Compton censoring model, averaged over both basis materials, was 258% greater than the ideal detector. If Compton events carried no spectral information, the Compton counting model would show similar variance. Instead, its basis material variance was 103% greater than the ideal detector, implying that Compton counts indeed carry significant spectral information. The Compton binning model had a basis material variance 60% greater than the ideal detector. The Compton binning model was not affected by a 5 keV minimum energy threshold, but the variance increased from 60% to 107% when charge sharing was included and to 78% with coarser energy bins. For optimized CNR images, the average variance was 149%, 12%, and 10% higher than the ideal detector for the Compton censoring, counting, and binning models, reinforcing the hypothesis that Compton counts are useful for detection tasks and that precise energy assignments are not necessary. CONCLUSIONS Substantial spectral information remains after Compton scattering events in silicon PCDs.
Collapse
Affiliation(s)
- Scott S Hsieh
- Department of Radiology, Mayo Clinic, Rochester, Minnesota, USA
| | - Katsuyuki Taguchi
- Departments of Radiology and Bioengineering, Johns Hopkins University, Baltimore, Maryland, USA
| |
Collapse
|
2
|
Hsieh SS. Possible improvements in effective fill factor using X-ray fluorescent interpixel reflectors. Med Phys 2024; 51:1617-1625. [PMID: 38259109 DOI: 10.1002/mp.16944] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Revised: 12/27/2023] [Accepted: 01/02/2024] [Indexed: 01/24/2024] Open
Abstract
BACKGROUND The spatial resolution of energy-integrating diagnostic CT scanners is limited by interpixel reflectors on the detector, which optically isolate pixels but create dead space. Because the width of the reflector cannot easily be decreased, fill factor diminishes as resolution increases. PURPOSE We propose loading (or mixing) a high-Z element into the reflectors, causing the reflectors to be X-ray fluorescent. Re-emitted characteristic X-rays could be detected in adjacent pixels, increasing the effective fill factor and compensating for fill factor loss with higher-resolution detectors. The purpose of this work is to understand the physical principles of this approach and to analyze its effectiveness using Monte Carlo simulations. METHODS Detector pixels were modeled using the GEANT4 Monte Carlo package. The width of the reflector was kept constant at 0.1 mm throughout, and we considered pixel pitches between 0.5 and 1 mm. The pixelated scintillator material was gadolinium oxysulfide, 3 mm thick. The baseline reflector material was chosen to be acrylic, and varying concentrations of a high-Z element were loaded into the material. We assumed that the optical characteristics of pixels were ideal (no absorption within pixels, perfect reflection at boundaries). The detector was irradiated uniformly with 10,000 X-ray photons to estimate its spectral response. The figure of merit was the variance of the detector signal at zero frequency normalized to that of an ideal single-bin photon-counting detector with 100% fill factor. Sensitivity analyses were conducted to understand the effect of varying the high-Z element concentration and the spectrum. RESULTS Initial simulations suggested that a k-edge near 50 keV would be ideal. Gd was therefore selected as the high-Z material. The relative variances for a conventional energy integrating detector without Gd at 1 mm pixel pitch (81% fill factor) and 0.5 mm pixel pitch (64% fill factor) were 1.38 and 1.74, compared to 1.00 for an ideal photon counting detector, implying a 26% variance penalty for 0.5 mm pitch. When 1 g/cm3 Gd was loaded into the interpixel reflector, the relative variance improved to 1.27 and 1.43, respectively, implying that the variance penalty for including Gd together with 0.5 mm pitch is only 4%. Performance was nearly maximized at 1.0 g/cm3 of Gd, but a concentration of 0.5 g/cm3 of Gd showed most of the benefit. Improvements depend weakly on kV, with lower kV associated with higher improvements. An external anti-scatter grid was not modeled in our simulations and would reduce the expected benefit, depending greatly on the pitch and dimensionality of the anti-scatter grid. CONCLUSIONS The losses in fill factor associated with smaller pixel pitch can be reduced if Gd or a similar element could be loaded into the interpixel reflector. These improvements in noise efficiency are yet to be verified experimentally.
Collapse
Affiliation(s)
- Scott S Hsieh
- Department of Radiology, Mayo Clinic, Rochester, Minnesota, USA
| |
Collapse
|
3
|
Richtsmeier D, Rodesch PA, Iniewski K, Bazalova-Carter M. Material decomposition with a prototype photon-counting detector CT system: expanding a stoichiometric dual-energy CT method via energy bin optimization and K-edge imaging. Phys Med Biol 2024; 69:055001. [PMID: 38306974 DOI: 10.1088/1361-6560/ad25c8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2023] [Accepted: 02/01/2024] [Indexed: 02/04/2024]
Abstract
Objective.Computed tomography (CT) has advanced since its inception, with breakthroughs such as dual-energy CT (DECT), which extracts additional information by acquiring two sets of data at different energies. As high-flux photon-counting detectors (PCDs) become available, PCD-CT is also becoming a reality. PCD-CT can acquire multi-energy data sets in a single scan by spectrally binning the incident x-ray beam. With this, K-edge imaging becomes possible, allowing high atomic number (high-Z) contrast materials to be distinguished and quantified. In this study, we demonstrated that DECT methods can be converted to PCD-CT systems by extending the method of Bourqueet al(2014). We optimized the energy bins of the PCD for this purpose and expanded the capabilities by employing K-edge subtraction imaging to separate a high-atomic number contrast material.Approach.The method decomposes materials into their effective atomic number (Zeff) and electron density relative to water (ρe). The model was calibrated and evaluated using tissue-equivalent materials from the RMI Gammex electron density phantom with knownρevalues and elemental compositions. TheoreticalZeffvalues were found for the appropriate energy ranges using the elemental composition of the materials.Zeffvaried slightly with energy but was considered a systematic error. Anex vivobovine tissue sample was decomposed to evaluate the model further and was injected with gold chloride to demonstrate the separation of a K-edge contrast agent.Main results.The mean root mean squared percent errors on the extractedZeffandρefor PCD-CT were 0.76% and 0.72%, respectively and 1.77% and 1.98% for DECT. The tissue types in theex vivobovine tissue sample were also correctly identified after decomposition. Additionally, gold chloride was separated from theex vivotissue sample with K-edge imaging.Significance.PCD-CT offers the ability to employ DECT material decomposition methods, along with providing additional capabilities such as K-edge imaging.
Collapse
Affiliation(s)
- Devon Richtsmeier
- Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Victoria, British Columbia V8P 5C2, Canada
| | - Pierre-Antoine Rodesch
- Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Victoria, British Columbia V8P 5C2, Canada
| | - Kris Iniewski
- Redlen Techologies, 1763 Sean Heights, Saanichton, British Columbia V8M 1X6, Canada
| | - Magdalena Bazalova-Carter
- Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Victoria, British Columbia V8P 5C2, Canada
| |
Collapse
|
4
|
Taguchi K, Hsieh SS. Direct energy binning for photon counting detectors: Simulation study. Med Phys 2024; 51:70-79. [PMID: 38011545 PMCID: PMC10842195 DOI: 10.1002/mp.16841] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 08/14/2023] [Accepted: 10/01/2023] [Indexed: 11/29/2023] Open
Abstract
BACKGROUND Photon counting detectors (PCDs) for x-ray computed tomography (CT) face spectral distortion from pulse pileup and charge sharing. The photon counting scheme used by many PCDs is threshold-subtract (TS) with pulse height analysis (PHA), where each counter counts up-crossing events when pulses exceed an energy threshold. PCD data are not Poisson-distributed due to charge sharing and pulse pileup, but the counting statistics have never been studied yet. PURPOSE The objectives of this study were (1) to propose a modified photon counting scheme, direct energy binning (DB), that is expected to be robust against pulse pileup; (2) to assess the performance of DB compared to TS; and (3) to evaluate its counting statistics. METHODS With DB scheme, counter k starts a timer upon an up-crossing event of energy threshold k, and adds a count only if the next higher energy threshold (k+1) was not crossed within a short time window (hence, the pulse peak belongs to the energy bin k). We used Monte Carlo (MC) simulation and assessed count-rate curves and count-rate-dependent spectral imaging task performance for conventional CT imaging as well as water thickness estimation, water-bone material decomposition, and K-edge imaging with tungsten as the K-edge material. We also assessed count-rate-dependent measurement statistics such as expectation, variance, and covariance of total counts as well as energy bin outputs. The agreement with counting statistics models was also evaluated. RESULTS The DB scheme improved the count-rate curve, that is, mean measured counts as a function of input count-rate, and peaked with 59% higher count-rate capability than the TS scheme (3.5 × 108 counts per second (cps)/mm2 versus 2.3 × 108 cps/mm2 ). The Cramér-Rao lower bounds (CRLB) of the variance of basis line integrals estimation for DB was better than those for TS by 2% for the conventional CT imaging, 30% for water-bone material decomposition, and 32% for K-edge imaging at 1000 mA (at 7.3 × 107 cps/sub-pixel after charge sharing). When count-rates were lower, PCD data statistics were dominated by charge sharing: the variance of total counts and lower energy bins was larger than the mean counts; the covariance of bin data was positive and non-zero. When count-rates were higher, PCD data statistics were dominated by pulse pileup: the variance of data was lower than the mean; the covariance of bin data was negative. The transition between the two regimes occurred smoothly, and pulse pileup dominated the statistics ≥400 mA (when the count-rate after charge sharing was 2.9 × 107 cps/sub-pixel and the probability of count-loss for DB was 37%). Both DB and TS had good agreement with Yu-Fessler's models of total counts; however, DB had a better agreement with Wang's variance and covariance models for energy bin data than TS did. CONCLUSIONS The proposed DB scheme had several advantages over TS. At low to moderate flux, DB could improve the resilience of PCDs to pulse pileup. Counting statistics deviated from the Poisson distribution due to charge sharing for lower count-rate conditions and pulse pileup for higher count-rate conditions.
Collapse
Affiliation(s)
- Katsuyuki Taguchi
- The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
| | - Scott S Hsieh
- Department of Radiology, Mayo Clinic, Rochester, Minnesota, USA
| |
Collapse
|
5
|
Lahoud E, Zaffrani Y, Yagil Y. Photon counting x-ray detectors as scatter probes. Med Phys 2024; 51:93-102. [PMID: 38043090 DOI: 10.1002/mp.16871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Revised: 10/16/2023] [Accepted: 11/17/2023] [Indexed: 12/05/2023] Open
Abstract
BACKGROUND Direct conversion x-ray Photon Counting Detectors (PCD) are posed to play a vital role in future medical imaging devices such as Computed Tomography (CT) scanners. PCD are expected to improve current CT technology on several fronts, such as resolution, dose utilization, and spectral performance. However, they are not readily expected to improve the handling of object scatter, one of the major sources of image artifacts in CT technology. PURPOSE We explore a potential method for obtaining in-situ object scatter estimation using the same PCD array used in the x-ray imaging system, such as in computed tomography. This unexpected benefit of using PCD has the potential to improve the image quality by providing better input into the scatter estimation and correction algorithms used in image reconstruction. METHODS In CT scanners the primary method for rejecting scatter signal originating from the scanned object relies on placing anti-scatter grids (ASG) close to the detector plane. This remains the case when transitioning to using PCD instead of energy integration detectors in CT. However, the combination of PCD and ASG opens a possibility to use some of the unique properties of PCD, namely, very low noise and coincidence counters to obtain, in addition to the attenuation data, a simultaneous and instantaneous estimate of the scatter signal reaching every detector element. When a small air gap is introduced between the ASG and the detector surface, the scatter radiation with large angular distribution has a greater probability of producing charge sharing events that can be detected by a coincidence counter. In this work we demonstrate the feasibility of such an approach in a tabletop experiment using PCD detector that lacks coincidence counting capability, instead we use the spectral signature of split charge events as proxy to coincidence counting. For this purpose, we first demonstrate the spectral impact of ASG misalignment using the same experimental setup. In addition, we perform a separate tabletop scattering experiment from a narrow column of water that demonstrates another potential use of the low noise capabilities of PCDs. RESULTS We measured and quantified the high sensitivity of the spectral response to ASG alignment on the PCD detector pixel array, we found that the probability of energy misregistration of 60 keV photons can increase by up to a factor of 3 when the ASG is poorly aligned. We then leveraged these results to obtain an estimate on the expected increase in coincidence counts for a wide range of scatter-to-primary (SPR) ratio and find a good match with expectations from a geometric modeling of the system, where the expected increase in coincidences was of the order of the SPR. Finally, the low noise detector also allowed us to measure the real space scatter signal associated with the coherent molecular form factor of water, revealing the ring-shaped scatter signal with an energy dependent distribution that was well captured by calculation. CONCLUSIONS The advent of PCD detectors and their imminent use in commercial CT scanners opens new and exciting possibilities for utilizing PCD detectors in unexpected ways. In this proof-of-concept study, we showed how charge sharing, a spectral information degrading effect, can instead be used to obtain in-situ scatter estimation. We also demonstrated the PCD ability to perform extremely sensitive measurements using affordable benchtop setup for investigations normally reserved for synchrotron facilities.
Collapse
|
6
|
章 浩, 李 树, 刘 颖, 路 鹤. [A comprehensive review on photon-counting computed tomography: Principles, technical hurdles and analysis of clinical applications]. SHENG WU YI XUE GONG CHENG XUE ZA ZHI = JOURNAL OF BIOMEDICAL ENGINEERING = SHENGWU YIXUE GONGCHENGXUE ZAZHI 2023; 40:1012-1018. [PMID: 37879932 PMCID: PMC10600420 DOI: 10.7507/1001-5515.202305015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Revised: 08/30/2023] [Indexed: 10/27/2023]
Abstract
In recent years, photon-counting computed tomography (PCD-CT) based on photon-counting detectors (PCDs) has become increasingly utilized in clinical practice. Compared with conventional CT, PCD-CT has the potential to achieve micron-level spatial resolution, lower radiation dose, negligible electronic noise, multi-energy imaging, and material identification, etc. This advancement facilitates the promotion of ultra-low dose scans in clinical scenarios, potentially detecting minimal and hidden lesions, thus significantly improving image quality. However, the current state of the art is limited and issues such as charge sharing, pulse pileup, K-escape and count rate drift remain unresolved. These issues could lead to a decrease in image resolution and energy resolution, while an increasing in image noise and ring artifact and so on. This article systematically reviewed the physical principles of PCD-CT, and outlined the structural differences between PCDs and energy integration detectors (EIDs), and the current challenges in the development of PCD-CT. In addition, the advantages and disadvantages of three detector materials were analysed. Then, the clinical benefits of PCD-CT were presented through the clinical application of PCD-CT in the three diseases with the highest mortality rate in China (cardiovascular disease, tumour and respiratory disease). The overall aim of the article is to comprehensively assist medical professionals in understanding the technological innovations and current technical limitations of PCD-CT, while highlighting the urgent problems that PCD-CT needs to address in the coming years.
Collapse
Affiliation(s)
- 浩伟 章
- 上海理工大学 健康科学与工程学院(上海 200093)School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, P. R. China
| | - 树晗 李
- 上海理工大学 健康科学与工程学院(上海 200093)School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, P. R. China
| | - 颖 刘
- 上海理工大学 健康科学与工程学院(上海 200093)School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, P. R. China
| | - 鹤晴 路
- 上海理工大学 健康科学与工程学院(上海 200093)School of Health Science and Engineering, University of Shanghai for Science and Technology, Shanghai 200093, P. R. China
| |
Collapse
|
7
|
Yang M, Wohlfahrt P, Shen C, Bouchard H. Dual- and multi-energy CT for particle stopping-power estimation: current state, challenges and potential. Phys Med Biol 2023; 68. [PMID: 36595276 DOI: 10.1088/1361-6560/acabfa] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2022] [Accepted: 12/15/2022] [Indexed: 12/23/2022]
Abstract
Range uncertainty has been a key factor preventing particle radiotherapy from reaching its full physical potential. One of the main contributing sources is the uncertainty in estimating particle stopping power (ρs) within patients. Currently, theρsdistribution in a patient is derived from a single-energy CT (SECT) scan acquired for treatment planning by converting CT number expressed in Hounsfield units (HU) of each voxel toρsusing a Hounsfield look-up table (HLUT), also known as the CT calibration curve. HU andρsshare a linear relationship with electron density but differ in their additional dependence on elemental composition through different physical properties, i.e. effective atomic number and mean excitation energy, respectively. Because of that, the HLUT approach is particularly sensitive to differences in elemental composition between real human tissues and tissue surrogates as well as tissue variations within and among individual patients. The use of dual-energy CT (DECT) forρsprediction has been shown to be effective in reducing the uncertainty inρsestimation compared to SECT. The acquisition of CT data over different x-ray spectra yields additional information on the material elemental composition. Recently, multi-energy CT (MECT) has been explored to deduct material-specific information with higher dimensionality, which has the potential to further improve the accuracy ofρsestimation. Even though various DECT and MECT methods have been proposed and evaluated over the years, these approaches are still only scarcely implemented in routine clinical practice. In this topical review, we aim at accelerating this translation process by providing: (1) a comprehensive review of the existing DECT/MECT methods forρsestimation with their respective strengths and weaknesses; (2) a general review of uncertainties associated with DECT/MECT methods; (3) a general review of different aspects related to clinical implementation of DECT/MECT methods; (4) other potential advanced DECT/MECT applications beyondρsestimation.
Collapse
Affiliation(s)
- Ming Yang
- The University of Texas MD Anderson Cancer Center, Department of Radiation Physics, 1515 Holcombe Blvd Houston, TX 77030, United States of America
| | - Patrick Wohlfahrt
- Massachusetts General Hospital and Harvard Medical School, Department of Radiation Oncology, Boston, MA 02115, United States of America
| | - Chenyang Shen
- University of Texas Southwestern Medical Center, Department of Radiation Oncology, 2280 Inwood Rd Dallas, TX 75235, United States of America
| | - Hugo Bouchard
- Département de physique, Université de Montréal, Complexe des sciences, 1375 Avenue Thérèse-Lavoie-Roux, Montréal, Québec H2V0B3, Canada.,Centre de recherche du Centre hospitalier de l'Université de Montréal, 900 Rue Saint-Denis, Montréal, Québec, H2X 0A9, Canada.,Département de radio-oncologie, Centre hospitalier de l'Université de Montréal, 1051 Rue Sanguinet, Montréal, Québec H2X 3E4, Canada
| |
Collapse
|
8
|
Wu X, Shields A, Nagesh SVS, Bednarek DR, Rudin S. Comparison of quantitative imaging characteristics between a new, larger-FOV 1000 fps high-speed angiographic (HSA) photon-counting detector (Aries) with a smaller HSA detector (Actaeon). PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2022; 12031:120310K. [PMID: 35982768 PMCID: PMC9385188 DOI: 10.1117/12.2611538] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
High Speed Angiography (HSA) requires imaging detectors with both high-temporal and high-spatial resolution. Both the Aries and Acteon detectors by Direct Conversion (Stockholm, Sweden) are CdTe direct photon-counting detectors (PCD) that have acquisition frame rates of up to 1000-fps and a 100-micrometer pixel pitch; however, the new Aries detector offers a larger field of view (512 × 768 pixels) compared to the smaller Actaeon detector (256 × 256 pixels). An expanded field of view is required for imaging of larger vasculature, thus the Aries offers this advantage. Evaluations were performed of both detectors under Anti-Coincidence Circuitry (ACC-ON) mode, which corrects for charge sharing between pixels. Initial evaluations of instrumentation noise and detector energy-threshold calibration using Am-241 gamma spectroscopy were performed for the new Aries detector. Linearity was also evaluated for the Aries for each of the 12 individual modules that compose the detector field to check for homogeneity in response to exposure throughout the detector. Finally, Normalized Noise Power Spectrum (NNPS), Modulation Transfer Function (MTF) and Detective Quantum Efficiency (DQE) were then compared between the Aries and Actaeon detectors at two different exposures and detector energy thresholds. The detectors are linear up to approximately 1000 μR and have no instrumentation noise above a threshold of 15 keV. As expected, the MTF's and DQE's are similar between the Aries and Actaeon detectors, and there are thus no tradeoff's in image quality to achieve the larger FOV.
Collapse
Affiliation(s)
- X Wu
- Canon Stroke and Vascular Research Center, University at Buffalo, Buffalo, NY
| | - A Shields
- Canon Stroke and Vascular Research Center, University at Buffalo, Buffalo, NY
| | - S V Setlur Nagesh
- Canon Stroke and Vascular Research Center, University at Buffalo, Buffalo, NY
| | - D R Bednarek
- Canon Stroke and Vascular Research Center, University at Buffalo, Buffalo, NY
| | - S Rudin
- Canon Stroke and Vascular Research Center, University at Buffalo, Buffalo, NY
| |
Collapse
|
9
|
Taguchi K, Iwanczyk JS. Assessment of multi-energy inter-pixel coincidence counters for photon-counting detectors at the presence of charge sharing and pulse pileup: A simulation study. Med Phys 2021; 48:4909-4925. [PMID: 34287966 PMCID: PMC9942613 DOI: 10.1002/mp.15112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2020] [Revised: 06/11/2021] [Accepted: 06/25/2021] [Indexed: 11/08/2022] Open
Abstract
PURPOSE Spectral distortion due to charge sharing (CS) and pulse pileup (PP) in photon-counting detectors (PCDs) degrades the quality of PCD data. We recently proposed multi-energy inter-pixel coincidence counters (MEICC) that provided spectral cross-talk information related to CS. When PP was absent, the normalized Cramér-Rao lower bounds (nCRLBs) of 225-µm pixel PCDs with MEICC was comparable to those of 450-µm pixel PCD without MEICC. The aim of this study was to assess the performance of PCDs with MEICC in the presence of both CS and PP using computer simulations. METHODS An in-house Monte Carlo program was modified to incorporate the following four temporal elements: (1) A pulse shape with a pulse duration of 20 ns, (2) delays of up to 10 ns in anode arrival times when photons were incident on pixel boundaries, (3) offsets proportional to a vertical separation between the primary and secondary charge clouds at the rate of ±4 ns per ±100 µm, and (4) a stochastic fluctuation of anode arrival times for all of the charge clouds with a standard deviation of 2 ns. We assessed the performance of five PCDs, (a)-(f), for three spectral tasks, (A)-(C): (a) The conventional PCD, (b) a PCD with MEICC, (c) a PCD with one coincidence counter (1CC), (d) a PCD with a 3 × 3 analog charge summing scheme (ACS), and (e) a PCD with a 3 × 3 digital count summing scheme (DCS); (A) conventional CT imaging with water (i.e., linear attenuation coefficient maps), (B) water-bone material decomposition, and (C) K-edge imaging with tungsten. The tube current was changed from 1 mA to 1000 mA and the nCRLB was assessed. RESULTS The recorded count rate curves were fitted by the non-paralyzable detection model with the effective deadtime parameter. The best fit was achieved by 25.8 ns for the conventional PCD, 18.6 ns for MEICC and 1CC, 140.5 ns for ACS, and 209.0 ns for DCS. The nCRLBs were strongly dependent on count rates. MEICC provided the best nCRLBs for all of the imaging tasks over the count rate range investigated except for a few conditions such as K-edge imaging at 1 mA. PP decreased the merit of MEICC over the conventional PCD in addressing CS. Nonetheless, MEICC consistently provided better nCRLBs than the conventional PCD did. The nCRLBs of MEICC were in the range of 49-58% of those of the conventional PCD for K-edge imaging, 45-76% for water-bone material decomposition, and 81-88% for the conventional CT imaging (i.e., linear attenuation coefficient maps). ACS provided better nCRLBs than the conventional PCD did only when the effect of PP was minor (e.g., when the counting efficiency of the conventional PCD was higher than 0.95 with the tube current of up to 100 mA). CONCLUSION Besides a few cases, MEICC provides the best nCRLBs for all of the tasks at all of the count rates. ACS and DCS provide better nCRLBs than the conventional PCD does only when count rates are very low.
Collapse
Affiliation(s)
- Katsuyuki Taguchi
- The Russell H. Morgan Department of Radiology and Radiological Science, Johns Hopkins University School of Medicine, Baltimore, Maryland 21287,Corresponding author.. 601 North Caroline Street, JHOC 4253, Baltimore, Maryland 21287, U.S.A., 443-287-2425 (telephone), 410-614-1060 (facsimile)
| | | |
Collapse
|
10
|
Ren L, Allmendinger T, Halaweish A, Schmidt B, Flohr T, McCollough CH, Yu L. Energy-integrating-detector multi-energy CT: Implementation and a phantom study. Med Phys 2021; 48:4857-4871. [PMID: 33988849 PMCID: PMC8455428 DOI: 10.1002/mp.14943] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 04/09/2021] [Accepted: 04/30/2021] [Indexed: 11/08/2022] Open
Abstract
PURPOSE Multi-energy computed tomography (MECT) has a great potential to enable many novel clinical applications such as simultaneous multi-contrast imaging. The purpose of this study was to implement triple-beam MECT on a traditional energy-integrating-detector (EID) CT platform (EID-MECT). METHODS This was accomplished by mounting a z-axis split-filter (0.05 mm Au, 0.6 mm Sn) on Tube A of a dual-source EID CT scanner. With the two split x-ray beams from Tube A and the third beam from Tube B, three beams with different x-ray spectra can be simultaneously acquired. With Tube B operated at 70 or 80 kV and Tube A at 120 or 140 kV, four different triple-beam configurations were calibrated for MECT measurements: 70/Au120/Sn120, 80/Au120/Sn120, 70/Au140/Sn140, and 80/Au140/Sn140 kV. Iodine (I), gadolinium (Gd), bismuth (Bi) samples, and their mixtures were prepared for 2 three-material-decomposition tasks and 1 four-material-decomposition task. For each task, samples were placed in a water phantom and scanned using each of the four triple-beam configurations. For comparison, the same phantom was also scanned using three other dual-energy CT (DECT) or MECT technologies: twin-beam DECT (TB-DECT), dual-source DECT (DS-DECT), and photon-counting-detector CT (PCD-CT), all with optimal x-ray spectrum settings and at equal volume CT dose index (CTDIvol). The phantom for four-material decomposition (I/Gd/Bi/Water imaging) was scanned using the PCD-CT only (140 kV with 25, 50, 75, and 90 keV). Image-based material decomposition was performed to acquire material-specific images, on which the mean basis material concentrations and noise levels were measured and compared across all triple-beam configurations in EID-MECT and various DECT/MECT systems. RESULTS The optimal triple-beam configuration was task-dependent with 70/Au120/Sn120, 70/Au140/Sn140, and 70/Au120/Sn120 kV for I/Gd/Water, I/Bi/Water, and I/Gd/Bi/Water material decomposition tasks, respectively. At equal radiation dose level, EID-MECT provided comparable or better quantification accuracy in material-specific images for all three material decomposition tasks, compared to EID-based DECT and PCD-CT systems. In terms of noise level comparison, EID-MECT-derived material-specific images showed lower noise levels than TB-DECT and DS-DECT, but slightly higher than that from PCD-CT in I/Gd/Water imaging. For I/Bi/Water imaging, EID-MECT showed a comparable noise level to DS-DECT, and a much lower noise level than TB-DECT and PCD-CT in all material-specific images. For the four-material decomposition task involving I/Gd/Bi/Water, the bismuth-specific image derived from EID-MECT was slightly noisier, but both iodine- and gadolinium-specific images showed much lower noise levels in comparison to PCD-CT. CONCLUSIONS For the first time, an EID-based MECT system that can simultaneously acquire three x-ray spectra measurements was implemented on a clinical scanner, which demonstrated comparable or better imaging performance than existing DECT and MECT systems.
Collapse
Affiliation(s)
- Liqiang Ren
- Department of Radiology, Mayo Clinic, Rochester, MN, US
| | | | | | | | | | | | - Lifeng Yu
- Department of Radiology, Mayo Clinic, Rochester, MN, US
| |
Collapse
|
11
|
Feng M, Ji X, Zhang R, Treb K, Dingle AM, Li K. An experimental method to correct low-frequency concentric artifacts in photon counting CT. Phys Med Biol 2021; 66. [PMID: 34315142 DOI: 10.1088/1361-6560/ac1833] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Accepted: 07/27/2021] [Indexed: 11/12/2022]
Abstract
Large-area photon counting detectors (PCDs) are usually built by tiling multiple semiconductor panels that often have slightly different spectral responses to input x-rays. As a result of this spectral inconsistency, experimental PCD-CT images of large, human-sized objects may show high-frequency ring artifacts and low-frequency band artifacts. Due to the much larger width of the bands compared with the rings, the concentric artifact problem in PCD-CT images of human-sized objects cannot be adequately addressed by conventional CT ring correction methods. This work presents an experimental method to correct the concentric artifacts in PCD-CT. The method is applicable to not only energy-discriminating PCDs with multiple bins but also PCDs with only a single threshold controller. Its principle is similar to the two-step beam hardening correction method, except that the proposed method uses pixel-specific polynomial functions to address the spectral inconsistency problem across the detector plane. The pixel-specific polynomial coefficients were experimentally calibrated using 15 acrylic sheets and 6 aluminum sheets of known thicknesses. The pixel-specific polynomial functions were used to convert the measured PCD-CT projection data to acrylic- and aluminum-equivalent thicknesses that are energy-independent. The proposed method was experimentally evaluated using a human cadaver head and multiple physical phantoms: two of them contain iodine and one phantom contains dual K-edge contrast materials (gadolinium and iodine). The results show that the proposed method can effectively remove the low-frequency concentric artifacts in PCD-CT images while reducing beam hardening artifacts. In contrast, the conventional CT ring correction algorithm did not adequately address the low-frequency band artifacts. Compared with the direct material decomposition-based correction method, the proposed method not only relaxes the requirement of multi-energy bins but also generates images with lower noise and fewer concentric artifacts.
Collapse
Affiliation(s)
- Mang Feng
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Xu Ji
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Ran Zhang
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Kevin Treb
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Aaron M Dingle
- Department of Surgery, University of Wisconsin-Madison, WI 53792, United States of America
| | - Ke Li
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America.,Department of Radiology, University of Wisconsin-Madison, 600 Highland Avenue, Madison, WI 53792, United States of America
| |
Collapse
|
12
|
Hsieh SS, Iniewski K. Improving digital charge sharing compensation in photon counting detectors with a low-threshold comparator. Med Phys 2021; 48:5819-5829. [PMID: 34370321 DOI: 10.1002/mp.15154] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2020] [Revised: 05/30/2021] [Accepted: 07/18/2021] [Indexed: 11/06/2022] Open
Abstract
PURPOSE Charge sharing is a major non-ideality in photon counting detectors (PCDs) and can increase variance in material decomposition images. Analog charge summing (ACS) is an effective mechanism for charge sharing compensation (CSC), but is complex to implement and may limit the maximum count rate of the PCD. Digital CSC mechanisms such as digital count summing (DCS) may be simpler to implement; however, earlier simulation studies suggest that digital CSC only provides half the benefit of ACS. We propose including an additional low-threshold comparator (LTC) underneath the noise floor of the PCD to improve the effectiveness of digital CSC. METHODS We simulated a PCD with four or eight equally spaced energy bins. X-ray photons arrived on the PCD following a Poisson distribution, and charge was allocated to PCD pixels following Monte Carlo techniques. Gaussian electronic noise was added with standard deviation of 2 keV and the signals were processed with four CSC schemes: no CSC, ACS, DCS, and DCS with LTC. The energy bins were placed from 25 to 100 keV at 25 keV intervals (for four bins) or from 25 to 112.5 keV at 12.5 keV intervals (for eight bins), and the LTC threshold was placed at 8 keV in both cases. The binned counts were transformed into estimates of water and iodine material thickness using a linear estimator that was fitted to the data. Our simulations were performed in the low-flux limit without any pileup, assuming a 120 kVp spectrum, 25 cm water object, and 0.3 mm PCD pixel size. RESULTS All CSC schemes decreased variance in basis material decomposition. In the four-bin PCD, the relative dose efficiencies (inverse of the variance) for iodine material decomposition were 1.0, 2.4, 3.2, and 4.3 for a PCD without CSC, DCS without LTC, DCS with LTC, and ACS, respectively. In the eight-bin PCD, the relative dose efficiencies were 1.1, 2.5, 3.1, and 4.8, respectively. In a sensitivity analysis, electronic noise had a stronger deleterious effect on ACS than DCS. In simulated visual images, DCS and ACS improved high frequency contrast in material decomposition images. CONCLUSIONS Introducing an LTC may reduce the performance differential between DCS and ACS. These findings have been derived from simulation studies only and have not been validated experimentally.
Collapse
Affiliation(s)
- Scott Sigao Hsieh
- Department of Radiology, Mayo Clinic Rochester, Rochester, Minnesota, USA
| | - Kris Iniewski
- Redlen Technologies, Saanichton, British Columbia, Canada
| |
Collapse
|
13
|
Wang AS, Pelc NJ. Spectral Photon Counting CT: Imaging Algorithms and Performance Assessment. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2021; 5:453-464. [PMID: 35419500 PMCID: PMC9000208 DOI: 10.1109/trpms.2020.3007380] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Adam S Wang
- Departments of Radiology and, by courtesy, Electrical Engineering, Stanford University, Stanford, CA 94305 USA
| | - Norbert J Pelc
- Department of Radiology, Stanford University, Stanford, CA 94305 USA
| |
Collapse
|
14
|
Ji X, Treb K, Li K. Anomalous edge response of cadmium telluride-based photon counting detectors jointly caused by high-flux radiation and inter-pixel communication. Phys Med Biol 2021; 66:10.1088/1361-6560/abf1fe. [PMID: 33765661 PMCID: PMC8086811 DOI: 10.1088/1361-6560/abf1fe] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Accepted: 03/25/2021] [Indexed: 12/20/2022]
Abstract
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.
Collapse
Affiliation(s)
- Xu Ji
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Kevin Treb
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
| | - Ke Li
- Department of Medical Physics, University of Wisconsin-Madison, 1111 Highland Avenue, Madison, WI 53705, United States of America
- Department of Radiology, University of Wisconsin-Madison, 600 Highland Avenue, Madison, WI 53792, United States of America
| |
Collapse
|
15
|
Affiliation(s)
- Amir Pourmorteza
- From the Department of Radiology and Imaging Sciences, Emory University, 1364 Clifton Rd NE, Atlanta, GA 30322
| |
Collapse
|
16
|
Pickford Scienti OLP, Bamber JC, Darambara DG. CdTe Based Energy Resolving, X-ray Photon Counting Detector Performance Assessment: The Effects of Charge Sharing Correction Algorithm Choice. SENSORS (BASEL, SWITZERLAND) 2020; 20:E6093. [PMID: 33120903 PMCID: PMC7663711 DOI: 10.3390/s20216093] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Revised: 10/23/2020] [Accepted: 10/24/2020] [Indexed: 11/17/2022]
Abstract
Most modern energy resolving, photon counting detectors employ small (sub 1 mm) pixels for high spatial resolution and low per pixel count rate requirements. These small pixels can suffer from a range of charge sharing effects (CSEs) that degrade both spectral analysis and imaging metrics. A range of charge sharing correction algorithms (CSCAs) have been proposed and validated by different groups to reduce CSEs, however their performance is often compared solely to the same system when no such corrections are made. In this paper, a combination of Monte Carlo and finite element methods are used to compare six different CSCAs with the case where no CSCA is employed, with respect to four different metrics: absolute detection efficiency, photopeak detection efficiency, relative coincidence counts, and binned spectral efficiency. The performance of the various CSCAs is explored when running on systems with pixel pitches ranging from 100 µm to 600µm, in 50 µm increments, and fluxes from 106 to 108 photons mm-2 s-1 are considered. Novel mechanistic explanations for the difference in performance of the various CSCAs are proposed and supported. This work represents a subset of a larger project in which pixel pitch, thickness, flux, and CSCA are all varied systematically.
Collapse
Affiliation(s)
- Oliver L. P. Pickford Scienti
- Joint Department of Physics, Institute of Cancer Research and Royal Marsden NHS Foundation Trust, London SM2 5NG, UK; (J.C.B.); (D.G.D.)
| | | | | |
Collapse
|
17
|
Taguchi K. Assessment of Multienergy Interpixel Coincidence Counters (MEICC) for Charge Sharing Correction or Compensation for Photon Counting Detectors With Boxcar Signals. IEEE TRANSACTIONS ON RADIATION AND PLASMA MEDICAL SCIENCES 2020; 5:465-475. [PMID: 34250325 DOI: 10.1109/trpms.2020.3003251] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
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.
Collapse
Affiliation(s)
- Katsuyuki Taguchi
- Russell H. Morgan Department of Radiology and Radiological Science, the Johns Hopkins University School of Medicine, Baltimore, MD
| |
Collapse
|