Comparison of ghosting effects for three commercial a-Si EPIDs

AAPM Task Group Report 307: Use of EPIDs for Patient-Specific IMRT and VMAT QA

PURPOSE

Electronic portal imaging devices (EPIDs) have been widely utilized for patient-specific quality assurance (PSQA) and their use for transit dosimetry applications is emerging. Yet there are no specific guidelines on the potential uses, limitations, and correct utilization of EPIDs for these purposes. The American Association of Physicists in Medicine (AAPM) Task Group 307 (TG-307) provides a comprehensive review of the physics, modeling, algorithms and clinical experience with EPID-based pre-treatment and transit dosimetry techniques. This review also includes the limitations and challenges in the clinical implementation of EPIDs, including recommendations for commissioning, calibration and validation, routine QA, tolerance levels for gamma analysis and risk-based analysis.

METHODS

Characteristics of the currently available EPID systems and EPID-based PSQA techniques are reviewed. The details of the physics, modeling, and algorithms for both pre-treatment and transit dosimetry methods are discussed, including clinical experience with different EPID dosimetry systems. Commissioning, calibration, and validation, tolerance levels and recommended tests, are reviewed, and analyzed. Risk-based analysis for EPID dosimetry is also addressed.

RESULTS

Clinical experience, commissioning methods and tolerances for EPID-based PSQA system are described for pre-treatment and transit dosimetry applications. The sensitivity, specificity, and clinical results for EPID dosimetry techniques are presented as well as examples of patient-related and machine-related error detection by these dosimetry solutions. Limitations and challenges in clinical implementation of EPIDs for dosimetric purposes are discussed and acceptance and rejection criteria are outlined. Potential causes of and evaluations of pre-treatment and transit dosimetry failures are discussed. Guidelines and recommendations developed in this report are based on the extensive published data on EPID QA along with the clinical experience of the TG-307 members.

CONCLUSION

TG-307 focused on the commercially available EPID-based dosimetric tools and provides guidance for medical physicists in the clinical implementation of EPID-based patient-specific pre-treatment and transit dosimetry QA solutions including intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) treatments.

Characterization of Gas Electron Multiplier-based detector for external beam radiation therapy dosimetry

PURPOSE

Electronic portal imaging devices (EPIDs) are commonly installed on modern linear accelerators (LINACs) and are convenient for imaging and, potentially, dosimetry. However, owing to their construction with metal and scintillating layers of high atomic number, they exhibit nonwater-equivalent response and oversensitivity to low-energy photons. Therefore, EPIDs are not ideal for dosimetry purposes. Additionally, nonlinearities due to the combined use of scintillators and photodiodes have been reported. Here, an EPID which employs a variable gain Gas Electron Multiplier (GEM) and direct detection of electrons is introduced. To investigate its dosimetric performance, measurements characterizing the novel EPID are performed and compared with measurements from ionization chambers and conventional EPIDs.

METHODS

Linearity, dose rate dependence, field size dependence, off-axis response, and transmission response were measured for all available energy settings (6, 10, 6 MV Flattening Filter Free (FFF) and 10 MVFFF) using three different detector gain settings. Additionally, an evaluation of the ghosting and image lag of the panel was completed. Reference ionization chamber measurements were performed for the off-axis and transmission response and existing data for conventional EPIDs and ionization chambers from equivalent measurements were used for comparison of the field size dependence. Elsewhere, values from the linac monitoring chambers were used.

RESULTS

In the range from 10 to 1000 Monitor Units (MU), the detector was linear within 1% for all combinations of gain settings and energies. The dose rate dependence was also within 1% for all energies and for two out of three gain settings. Regarding field size dependency, the ratio of ionization chamber and panel values was 0.94 and 0.98 for the conventional EPID and GEMini respectively, at 20 × 20 cm² and 10 MV. For 6 MV, 6 MVFFF, and 10MVFFF these ratios were 0.97, 0.98, and 0.99 for the GEMini, and 0.95, 0.97, and 0.97 for the conventional EPID. Similar performance between the GEMini and conventional EPID is observed for field sizes smaller than 10 × 10 cm² . The transmission response was within 5% for all energies for thicknesses up to 30 cm, compared to 10-20% for a conventional EPID. The off-axis response for shifts up to 16 cm was within 1% and 3% for 6 MV and 10 MV, with and without phantom. The rise and fall of the signal from the detector correspond well to monitor chamber measurements indicating little ghosting and image lag, regardless of gain setting.

CONCLUSION

The GEM EPID exhibits dose rate dependence and linearity within 1%, and negligible ghosting and image lag. In this regard, it performs particularly well using 50 and 250 V of gain, and either could be chosen. For higher sensitivity, 250 V is the recommended base gain setting, although other applications may warrant different gains. For most tests performed in this study, the GEM EPID demonstrates a more water-equivalent response than conventional EPIDs making GEMs a viable technology for dosimetry in radiation therapy.

Dosimetric performance of continuous EPID imaging in stereotactic treatment conditions

PURPOSE

This study aims at investigating the dosimetric characteristics of a Varian aS1000 EPID, focusing on its continuous acquisition mode under the challenging conditions that can be met in stereotactic radiotherapy verification.

METHODS

An aS1000 EPID installed on a Varian TrueBeamSTx was irradiated with 6 and 10 MV unflattened and flattened photon beams. In order to avoid detector saturation, the source-to-detector distance (SDD) was set to 150 or 180 cm depending on the dose rate. EPID image sets were acquired in continuous mode (CM) and also in the commonly used integrated mode (IM) for comparison, to evaluate dose linearity (including dose rate dependence), repeatability, reproducibility, stability, ghosting effect and field size dependence.

RESULTS

CM response linearity was found to be within 0.8% of IM and independent of dose rate. Response repeatability was slightly better for IM and FF beams, being in all cases within 0.9%. Reproducibility was within 0.6% for both modes and all beam qualities. Response stability between continuous frames varied within 1% for dynamic and static irradiations and for all the beam qualities, showing its independence from these parameters. Ghosting effect was not significant, being comparable to signal variations between continuous frames (±1%). Field size dependence in both modes agreed within 1%.

