A framework for performance characterization of energy-resolving photon-counting detectors

PURPOSE

Photon-counting, energy-resolving detectors are subject to intense research interest, and there is a need for a general framework for performance assessment of these detectors. The commonly used linear-systems theory framework, which measures detector performance in terms of noise-equivalent quanta (NEQ) and detective quantum efficiency (DQE) is widely used for characterizing conventional x-ray detectors but does not take energy-resolving capabilities into account. The purpose of this work is to extend this framework to encompass energy-resolving photon-counting detectors and elucidate how the imperfect energy response and other imperfections in real-world detectors affect imaging performance, both for feature detection and for material quantification tasks.

METHOD

We generalize NEQ and DQE to matrix-valued quantities as functions of spatial frequency, and show how these matrices can be calculated from simple Monte Carlo simulations. To demonstrate how the new metrics can be interpreted, we compute them for simplified models of fluorescence and Compton scatter in a photon-counting detector and for a Monte Carlo model of a CdTe detector with 0.5 × 0.5 mm 2 pixels.

RESULTS

Our results show that the ideal-linear-observer performance for any detection or material quantification task can be calculated from the proposed generalized NEQ and DQE metrics. We also demonstrate that the proposed NEQ metric is closely related to a generalized version of the Cramér-Rao lower bound commonly used for assessing material quantification performance. Off-diagonal elements in the NEQ and DQE matrices are shown to be related to loss of energy information due to imperfect energy resolution. The Monte Carlo model of the CdTe detector predicts a zero-frequency dose efficiency relative to an ideal detector of 0.86 and 0.65 for detecting water and bone, respectively. When the task instead is to quantify these materials, the corresponding values are 0.34 for water and 0.26 for bone.

CONCLUSIONS

We have developed a framework for assessing the performance of photon-counting energy-resolving detectors and shown that the matrix-valued NEQ and DQE metrics contain sufficient information for calculating the dose efficiency for both detection and quantification tasks, the task having any spatial and energy dependence. This framework will be beneficial for the development and optimization of photon-counting x-ray detectors.

Theoretical Framework for the Dual-Energy Cone-Beam CT Noise-Power Spectrum, NEQ, and Tasked-Based Detectability Index

The optimization of dual-energy computed tomography (DE-CT) is challenged by the lack of a theoretical foundation for image quality. This work reports a cascaded systems analysis model that was used to derive signal and noise propagation in DE-CBCT in prevalent Fourier metrics such as the noise-power spectrum (NPS) and noise-equivalent quanta (NEQ). The model was validated in comparison to measurements of the 3D NPS and NEQ in DE-CBCT images acquired using an experimental imaging bench. Task-based detectability index was derived using DE-NPS and NEQ as an objective function in optimizing DE imaging parameters such as the dose allocation factor (D_A) and kVp pair. The resulting dose allocation optimization is in agreement with the practice of assigning more dose to the high-energy image (D_A < 0.5), and the model provides a quantitative basis for examining the optimal dose allocation as a function of total dose, kVp pair, the presence of electronics noise, and the imaging task. An example optimization is shown for a breast tumor detection task. Using DE decomposition to cancel fibroglandular tissue (rendering a DE-CBCT image of breast tumor against an adipose tissue background) and assuming a total dose of 15mGy, the optimal kVp pair is identified at [45, 105]kVp with D_A=0.46. The model is sufficiently general for applications beyond this example, demonstrating utility in the optimization in a broad range of imaging parameters. The model provides a new, valuable framework for understanding the theoretical limits of DE-CBCT imaging performance and maximizing image quality while minimizing radiation dose.

Task-based modeling and optimization of a cone-beam CT scanner for musculoskeletal imaging

PURPOSE

This work applies a cascaded systems model for cone-beam CT imaging performance to the design and optimization of a system for musculoskeletal extremity imaging. The model provides a quantitative guide to the selection of system geometry, source and detector components, acquisition techniques, and reconstruction parameters.

METHODS

