Daniel Morf

Rapid Monte Carlo simulation of detector DQE(f)

PURPOSE Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10(7) - 10(9) detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS. METHODS The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels. RESULTS On the order of only 10-100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10(7) reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10(6) CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU. CONCLUSIONS The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.

Tracking latency in image-based dynamic MLC tracking with direct image access

Abstract Purpose. Target tracking is a promising method for motion compensation in radiotherapy. For image-based dynamic multileaf collimator (DMLC) tracking, latency has been shown to be the main contributor to geometrical errors in tracking of respiratory motion, specifically due to slow transfer of image data from the image acquisition system to the tracking system via image file storage on a hard disk. The purpose of the current study was to integrate direct image access with a DMLC tracking system and to quantify the tracking latency of the integrated system for both kV and MV image-based tracking. Method. A DMLC tracking system integrated with a linear accelerator was used for tracking of a motion phantom with an embedded tungsten marker. Real-time target localization was based on x-ray images acquired either with a portal imager or a kV imager mounted orthogonal to the treatment beam. Images were processed directly without intermediate disk access. Continuous portal images and system log files were stored during treatment delivery for detailed offline analysis of the tracking latency. Results. The mean tracking system latency for kV and MV image-based tracking as function of the imaging interval ΔTimage increased linearly with ΔTimage as 148 ms + 0.58 * ΔTimage (kV) and 162 ms + 1.1 * ΔTimage (MV). The latency contribution from image acquisition and image transfer for kV image-based tracking was independent on ΔTimage at 103 ± 14 ms. For MV-based tracking, it increased with ΔTimage as 124 ms + 0.44 * ΔTimage. For ΔTimage = 200 ms (5 Hz imaging), the total latency was reduced from 550 ms to 264 ms for kV image-based tracking and from 500 ms to 382 ms for MV image-based tracking as compared to the previously used indirect image transfer via image file storage on a hard disk. Conclusion. kV and MV image-based DMLC tracking was successfully integrated with direct image access. It resulted in substantial tracking latency reductions compared with image-based tracking without direct image access.

Development and testing of an improved dosimetry system using a backscatter shielded electronic portal imaging device

PURPOSE To investigate the properties of a modified backscatter shielded electronic portal imaging device (BSS-EPID) and to develop a dose model to convert BSS-EPID images to dose in water as part of an improved system for dosimetry using EPIDs. METHODS The effectiveness of the shielding of the BSS-EPID was studied by comparing images measured with the BSS-EPID mounted on the support arm to images measured with the BSS-EPID removed from the support arm. A dose model was developed and optimized to reconstruct dose in water at different depths from measured BSS-EPID images. The accuracy of the dose model was studied using BSS-EPID images of 28 IMRT fields to reconstruct dose in water at depths of 2, 5, 10, and 20 cm and comparing to measured dose in water from a two-dimensional diode array at the same depths. The ability of the BSS-EPID system to operate independently of detector position was demonstrated by comparing the dose reconstruction of a 10 × 10 cm(2) field using different detector offsets to that measured by a two-dimensional diode array. RESULTS The shielding of the BSS-EPID was found to be effective, with more than 99% of pixels showing less than 0.5% change due to the presence of the support arm and at most a 0.2% effect on the central axis for 2 × 2 cm(2) fields to fully open 30 × 40 cm(2) images. The dose model was shown to accurately reconstruct measurements of dose in water using BSS-EPID images with average γ pass rates (2%, 2 mm criteria) of 92.5%, 98.7%, 97.4%, and 97.2% at depths of 2, 5, 10, and 20 cm, respectively, when compared to two-dimensional diode array measurements. When using 3%, 3 mm γ criteria, the average pass rate was greater than 97% at all depths. Reconstructed dose in water for a 10 × 10 cm(2) field measured with detector offsets as large as 10 cm agreed with each other and two-dimensional diode array measurements within 0.9%. CONCLUSIONS The modified BSS-EPID and associated dose model provide an improved system for dosimetry measurements using EPIDs. Several important limitations of the current hardware and software are addressed by this system.

A novel EPID design for enhanced contrast and detective quantum efficiency