CONCLUSIONS

The dosimetric response of the aS1000 EPID in CM with FFF beams and high dose rates is comparable to that in IM and for flattened beams provided that the appropriate SDD is used. aS1000 EPID in continuous acquisition mode is therefore suitable for stereotactic applications.

Remote dosimetric auditing of clinical trials: The need for vendor specific models to convert images to dose

INTRODUCTION

A previous pilot study has demonstrated the feasibility of a novel image-based approach for remote dosimetric auditing of clinical trials. The approach uses a model to convert in-air acquired intensity modulated radiotherapy (IMRT) images to delivered dose inside a virtual phantom. The model was developed using images from an electronic portal imaging device (EPID) on a Varian linear accelerator. It was tuned using beam profiles and field size factors (FSFs) of a series of square fields measured in water tank. This work investigates the need for vendor specific conversion models for image-based auditing. The EPID measured profile and FSF data for Varian (vendor 1) and Elekta (vendor 2) systems are compared along with the performance of the existing Varian model (VM) and a new Elekta model (EM) for a series of audit IMRT fields measured on vendor 2 systems.

MATERIALS AND METHODS

The EPID measured beam profile and FSF data were studied for the two vendors to quantify and understand their relevant dosimetric differences. Then, an EM was developed converting EPID to dose in the virtual water phantom using a vendor 2 water tank data and images from corresponding EPID. The VM and EM were compared for predicting vendor 2 measured dose in water tank. Then, the performance of the new EM was compared to the VM for auditing of 54 IMRT fields from four vendor 2 facilities. Statistical significance of using vendor specific models was determined.

RESULTS

Observed dosimetry differences between the two vendors suggested developing an EM would be beneficial. The EM performed better than VM for vendor 2 square and IMRT fields. The IMRT audit gamma pass rates were (99.8 ± 0.5)%, (98.6 ± 2.3)% and (97.0 ± 3.0)% at respectively 3%/3 mm, 3%/2 mm and 2%/2 mm with improvements at most fields compared with using the VM. For the pilot audit, the difference between gamma results of the two vendors was reduced when using vendor specific models (VM: P < 0.0001, vendor specific models: P = 0.0025).

CONCLUSION

A new model was derived to convert images from vendor 2 EPIDs to dose for remote auditing vendor 2 deliveries. Using vendor specific models is recommended to remotely audit systems from different vendors, however, the improvements found were not major.

Characterization of the a-Si EPID in the unity MR-linac for dosimetric applications

Electronic portal imaging devices (EPIDs) are frequently used in external beam radiation therapy for dose verification purposes. The aim of this study was to investigate the dose-response characteristics of the EPID in the Unity MR-linac (Elekta AB, Stockholm, Sweden) relevant for dosimetric applications under clinical conditions. EPID images and ionization chamber (IC) measurements were used to study the effects of the magnetic field, the scatter generated in the MR housing reaching the EPID, and inhomogeneous attenuation from the MR housing. Dose linearity and dose rate dependencies were also determined. The magnetic field strength at EPID level did not exceed 10 mT, and dose linearity and dose rate dependencies proved to be comparable to that on a conventional linac. Profiles of fields, delivered with and without the magnetic field, were indistinguishable. The EPID center had an offset of 5.6 cm in the longitudinal direction, compared to the beam central axis, meaning that large fields in this direction will partially fall outside the detector area and not be suitable for verification. Beam attenuation by the MRI scanner and the table is gantry angle dependent, presenting a minimum attenuation of 67% relative to the 90° measurement. Repeatability, observed over two months, was within 0.5% (1 SD). In order to use the EPID for dosimetric applications in the MR-linac, challenges related to the EPID position, scatter from the MR housing, and the inhomogeneous, gantry angle-dependent attenuation of the beam will need to be solved.

Online pretreatment verification of high-dose rate brachytherapy using an imaging panel

Brachytherapy is employed to treat a wide variety of cancers. However, an accurate treatment verification method is currently not available. This study describes a pre-treatment verification system that uses an imaging panel (IP) to verify important aspects of the treatment plan. A detailed modelling of the IP was only possible with an extensive calibration performed using a robotic arm. Irradiations were performed with a high dose rate (HDR) ¹⁹²Ir source within a water phantom. An empirical fit was applied to measure the distance between the source and the detector so 3D Cartesian coordinates of the dwell positions can be obtained using a single panel. The IP acquires 7.14 fps to verify the dwell times, dwell positions and air kerma strength (Sk). A gynecological applicator was used to create a treatment plan that was registered with a CT image of the water phantom used during the experiments for verification purposes. Errors (shifts, exchanged connections and wrong dwell times) were simulated to verify the proposed verification system. Cartesian source positions (panel measurement plane) have a standard deviation of about 0.02 cm. The measured distance between the source and the panel (z-coordinate) have a standard deviation up to 0.16 cm and maximum absolute error of  ≈0.6 cm if the signal is close to sensitive limit of the panel. The average response of the panel is very linear with Sk. Therefore, Sk measurements can be performed with relatively small errors. The measured dwell times show a maximum error of 0.2 s which is consistent with the acquisition rate of the panel. All simulated errors were clearly identified by the proposed system. The use of IPs is not common in brachytherapy, however, it provides considerable advantages. It was demonstrated that the IP can accurately measure Sk, dwell times and dwell positions.

Normalize the response of EPID in pursuit of linear accelerator dosimetry standardization

