Deriving depth-dependent light escape efficiency and optical Swank factor from measured pulse height spectra of scintillators

Protocol optimization for functional cardiac CT imaging using noise emulation in the raw data domain

BACKGROUND

Four-dimensional (4D) wide coverage computed tomography (CT) is an effective imaging modality for measuring the mechanical function of the myocardium. However, repeated CT measurement across a number of heartbeats is still a concern.

PURPOSE

A projection-domain noise emulation method is presented to generate accurate low-dose (mA modulated) 4D cardiac CT scans from high-dose scans, enabling protocol optimization to deliver sufficient image quality for functional cardiac analysis while using a dose level that is as low as reasonably achievable (ALARA).

METHODS

Given a targeted low-dose mA modulation curve, the proposed noise emulation method injects both quantum and electronic noise of proper magnitude and correlation to the high-dose data in projection domain. A spatially varying (i.e., channel-dependent) detector gain term as well as its calibration method were proposed to further improve the noise emulation accuracy. To determine the ALARA dose threshold, a straightforward projection domain image quality (IQ) metric was proposed that is based on the number of projection rays that do not fall under the non-linear region of the detector response. Experiments were performed to validate the noise emulation method with both phantom and clinical data in terms of visual similarity, contrast-to-noise ratio (CNR), and noise-power spectrum (NPS).

RESULTS

For both phantom and clinical data, the low-dose emulated images exhibited similar noise magnitude (CNR difference within 2%), artifacts, and texture to that of the real low-dose images. The proposed channel-dependent detector gain term resulted in additional increase in emulation accuracy. Using the proposed IQ metric, recommended kVp and mA settings were calculated for low dose 4D Cardiac CT acquisitions for patients of different sizes.

CONCLUSIONS

A detailed method to estimate system-dependent parameters for a raw-data based low dose emulation framework was described. The method produced realistic noise levels, artifacts, and texture with phantom and clinical studies. The proposed low-dose emulation method can be used to prospectively select patient-specific minimal-dose protocols for functional cardiac CT.

Protocol Optimization for Functional Cardiac CT Imaging Using Noise Emulation in the Raw Data Domain

Background

Four-dimensional (4D) wide coverage computed tomography (CT) is an effective imaging modality for measuring the mechanical function of the myocardium. However, repeated CT measurement across a number of heartbeats is still a concern.

Purpose

A projection-domain noise emulation method is presented to generate accurate low-dose (mA modulated) 4D cardiac CT scans from high-dose scans, enabling protocol optimization to deliver sufficient image quality for functional cardiac analysis while using a dose level that is as low as reasonably achievable (ALARA).

Methods

Given a targeted low-dose mA modulation curve, the proposed noise emulation method injects both quantum and electronic noise of proper magnitude and correlation to the high-dose data in projection domain. A spatially varying (i.e., channel-dependent) detector gain term as well as its calibration method were proposed to further improve the noise emulation accuracy. To determine the ALARA dose threshold, a straightforward projection domain image quality (IQ) metric was proposed that is based on the number of projection rays that do not fall under the non-linear region of the detector response. Experiments were performed to validate the noise emulation method with both phantom and clinical data in terms of visual similarity, contrast-to-noise ratio (CNR), and noise-power spectrum (NPS).

Results

For both phantom and clinical data, the low-dose emulated images exhibited similar noise magnitude (CNR difference within 2%), artifacts, and texture to that of the real low-dose images. The proposed channel-dependent detector gain term resulted in additional increase in emulation accuracy. Using the proposed IQ metric, recommended kVp and mA settings were calculated for low dose 4D Cardiac CT acquisitions for patients of different sizes.

Conclusions

A detailed method to estimate system-dependent parameters for a raw-data based low dose emulation framework was described. The method produced realistic noise levels, artifacts, and texture with phantom and clinical studies. The proposed low-dose emulation method can be used to prospectively select patient-specific minimal-dose protocols for functional cardiac CT.