Beams-eye-view imaging applications such as real-time soft-tissue motion estimation are hindered by the inherently low image contrast of electronic portal imaging devices (EPID) currently available for clinical use. We introduce and characterize a novel EPID design that provides substantially increased detective quantum efficiency (DQE), contrast-to-noise ratio (CNR) and sensitivity without degradation in spatial resolution. The prototype design features a stack of four conventional EPID layers combined with low noise integrated readout electronics. Each layer consists of a copper plate, a scintillator ([Formula: see text]) and a photodiode/TFT-switch (aSi:H). We characterize the prototypes signal response to a 6 MV photon beam in terms of modulation transfer function (MTF), DQE and CNR. The presampled MTF is estimated using a slanted slit technique, the DQE is calculated from measured normalized noise power spectra (nNPS) and the MTF and CNR is estimated using a Las Vegas contrast phantom. The prototype has been designed and built to be interchangeable with the current clinical EPID on the Varian TrueBeam platform (AS-1200) in terms of size and data output specifications. Performance evaluation is conducted in absolute values as well as in relative terms using the Varian AS-1200 EPID as a reference detector. A fivefold increase of DQE(0) to about 6.7% was observed by using the four-layered design versus the AS-1200 reference detector. No substantial differences are observed between each layers individual MTF and the one for all four layers operating combined indicating that defocusing due to beam divergence is negligible. Also, using four layers instead of one increases the signal to noise ratio by a factor of 1.7.

Simultaneous MV-kV imaging for intrafractional motion management during volumetric-modulated arc therapy delivery*

The purpose of this study was to evaluate the accuracy and clinical feasibility of a motion monitoring method employing simultaneously acquired MV and kV images during volumetric‐modulated arc therapy (VMAT). Short‐arc digital tomosynthesis (SA‐DTS) is used to improve the quality of the MV images that are then combined with orthogonally acquired kV images to assess 3D motion. An anthropomorphic phantom with implanted gold seeds was used to assess accuracy of the method under static, typical prostatic, and respiratory motion scenarios. Automatic registration of kV images and single MV frames or MV SA‐DTS reconstructed with arc lengths from 2° to 7° with the appropriate reference fiducial template images was performed using special purpose‐built software. Clinical feasibility was evaluated by retrospectively analyzing images acquired over four or five sessions for each of three patients undergoing hypofractionated prostate radiotherapy. The standard deviation of the registration error in phantom using MV SA‐DTS was similar to single MV images for the static and prostate motion scenarios (σ=0.25 mm). Under respiratory motion conditions, the standard deviation of the registration error increased to 0.7 mm and 1.7 mm for single MV and MV SA‐DTS, respectively. Registration failures were observed with the respiratory scenario only and were due to motion‐induced fiducial blurring. For the three patients studied, the mean and standard deviation of the difference between automatic registration using 4° MV SA‐DTS and manual registration using single MV images results was 0.07±0.52 mm. The MV SA‐DTS results in patients were, on average, superior to single‐frame MV by nearly 1 mm — significantly more than what was observed in phantom. The best MV SA‐DTS results were observed with arc lengths of 3° to 4°. Registration failures in patients using MV SA‐DTS were primarily due to blockage of the gold seeds by the MLC. The failure rate varied from 2% to 16%. Combined MV SA‐DTS and kV imaging is feasible for intratreatment motion monitoring during VMAT of anatomic sites where limited motion is expected, and improves registration accuracy compared to single MV/kV frames. To create a clinically robust technique, further improvements to ensure visualization of fiducials at the desired control points without degradation of the treatment plan are needed. PACS number(s): 87.55.km, 87.55.N‐

Dosimetric properties of an amorphous silicon EPID for verification of modulated electron radiotherapy