Normalize the response of electronic portal imaging device (EPID) is the first step toward an EPID‐based standardization of Linear Accelerator (linac) dosimetry quality assurance. In this study, we described an approach to generate two‐dimensional (2D) pixel sensitivity maps (PSM) for EPIDs response normalization utilizing an alternative beam and dark‐field (ABDF) image acquisition technique and large overlapping field irradiations. The automated image acquisition was performed by XML‐controlled machine operation and the PSM was generated based on a recursive calculation algorithm for Varian linacs equipped with aS1000 and aS1200 imager panels. Cross‐comparisons of normalized beam profiles and 1.5%/1.5 mm 1D Gamma analysis was adopted to quantify the improvement of beam profile matching before and after PSM corrections. PSMs were derived for both photon (6, 10, 15 MV) and electron (6, 20 MeV) beams via proposed method. The PSM‐corrected images reproduced a horn‐shaped profile for photon beams and a relative uniform profiles for electrons. For dosimetrically matched linacs equipped with aS1000 panels, PSM‐corrected images showed increased 1D‐Gamma passing rates for all energies, with an average 10.5% improvement for crossline and 37% for inline beam profiles. Similar improvements in the phantom study were observed with a maximum improvement of 32% for 15 MV and 22% for 20 MeV. The PSM value showed no significant change for all energies over a 3‐month period. In conclusion, the proposed approach correct EPID response for both aS1000 and aS1200 panels. This strategy enables the possibility to standardize linac dosimetry QA and to benchmark linac performance utilizing EPID as the common detector.

EPID-based dosimetry to verify IMRT planar dose distribution for the aS1200 EPID and FFF beams

We proposed to perform a basic dosimetry commissioning on a new imager system, the Varian aS1200 electronic portal imaging device (EPID) and TrueBeam 2.0 linear accelerator for flattened (FF) and flattening filter‐free (FFF) beams, then to develop an image‐based quality assurance (QA) model for verification of the system delivery accuracy for intensity‐modulated radiation therapy (IMRT) treatments. For dosimetry testing, linearity of dose response with MU, imager lag, and effectiveness of backscatter shielding were investigated. Then, an image‐based model was developed to convert images to planar dose onto a virtual water phantom. The model parameters were identified using energy fluence of the Acuros treatment planning system (TPS) and, reference dose profiles and output factors measured at depths of 5, 10, 15, and 20 cm in water phantom for square fields. To validate the model, its calculated dose was compared to measured dose from MapCHECK 2 diode arrays for 36 IMRT fields at 10 cm depth delivered with 6X, 6XFFF, 10X, and 10XFFF energies. An in‐house gamma function was used to compare planar doses pixel‐by‐pixel. Finally, the method was applied to the same IMRT fields to verify their pretreatment delivery dose compared with Eclipse TPS dose. For the EPID commissioning, dose linearity was within 0.4% above 5 MU and ∼1% above 2 MU, measured lag was smaller than the previous EPIDs, and profile symmetry was improved. The model was validated with mean gamma pass rates (standard deviation) of 99.0% (0.4%), 99.5% (0.6%), 99.3% (0.4%), and 98.0% (0.8%) at 3%/3 mm for respectively 6X, 6XFFF, 10X, and 10XFFF beams. Using the same comparison criteria, the beam deliveries were verified with mean pass rates of 100% (0.0%), 99.6% (0.3%), 99.9% (0.1%), and 98.7% (1.4%). Improvements were observed in dosimetric response of the aS1200 imager compared to previous EPID models, and the model was successfully developed for the new system and delivery energies of 6 and 10 MV, FF, and FFF modes.

PACS number(s): 87.53.Oq, 87.53.Xd

On flattening filter-free portal dosimetry

Varian introduced (in 2010) the option of removing the flattening filter (FF) in their C‐Arm linacs for intensity‐modulated treatments. This mode, called flattening filter‐free (FFF), offers the advantage of a greater dose rate. Varian's “Portal Dosimetry” is an electronic portal imager device (EPID)‐based tool for IMRT verification. This tool lacks the capability of verifying flattening filter‐free (FFF) modes due to saturation and lack of an image prediction algorithm. (Note: the latest versions of this software and EPID correct these issues.) The objective of the present study is to research the feasibility of said verifications (with the older versions of the software and EPID). By placing the EPID at a greater distance, the images can be acquired without saturation, yielding a linearity similar to the flattened mode. For the image prediction, a method was optimized based on the clinically used algorithm (analytical anisotropic algorithm (AAA)) over a homogeneous phantom. The depth inside the phantom and its electronic density were tailored. An application was developed to allow the conversion of a dose plane (in DICOM format) to Varian's custom format for Portal Dosimetry. The proposed method was used for the verification of test and clinical fields for the three qualities used in our institution for IMRT: 6X, 6FFF and 10FFF. The method developed yielded a positive verification (more than 95% of the points pass a 2%/2 mm gamma) for both the clinical and test fields. This method was also capable of “predicting” static and wedged fields. A workflow for the verification of FFF fields was developed. This method relies on the clinical algorithm used for dose calculation and is able to verify the FFF modes, as well as being useful for machine quality assurance. The procedure described does not require new hardware. This method could be used as a verification of Varian's Portal Dose Image Prediction.

PACS number(s): 87.53.Kn, 87.55.T‐, 87.56.bd, 87.59.‐e

Development of an accurate EPID-based output measurement and dosimetric verification tool for electron beam therapy

PURPOSE

To develop an efficient and robust tool for output measurement and absolute dose verification of electron beam therapy by using a high spatial-resolution and high frame-rate amorphous silicon flat panel electronic portal imaging device (EPID).

METHODS

The dosimetric characteristics of the EPID, including saturation, linearity, and ghosting effect, were first investigated on a Varian Clinac 21EX accelerator. The response kernels of the individual pixels of the EPID to all available electron energies (6, 9, 12, 16, and 20 MeV) were calculated by using Monte Carlo (MC) simulations, which formed the basis to deconvolve an EPID raw images to the incident electron fluence map. The two-dimensional (2D) dose distribution at reference depths in water was obtained by using the constructed fluence map with a MC simulated pencil beam kernel with consideration of the geometric and structural information of the EPID. Output factor measurements were carried out with the EPID at a nominal source-surface distance of 100 cm for 2 × 2, 3 × 3, 6 × 6, 10 × 10, and 15 × 15 cm(2) fields for all available electron energies, and the results were compared with that measured in a solid water phantom using film and a Farmer-type ion chamber. The dose distributions at a reference depth specific to each energy and the flatness and symmetry of the 10 × 10 cm(2) electron beam were also measured using EPID, and the results were compared with ion chamber array and water scan measurements. Finally, three patient cases with various field sizes and irregular cutout shapes were also investigated.