Feasibility of 4D HDR brachytherapy source tracking using x-ray tomosynthesis: Monte Carlo investigation

PURPOSE

High dose rate (HDR) brachytherapy rapidly delivers dose to targets with steep dose gradients. This treatment method must adhere to prescribed treatment plans with high spatiotemporal accuracy and precision, as failure to do so may degrade clinical outcomes. One approach to achieving this goal is to develop imaging techniques to track HDR sources in vivo in reference to surrounding anatomy. This work investigates the feasibility of using an isocentric C-arm x-ray imager and tomosynthesis methods to track Ir-192 HDR brachytherapy sources in vivo over time (4D).

METHODS

A tomosynthesis imaging workflow was proposed and its achievable source detectability, localization accuracy, and spatiotemporal resolution were investigated in silico. An anthropomorphic female XCAT phantom was modified to include a vaginal cylinder applicator and Ir-192 HDR source (0.5 × 0.5 × 5.0 mm3 ), and the workflow was carried out using the MC-GPU Monte Carlo image simulation platform. Source detectability was characterized using the reconstructed source signal-difference-to-noise-ratio (SDNR), localization accuracy by the absolute 3D error in its measured centroid location, and spatiotemporal resolution by the full-width-at-half-maximum (FWHM) of line profiles through the source in each spatial dimension considering a maximum C-arm angular velocity of 30° per second. The dependence of these parameters on acquisition angular range (θtot = 0°-90°), number of views, angular increment between views (Δθ = 0°-15°), and volumetric constraints imposed in reconstruction was evaluated. Organ voxel doses were tallied to derive the workflow's attributable effective dose.

RESULTS

The HDR source was readily detected and its centroid was accurately localized with the proposed workflow and method (SDNR: 10-40, 3D error: 0-0.144 mm). Tradeoffs were demonstrated for various combinations of image acquisition parameters; namely, increasing the tomosynthesis acquisition angular range improved resolution in the depth-encoded direction, for example from 2.5 mm to 1.2 mm between θtot = 30o and θtot = 90o , at the cost of increasing acquisition time from 1 to 3 s. The best-performing acquisition parameters (θtot = 90o , Δθ = 1°) yielded no centroid localization error, and achieved submillimeter source resolution (0.57 × 1.21 × 5.04 mm3 apparent source dimensions, FWHM). The total effective dose for the workflow was 263 µSv for its required pre-treatment imaging component and 7.59 µSv per mid-treatment acquisition thereafter, which is comparable to common diagnostic radiology exams.

CONCLUSIONS

A system and method for tracking HDR brachytherapy sources in vivo using C-arm tomosynthesis was proposed and its performance investigated in silico. Tradeoffs in source conspicuity, localization accuracy, spatiotemporal resolution, and dose were determined. The results suggest this approach is feasible for localizing an Ir-192 HDR source in vivo with submillimeter spatial resolution, 1-3 second temporal resolution and minimal additional dose burden.

Computational models of direct and indirect X-ray breast imaging detectors for in silico trials

BACKGROUND

To facilitate in silico studies that investigate digital mammography (DM) and breast tomosynthesis (DBT), models replicating the variety in imaging performance of the DM and DBT systems, observed across manufacturers are needed.

PURPOSE

The main purpose of this work is to develop generic physics models for direct and indirect detector technology used in commercially available systems, with the goal of making them available open source to manufacturers to further tweak and develop the exact in silico replicas of their systems.

METHODS

