Modeling the temporal-spatial nature of the readout of an electronic portal imaging device (EPID).

PURPOSE: In real-time electronic portal imaging device (EPID) dosimetry applications where on-treatment measured transmission images are compared to an ideal predicted image, ideally a tight tolerance should be set on the quantitative image comparison in order to detect a wide variety of possible delivery errors. However, this is currently not possible due to the appearance of banding artifacts in individual frames of the measured EPID image sequences. The purpose of this work was to investigate simulating banding artifacts in our cine-EPID predicted image sequences to improve matching of individual image frames to the acquired image sequence. Increased sensitivity of this method to potential treatment delivery errors would represent an improvement in patient safety and treatment accuracy. METHODS: A circuit board was designed and built to capture the target current (TARG-I) and forward power signals produced by the linac to help model the discrete beam-formation process of the linac. To simulate the temporal-spatial nature of the EPID readout, a moving read out mask was applied with the timing of the application of the readout mask synchronized to the TARG-I pulses. Since identifying the timing of the first TARG-I pulse affected the location of the banding artifacts throughout the image sequence, and furthermore the first several TARG-I pulses at the beginning of "beam on" are not at full height yet (i.e., dose rate is ramping up), the forward-power signal was also used to assist in reliable detection of the first radiation pulse of the beam delivery. The predicted EPID cine-image sequence obtained using a comprehensive physics-based model was modified to incorporate the discrete nature of the EPID frame readout. This modified banding predicted EPID (MBP-EPID) image sequence was then compared to its corresponding measured EPID cine-image sequence on a frame-by-frame basis. The EPID was mounted on a Clinac 2100ix linac (Varian Medical Systems, Palo Alto, CA). The field size was set to 21.4  ×  28.6 cm2 with no MLC modulation, beam energy of 6 MV, dose rate of 600 MU/min, and 700 MU were delivered for each clockwise (CW) and counter-clockwise (CCW) arc. No phantoms were placed in the beam. RESULTS: The dose rate ramp up effect was observed at the beginning irradiations, and the identification and timing of the radiation pulses, even during the dose rate ramp up, were able to be quantified using the TARG-I and forward power signals. The approach of capturing individual dose pulses and synchronizing with the mask image applied to the original predicted EPID image sequence was demonstrated to model the actual EPID readout. The MBP-EPID image sequences closely reproduced the location and magnitude of the banding features observed in the acquired (i.e., measured) image sequence, for all test irradiations examined here. CONCLUSIONS: The banding artifacts observed in the measured EPID cine-frame sequences were reproduced in the predicted EPID cine-frames by simulating the discrete temporal-spatial nature of the EPID read out. The MBP-EPID images showed good agreement qualitatively to the corresponding measured EPID frame sequence of a simple square test field, without any phantom in the beam. This approach will lead to improved image comparison tolerances for real-time patient dosimetry applications.

In vivo dosimetry for lung radiotherapy including SBRT.

SBRT for lung cancer is being rapidly adopted as a treatment option in modern radiotherapy centres. This treatment is one of the most complex in common clinical use, requiring significant expertise and resources. It delivers a high dose per fraction (typically ∼6-30Gy/fraction) over few fractions. The complexity and high dose delivered in only a few fractions make powerful arguments for the application of in vivo dosimetry methods for these treatments to enhance patient safety. In vivo dosimetry is a group of techniques with a common objective - to estimate the dose delivered to the patient through a direct measurement of the treatment beam(s). In particular, methods employing an electronic portal imaging device have been intensely investigated over the past two decades. Treatment verification using in vivo dosimetry approaches has been shown to identify errors that would have been missed with other common quality assurance methods. With the addition of in vivo dosimetry to verify treatments, medical physicists and clinicians have a higher degree of confidence that the dose has been delivered to the patient as intended. In this review, the technical aspects and challenges of in vivo dosimetry for lung SBRT will be presented, focusing on transit dosimetry applications using electronic portal imaging devices (EPIDs). Currently available solutions will be discussed and published clinical experiences, which are very limited to date, will be highlighted.

Clinical Implementation of a Model-Based In Vivo Dose Verification System for Stereotactic Body Radiation Therapy-Volumetric Modulated Arc Therapy Treatments Using the Electronic Portal Imaging Device.