RESULTS

EPID-measured dose changed linearly with the monitor units and showed little ghosting effect for dose rate up to 600 MU/min. The flatness and symmetry measured with the EPID were found to be consistent with ion chamber array and water scan measurements. The EPID-measured output factors for standard square fields of 2 × 2, 3 × 3, 6 × 6, 10 × 10, 15 × 15 cm(2) agreed with film and ion chamber measurements. The average discrepancy between EPID and ion chamber/film measurements was 0.81% ± 0.60% (SD) and 1.34% ± 0.75%, respectively. For the three clinical cases, the difference in output between the EPID- and ion chamber array measured values was found to be 1.13% ± 0.11%, 0.54% ± 0.10%, and 0.74% ± 0.11%, respectively. Furthermore, the γ-index analysis showed an excellent agreement between the EPID- and ion chamber array measured dose distributions: 100% of the pixels passed the criteria of 3%/3 mm. When the γ-index was set to be 2%/2 mm, the pass rate was found to be 99.0% ± 0.07%, 98.2% ± 0.14%, and 100% for the three cases.

CONCLUSIONS

The EPID dosimetry system developed in this work provides an accurate and reliable tool for routine output measurement and dosimetric verification of electron beam therapy. Coupled with its portability and ease of use, the proposed system promises to replace the current film-based approach for fast and reliable assessment of small and irregular electron field dosimetry.

Use of electronic portal imaging devices for electron treatment verification

This study aims to help broaden the use of electronic portal imaging devices (EPIDs) for pre-treatment patient positioning verification, from photon-beam radiotherapy to photon- and electron-beam radiotherapy, by proposing and testing a method for acquiring clinically-useful EPID images of patient anatomy using electron beams, with a view to enabling and encouraging further research in this area. EPID images used in this study were acquired using all available beams from a linac configured to deliver electron beams with nominal energies of 6, 9, 12, 16 and 20 MeV, as well as photon beams with nominal energies of 6 and 10 MV. A widely-available heterogeneous, approximately-humanoid, thorax phantom was used, to provide an indication of the contrast and noise produced when imaging different types of tissue with comparatively realistic thicknesses. The acquired images were automatically calibrated, corrected for the effects of variations in the sensitivity of individual photodiodes, using a flood field image. For electron beam imaging, flood field EPID calibration images were acquired with and without the placement of blocks of water-equivalent plastic (with thicknesses approximately equal to the practical range of electrons in the plastic) placed upstream of the EPID, to filter out the primary electron beam, leaving only the bremsstrahlung photon signal. While the electron beam images acquired using a standard (unfiltered) flood field calibration were observed to be noisy and difficult to interpret, the electron beam images acquired using the filtered flood field calibration showed tissues and bony anatomy with levels of contrast and noise that were similar to the contrast and noise levels seen in the clinically acceptable photon beam EPID images. The best electron beam imaging results (highest contrast, signal-to-noise and contrast-to-noise ratios) were achieved when the images were acquired using the higher energy electron beams (16 and 20 MeV) when the EPID was calibrated using an intermediate (12 MeV) electron beam energy. These results demonstrate the feasibility of acquiring clinically-useful EPID images of patient anatomy using electron beams and suggest important avenues for future investigation, thus enabling and encouraging further research in this area. There is manifest potential for the EPID imaging method proposed in this work to lead to the clinical use of electron beam imaging for geometric verification of electron treatments in the future.

Dosimetric properties of fluoroscopic EPID for transit dosimetry

The aim of this work was to evaluate dose response of fluoroscopic EPID for transit dosimetry applications. Properties studied included warm up time, build-up thickness evaluation, dose history, linearity, stability, and short and long-term reproducibility of EPID response, as well as field size dependence.

Pixel value matrices of electronic portal images in DICOM format were analysed in central and 8 off axis points using customised written codes in Matlab. In order to do this, nine 26×26 pixel matrices were selected as regions of interest, the regions represented by these arrays were 1×1 and 0·65×0·65 cm2 at the EPID and isocentre level, respectively.

Necessary warm up time for stable operation of EPID is 30 minutes, and there is no need for extra build-up layer to increase the dose response. Linearity tests indicate charged coupled device camera of EPID saturates at 50 cGy level, and does not have linear relationship with dose. Reproducibility and stability of the measurements were excellent and the detector showed same signal with a maximum deviation of <0·3% both in short and long terms. Results of dosimetric evaluation have shown the TheraView fluoroscopic EPID can be used for transit dosimetry purposes.

The impact of cine EPID image acquisition frame rate on markerless soft-tissue tracking

PURPOSE

Although reduction of the cine electronic portal imaging device (EPID) acquisition frame rate through multiple frame averaging may reduce hardware memory burden and decrease image noise, it can hinder the continuity of soft-tissue motion leading to poor autotracking results. The impact of motion blurring and image noise on the tracking performance was investigated.

METHODS

Phantom and patient images were acquired at a frame rate of 12.87 Hz with an amorphous silicon portal imager (AS1000, Varian Medical Systems, Palo Alto, CA). The maximum frame rate of 12.87 Hz is imposed by the EPID. Low frame rate images were obtained by continuous frame averaging. A previously validated tracking algorithm was employed for autotracking. The difference between the programmed and autotracked positions of a Las Vegas phantom moving in the superior-inferior direction defined the tracking error (δ). Motion blurring was assessed by measuring the area change of the circle with the greatest depth. Additionally, lung tumors on 1747 frames acquired at 11 field angles from four radiotherapy patients are manually and automatically tracked with varying frame averaging. δ was defined by the position difference of the two tracking methods. Image noise was defined as the standard deviation of the background intensity. Motion blurring and image noise are correlated with δ using Pearson correlation coefficient (R).

RESULTS