We recently reported on an in silico version of the SIEMENS Mammomat Inspiration DM/DBT system using an open-source GPU-accelerated Monte Carlo x-ray imaging simulation code (MC-GPU). We build on the previous version of the MC-GPU codes to mimic the imaging performances of two other Food and Drug Administration (FDA)-approved DM/DBT systems, such as Hologic Selenia Dimensions (HSD) and the General Electric Senographe Pristina (GSP) systems. In this work, we developed a hybrid technique to model the optical spread and signal crosstalk observed in the GSP and HSD systems. MC simulations are used to track each x-ray photon till its first interaction within the x-ray detector. On the other hand, the signal spread in the x-ray detectors is modeled using previously developed analytical equations. This approach allows us to preserve the modeling accuracy offered by MC methods in the patient body, while speeding up secondary carrier transport (either electron-hole pairs or optical photons) using analytical equations in the detector. The analytical optical spread model for the indirect detector includes the depth-dependent spread and collection of optical photons and relies on a pre-computed set of point response functions that describe the optical spread as a function of depth. To understand the capabilities of the computational x-ray detector models, we compared image quality metrics like modulation transfer function (MTF), normalized noise power spectrum (NNPS), and detective quantum efficiency (DQE), simulated with our models against measured data. Please note that the purpose of these comparisons with measured data would be to gauge if the model developed as part of this work could replicate commercially used direct and indirect technology in general and not to achieve perfect fits with measured data.

RESULTS

We found that the simulated image quality metrics such as MTF, NNPS, and DQE were in reasonable agreement with experimental data. To demonstrate the imaging performance of the three DM/DBT systems, we integrated the detector models with the VICTRE pipeline and simulated DM images of a fatty breast model containing a spiculated mass and a calcium oxalate cluster. In general, we found that the images generated using the indirect model appeared more blurred with a different noise texture and contrast as compared to the systems with direct detectors.

CONCLUSIONS

We have presented computational models of three commercially available FDA-approved DM/DBT systems, which implement both direct and indirect detector technology. The updated versions of the MC-GPU codes that can be used to replicate three systems are available in open source format through GitHub.

High-resolution model-based material decomposition in dual-layer flat-panel CBCT

PURPOSE

Spectral CT uses energy-dependent measurements that enable material discrimination in addition to reconstruction of structural information. Flat-panel detectors (FPDs) have been widely used in dedicated and interventional systems to deliver high spatial resolution, volumetric cone-beam CT (CBCT) in compact and OR-friendly designs. In this work, we derive a model-based method that facilitates high-resolution material decomposition in a spectral CBCT system equipped with a prototype dual-layer FPD. Through high-fidelity modeling of multilayer detector, we seek to avoid resolution loss that is present in more traditional processing and decomposition approaches.

METHOD

A physical model for spectral measurements in dual-layer flat-panel CBCT is developed including layer-dependent differences in system geometry, spectral sensitivities, and detector blur (e.g., due to varied scintillator thicknesses). This forward model is integrated into a model-based material decomposition (MBMD) method based on minimization of a penalized weighted least-squared (PWLS) objective function. The noise and resolution performance of this approach was compared with traditional projection-domain decomposition (PDD) and image-domain decomposition (IDD) approaches as well as one-step MBMD with lower-fidelity models that use approximated geometry, projection interpolation, or an idealized system geometry without system blur model. Physical studies using high-resolution three-dimensional (3D)-printed water-iodine phantoms were conducted to demonstrate the high-resolution imaging performance of the compared decomposition methods in iodine basis images and synthetic monoenergetic images.

RESULTS

Physical experiments demonstrate that the MBMD methods incorporating an accurate geometry model can yield higher spatial resolution iodine basis images and synthetic monoenergetic images than PDD and IDD results at the same noise level. MBMD with blur modeling can further improve the spatial-resolution compared with the decomposition results obtained with IDD, PDD, and MBMD methods with lower-fidelity models. Using the MBMD without or with blur model can increase the absolute modulation at 1.75 lp/mm by 10% and 22% compared with IDD at the same noise level.

CONCLUSION

The proposed model-based material decomposition method for a dual-layer flat-panel CBCT system has demonstrated an ability to extend high-resolution performance through sophisticated detector modeling including the layer-dependent blur. The proposed work has the potential to not only facilitate high-resolution spectral CT in interventional and dedicated CBCT systems, but may also provide the opportunity to evaluate different flat-panel design trade-offs including multilayer FPDs with mismatched geometries, scintillator thicknesses, and spectral sensitivities.