PURPOSE To investigate the dosimetric properties of an electronic portal imaging device (EPID) for electron beam detection and to evaluate its potential for quality assurance (QA) of modulated electron radiotherapy (MERT). METHODS A commercially available EPID was used to detect electron beams shaped by a photon multileaf collimator (MLC) at a source-surface distance of 70 cm. The fundamental dosimetric properties such as reproducibility, dose linearity, field size response, energy response, and saturation were investigated for electron beams. A new method to acquire the flood-field for the EPID calibration was tested. For validation purpose, profiles of open fields and various MLC fields (square and irregular) were measured with a diode in water and compared to the EPID measurements. Finally, in order to use the EPID for QA of MERT delivery, a method was developed to reconstruct EPID two-dimensional (2D) dose distributions in a water-equivalent depth of 1.5 cm. Comparisons were performed with film measurement for static and dynamic monoenergy fields as well as for multienergy fields composed by several segments of different electron energies. RESULTS The advantageous EPID dosimetric properties already known for photons as reproducibility, linearity with dose, and dose rate were found to be identical for electron detection. The flood-field calibration method was proven to be effective and the EPID was capable to accurately reproduce the dose measured in water at 1.0 cm depth for 6 MeV, 1.3 cm for 9 MeV, and 1.5 cm for 12, 15, and 18 MeV. The deviations between the output factors measured with EPID and in water at these depths were within ±1.2% for all the energies with a mean deviation of 0.1%. The average gamma pass rate (criteria: 1.5%, 1.5 mm) for profile comparison between EPID and measurements in water was better than 99% for all the energies considered in this study. When comparing the reconstructed EPID 2D dose distributions at 1.5 cm depth to film measurements, the gamma pass rate (criteria: 2%, 2 mm) was better than 97% for all the tested cases. CONCLUSIONS This study demonstrates the high potential of the EPID for electron dosimetry, and in particular, confirms the possibility to use it as an efficient verification tool for MERT delivery.PURPOSE To investigate the dosimetric properties of an electronic portal imaging device (EPID) for electron beam detection and to evaluate its potential for quality assurance (QA) of modulated electron radiotherapy (MERT). METHODS A commercially available EPID was used to detect electron beams shaped by a photon multileaf collimator (MLC) at a source-surface distance of 70 cm. The fundamental dosimetric properties such as reproducibility, dose linearity, field size response, energy response, and saturation were investigated for electron beams. A new method to acquire the flood-field for the EPID calibration was tested. For validation purpose, profiles of open fields and various MLC fields (square and irregular) were measured with a diode in water and compared to the EPID measurements. Finally, in order to use the EPID for QA of MERT delivery, a method was developed to reconstruct EPID two-dimensional (2D) dose distributions in a water-equivalent depth of 1.5 cm. Comparisons were performed with film measurement for static and dynamic monoenergy fields as well as for multienergy fields composed by several segments of different electron energies. RESULTS The advantageous EPID dosimetric properties already known for photons as reproducibility, linearity with dose, and dose rate were found to be identical for electron detection. The flood-field calibration method was proven to be effective and the EPID was capable to accurately reproduce the dose measured in water at 1.0 cm depth for 6 MeV, 1.3 cm for 9 MeV, and 1.5 cm for 12, 15, and 18 MeV. The deviations between the output factors measured with EPID and in water at these depths were within ±1.2% for all the energies with a mean deviation of 0.1%. The average gamma pass rate (criteria: 1.5%, 1.5 mm) for profile comparison between EPID and measurements in water was better than 99% for all the energies considered in this study. When comparing the reconstructed EPID 2D dose distributions at 1.5 cm depth to film measurements, the gamma pass rate (criteria: 2%, 2 mm) was better than 97% for all the tested cases. CONCLUSIONS This study demonstrates the high potential of the EPID for electron dosimetry, and in particular, confirms the possibility to use it as an efficient verification tool for MERT delivery.

Low-dose 2.5 MV cone-beam computed tomography with thick CsI flat-panel imager

Most of the treatment units, both new and old models, are equipped with a megavoltage portal imager but its use for volumetric imaging is limited. This is mainly due to the poor image quality produced by the high‐energy treatment beam (>6 MV). A linac at our center is equipped with a prototype 2.5 MV imaging beam. This study evaluates the feasibility of low‐dose megavoltage cone‐beam imaging with the 2.5 MV beam and a thick cesium iodide detector, which is a high‐efficiency imager. Basic imaging properties such as spatial resolution and modulation transfer function were assessed for the 2.5 MV prototype imaging system. For image quality and imaging dose, a series of megavoltage cone‐beam scans were acquired for the head, thorax, and pelvis of an anthropomorphic phantom and were compared to kilovoltage cone‐beam and 6X megavoltage cone‐beam images. To demonstrate the advantage of MV imaging, a phantom with metallic inserts was scanned and the image quality was compared to CT and kilovoltage cone‐beam scans. With a lower energy beam and higher detector efficiency, the 2.5 MV imaging system generally yields better image quality than does the 6 MV imaging system with the conventional MV imager. In particular, with the anthropomorphic phantom studies, the contrast to noise of bone to tissue is generally improved in the 2.5 MV images compared to 6 MV. With an image quality sufficient for bony alignment, the imaging dose for 2.5 MV cone‐beam images is 2.4−3.4 MU compared to 26 MU in 6 MV cone‐beam scans for the head, thorax, and pelvis regions of the phantom. Unlike kilovoltage cone‐beam, the 2.5 MV imaging system does not suffer from high‐Z image artifacts. This can be very useful for treatment planning in cases where high‐Z prostheses are present. PACS number(s): 87.57.Q‐Most of the treatment units, both new and old models, are equipped with a megavoltage portal imager but its use for volumetric imaging is limited. This is mainly due to the poor image quality produced by the high-energy treatment beam (>6 MV). A linac at our center is equipped with a prototype 2.5 MV imaging beam. This study evaluates the feasibility of low-dose megavoltage cone-beam imaging with the 2.5 MV beam and a thick cesium iodide detector, which is a high-efficiency imager. Basic imaging properties such as spatial resolution and modulation transfer function were assessed for the 2.5 MV prototype imaging system. For image quality and imaging dose, a series of megavoltage cone-beam scans were acquired for the head, thorax, and pelvis of an anthropomorphic phantom and were compared to kilovoltage cone-beam and 6X megavoltage cone-beam images. To demonstrate the advantage of MV imaging, a phantom with metallic inserts was scanned and the image quality was compared to CT and kilovoltage cone-beam scans. With a lower energy beam and higher detector efficiency, the 2.5 MV imaging system generally yields better image quality than does the 6 MV imaging system with the conventional MV imager. In particular, with the anthropomorphic phantom studies, the contrast to noise of bone to tissue is generally improved in the 2.5 MV images compared to 6 MV. With an image quality sufficient for bony alignment, the imaging dose for 2.5 MV cone-beam images is 2.4-3.4 MU compared to 26 MU in 6 MV cone-beam scans for the head, thorax, and pelvis regions of the phantom. Unlike kilovoltage cone-beam, the 2.5 MV imaging system does not suffer from high-Z image artifacts. This can be very useful for treatment planning in cases where high-Z prostheses are present. PACS number(s): 87.57.Q.