For both phantom and patient studies, the autotracking errors increased at frame rates lower than 4.29 Hz. Above 4.29 Hz, changes in errors were negligible withδ < 1.60 mm. Motion blurring and image noise were observed to increase and decrease with frame averaging, respectively. Motion blurring and tracking errors were significantly correlated for the phantom (R = 0.94) and patient studies (R = 0.72). Moderate to poor correlation was found between image noise and tracking error with R -0.58 and -0.19 for both studies, respectively.

CONCLUSIONS

Cine EPID image acquisition at the frame rate of at least 4.29 Hz is recommended. Motion blurring in the images with frame rates below 4.29 Hz can significantly reduce the accuracy of autotracking.

Dose calibration of EPIDs for segmented IMRT dosimetry

The purpose of this study was to investigate the dose response of amorphous silicon (a‐Si) electronic portal imaging devices (EPIDs) under different acquisition settings for both open jaw defined fields and segmented intensity‐modulated radiation therapy (IMRT) fields. Four different EPIDs were used. Two Siemens and one Elekta plus a standalone Perkin Elmer research EPID. Each was operated with different acquisition systems and settings. Dose response linearity was measured for open static jaw defined fields and ‘simple’ segmented IMRT fields for a range of equipment and system settings. Six ‘simple’ segmented IMRT fields were used. The segments of each IMRT field were fixed at 10×10cm2 field size with equal MU per segment, each field having a total of 20 MU. Simultaneous measurements with an ionization chamber array (ICA) and EPID were performed to separate beam and detector response characteristics. Three different pixel calibration methods were demonstrated and compared for an example ‘clinical IMRT field’. The dose response with the Elekta EPID for ‘simple’ segmented IMRT fields versus static fields agreed to within 2.5% for monitor unit (MU)≥2. The dose response for the Siemens systems was difficult to interpret due to the poor reproducibility for segmented delivery, at MU≤5, which was not observed with the standalone research EPID nor ICA on the same machine. The dose response measured under different acquisition settings and different linac/EPID combinations matched closely (≤1%), except for the Siemens EPID. Clinical IMRT EPID dosimetry implemented with the different pixel‐to‐dose calibration methods indicated that calibration at 20 MU provides equivalent results to implementing a ghosting correction model. The nonlinear dose response was consistent across both clinical EPIDs and the standalone research EPID, with the exception of the poor reproducibility seen with Siemens EPID images of IMRT fields. The nonlinear dose response was relatively insensitive to acquisition settings and appears to be primarily due to gain ghosting effects. No additional ghosting correction factor is necessary when the pixel‐to‐dose calibration factor at small MU calibration method is used.

PACS numbers: 87.53.Bn, 87.55.Qr, 87.56.Fc, 87.57.uq

Time dependent pre-treatment EPID dosimetry for standard and FFF VMAT

Methods to calibrate Megavoltage electronic portal imaging devices (EPIDs) for dosimetry have been previously documented for dynamic treatments such as intensity modulated radiotherapy (IMRT) using flattened beams and typically using integrated fields. While these methods verify the accumulated field shape and dose, the dose rate and differential fields remain unverified. The aim of this work is to provide an accurate calibration model for time dependent pre-treatment dose verification using amorphous silicon (a-Si) EPIDs in volumetric modulated arc therapy (VMAT) for both flattened and flattening filter free (FFF) beams. A general calibration model was created using a Varian TrueBeam accelerator, equipped with an aS1000 EPID, for each photon spectrum 6 MV, 10 MV, 6 MV-FFF, 10 MV-FFF. As planned VMAT treatments use control points (CPs) for optimization, measured images are separated into corresponding time intervals for direct comparison with predictions. The accuracy of the calibration model was determined for a range of treatment conditions. Measured and predicted CP dose images were compared using a time dependent gamma evaluation using criteria (3%, 3 mm, 0.5 sec). Time dependent pre-treatment dose verification is possible without an additional measurement device or phantom, using the on-board EPID. Sufficient data is present in trajectory log files and EPID frame headers to reliably synchronize and resample portal images. For the VMAT plans tested, significantly more deviation is observed when analysed in a time dependent manner for FFF and non-FFF plans than when analysed using only the integrated field. We show EPID-based pre-treatment dose verification can be performed on a CP basis for VMAT plans. This model can measure pre-treatment doses for both flattened and unflattened beams in a time dependent manner which highlights deviations that are missed in integrated field verifications.

Image quality improvements of electronic portal imaging devices by multi-level gain calibration and temperature correction

Amorphous silicon (aSi:H) flat panel detectors are prevalent in radiotherapy for megavoltage imaging tasks. Any clinical and dosimetrical application requires a well-defined dose response of the system to achieve meaningful results. Due to radiation damages, panels deteriorate and the linearity of pixel response to dose as well as the stability with regard to changing operating temperatures get worse with time. Using a single level gain correction can lead to an error of about 23% when irradiating a flood field image with 100 MU min(-1) on an old detector. A multi-level gain (MLG) correction is introduced, emending the nonlinearities and subpanel-related artifacts caused by insufficient radiation hardness of amplifiers in the read-out electronics. With rising temperature, offset values typically increase (up to 300 gray values) while the response at higher dose values per frame remain constant for a majority of pixels. To account for temperature-related image artifacts, two additional temperature correction methods have been developed. MLG in combination with temperature corrections can re-establish the aSi:H image quality to the performance required by reliable medical verification tools. Furthermore, the life span and recalibration intervals of these costly devices can be prolonged decisively.

2D lag and signal nonlinearity correction in an amorphous silicon EPID and their impact on pretreatment dosimetric verification

PURPOSE

This investigation provides measurements of signal lag and nonlinearity separately for the Varian aS500 electronic portal imaging device (EPID), and an algorithm to correct for these effects in 2D; their potential impact on intensity modulated radiation therapy (IMRT) verification is also investigated. The authors quantify lag, as a function of both delivered monitor units (MU) and time, by using a range of MUs delivered at a clinically used rate of 400 MU∕min. Explicit cumulative lag curves are thus determined for a range of MUs and times between the end of irradiation and the end of image acquisition. Signal nonlinearity is also investigated as a function of total MUs delivered. The family of cumulative lag curves and signal nonlinearity are then used to determine their effects on dynamic multileaf collimator (MLC) (IMRT) deliveries, and to correct for theses effects in 2D.

