Spatial frequency‐dependent pulse‐height spectrum and method for analyzing detector DQE( f ) from ensembles of single X‐ray images

Scott Dow,Adrian Howansky,Anthony R. Lubinsky,Wei Zhao
DOI: https://doi.org/10.1002/mp.15344
IF: 4.506
2021-11-29
Medical Physics
Abstract:PURPOSE: Scintillators and photoconductors used in energy integrating detectors (EIDs) have inherent variations in their imaging response to single-detected X-rays due to variations in X-ray energy deposition and secondary quanta generation and transport, which degrades DQE(f). The imaging response of X-ray scintillators to single X-rays may be recorded and studied using single X-ray imaging (SXI) experiments; however, no method currently exists for relating SXI experimental results to EID DQE(f). This work proposes a general analytical framework for computing and analyzing the DQE(f) performance of EIDs from single X-ray image ensembles using a spatial frequency-dependent pulse-height spectrum.METHODS: A spatial frequency (f)-dependent gain, <math> <semantics> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\tilde{g}(f)$</annotation></semantics> </math> , is defined as the Fourier transform of the imaging response of an EID to a single-detected X-ray. A f-dependent pulse-height spectrum, <math> <semantics><mrow><mo>Pr</mo> <mo>[</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> <mo>]</mo></mrow> <annotation>$\Pr [\tilde{g}(f)]$</annotation></semantics> </math> , is defined as the 2D probability density function of <math> <semantics> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\tilde{g}(f)$</annotation></semantics> </math> over the complex plane. <math> <semantics><mrow><mo>Pr</mo> <mo>[</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> <mo>]</mo></mrow> <annotation>$\Pr [\tilde{g}(f)]$</annotation></semantics> </math> is used to define a f-dependent Swank factor, A<sub>S</sub> (f), which fully characterizes the DQE(f) degradation due to single X-ray noise. A<sub>S</sub> (f) is analyzed in terms of its degradation due to Swank noise, variations in the frequency-dependent attenuation of <math> <semantics> <mrow><mrow><mo>|</mo></mrow> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <mrow><mo>|</mo></mrow> </mrow> <annotation>$| {\tilde{g}(f)} |$</annotation></semantics> </math> , and noise in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> which occurs due to variations in the asymmetry in each single X-ray's imaging response. Three example imaging systems are simulated to demonstrate the impact of depth-dependent variation in <math> <semantics> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\tilde{g}(f)$</annotation></semantics> </math> , remote energy deposition, and a finite number of secondary quanta, on <math> <semantics><mrow><mo>Pr</mo> <mo>[</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> <mo>]</mo></mrow> <annotation>$\Pr [\tilde{g}(f)]$</annotation></semantics> </math> , A<sub>S</sub> (f), MTF(f), and NPS(f)/NPS(0), which are computed from ensembles of single X-ray images. The same is also demonstrated by simulating a realistic imaging system; that is, a Gd<sub>2</sub> O<sub>2</sub> S-based EID. Using the latter imaging system, the convergence of A<sub>S</sub> (f) estimates is investigated as a function of the number of detected X-rays per ensemble.RESULTS: Depth-dependent <math> <semantics> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\tilde{g}(f)$</annotation></semantics> </math> variation resulted in A<sub>S</sub> (f) degradation exclusively due to depth-dependent optical Swank noise and the Lubberts effect. Conversely, the majority of A<sub>S</sub> (f) degradation caused by remote energy deposition and finite secondary quanta occurred due to variations in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> . When using input X-ray energies below the K-edge of Gd, variations in the frequency-dependent attenuation of <math> <semantics> <mrow><mrow><mo>|</mo></mrow> <mrow><mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <mrow><mo>|</mo></mrow> </mrow> <annotation>$| {\tilde{g}(f)} |$</annotation></semantics> </math> accounted for the majority of A<sub>S</sub> (f) degradation in the GOS-based EID, and very little Swank noise and variations in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> were observed. Above the K-edge, however, A<sub>S</sub> (f) degradation due to Swank noise and variations in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> greatly increased. The convergence of A<sub>S</sub> (f) was limited by variation in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> ; imaging systems with more variation in <math> <semantics><mrow><mo>arg</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> </mrow> <annotation>$\arg \tilde{g}(f)$</annotation></semantics> </math> required more detected X-rays per ensemble.CONCLUSIONS: An analytical framework is proposed that generalizes the pulse-height spectrum and Swank factor to arbitrary f. The impact of single X-ray noise sources, such as the Lubberts effect, remote energy deposition, and finite secondary quanta on detector performance, may be represented using <math> <semantics><mrow><mo>Pr</mo> <mo>[</mo> <mover><mi>g</mi> <mo>∼</mo></mover> <mrow><mo>(</mo> <mi>f</mi> <mo>)</mo></mrow> <mo>]</mo></mrow> <annotation>$\Pr [\tilde{g}(f)]$</annotation></semantics> </math> , and quantified using A<sub>S</sub> (f). The approach may be used to compute MTF(f), NPS(f), and DQE(f) from ensembles of single X-ray images and provides an additional tool to analyze proposed EID designs.
radiology, nuclear medicine & medical imaging
What problem does this paper attempt to address?