Comparison of CsI:Tl and Gd2 O2 S:Tb indirect flat panel detector x-ray imaging performance in front- and back-irradiation geometries

PURPOSE

The detective quantum efficiency (DQE) of indirect flat panel detectors (I-FPDs) is limited at higher x-ray energies (e.g., 100-140 kVp) by low absorption in their scintillating x-ray conversion layer. While increasing the thickness of the scintillator can improve its x-ray absorption efficiency, this approach is potentially limited by reduced spatial resolution and increased noise due to depth dependence in the scintillator's response to x rays. One strategy proposed to mitigate these deleterious effects is to irradiate the scintillator through the pixel sensor in a "back-irradiation" geometry. This work directly evaluates the impact of irradiation geometry on the inherent imaging performance of I-FPDs composed with columnar CsI:Tl and powder Gd2 O2 S:Tb (GOS) scintillators.

METHODS

A "bidirectional" FPD was constructed which allows scintillator samples to be interchangeably coupled with the detector's active matrix to compose an I-FPD. Radio-translucent windows in the detector's housing permit imaging in both "front-irradiation" (FI) and "back-irradiation" (BI) geometries. This test device was used to evaluate the impact of irradiation geometry on the x-ray sensitivity, modulation transfer function (MTF), noise power spectrum (NPS), and DQE of four I-FPDs composed using columnar CsI:Tl scintillators of varying thickness (600-1000 µm) and optical backing, and a Fast Back GOS screen. All experiments used an RQA9 x-ray beam.

RESULTS

Each I-FPD's x-ray sensitivity, MTF, and DQE was greater or equal in BI geometry than in FI. The I-FPD composed with CsI:Tl (1 mm) and an optically absorptive backing had the largest variation in sensitivity (17%) between FI and BI geometries. The detector composed with GOS had the largest improvement in limiting resolution (31%). Irradiation geometry had little impact on MTF(f) and DQE(f) measurements near zero frequency, however, the difference between FI and BI measurements generally increased with spatial frequency. The CsI:Tl scintillator with optically absorptive backing (1 mm) in BI geometry had the highest spatial resolution and DQE over all frequencies.

CONCLUSIONS

Back irradiation may improve the inherent x-ray imaging performance of I-FPDs composed with CsI:Tl and GOS scintillators. This approach can be leveraged to improve tradeoffs between detector dose efficiency, spatial resolution and noise for higher energy x-ray imaging.

Back-irradiated and dual-screen sandwich detector configurations for radiography

Recent advances in thin film transistor array technology have enabled the possibility of "back-irradiated" (BI) indirect active-matrix flat-panel imagers (AMFPIs), in which x-rays first expose the scintillator through the optical sensor, and "dual-screen" AMFPIs, in which two scintillating screens are sandwiched around a bidirectional active matrix. We developed a theoretical treatment of such detectors. The theory is used to investigate possible imaging performance improvements over conventional "front-irradiation" (FI) AMFPIs, where the active matrix is opposite the x-ray entrance surface. Simple expressions for the modulation transfer function, normalized noise power spectrum, Swank factor ( A s ), Lubberts function L ( f ) , and spatial frequency-dependent detective quantum efficiency DQE ( f ) are derived and used to compute these quantities for a variety of FI, BI, and dual-screen detector configurations. DQE ( f ) is used as the figure of merit for optimizing and comparing the performance of the various configurations. Large performance improvements over FI single-screen systems are found possible with BI. Further improvements are found possible with dual-screen configurations. The ratio of the thicknesses of the two screens that optimizes DQE is generally asymmetric, with the thinner screen facing the incident flux. The optimum ratio depends on the x-ray attenuation length in the screen.