METHODS

Images acquired with an aS500 EPID and Varis Portal-Vision software were used to quantify detector lag and signal-nonlinearity. For the signal lag investigation, Portal-Vision's service monitor was used to acquire EPID images at a rate of 8 frames/s. The images were acquired during irradiation and 66 s thereafter, by inhibiting the M-holdoff-In signal of the Linac for a range of 4.5-198.5 MUs. Relative cumulative lag was calculated by integrating the EPID signal for a time after beam-off, and normalizing this to the integrated EPID signal accumulated during radiation. Signal nonlinearity was studied by acquiring 10 × 10 cm(2) open-field EPID images in "integrated image" mode for a range of 2-500 MUs, and normalized to the 100 MU case. All data were incorporated into in-house written software to create a 2D correction map for these effects, using the field's MLC file and a field-specific calculated 2D "time-map," which keeps track of the time elapsed from the last fluence delivered at each given point in the image to the end of the beam delivery.

RESULTS

Relative cumulative lag curves reveal that the lag alone can deviate the EPID's perceived dose by as large as 6% (1 MU delivery, 60 s postirradiation). For signal nonlinearity relative to 100 MU, EPID signals per MU of 0.84 and 1.01 were observed for 2 and 500 MUs, respectively. Correction maps were applied to a 1 cm sweeping-window 14 × 14 cm(2) field and clinical head-and-neck IMRT field. A mean correction of 1.028 was implemented in the head-and-neck field, which significantly reduced lag-related asymmetries in the EPID images, and restored linearity to the EPID imager's dose response. Corrections made to the sweeping-field showed good agreement with the treatment planning system-predicted field, yielding an average percent difference of 0.05% ± 0.91%, compared to the -1.32% ± 1.02% before corrections, or 1.75% ± 1.04% when only a signal nonlinearity correction is made.

CONCLUSIONS

Lag and signal-nonlinearity have been quantified for an aS500 EPID imager, and an effective 2D correction method has been developed which effectively removes nonlinearity and lag effects. Both of these effects were shown to negatively impact IMRT verifications. Especially fields that involve prolonged irradiation and small overall MUs should be corrected for in 2D.

A simple approach for EPID dosimetric calibration to overcome the effect of image-lag and ghosting

EPID dosimetry has known drawbacks. The main issue is that a measurable residual signal is observed after the end of irradiation for prolonged periods of time, thus making measurement difficult. We present a detailed analysis of EPID response and suggest a simple, yet accurate approach for calibration that avoids the complexity of incorporating ghosting and image-lag using the maximum integrated signal instead of the total integrated signal. This approach is linear with dose and independent of dose rate.

Calculation of exit dose for conformal and dynamically-wedged fields, based on water-equivalent path length measured with an amorphous silicon electronic portal imaging device

In this study, we use the quadratic calibration method (QCM), in which an EPID image is converted into a matrix of equivalent path lengths (EPLs) and, therefore, exit doses, so as to model doses in conformal and enhanced dynamic wedge (EDW) fields. The QCM involves acquiring series of EPID images at a reference field size for different thicknesses of homogeneous solid water blocks. From these, a set of coefficients is established that is used to compute the EPL of any other irradiated material. To determine the EPL, the irradiated area must be known in order to establish the appropriate scatter correction. A method was devised for the automatic calculation of areas from the EPID image that facilitated the calculation of EPL for any field and exit dose. For EDW fields, the fitting coefficients were modified by utilizing the linac manufacturer's golden segmented treatment tables (GSTT) methodology and MU fraction model. The nonlinear response of the EPL with lower monitor units (MUs) was investigated and slight modification of the algorithm performed to account for this. The method permits 2D dose distributions at the exit of phantom or patient to be generated by relating the EPL with an appropriate depth dose table. The results indicate that the inclusion of MU correction improved the EPL determination. The irradiated field areas can be accurately determined from EPID images to within ± 1% uncertainty. Cross-plane profiles and 2D dose distributions of EPID predicted doses were compared with those calculated with the Eclipse treatment planning system (TPS) and those measured directly with MapCHECK 2 device. Comparison of the 2D EPID dose maps to those from TPS and MapCHECK shows that more than 90% of all points passed the gamma index acceptance criteria of 3% dose difference and 3 mm distance to agreement (DTA), for both conformal and EDW study cases. We conclude that the EPID QCM is an accurate and convenient method for in vivo dosimetry and may, therefore, complement existing techniques.

A generalized calibration procedure for in vivo transit dosimetry using siemens electronic portal imaging devices

A practical and accurate generalized in vivo dosimetry procedure has been implemented for Siemens linacs supplying 6, 10, and 15 MV photon beams, equipped with aSi electronic portal imaging devices (EPIDs). The in vivo dosimetry method makes use of correlation ratios between EPID transit signal, s (t) (0) (TPR,w,L), and phantom mid-plane dose, D (0)(TPR,w,L), as functions of phantom thickness, w, square field dimensions, L, and tissue-phantom ratio TPR(20,10). The s (t) (0) (TPR,w,L) and D (0)(TPR,w,L) values were defined to be independent of the EPID sensitivity and monitor unit calibration, while their dependence on TPR(20,10) was investigated to determine a set of generalized correlation ratios to be used for beams with TPR(20,10) falling in the examined range. This way, other radiotherapy centers can use the method with no need to locally perform the whole set of measurements in solid water phantoms, required to implement it. Tolerance levels for 3D conformal treatments, ranging between ±5 and ±6% according to tumor type and location, were estimated for comparison purposes between reconstructed isocenter dose, D (iso), and treatment planning system (TPS) computed dose D (iso,TPS). Finally a dedicated software, interfaceable with record and verify (R&V) systems used in the centers, was developed to obtain in vivo dosimetry results in less than 2 min after beam delivery.

Monte Carlo-based adaptive EPID dose kernel accounting for different field size responses of imagers