PURPOSE: To report findings from an in vivo dosimetry program implemented for all stereotactic body radiation therapy patients over a 31-month period and discuss the value and challenges of utilizing in vivo electronic portal imaging device (EPID) dosimetry clinically. METHODS AND MATERIALS: From December 2013 to July 2016, 117 stereotactic body radiation therapy-volumetric modulated arc therapy patients (100 lung, 15 spine, and 2 liver) underwent 602 EPID-based in vivo dose verification events. A developed model-based dose reconstruction algorithm calculates the 3-dimensional dose distribution to the patient by back-projecting the primary fluence measured by the EPID during treatment. The EPID frame-averaging was optimized in June 2015. For each treatment, a 3%/3-mm γ comparison between our EPID-derived dose and the Eclipse AcurosXB-predicted dose to the planning target volume (PTV) and the ≥20% isodose volume were performed. Alert levels were defined as γ pass rates <85% (lung and liver) and <80% (spine). Investigations were carried out for all fractions exceeding the alert level and were classified as follows: EPID-related, algorithmic, patient setup, anatomic change, or unknown/unidentified errors. RESULTS: The percentages of fractions exceeding the alert levels were 22.6% for lung before frame-average optimization and 8.0% for lung, 20.0% for spine, and 10.0% for liver after frame-average optimization. Overall, mean (± standard deviation) planning target volume γ pass rates were 90.7% ± 9.2%, 87.0% ± 9.3%, and 91.2% ± 3.4% for the lung, spine, and liver patients, respectively. CONCLUSIONS: Results from the clinical implementation of our model-based in vivo dose verification method using on-treatment EPID images is reported. The method is demonstrated to be valuable for routine clinical use for verifying delivered dose as well as for detecting errors.

Validation of a method for in vivo 3D dose reconstruction for IMRT and VMAT treatments using on-treatment EPID images and a model-based forward-calculation algorithm.

PURPOSE: Radiation treatments are trending toward delivering higher doses per fraction under stereotactic radiosurgery and hypofractionated treatment regimens. There is a need for accurate 3D in vivo patient dose verification using electronic portal imaging device (EPID) measurements. This work presents a model-based technique to compute full three-dimensional patient dose reconstructed from on-treatment EPID portal images (i.e., transmission images). METHODS: EPID dose is converted to incident fluence entering the patient using a series of steps which include converting measured EPID dose to fluence at the detector plane and then back-projecting the primary source component of the EPID fluence upstream of the patient. Incident fluence is then recombined with predicted extra-focal fluence and used to calculate 3D patient dose via a collapsed-cone convolution method. This method is implemented in an iterative manner, although in practice it provides accurate results in a single iteration. The robustness of the dose reconstruction technique is demonstrated with several simple slab phantom and nine anthropomorphic phantom cases. Prostate, head and neck, and lung treatments are all included as well as a range of delivery techniques including VMAT and dynamic intensity modulated radiation therapy (IMRT). RESULTS: Results indicate that the patient dose reconstruction algorithm compares well with treatment planning system computed doses for controlled test situations. For simple phantom and square field tests, agreement was excellent with a 2%/2 mm 3D chi pass rate ≥98.9%. On anthropomorphic phantoms, the 2%/2 mm 3D chi pass rates ranged from 79.9% to 99.9% in the planning target volume (PTV) region and 96.5% to 100% in the low dose region (>20% of prescription, excluding PTV and skin build-up region). CONCLUSIONS: An algorithm to reconstruct delivered patient 3D doses from EPID exit dosimetry measurements was presented. The method was applied to phantom and patient data sets, as well as for dynamic IMRT and VMAT delivery techniques. Results indicate that the EPID dose reconstruction algorithm presented in this work is suitable for clinical implementation.

An investigation of gantry angle data accuracy for cine-mode EPID images acquired during arc IMRT.

EPID images acquired in cine mode during arc therapy have inaccurate gantry angles recorded in their image headers. In this work, methods were developed to assess the accuracy of the gantry potentiometer for linear accelerators. As well, assessments of the accuracy of other, more accessible, sources of gantry angle information (i.e., treatment log files, analysis of EPID image headers) were investigated. The methods used in this study are generally applicable to any linear accelerator unit, and have been demonstrated here with Clinac/Trilogy systems. Gantry angle data were simultaneously acquired using three methods: i) a direct gantry potentiometer measurement, ii) an incremental rotary encoder, and iii) a custom-made radiographic gantry-angle phantom which produced unique wire intersections as a function of gantry angle. All methods were compared to gantry angle data from the EPID image header and the linac MLC DynaLog file. The encoder and gantry-angle phantom were used to validate the accuracy of the linac's potentiometer. The EPID image header gantry angles and the DynaLog file gantry angles were compared to the potentiometer. The encoder and gantry-angle phantom mean angle differences with the potentiometer were 0.13° ± 0.14° and 0.10°± 0.30°, respectively. The EPID image header angles analyzed in this study were within ± 1° of the potentiometer angles only 35% of the time. In some cases, EPID image header gantry angles disagreed by as much as 3° with the potentiometer. A time delay in frame acquisition was determined using the continuous acquisition mode of the EPID. After correcting for this time delay, 75% of the header angles, on average, were within ± 1° of the true gantry angle, compared to an average of only 35% without the correction. Applying a boxcar smoothing filter to the corrected gantry angles further improved the accuracy of the header-derived gantry angles to within ± 1° for almost all images (99.4%). An angle accuracy of 0.11° ± 0.04° was determined using a point-by-point comparison of the gantry angle data in the MLC DynaLog file and the potentiometer data. These simple correction methods can be easily applied to individual treatment EPID images in order to more accurately define the gantry angle.