An apparatus and method for directly measuring the depth-dependent gain and spatial resolution of turbid scintillators

PURPOSE

Turbid (powder or columnar-structured) scintillators are widely used in indirect flat panel detectors (I-FPDs) for scientific, industrial, and medical radiography. Light diffusion and absorption within these scintillators is expected to cause depth-dependent variations in their x ray conversion gain and spatial blur. These variations degrade the detective quantum efficiency of I-FPDs at all spatial frequencies. Despite their importance, there are currently no established methods for directly measuring scintillator depth effects. This work develops the instrumentation and methods to achieve this capability.

METHODS

An ultra-high-sensitivity camera was assembled for imaging single x ray interactions in two commercial Gd₂ O₂ S:Tb (GOS) screens (Lanex Regular and Fast Back, Eastman Kodak Company). X ray interactions were localized to known depths in the screens using a slit beam of parallel synchrotron radiation (32 keV), with beam width (~20 μm) much narrower than the screen thickness. Depth-localized x ray interaction images were acquired in 30 μm depth-intervals, and analyzed to measure each scintillator's depth-dependent average gain g ¯ ( z ) and modulation transfer function MTF(z,f). These measurements were used to calculate each screen's expected MTF(f) in an energy-integrating detector (e.g., I-FPD). Calculations were compared to presampling MTF measurements made by coupling each screen to a high-resolution CMOS image sensor (48 μm pixel) and using the slanted-edge method.

RESULTS

Both g ¯ ( z ) and MTF(z,f) continuously increased as interactions occurred closer to each screen's sensor-coupled surface. The Regular yielded 1351 ± 66 and 2117 ± 54 photons per absorbed x ray (42-66 keV^-1 ) in interactions occurring furthest from and nearest to the image sensor, while the Fast Back yielded 833 ± 22 and 1910 ± 39 photons (26-60 keV^-1 ). At f = 1 mm^-1 , MTF(z,f) varied between 0.63 and 0.78 in the Regular and 0.30-0.76 in the Fast Back. Calculations of presampling MTF(f) using g ¯ ( z ) and MTF(z,f) showed excellent agreement with slanted-edge measurements.

CONCLUSIONS

The developed instrument and method enable direct measurements of the depth-dependent gain and spatial resolution of turbid scintillators. This knowledge can be used to predict, understand, and potentially improve I-FPD imaging performance.

Modeling and evaluation of a high-resolution CMOS detector for cone-beam CT of the extremities

PURPOSE

Quantitative assessment of trabecular bone microarchitecture in extremity cone-beam CT (CBCT) would benefit from the high spatial resolution, low electronic noise, and fast scan time provided by complementary metal-oxide semiconductor (CMOS) x-ray detectors. We investigate the performance of CMOS sensors in extremity CBCT, in particular with respect to potential advantages of thin (<0.7 mm) scintillators offering higher spatial resolution.

METHODS

A cascaded systems model of a CMOS x-ray detector incorporating the effects of CsI:Tl scintillator thickness was developed. Simulation studies were performed using nominal extremity CBCT acquisition protocols (90 kVp, 0.126 mAs/projection). A range of scintillator thickness (0.35-0.75 mm), pixel size (0.05-0.4 mm), focal spot size (0.05-0.7 mm), magnification (1.1-2.1), and dose (15-40 mGy) was considered. The detectability index was evaluated for both CMOS and a-Si:H flat-panel detector (FPD) configurations for a range of imaging tasks emphasizing spatial frequencies associated with feature size aobj. Experimental validation was performed on a CBCT test bench in the geometry of a compact orthopedic CBCT system (SAD = 43.1 cm, SDD = 56.0 cm, matching that of the Carestream OnSight 3D system). The test-bench studies involved a 0.3 mm focal spot x-ray source and two CMOS detectors (Dalsa Xineos-3030HR, 0.099 mm pixel pitch) - one with the standard CsI:Tl thickness of 0.7 mm (C700) and one with a custom 0.4 mm thick scintillator (C400). Measurements of modulation transfer function (MTF), detective quantum efficiency (DQE), and CBCT scans of a cadaveric knee (15 mGy) were obtained for each detector.