The aim of this study is to present an efficient method to generate imager-specific Monte Carlo (MC)-based dose kernels for amorphous silicon-based electronic portal image device dose prediction and determine the effective backscattering thicknesses for such imagers. EPID field size-dependent responses were measured for five matched Varian accelerators from three institutions with 6 MV beams at the source to detector distance (SDD) of 105 cm. For two imagers, measurements were made with and without the imager mounted on the robotic supporting arm. Monoenergetic energy deposition kernels with 0-2.5 cm of water backscattering thicknesses were simultaneously computed by MC to a high precision. For each imager, the backscattering thickness required to match measured field size responses was determined. The monoenergetic kernel method was validated by comparing measured and predicted field size responses at 150 cm SDD, 10 x 10 cm2 multileaf collimator (MLC) sliding window fields created with 5, 10, 20, and 50 mm gaps, and a head-and-neck (H&N) intensity modulated radiation therapy (IMRT) patient field. Field size responses for the five different imagers deviated by up to 1.3%. When imagers were removed from the robotic arms, response deviations were reduced to 0.2%. All imager field size responses were captured by using between 1.0 and 1.6 cm backscatter. The predicted field size responses by the imager-specific kernels matched measurements for all involved imagers with the maximal deviation of 0.34%. The maximal deviation between the predicted and measured field size responses at 150 cm SDD is 0.39%. The maximal deviation between the predicted and measured MLC sliding window fields is 0.39%. For the patient field, gamma analysis yielded that 99.0% of the pixels have gamma < 1 by the 2%, 2 mm criteria with a 3% dose threshold. Tunable imager-specific kernels can be generated rapidly and accurately in a single MC simulation. The resultant kernels are imager position independent and are able to predict fields with varied incident energy spectra and a H&N IMRT patient field. The proposed adaptive EPID dose kernel method provides the necessary infrastructure to build reliable and accurate portal dosimetry systems.

The use of EPID-measured leaf sequence files for IMRT dose reconstruction in adaptive radiation therapy

For intensity modulated radiation treatment (IMRT) dose reconstruction, multileaf collimator (MLC) log files have been shown applicable for deriving delivered fluence maps. However, MLC log files are dependent on the accuracy of leaf calibration and only available from one linear accelerator manufacturer. This paper presents a proof of feasibility and principles in (1) using an amorphous silicon electronic portal imaging device (aSi-EPID) to capture the MLC segments during an IMRT delivery and (2) reconstituting a leaf sequence (LS) file based on the leaf end positions calculated from the MLC segments and their associated fractional monitor units. These EPID-measured LS files are then used to derive delivered fluence maps for dose reconstruction. The developed approach was tested on a pelvic phantom treated with a typical prostate IMRT plan. The delivered fluence maps, which were derived from the EPID-measured LS files, showed slight differences in the intensity levels compared with the corresponding planned ones. The dose distribution calculated with the delivered fluence maps showed a discernible difference in the high dose region when compared to that calculated with the planned fluence maps. The maximum dose in the former distribution was also 2.5% less than that in the latter one. The EPID-measured LS file can serve the same purpose as a MLC log files does for the derivation of the delivered fluence map and yet is independent of the leaf calibration. The approach also allows users who do not have access to MLC log files to probe the actual IMRT delivery and translate the information gained for dose reconstruction in adaptive radiation therapy.

An empirical model of electronic portal imager response implemented within a commercial treatment planning system for verification of intensity-modulated radiation therapy fields

Quality assurance (QA) of an intensity‐modulated radiation therapy (IMRT) plan is more complex than that of a conventional plan. To improve the efficiency of QA, electronic portal imaging devices (EPIDs) can be used. The major objective of the present work was to use a commercial treatment planning system to model EPID response for the purpose of pre‐treatment IMRT dose verification.

Images were acquired with an amorphous silicon flat panel portal imager (aS500: Varian Medical Systems, Palo Alto, CA) directly irradiated with a 6‐MV photon beam from a Clinac 21EX linear accelerator (Varian Medical Systems). Portal images were acquired for a variety of rectangular fields, from which profiles and relative output factors were extracted. A dedicated machine model was created using the physics tools of the Pinnacle³ (Philips Medical Systems, Madison, WI) treatment planning system to model the data. Starting with the known photon spectrum and assuming an effective depth of 7 cm, machine model parameters were adjusted to best fit measured profile and output factors. The machine parameters of a second model, which assumed a 0.8 MeV monoenergetic photon spectrum and an effective depth in water of 3 cm, were also optimized. The second EPID machine model was used to calculate planar dose maps of simple geometric IMRT fields as well as a 9‐field IMRT plan developed for clinical trials credentialing purposes.

The choice of energy and depth for an EPID machine model influenced the best achievable fit of the optimized machine model to the measured data. When both energy and depth were reduced by a significant amount, a better overall fit was achieved. In either case, the secondary source size and strength could be adjusted to give reasonable agreement with measured data. The gamma evaluation method was used to compare planar dose maps calculated using the second EPID machine model with the EPID images of small IMRT fields. In each case, more than 95% of points fell within 3% of the maximum dose or 3 mm distance to agreement. These results are slightly poorer than those obtained using an ion chamber array, which confirms agreement to within 2% of the maximum dose or 2 mm distance to agreement for all points within these fields.

PACS numbers: 87.55Qr, 87.56.Fc

A literature review of electronic portal imaging for radiotherapy dosimetry