A Monte Carlo study of the impact of phosphor optical properties on EPID imaging performance

We have developed a Monte Carlo computational model of a clinically employed electronic portal imaging device (EPID), and demonstrated the impact of phosphor optical properties on the imaging performance. The EPID model was built with Geant4 application for tomographic emission. Both radiative and optical transport were included in the model. Modulation transfer function (MTF), normalized noise-power spectrum times the incident x-ray fluence (qNNPS), and detective quantum efficiency (DQE) were calculated for simulated and measured data, and their agreement was quantified by the normalized root-mean-square error (NRMSE). MTF was computed using a 100 µm wide slit tilted by 1.5° and qNNPS was estimated using the Fujita-Lubberts-Swank method. DQE was calculated from MTF and qNNPS data. The NRMSE value was 0.0467 for MTF, 0.0217 for qNNPS, and 0.0885 for DQE, showing good agreement between measurement and simulation. Five major optical properties, phosphor grain size, phosphor thickness, phosphor refractive index, binder refractive index, and packing ratio were tested for their influence on the qNNPS, MTF, and DQE(0) of the model. Generally, the effect on the qNNPS is greater than MTF, and no impact on DQE(0), except from phosphor thickness, was observed. Multiple applications, such as imager design optimization and investigations of the dosimetric performance, are expected to benefit from the validated model.

Multi-layer imager design for mega-voltage spectral imaging

The architecture of multi-layer imagers (MLIs) can be exploited to provide megavoltage spectral imaging (MVSPI) for specific imaging tasks. In the current work, we investigated bone suppression and gold fiducial contrast enhancement as two clinical tasks which could be improved with spectral imaging. A method based on analytical calculations that enables rapid investigation of MLI component materials and thicknesses was developed and validated against Monte Carlo computations. The figure of merit for task-specific imaging performance was the contrast-to-noise ratio (CNR) of the gold fiducial when the CNR of bone was equal to zero after a weighted subtraction of the signals obtained from each MLI layer. Results demonstrated a sharp increase in the CNR of gold when the build-up component or scintillation materials and thicknesses were modified. The potential for low-cost, prompt implementation of specific modifications (e.g. composition of the build-up component) could accelerate clinical translation of MVSPI.

Spectral imaging using clinical megavoltage beams and a novel multi-layer imager

We assess the feasibility of clinical megavoltage (MV) spectral imaging for material and bone separation with a novel multi-layer imager (MLI) prototype. The MLI provides higher detective quantum efficiency and lower noise than conventional electronic portal imagers. Simulated experiments were performed using a validated Monte Carlo model of the MLI to estimate energy absorption and energy separation between the MLI components. Material separation was evaluated experimentally using solid water and aluminum (Al), copper (Cu) and gold (Au) for 2.5 MV, 6 MV and 6 MV flattening filter free (FFF) clinical photon beams. An anthropomorphic phantom with implanted gold fiducials was utilized to further demonstrate bone/gold separation. Weighted subtraction imaging was employed for material and bone separation. The weighting factor (w) was iteratively estimated, with the optimal w value determined by minimization of the relative signal difference ([Formula: see text]) and signal-difference-to-noise ratio (SDNR) between material (or bone) and the background. Energy separation between layers of the MLI was mainly the result of beam hardening between components with an average energy separation between 34 and 47 keV depending on the x-ray beam energy. The minimum average energy of the detected spectrum in the phosphor layer was 123 keV in the top layer of the MLI with the 2.5 MV beam. The w values that minimized [Formula: see text] and SDNR for Al, Cu and Au were 0.89, 0.76 and 0.64 for 2.5 MV; for 6 MV FFF, w was 0.98, 0.93 and 0.77 respectively. Bone suppression in the anthropomorphic phantom resulted in improved visibility of the gold fiducials with the 2.5 MV beam. Optimization of the MLI design is required to achieve optimal separation at clinical MV beam energies.