Monte Carlo model of a prototype flat-panel detector for multi-energy applications in radiotherapy.

BACKGROUND: The incorporation of multi-energy capabilities into radiotherapy flat-panel detectors offers advantages including enhanced soft tissue visualization by reduction of signal from overlapping anatomy such as bone in 2D image projections; creation of virtual monoenergetic images for 3D contrast enhancement, metal artefact reduction and direct acquisition of relative electron density. A novel dual-layer on-board imager offering dual energy processing capabilities is being designed. As opposed to other dual-energy implementation techniques which require separate acquisition with two different x-ray spectra, the dual-layer detector design enables simultaneous acquisition of high and low energy images with a single exposure. A computational framework is required to optimize the design parameters and evaluate detector performance for specific clinical applications. PURPOSE: In this study, we report on the development of a Monte Carlo (MC) model of the imager including model validation. METHODS: The stack-up of the dual-layer imager (DLI) was implemented in GEANT4 Application for Tomographic Emission (GATE). The DLI model has an active area of 43×43 cm2 , with top and bottom Cesium Iodide (CsI) scintillators of 600 and 800 µm thickness, respectively. Measurement of spatial resolution and imaging of dedicated multi-material dual-energy (DE) phantoms were used to validate the model. The modulation transfer function (MTF) of the detector was calculated for a 120 kVp x-ray spectrum using a 0.5 mm thick tantalum edge rotated by 2.5o . For imaging validation, the DE phantom was imaged using a 140 kVp x-ray spectrum. For both validation simulations, corresponding measurements were done using an initial prototype of the imager. Agreement between simulations and measurement was assessed using normalized root mean square error (NRMSE) and 1D profile difference for the MTF and phantom images respectively. Further comparison between measurement and simulation was made using virtual monoenergetic images (VMIs) generated from basis material images derived using precomputed look-up tables. RESULTS: The MTF of the bottom layer of the dual-layer model shows values decreasing more quickly with spatial frequency, compared to the top layer, due to the thicker bottom scintillator thickness and scatter from the top layer. A comparison with measurement shows NRMSE of 0.013 and 0.015 as well as identical MTF50 of 0.8 mm1 and 1.0 mm1 for the top and bottom layer respectively. For the DE imaging of the DE-phantom, although a maximum deviation of 3.3% is observed for the 10 mm aluminum and Teflon inserts at the top layer, the agreement for all other inserts is less than 2.2% of the measured value at both layers. Material decomposition of simulated scatter-free DE images gives an average accuracy in PMMA and aluminum composition of 4.9% and 10.3% for 11-30 mm PMMA and 1-10 mm aluminum objects respectively. A comparison of decomposed values using scatter containing measured and simulated DE images shows good agreement within statistical uncertainty. CONCLUSION: Validation using both MTF and phantom imaging shows good agreement between simulation and measurements. With the present configuration of the digital prototype, the model can generate material decomposed images and virtual monoenergetic images.

Frequency-dependent optimal weighting approach for megavoltage multilayer imagers.

Multi-layer imaging (MLI) devices improve the detective quantum efficiency (DQE) while maintaining the spatial resolution of conventional mega-voltage (MV) x-ray detectors for applications in radiotherapy. To date, only MLIs with identical detector layers have been explored. However, it may be possible to instead use different scintillation materials in each layer to improve the final image quality. To this end, we developed and validated a method for optimally combining the individual images from each layer of MLI devices that are built with heterogeneous layers. Two configurations were modeled within the GATE Monte Carlo package by stacking different layers of a terbium doped gadolinium oxysulfide Gd2O2S:Tb (GOS) phosphor and a LKH-5 glass scintillator. Detector response was characterized in terms of the modulation transfer function (MTF), normalized noise power spectrum (NNPS) and DQE. Spatial frequency-dependent weighting factors were then analytically derived for each layer such that the total DQE of the summed combination image would be maximized across all spatial modes. The final image is obtained as the weighted sum of the sub-images from each layer. Optimal weighting factors that maximize the DQE were found to be the quotient of MTF and NNPS of each layer in the heterogeneous MLI detector. Results validated the improvement of the DQE across the entire frequency domain. For the LKH-5 slab configuration, DQE(0) increases between 2%-3% (absolute), while the corresponding improvement for the LKH-5 pixelated configuration was 7%. The performance of the weighting method was quantitatively evaluated with respect to spatial resolution, contrast-to-noise ratio (CNR) and signal-to-noise ratio (SNR) of simulated planar images of phantoms at 2.5 and 6 MV. The line pair phantom acquisition exhibited a twofold increase in CNR and SNR, however MTF was degraded at spatial frequencies greater than 0.2 lp mm-1. For the Las Vegas phantom, the weighting improved the CNR by around 30% depending on the contrast region while the SNR values are higher by a factor of 2.5. These results indicate that the imaging performance of MLI systems can be enhanced using the proposed frequency-dependent weighting scheme. The CNR and SNR of the weighted combined image are improved across all spatial scales independent of the detector combination or photon beam energy.