Electronic portal imaging devices (EPIDs) have been the preferred tools for verification of patient positioning for radiotherapy in recent decades. Since EPID images contain dose information, many groups have investigated their use for radiotherapy dose measurement. With the introduction of the amorphous-silicon EPIDs, the interest in EPID dosimetry has been accelerated because of the favourable characteristics such as fast image acquisition, high resolution, digital format, and potential for in vivo measurements and 3D dose verification. As a result, the number of publications dealing with EPID dosimetry has increased considerably over the past approximately 15 years. The purpose of this paper was to review the information provided in these publications. Information available in the literature included dosimetric characteristics and calibration procedures of various types of EPIDs, strategies to use EPIDs for dose verification, clinical approaches to EPID dosimetry, ranging from point dose to full 3D dose distribution verification, and current clinical experience. Quality control of a linear accelerator, pre-treatment dose verification and in vivo dosimetry using EPIDs are now routinely used in a growing number of clinics. The use of EPIDs for dosimetry purposes has matured and is now a reliable and accurate dose verification method that can be used in a large number of situations. Methods to integrate 3D in vivo dosimetry and image-guided radiotherapy (IGRT) procedures, such as the use of kV or MV cone-beam CT, are under development. It has been shown that EPID dosimetry can play an integral role in the total chain of verification procedures that are implemented in a radiotherapy department. It provides a safety net for simple to advanced treatments, as well as a full account of the dose delivered. Despite these favourable characteristics and the vast range of publications on the subject, there is still a lack of commercially available solutions for EPID dosimetry. As strategies evolve and commercial products become available, EPID dosimetry has the potential to become an accurate and efficient means of large-scale patient-specific IMRT dose verification for any radiotherapy department.

[MLC positioning checks and calibration with a portal imaging system]

For checking the leaf positions of a MLC (Multi Leaf Collimator) images are acquired with an EPID (Electronic Portal Imaging Device) and then evaluated with a programme (MLC check) developed in-house. During image acquisition a Perspex tray with two metal markers of known position (in the radiation field) is inserted in the satellite tray holder. After determination of the marker positions within the image coordinate system, the image can be transformed to the radiation field coordinate system. This allows the exact determination of the leaf tip positions relative to the radiation field. This evaluation can be applied to images of arbitrary field shapes, provided they were acquired in the same geometry (EPID position, gantry and collimator angles). The entire measurement procedure is based on images in bmp (Windows Bitmap) format, with a 1024 x 1024 matrix and a pixel depth of 24 bit (8 bit per color channel). A suitable marker plate can be manufactured without sophisticated workload; thus the method may be easily and cost effectively adapted at other locations.

Experimental investigation of the response of an amorphous silicon EPID to intensity modulated radiotherapy beams

The aim of this work was to experimentally determine the difference in response of an amorphous silicon (a-Si) electronic portal imaging device (EPID) to the open and multileaf collimator (MLC) transmitted beam components of intensity modulated radiation therapy (IMRT) beams. EPID dose response curves were measured for open and MLC transmitted (MLCtr) 10 x 10 cm2 beams at central axis and with off axis distance using a shifting field technique. The EPID signal was obtained by replacing the flood-field correction with a pixel sensitivity variation matrix correction. This signal, which includes energy-dependent response, was then compared to ion-chamber measurements. An EPID calibration method to remove the effect of beam energy variations on EPID response was developed for IMRT beams. This method uses the component of open and MLCtr fluence to an EPID pixel calculated from the MLC delivery file and applies separate radially dependent calibration factors for each component. The calibration procedure does not correct for scatter differences between ion chamber in water measurements and EPID response; these must be accounted for separately with a kernel-based approach or similar method. The EPID response at central axis for the open beam was found to be 1.28 +/- 0.03 of the response for the MLCtr beam, with the ratio increasing to 1.39 at 12.5 cm off axis. The EPID response to MLCtr radiation did not change with off-axis distance. Filtering the beam with copper plates to reduce the beam energy difference between open and MLCtr beams was investigated; however, these were not effective at reducing EPID response differences. The change in EPID response for uniform sliding window IMRT beams with MLCtr dose components from 0.3% to 69% was predicted to within 2.3% using the separate EPID response calibration factors for each dose component. A clinical IMRT image calibrated with this method differed by nearly 30% in high MLCtr regions from an image calibrated with an open beam calibration factor only. Accounting for the difference in EPID response to open and MLCtr radiation should improve IMRT dosimetry with a-Si EPIDs.

Off-axis dose response characteristics of an amorphous silicon electronic portal imaging device

Amorphous silicon (a-Si) electronic portal imaging devices (EPIDs) have typically been calibrated to dose at central axis (CAX). Division of acquired images by the flood-field (FF) image that corrects for pixel sensitivity variation as well as open field energy-dependent off-axis response variation should result in a flat EPID response over the entire matrix for the same field size. While the beam profile can be reintroduced to the image by an additional correction matrix, the CAX EPID response to dose calibration factor is assumed to apply to all pixels in the detector. The aim of this work was to investigate the dose response of the Varian aS500 amorphous silicon detector across the entire detector area. First it was established that the EPID response across the panel became stable (within approximately 0.2%) for MU settings greater than approximately 200 MU. The EPID was then FF calibrated with a high MU setting of approximately 400 for all subsequent experiments. Whole detector images with varying MU settings from 2-500 were then acquired for two dose rates (300 and 600 MU/min) for 6 MV photons for two EPIDs. The FF corrected EPID response was approximately flat or uniform across the detector for greater than 100 MU delivered (within 0.5%). However, the off-axis EPID response was greater than the CAX response for small MU irradiations, giving a raised EPID profile. Up to 5% increase in response at 20 cm off-axis compared to CAX was found for very small MU settings for one EPID, while it was within 2% for the second (newer) EPID. Off-axis response nonuniformities attributed to detector damage were also found for the older EPID. Similar results were obtained with the EPID at 18 MV energy and operating in asynchronous mode (acquisition not synchronized with beam pulses), however the profiles were flatter and more irregular for the small MU irradiations. By moving the detector laterally and repeating the experiments, the increase in response off-axis was found to depend on the pixel position relative to the beam CAX. When the beam was heavily filtered by a phantom the off-axis response variation was reduced markedly to within 0.5% for all MU settings. Independent measurements of off-axis point doses with ion chamber did not show any change in off-axis factor with MUs. Measurements of beam quality (TMR20-10) for MU settings of 2, 5, and 100 at central axis and at 15 cm off-axis could not explain the effect. The response change is unlikely to be significant for clinical IMRT verification with this imaging/acclerator system where MUs are of the order of 100-300, provided the detector does not exhibit radiation damage artifacts.