Gantry angle determination during arc IMRT: evaluation of a simple EPID-based technique and two commercial inclinometers.

The increasing popularity of intensity-modulated arc therapy (IMAT) treatments requires specifically designed linac quality assurance (QA) programs. Gantry angle is one of the parameters that has a major effect on the outcome of IMAT treatments since dose reconstruction for patient-specific QA relies on the gantry angle; therefore, it is essential to ensure its accuracy for correct delivery of the prescribed dose. In this study, a simple measurement method and algorithm are presented for QA of gantry angle measurements based on integrated EPID images acquired at distinct gantry angles and cine EPID images during an entire 360° arc. A comprehensive study was carried out to evaluate this method, as well as to evaluate two commercially available inclinometers (NG360 and IBA GAS supplied in conjunction with popular array dosimeters Delta4 and MatriXXEvolution, respectively) by comparing their simultaneous angle measurement results with the linac potentiometer readouts at five gantry speeds. In all tested measurement systems, the average differences with the reference angle data were less than 0.3° in static mode. In arc mode, at all tested gantry speeds the average difference was less than 0.1° for the IBA GAS and the proposed EPID-based method, and 0.6° for the NG360 after correction for the inherent systematic time delay of the inclinometer. The gantry rotation speed measured by the three independent systems had an average deviation of about 0.01°/s from the nominal gantry speed.

Detection and correction for EPID and gantry sag during arc delivery using cine EPID imaging.

PURPOSE: Electronic portal imaging devices (EPIDs) have been studied and used for pretreatment and in-vivo dosimetry applications for many years. The application of EPIDs for dosimetry in arc treatments requires accurate characterization of the mechanical sag of the EPID and gantry during rotation. Several studies have investigated the effects of gravity on the sag of these systems but each have limitations. In this study, an easy experiment setup and accurate algorithm have been introduced to characterize and correct for the effect of EPID and gantry sag during arc delivery. METHODS: Three metallic ball bearings were used as markers in the beam: two of them fixed to the gantry head and the third positioned at the isocenter. EPID images were acquired during a 360° gantry rotation in cine imaging mode. The markers were tracked in EPID images and a robust in-house developed MATLAB code was used to analyse the images and find the EPID sag in three directions as well as the EPID + gantry sag by comparison to the reference gantry zero image. The algorithm results were then tested against independent methods. The method was applied to compare the effect in clockwise and counter clockwise gantry rotations and different source-to-detector distances (SDDs). The results were monitored for one linear accelerator over a course of 15 months and six other linear-accelerators from two treatment centers were also investigated using this method. The generalized shift patterns were derived from the data and used in an image registration algorithm to correct for the effect of the mechanical sag in the system. The Gamma evaluation (3%, 3 mm) technique was used to investigate the improvement in alignment of cine EPID images of a fixed field, by comparing both individual images and the sum of images in a series with the reference gantry zero image. RESULTS: The mechanical sag during gantry rotation was dependent on the gantry angle and was larger in the in-plane direction, although the patterns were not identical for various linear-accelerators. The reproducibility of measurements was within 0.2 mm over a period of 15 months. The direction of gantry rotation and SDD did not affect the results by more than 0.3 mm. Results of independent tests agreed with the algorithm within the accuracy of the measurement tools. When comparing summed images, the percentage of points with Gamma index <1 increased from 85.4% to 94.1% after correcting for the EPID sag, and to 99.3% after correction for gantry + EPID sag. CONCLUSIONS: The measurement method and algorithms introduced in this study use cine-images, are highly accurate, simple, fast, and reproducible. It tests all gantry angles and provides a suitable automatic analysis and correction tool to improve EPID dosimetry and perform comprehensive linac QA for arc treatments.