GPU-accelerated Monte Carlo simulation of MV-CBCT.

Monte Carlo simulation (MCS) is one of the most accurate computation methods for dose calculation and image formation in radiation therapy. However, the high computational complexity and long execution time of MCS limits its broad use. In this paper, we present a novel strategy to accelerate MCS using a graphic processing unit (GPU), and we demonstrate the application in mega-voltage (MV) cone-beam computed tomography (CBCT) simulation. A new framework that generates a series of MV projections from a single simulation run is designed specifically for MV-CBCT acquisition. A Geant4-based GPU code for photon simulation is incorporated into the framework for the simulation of photon transport through a phantom volume. The FastEPID method, which accelerates the simulation of MV images, is modified and integrated into the framework. The proposed GPU-based simulation strategy was tested for its accuracy and efficiency in a Catphan 604 phantom and an anthropomorphic pelvis phantom with beam energies at 2.5 MV, 6 MV, and 6 MV FFF. In all cases, the proposed GPU-based simulation demonstrated great simulation accuracy and excellent agreement with measurement and CPU-based simulation in terms of reconstructed image qualities. The MV-CBCT simulation was accelerated by factors of roughly 900-2300 using an NVIDIA Tesla V100 GPU card against a 2.5 GHz AMD Opteron™ Processor 6380.

A rapid, accurate image simulation strategy for mega-voltage cone-beam computed tomography.

Intensive computation time is required to simulate images of electronic portal imaging device (EPID) using Monte Carlo (MC) technique, limiting the development of applications associated with EPID, such as mega-voltage cone-beam computed tomography (MV-CBCT). In this study, a fast, accurate simulation strategy for MV-CBCT utilizing the FastEPID technique has been developed and validated. During FastEPID simulation, photon detection was determined by pre-calculated photon energy deposition efficiency (η) and particle transport within the EPID was replaced with a pre-calculated optical photon spread function. This method is capable of reducing the time required for EPID image simulation by a factor of 90-140, without compromising image quality. MV-CBCT images reconstructed from the FastEPID simulated projections have been validated against measurement in terms of mean Hounsfield unit (HU), noise, and cupping artifact. These images were obtained with both a Catphan 604 phantom and an anthropomorphic pelvis phantom, under treatment beam energies of 2.5 MV, 6 MV, and 6 MV flattening filter free. The agreement between measurement and simulation was excellent in all cases. This novel strategy was capable of reducing the run time of a full scan simulation of MV-CBCT performed on a CPU cluster to a matter of hours, rather than weeks or months required by a conventional approach. Multiple applications associated with MV-CBCT (e.g. imager design optimization) are anticipated to gain from the implementation of this novel simulation strategy.

A novel method for fast image simulation of flat panel detectors.

We have developed a novel method for fast image simulation of flat panel detectors, based on the photon energy deposition efficiency and the optical spread function (OSF). The proposed method, FastEPID, determines the photon detection using photon energy deposition and replaces particle transport within the detector with precalculated OSFs. The FastEPID results are validated against experimental measurement and conventional Monte Carlo simulation in terms of modulation transfer function (MTF), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), contrast, and relative difference of pixel value, obtained with a slanted slit image, Las Vegas phantom, and anthropomorphic pelvis phantom. Excellent agreement is observed between simulation and measurement in all cases. Without degrading image quality, the FastEPID method is capable of reducing simulation time up to a factor of 150. Multiple applications, such as imager design optimization for planar and volumetric imaging, are expected to benefit from the implementation of the FastEPID method.