RESULTS

Optimal detectability for high-frequency tasks (feature size of ~0.06 mm, consistent with the size of trabeculae) was ~4× for the C700 CMOS detector compared to the a-Si:H FPD at nominal system geometry of extremity CBCT. This is due to ~5× lower electronic noise of a CMOS sensor, which enables input quantum-limited imaging at smaller pixel size. Optimal pixel size for high-frequency tasks was <0.1 mm for a CMOS, compared to ~0.14 mm for an a-Si:H FPD. For this fine pixel pitch, detectability of fine features could be improved by using a thinner scintillator to reduce light spread blur. A 22% increase in detectability of 0.06 mm features was found for the C400 configuration compared to C700. An improvement in the frequency at 50% modulation (f50 ) of MTF was measured, increasing from 1.8 lp/mm for C700 to 2.5 lp/mm for C400. The C400 configuration also achieved equivalent or better DQE as C700 for frequencies above ~2 mm-1 . Images of cadaver specimens confirmed improved visualization of trabeculae with the C400 sensor.

CONCLUSIONS

The small pixel size of CMOS detectors yields improved performance in high-resolution extremity CBCT compared to a-Si:H FPDs, particularly when coupled with a custom 0.4 mm thick scintillator. The results indicate that adoption of a CMOS detector in extremity CBCT can benefit applications in quantitative imaging of trabecular microstructure in humans.

Toward Scintillator High-Gain Avalanche Rushing Photoconductor Active Matrix Flat Panel Imager (SHARP-AMFPI): Initial fabrication and characterization

PURPOSE

We present the first prototype Scintillator High-Gain Avalanche Rushing Photoconductor Active Matrix Flat Panel Imager (SHARP-AMFPI). This detector includes a layer of avalanche amorphous Selenium (a-Se) (HARP) as the photoconductor in an indirect detector to amplify the signal and reduce the effects of electronic noise to obtain quantum noise-limited images for low-dose applications. It is the first time avalanche a-Se has been used in a solid-state imaging device and poses as a possible solution to eliminate the effects of electronic noise, which is crucial for low-dose imaging performance of AMFPI.

METHODS

We successfully deposited a solid-state HARP structure onto a 24 × 30 cm² array of thin-film transistors (TFT array) with a pixel pitch of 85 μm. The HARP layer consists of 16 μm of a-Se with a hole-blocking and electron-blocking layer to prevent charge injection from the high-voltage bias and pixel electrodes, respectively. An electric field (E_Se ) up to 105 V μm^-1 was applied across the a-Se layer without breakdown. A 150 μm thick-structured CsI:Tl scintillator was used to form SHARP-AMFPI. The x-ray imaging performance is characterized using a 30 kVp Mo/Mo beam. We evaluate the spatial resolution, noise power, and detective quantum efficiency at zero frequency of the system with and without avalanche gain. The results are analyzed using cascaded linear system model (CLSM).

RESULTS

An avalanche gain of 76 ± 5 was measured at E_Se = 105 V μm^-1 . We demonstrate that avalanche gain can amplify the signal to overcome electronic noise. As avalanche gain is increased, image quality improves for a constant (0.76 mR) exposure until electronic noise is overcome. Our system is currently limited by poor optical transparency of our high-voltage electrode and long integrating time which results in dark current noise. These two effects cause high-spatial frequency noise to dominate imaging performance.

CONCLUSIONS

We demonstrate the feasibility of a solid-state HARP x-ray imager and have fabricated the largest active area HARP sensor to date. Procedures to reduce secondary quantum and dark noise are outlined. Future work will improve optical coupling and charge transport which will allow for frequency DQE and temporal metrics to be obtained.