The model is based on cascaded systems analysis of the 3D noise-power spectrum (NPS) and noise-equivalent quanta (NEQ) combined with factors of system geometry (magnification, focal spot size, and scatter-to-primary ratio) and anatomical background clutter. The model was extended to task-based analysis of detectability index (d') for tasks ranging in contrast and frequency content, and d' was computed as a function of system magnification, detector pixel size, focal spot size, kVp, dose, electronic noise, voxel size, and reconstruction filter to examine trade-offs and optima among such factors in multivariate analysis. The model was tested quantitatively versus the measured NPS and qualitatively in cadaver images as a function of kVp, dose, pixel size, and reconstruction filter under conditions corresponding to the proposed scanner.

RESULTS

The analysis quantified trade-offs among factors of spatial resolution, noise, and dose. System magnification (M) was a critical design parameter with strong effect on spatial resolution, dose, and x-ray scatter, and a fairly robust optimum was identified at M ∼ 1.3 for the imaging tasks considered. The results suggested kVp selection in the range of ∼65-90 kVp, the lower end (65 kVp) maximizing subject contrast and the upper end maximizing NEQ (90 kVp). The analysis quantified fairly intuitive results-e.g., ∼0.1-0.2 mm pixel size (and a sharp reconstruction filter) optimal for high-frequency tasks (bone detail) compared to ∼0.4 mm pixel size (and a smooth reconstruction filter) for low-frequency (soft-tissue) tasks. This result suggests a specific protocol for 1 × 1 (full-resolution) projection data acquisition followed by full-resolution reconstruction with a sharp filter for high-frequency tasks along with 2 × 2 binning reconstruction with a smooth filter for low-frequency tasks. The analysis guided selection of specific source and detector components implemented on the proposed scanner. The analysis also quantified the potential benefits and points of diminishing return in focal spot size, reduced electronic noise, finer detector pixels, and low-dose limits of detectability. Theoretical results agreed quantitatively with the measured NPS and qualitatively with evaluation of cadaver images by a musculoskeletal radiologist.

CONCLUSIONS

A fairly comprehensive model for 3D imaging performance in cone-beam CT combines factors of quantum noise, system geometry, anatomical background, and imaging task. The analysis provided a valuable, quantitative guide to design, optimization, and technique selection for a musculoskeletal extremities imaging system under development.

Anatomical background and generalized detectability in tomosynthesis and cone-beam CT

PURPOSE

Anatomical background presents a major impediment to detectability in 2D radiography as well as 3D tomosynthesis and cone-beam CT (CBCT). This article incorporates theoretical and experimental analysis of anatomical background "noise" in cascaded systems analysis of 2D and 3D imaging performance to yield "generalized" metrics of noise-equivalent quanta (NEQ) and detectability index as a function of the orbital extent of the (circular arc) source-detector orbit.

METHODS

A physical phantom was designed based on principles of fractal self-similarity to exhibit power-law spectral density (kappa/Fbeta) comparable to various anatomical sites (e.g., breast and lung). Background power spectra [S(B)(F)] were computed as a function of source-detector orbital extent, including tomosynthesis (approximately 10 degrees -180 degrees) and CBCT (180 degrees + fan to 360 degrees) under two acquisition schemes: (1) Constant angular separation between projections (variable dose) and (2) constant total number of projections (constant dose). The resulting S(B) was incorporated in the generalized NEQ, and detectability index was computed from 3D cascaded systems analysis for a variety of imaging tasks.

RESULTS

The phantom yielded power-law spectra within the expected spatial frequency range, quantifying the dependence of clutter magnitude (kappa) and correlation (beta) with increasing tomosynthesis angle. Incorporation of S(B) in the 3D NEQ provided a useful framework for analyzing the tradeoffs among anatomical, quantum, and electronic noise with dose and orbital extent. Distinct implications are posed for breast and chest tomosynthesis imaging system design-applications varying significantly in kappa and beta, and imaging task and, therefore, in optimal selection of orbital extent, number of projections, and dose. For example, low-frequency tasks (e.g., soft-tissue masses or nodules) tend to benefit from larger orbital extent and more fully 3D tomographic imaging, whereas high-frequency tasks (e.g., microcalcifications) require careful, application-specific selection of orbital extent and number of projections to minimize negative effects of quantum and electronic noise.

CONCLUSIONS

The complex tradeoffs among anatomical background, quantum noise, and electronic noise in projection imaging, tomosynthesis, and CBCT can be described by generalized cascaded systems analysis, providing a useful framework for system design and optimization.

The Generalized NEQ and Detectability Index for Tomosynthesis and Cone-Beam CT: From Cascaded Systems Analysis to Human Observers

PURPOSE

In the early development of new imaging modalities - such as tomosynthesis and cone-beam CT (CBCT) - an accurate predictive model for imaging performance is particularly valuable in identifying the physical factors that govern image quality and guiding system optimization. In this work, a task-based cascaded systems model for detectability index is proposed that describes not only the signal and noise propagation in the 2D (projection) and 3D (reconstruction) imaging chain but also the influence of background anatomical noise. The extent to which generalized detectability index provides a valid metric for imaging performance was assessed through direct comparison to human observer experiments.

METHODS

Detectability index (d') was generalized to include anatomical background noise in the same manner as the generalized noise-equivalent quanta (NEQ) proposed by Barrett et al. (Proc. SPIE Med. Imaging, Vol. 1090, 1989). Anatomical background noise was measured from a custom phantom designed to present power-law spectral density comparable to various anatomical sites (e.g., breast and lung). Theoretical calculations of d' as a function of the source-detector orbital extent (θtot) was obtained from a 3D cascaded systems analysis model for tomosynthesis and cone-beam CT (CBCT). Four model observers were considered in the calculation of d': prewhitening (PW), non-prewhitening (NPW), prewhitening with eye filter and internal noise (PWE), and non-prewhitening with eye filter and internal noise (NPWE). Human observer performance was measured from 9AFC tests for a variety of idealized imaging tasks presented within a clutter phantom. Theoretical results (d') were converted to area under the ROC curve (Az ) and compared directly to human observer performance as a function of imaging task and orbital extent.

RESULTS

Theoretical results demonstrated reasonable correspondence with human observer response for all tasks across the continuum in θtot ranging from low-angle tomosynthesis (θtot ~10°) to CBCT (θtot ~180°). Both theoretical and experimental Az were found to increase with acquisition angle, consistent with increased rejection of out-of-plane clutter for larger tomosynthesis angle. Of the four theoretical model observers considered, the prewhitening models tended to overestimate real observer performance, while the non-prewhitening models demonstrated reasonable agreement.

CONCLUSIONS

Generalized detectability index was shown to provide a meaningful metric for imaging performance, helping to bridge the gap between real observer performance and prevalent Fourier-based metrics based in first principles of spatial-frequency-dependent NEQ and imaging task.