Rapid in vivo EPID image prediction using a combination of analytically calculated attenuation and AI predicted scatter

BACKGROUND

The electronic portal imaging device (EPID) can be used in vivo, to detect on-treatment errors by evaluating radiation exiting a patient. To detect deviations from the planning intent, image predictions need to be modeled based on the patient's anatomy and plan information. To date in vivo transit images have been predicted using Monte Carlo (MC) algorithms. A deep learning approach can make predictions faster than MC and only requires patient information for training.

PURPOSE

To test the feasibility and reliability of creating a deep-learning model with patient data for predicting in vivo EPID images for IMRT treatments.

METHODS

In our approach, the in vivo EPID image was separated into contributions from primary and scattered photons. A primary photon attenuation function was determined by measuring attenuation factors for various thicknesses of solid water. The scatter component of in vivo EPID images was estimated using a convolutional neural network (CNN). The CNN input was a 3-channel image comprised of the non-transit EPID image and ray tracing projections through a pretreatment CBCT. The predicted scatter component was added to the primary attenuation component to give the full predicted in vivo EPID image. We acquired 193 IMRT fields/images from 93 patients treated on the Varian Halcyon. Model training:validation:test dataset ratios were 133:20:40 images. Additional patient plans were delivered to anthropomorphic phantoms, yielding 75 images for further validation. We assessed model accuracy by comparing model-calculated and measured in vivo images with a gamma comparison.

RESULTS

Comparing the model-calculated and measured in vivo images gives a mean gamma pass rate for the training:validation:test datasets of 95.4%:94.1%:92.9% for 3%/3 mm and 98.4%:98.4%:96.8% for 5%/3 mm. For images delivered to phantom data sets the average gamma pass rate was 96.4% (3%/3 mm criteria). In all data sets, the lower passing rates of some images were due to CBCT artifacts and patient motion that occurred between the time of CBCT and treatment.  CONCLUSIONS: The developed deep-learning-based model can generate in vivo EPID images with a mean gamma pass rate greater than 92% (3%/3 mm criteria). This approach provides an alternative to MC prediction algorithms. Image predictions can be made in 30 ms on a standard GPU. In future work, image predictions from this model can be used to detect in vivo treatment errors and on-treatment changes in patient anatomy, providing an additional layer of patient-specific quality assurance.

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.

A feasibility study for in vivo treatment verification of IMRT using Monte Carlo dose calculation and deep learning-based modelling of EPID detector response

Background

This paper describes the development of a predicted electronic portal imaging device (EPID) transmission image (TI) using Monte Carlo (MC) and deep learning (DL). The measured and predicted TI were compared for two-dimensional in vivo radiotherapy treatment verification.

Methods

The plan CT was pre-processed and combined with solid water and then imported into PRIMO. The MC method was used to calculate the dose distribution of the combined CT. The U-net neural network-based deep learning model was trained to predict EPID TI based on the dose distribution of solid water calculated by PRIMO. The predicted TI was compared with the measured TI for two-dimensional in vivo treatment verification.

Results

The EPID TI of 1500 IMRT fields were acquired, among which 1200, 150, and 150 fields were used as the training set, the validation set, and the test set, respectively. A comparison of the predicted and measured TI was carried out using global gamma analyses of 3%/3 mm and 2%/2 mm (5% threshold) to validate the model's accuracy. The gamma pass rates were greater than 96.7% and 92.3%, and the mean gamma values were 0.21 and 0.32, respectively.

Conclusions

Our method facilitates the modelling process more easily and increases the calculation accuracy when using the MC algorithm to simulate the EPID response, and has potential to be used for in vivo treatment verification in the clinic.

A method for in vivo treatment verification of IMRT and VMAT based on electronic portal imaging device

Background

Intensity-modulated radiation therapy (IMRT) and volume-modulated arc therapy (VMAT) are rather complex treatment techniques and require patient-specific quality assurance procedures. Electronic portal imaging devices (EPID) are increasingly used in the verification of radiation therapy (RT). This work aims to develop a novel model to predict the EPID transmission image (TI) with fluence maps from the RT plan. The predicted TI is compared with the measured TI for in vivo treatment verification.

Methods

The fluence map was extracted from the RT plan and corrections of penumbra, response, global field output, attenuation, and scatter were applied before the TI was calculated. The parameters used in the model were calculated separately for central axis and off-axis points using a series of EPID measurement data. Our model was evaluated using a CIRS thorax phantom and 20 clinical plans (10 IMRT and 10 VMAT) optimized for head and neck, breast, and rectum treatments.

Results

Comparisons of the predicted and measured images were carried out using a global gamma analysis of 3%/2 mm (10% threshold) to validate the accuracy of the model. The gamma pass rates for IMRT and VMAT were greater than 97.2% and 94.5% at 3%/2 mm, respectively.

Conclusion

We have developed an accurate and straightforward EPID-based quality assurance model that can potentially be used for in vivo treatment verification of the IMRT and VMAT delivery.

Toward real-time verification for MLC tracking treatments using time-resolved EPID imaging

PURPOSE

In multileaf collimator (MLC) tracking, the MLC positions from the original treatment plan are continuously modified to account for intrafraction tumor motion. As the treatment is adapted in real time, there is additional risk of delivery errors which cannot be detected using traditional pretreatment dose verification. The purpose of this work is to develop a system for real-time geometric verification of MLC tracking treatments using an electronic portal imaging device (EPID).

METHODS

MLC tracking was utilized during volumetric modulated arc therapy (VMAT). During these deliveries, treatment beam images were taken at 9.57 frames per second using an EPID and frame grabber computer. MLC positions were extracted from each image frame and used to assess delivery accuracy using three geometric measures: the location, size, and shape of the radiation field. The EPID-measured field location was compared to the tumor motion measured by implanted electromagnetic markers. The size and shape of the beam were compared to the size and shape from the original treatment plan, respectively. This technique was validated by simulating errors in phantom test deliveries and by comparison between EPID measurements and treatment log files. The method was applied offline to images acquired during the LIGHT Stereotactic Ablative Body Radiotherapy (SABR) clinical trial, where MLC tracking was performed for 17 lung cancer patients. The EPID-based verification results were subsequently compared to post-treatment dose reconstruction.

RESULTS

Simulated field location errors were detected during phantom validation tests with an uncertainty of 0.28 mm (parallel to MLC motion) and 0.38 mm (perpendicular), expressed as a root-mean-square error (RMS_Error ). For simulated field size errors, the RMS_Error was 0.47 cm² and field shape changes were detected for random errors with standard deviation ≥ 2.5 mm. For clinical lung SABR deliveries, field location errors of 1.6 mm (parallel MLC motion) and 4.9 mm (perpendicular) were measured (expressed as a full-width-half-maximum). The mean and standard deviation of the errors in field size and shape were 0.0 ± 0.3 cm² and 0.3 ± 0.1 (expressed as a translation-invariant normalized RMS). No correlation was observed between geometric errors during each treatment fraction and dosimetric errors in the reconstructed dose to the target volume for this cohort of patients.

CONCLUSION

A system for real-time delivery verification has been developed for MLC tracking using time-resolved EPID imaging. The technique has been tested offline in phantom-based deliveries and clinical patient deliveries and was used to independently verify the geometric accuracy of the MLC during MLC tracking radiotherapy.

Analyses of integrated EPID images for on-treatment quality assurance to account for interfractional variations in volumetric modulated arc therapy

Purpose

To investigate the effects of interfractional variation, such as anatomical changes and setup errors, on dose delivery during treatment for prostate cancer (PC) and head and neck cancer (HNC) by courses of volumetric modulated arc therapy (VMAT) aided by on‐treatment electronic portal imaging device (EPID) images.

Methods

Seven patients with PC and 20 patients with HNC who had received VMAT participated in this study. After obtaining photon fluence at the position of the EPID for each treatment arc from on‐treatment integrated EPID images, we calculated the differences between the fluence for the first fraction and each subsequent fraction for each arc. The passing rates were investigated based on a tolerance level of 3% of the maximum fluence during the treatment courses and the correlations between the passing rates and anatomical changes.

Results

In PC, the median and lowest passing rates were 99.8% and 95.2%, respectively. No correlations between passing rates and interfractional variation were found. In HNC, the median passing rate of all fractions was 93.0%, and the lowest passing rate was 79.6% during the 35th fraction. Spearman’s correlation coefficients between the passing rates and changes in weight or neck volume were − 0.77 and − 0.74, respectively.

Conclusions

Analyses of the on‐treatment EPID images facilitates estimates of the interfractional anatomical variation in HNC patients during VMAT and thus improves assessments of the need for re‐planning or adaptive strategies and the timing thereof.

A method to verify sections of arc during intrafraction portal dosimetry for prostate VMAT

This study investigates the use of a running sum of images during segment-resolved intrafraction portal dosimetry for volumetric modulated arc therapy (VMAT), so as to alert the operator to an error before it becomes irremediable. At the time of treatment planning, predicted portal images were created for each segment of the VMAT arc, and at the time of delivery, intrafraction monitoring software polled the portal imager to read new images as they became available. The predicted and measured images were compared and displayed on a segment basis. In particular, a running sum of images from ten segments (a 'section') was investigated, with mean absolute difference between predicted and measured images being quantified. Images for 13 prostate patients were used to identify appropriate tolerance values for this statistic. Errors in monitor units of 2%-10%, field size of 2-10 mm, field position of 2-10 mm and path length of 10-50 mm were deliberately introduced into the treatment plans and delivered to a water-equivalent phantom and the sensitivity of the method to these errors was investigated. Gross errors were also considered for one case. The patient images show considerable variability from segment to segment, but when using a section of the arc the variability is reduced, so that the maximum value of mean absolute difference between predicted and measured images is reduced to below 12%, after excluding the first 10% of segments. This tolerance level is also found to be applicable for delivery of the plans to a water-equivalent phantom. Using this as a tolerance level for the error plans, a 10% increase in monitor units is detected, 4 mm increase or shift in multileaf collimator settings can be detected, and an air gap of dimensions 40 mm  ×  50 mm is detected. Gross errors can also be detected instantly after the first 10% of segments. The running difference between predicted and measured images over ten segments is able to identify errors at specific regions of the arc, as well as in the overall treatment.

Evaluating the sensitivity of Halcyon's automatic transit image acquisition for treatment error detection: A phantom study using static IMRT

PURPOSE

The Varian Halcyon™ electronic portal imaging detector is always in-line with the beam and automatically acquires transit images for every patient with full-field coverage. These images could be used for "every patient, every monitor unit" quality assurance (QA) and eventually adaptive radiotherapy. This study evaluated the imager's sensitivity to potential clinical errors and day-to-day variations from clinical exit images.

METHODS

Open and modulated fields were delivered for each potential error. To evaluate output changes, monitor units were scaled by 2%-10% and delivered to solid water slabs and a homogeneous CIRS phantom. To mimic weight changes, 0.5-5.0 cm of buildup was added to the solid water. To evaluate positioning changes, a homogeneous and heterogeneous CIRS phantom were shifted 2-10 cm and 0.2-1.5 cm, respectively. For each test, mean relative differences (MRDs) and standard deviations in the pixel-difference histograms (σ_RD ) between test and baseline images were calculated. Lateral shift magnitudes were calculated using cross-correlation and edge-detection filtration. To assess patient variations, MRD and σ_RD were calculated from six prostate patients' daily exit images and compared between fractions with and without gas present.

RESULTS

MRDs responded linearly to output and buildup changes with a standard deviation of 0.3%, implying a 1% output change and 0.2 cm changes in buildup could be detected with 2.5σ confidence. Shifting the homogenous phantom laterally resulted in detectable MRD and σ_RD changes, and the cross-correlation function calculated the shift to within 0.5 mm for the heterogeneous phantom. MRD and σ_RD values were significantly associated with the presence of gas for five of the six patients.

CONCLUSIONS

Rapid analyses of automatically acquired Halcyon™ exit images could detect mid-treatment changes with high sensitivity, though appropriate thresholds will need to be set. This study presents the first steps toward developing effortless image evaluation for all aspects of every patient's treatment.

A novel approach to EPID-based 3D volumetric dosimetry for IMRT and VMAT QA

Intensity modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) are relatively complex treatment delivery techniques and require quality assurance (QA) procedures. Pre-treatment dosimetric verification represents a fundamental QA procedure in daily clinical routine in radiation therapy. The purpose of this study is to develop an EPID-based approach to reconstruct a 3D dose distribution as imparted to a virtual cylindrical water phantom to be used for plan-specific pre-treatment dosimetric verification for IMRT and VMAT plans. For each depth, the planar 2D dose distributions acquired in air were back-projected and convolved by depth-specific scatter and attenuation kernels. The kernels were obtained by making use of scatter and attenuation models to iteratively estimate the parameters from a set of reference measurements. The derived parameters served as a look-up table for reconstruction of arbitrary measurements. The summation of the reconstructed 3D dose distributions resulted in the integrated 3D dose distribution of the treatment delivery. The accuracy of the proposed approach was validated in clinical IMRT and VMAT plans by means of gamma evaluation, comparing the reconstructed 3D dose distributions with Octavius measurement. The comparison was carried out using (3%, 3 mm) criteria scoring 99% and 96% passing rates for IMRT and VMAT, respectively. An accuracy comparable to the one of the commercial device for 3D volumetric dosimetry was demonstrated. In addition, five IMRT and five VMAT were validated against the 3D dose calculation performed by the TPS in a water phantom using the same passing rate criteria. The median passing rates within the ten treatment plans was 97.3%, whereas the lowest was 95%. Besides, the reconstructed 3D distribution is obtained without predictions relying on forward dose calculation and without external phantom or dosimetric devices. Thus, the approach provides a fully automated, fast and easy QA procedure for plan-specific pre-treatment dosimetric verification.

Comparison of forward- and back-projection in vivo EPID dosimetry for VMAT treatment of the prostate

In the forward-projection method of portal dosimetry for volumetric modulated arc therapy (VMAT), the integrated signal at the electronic portal imaging device (EPID) is predicted at the time of treatment planning, against which the measured integrated image is compared. In the back-projection method, the measured signal at each gantry angle is back-projected through the patient CT scan to give a measure of total dose to the patient. This study aims to investigate the practical agreement between the two types of EPID dosimetry for prostate radiotherapy. The AutoBeam treatment planning system produced VMAT plans together with corresponding predicted portal images, and a total of 46 sets of gantry-resolved portal images were acquired in 13 patients using an iViewGT portal imager. For the forward-projection method, each acquisition of gantry-resolved images was combined into a single integrated image and compared with the predicted image. For the back-projection method, iViewDose was used to calculate the dose distribution in the patient for comparison with the planned dose. A gamma index for 3% and 3 mm was used for both methods. The results were investigated by delivering the same plans to a phantom and repeating some of the deliveries with deliberately introduced errors. The strongest agreement between forward- and back-projection methods is seen in the isocentric intensity/dose difference, with moderate agreement in the mean gamma. The strongest correlation is observed within a given patient, with less correlation between patients, the latter representing the accuracy of prediction of the two methods. The error study shows that each of the two methods has its own distinct sensitivity to errors, but that overall the response is similar. The forward- and back-projection EPID dosimetry methods show moderate agreement in this series of prostate VMAT patients, indicating that both methods can contribute to the verification of dose delivered to the patient.

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.

Commissioning and quality assurance for VMAT delivery systems: An efficient time-resolved system using real-time EPID imaging

PURPOSE

An ideal commissioning and quality assurance (QA) program for Volumetric Modulated Arc Therapy (VMAT) delivery systems should assess the performance of each individual dynamic component as a function of gantry angle. Procedures within such a program should also be time-efficient, independent of the delivery system and be sensitive to all types of errors. The purpose of this work is to develop a system for automated time-resolved commissioning and QA of VMAT control systems which meets these criteria.

METHODS

The procedures developed within this work rely solely on images obtained, using an electronic portal imaging device (EPID) without the presence of a phantom. During the delivery of specially designed VMAT test plans, EPID frames were acquired at 9.5 Hz, using a frame grabber. The set of test plans was developed to individually assess the performance of the dose delivery and multileaf collimator (MLC) control systems under varying levels of delivery complexities. An in-house software tool was developed to automatically extract features from the EPID images and evaluate the following characteristics as a function of gantry angle: dose delivery accuracy, dose rate constancy, beam profile constancy, gantry speed constancy, dynamic MLC positioning accuracy, MLC speed and acceleration constancy, and synchronization between gantry angle, MLC positioning and dose rate. Machine log files were also acquired during each delivery and subsequently compared to information extracted from EPID image frames.

RESULTS

The largest difference between measured and planned dose at any gantry angle was 0.8% which correlated with rapid changes in dose rate and gantry speed. For all other test plans, the dose delivered was within 0.25% of the planned dose for all gantry angles. Profile constancy was not found to vary with gantry angle for tests where gantry speed and dose rate were constant, however, for tests with varying dose rate and gantry speed, segments with lower dose rate and higher gantry speed exhibited less profile stability. MLC positional accuracy was not observed to be dependent on the degree of interdigitation. MLC speed was measured for each individual leaf and slower leaf speeds were shown to be compensated for by lower dose rates. The test procedures were found to be sensitive to 1 mm systematic MLC errors, 1 mm random MLC errors, 0.4 mm MLC gap errors and synchronization errors between the MLC, dose rate and gantry angle controls systems of 1°. In general, parameters measured by both EPID and log files agreed with the plan, however, a greater average departure from the plan was evidenced by the EPID measurements.

CONCLUSION

QA test plans and analysis methods have been developed to assess the performance of each dynamic component of VMAT deliveries individually and as a function of gantry angle. This methodology relies solely on time-resolved EPID imaging without the presence of a phantom and has been shown to be sensitive to a range of delivery errors. The procedures developed in this work are both comprehensive and time-efficient and can be used for streamlined commissioning and QA of VMAT delivery systems.

A Swiss cheese error detection method for real-time EPID-based quality assurance and error prevention

PURPOSE

To develop a robust and efficient process that detects relevant dose errors (dose errors of ≥5%) in external beam radiation therapy and directly indicates the origin of the error. The process is illustrated in the context of electronic portal imaging device (EPID)-based angle-resolved volumetric-modulated arc therapy (VMAT) quality assurance (QA), particularly as would be implemented in a real-time monitoring program.

METHODS

A Swiss cheese error detection (SCED) method was created as a paradigm for a cine EPID-based during-treatment QA. For VMAT, the method compares a treatment plan-based reference set of EPID images with images acquired over each 2° gantry angle interval. The process utilizes a sequence of independent consecutively executed error detection tests: an aperture check that verifies in-field radiation delivery and ensures no out-of-field radiation; output normalization checks at two different stages; global image alignment check to examine if rotation, scaling, and translation are within tolerances; pixel intensity check containing the standard gamma evaluation (3%, 3 mm) and pixel intensity deviation checks including and excluding high dose gradient regions. Tolerances for each check were determined. To test the SCED method, 12 different types of errors were selected to modify the original plan. A series of angle-resolved predicted EPID images were artificially generated for each test case, resulting in a sequence of precalculated frames for each modified treatment plan. The SCED method was applied multiple times for each test case to assess the ability to detect introduced plan variations. To compare the performance of the SCED process with that of a standard gamma analysis, both error detection methods were applied to the generated test cases with realistic noise variations.

RESULTS

Averaged over ten test runs, 95.1% of all plan variations that resulted in relevant patient dose errors were detected within 2° and 100% within 14° (<4% of patient dose delivery). Including cases that led to slightly modified but clinically equivalent plans, 89.1% were detected by the SCED method within 2°. Based on the type of check that detected the error, determination of error sources was achieved. With noise ranging from no random noise to four times the established noise value, the averaged relevant dose error detection rate of the SCED method was between 94.0% and 95.8% and that of gamma between 82.8% and 89.8%.

CONCLUSIONS

An EPID-frame-based error detection process for VMAT deliveries was successfully designed and tested via simulations. The SCED method was inspected for robustness with realistic noise variations, demonstrating that it has the potential to detect a large majority of relevant dose errors. Compared to a typical (3%, 3 mm) gamma analysis, the SCED method produced a higher detection rate for all introduced dose errors, identified errors in an earlier stage, displayed a higher robustness to noise variations, and indicated the error source.

The dosimetric impact of control point spacing for sliding gap MLC fields

Dynamic sliding gap multileaf collimator (MLC) fields are used to model MLC properties within the treatment planning system (TPS) for dynamic treatments. One of the key MLC properties in the Eclipse TPS is the dosimetric leaf gap (DLG) and precise determination of this parameter is paramount to ensuring accurate dose delivery. In this investigation, we report on how the spacing between control points (CPs) for sliding gap fields impacts the dose delivery, MLC positioning accuracy, and measurement of the DLG. The central axis dose was measured for sliding gap MLC fields with gap widths ranging from 2 to 40 mm. It was found that for deliveries containing two CPs, the central axis dose was underestimated by the TPS for all gap widths, with the maximum difference being 8% for a 2 mm gap field. For the same sliding gap fields containing 50 CPs, the measured dose was always within ±2% of the TPS dose. By directly measuring the MLC trajectories we show that this dose difference is due to a systematic MLC gap error for fields containing two CPs, and that the cause of this error is due to the leaf position offset table which is incorrectly applied when the spacing between CPs is too large. This MLC gap error resulted in an increase in the measured DLG of 0.5 mm for both 6 MV and 10 MV, when using fields with 2 CPs compared to 50 CPs. Furthermore, this change in DLG was shown to decrease the mean TPS‐calculated dose to the target volume by 2.6% for a clinical IMRT test plan. This work has shown that systematic MLC positioning errors occur for sliding gap MLC fields containing two CPs and that using these fields to model critical TPS parameters, such as the DLG, may result in clinically significant systematic dose calculation errors during subsequent dynamic MLC treatments.

PACS number(s): 87.56.nk

An EPID-based system for gantry-resolved MLC quality assurance for VMAT

Multileaf collimator (MLC) positions should be precisely and independently mea-sured as a function of gantry angle as part of a comprehensive quality assurance (QA) program for volumetric-modulated arc therapy (VMAT). It is also ideal that such a QA program has the ability to relate MLC positional accuracy to patient-specific dosimetry in order to determine the clinical significance of any detected MLC errors. In this work we propose a method to verify individual MLC trajectories during VMAT deliveries for use as a routine linear accelerator QA tool. We also extend this method to reconstruct the 3D patient dose in the treatment planning sys-tem based on the measured MLC trajectories and the original DICOM plan file. The method relies on extracting MLC positions from EPID images acquired at 8.41fps during clinical VMAT deliveries. A gantry angle is automatically tagged to each image in order to obtain the MLC trajectories as a function of gantry angle. This analysis was performed for six clinical VMAT plans acquired at monthly intervals for three months. The measured trajectories for each delivery were compared to the MLC positions from the DICOM plan file. The maximum mean error detected was 0.07 mm and a maximum root-mean-square error was 0.8 mm for any leaf of any delivery. The sensitivity of this system was characterized by introducing random and systematic MLC errors into the test plans. It was demonstrated that the system is capable of detecting random and systematic errors on the range of 1-2mm and single leaf calibration errors of 0.5 mm. The methodology developed in the work has potential to be used for efficient routine linear accelerator MLC QA and pretreatment patient-specific QA and has the ability to relate measured MLC positional errors to 3D dosimetric errors within a patient volume.

Investigation of a real-time EPID-based patient dose monitoring safety system using site-specific control limits

Purpose

The aim of this study is to investigate the performance and limitations of a real-time transit electronic portal imaging device (EPID) dosimetry system for error detection during dynamic intensity modulated radiation therapy (IMRT) treatment delivery. Sites studied are prostate, head and neck (HN), and rectal cancer treatments.

Methods

The system compares measured cumulative transit EPID image frames with predicted cumulative image frames in real-time during treatment using a χ comparison with 4 %, 4 mm criteria. The treatment site-specific thresholds (prostate, HN and rectum IMRT) were determined using initial data collected from 137 patients (274 measured treatment fractions) and a statistical process control methodology. These thresholds were then applied to data from 15 selected patients including 5 prostate, 5 HN, and 5 rectum IMRT treatments for system evaluation and classification of error sources.

Results

Clinical demonstration of real-time transit EPID dosimetry in IMRT was presented. For error simulation, the system could detect gross errors (i.e. wrong patient, wrong plan, wrong gantry angle) immediately after EPID stabilisation; 2 seconds after the start of treatment. The average rate of error detection was 7.0 % (prostate = 5.6 %, HN= 8.7 % and rectum = 6.7 %). The detected errors were classified as either clinical in origin (e.g. patient anatomical changes), or non-clinical in origin (e.g. detection system errors). Classified errors were 3.2 % clinical and 3.9 % non-clinical.

Conclusion

An EPID-based real-time error detection method for treatment verification during dynamic IMRT has been developed and tested for its performance and limitations. The system is able to detect gross errors in real-time, however improvement in system robustness is required to reduce the non-clinical sources of error detection.

Time-Resolved Versus Integrated Transit Planar Dosimetry for Volumetric Modulated Arc Therapy

Purpose:

It is desirable that dosimetric deviations during radiation treatments are detected. Integrated transit planar dosimetry is commonly used to evaluate external beam treatments such as volumetric-modulated arc therapy. This work focuses on patient geometry changes which result in differences between the planned and the delivered radiation dose. Integrated transit planar dosimetry will average out some deviations. Novel time-resolved transit planar dosimetry compares the delivered dose of volumetric-modulated arc therapy to the planned dose at various time points. Four patient cases are shown where time-resolved transit planar dosimetry detects patient geometry changes during treatment.

Methods:

A control point to control point comparison between the planned dose and the treatment dose of volumetric-modulated arc therapy beams is calculated using the planning computed tomography and the kV cone-beam computed tomography of the day and evaluated with a time-resolved γ function. Results were computed for 4 patients treated with volumetric-modulated arc therapy, each showing an anatomical change: pleural effusion, rectal gas pockets, and tumor regression.

Results:

In all cases, the geometrical change was detected by time-resolved transit planar dosimetry, whereas integrated transit planar dosimetry showed minor or no indication of the dose discrepancy. Both tumor regression cases were detected earlier in the treatment with time-resolved planar dosimetry in comparison to integrated transit planar dosimetry. The pleural effusion and the gas pocket were detected exclusively with time-resolved transit planar dosimetry.

Conclusions:

Clinical cases were presented in this proof-of-principle study in which integrated transit planar dosimetry did not detect dosimetrically relevant deviations to the same extent time-resolved transit planar dosimetry was able to. Time-resolved transit planar dosimetry also provides results that can be presented as a function of arc delivery angle allowing easier interpretation compared to integrated transit planar dosimetry.

Investigation of the mechanical performance of Siemens linacs components during arc: gantry, MLC, and electronic portal imaging device

Background

In radiotherapy treatments, it is crucial to monitor the performance of linac components including gantry, collimation system, and electronic portal imaging device (EPID) during arc deliveries. In this study, a simple EPID-based measurement method is suggested in conjunction with an algorithm to investigate the stability of these systems at various gantry angles with the aim of evaluating machine-related errors in treatments.

Methods

The EPID sag, gantry sag, changes in source-to-detector distance (SDD), EPID and collimator skewness, EPID tilt, and the sag in leaf bank assembly due to linac rotation were separately investigated by acquisition of 37 EPID images of a simple phantom with five ball bearings at various gantry angles. A fast and robust software package was developed for automated analysis of image data. Three Siemens linacs were investigated.

Results

The average EPID sag was within 1 mm for all tested linacs. Two machines showed >1 mm gantry sag. Changes in the SDD values were within 7.5 mm. EPID skewness and tilt values were <1° in all machines. The maximum sag in leaf bank assembly was <1 mm.

Conclusion

The method and software developed in this study provide a simple tool for effective investigation of the behavior of Siemens linac components with gantry rotation. Such a comprehensive study has been performed for the first time on Siemens machines.

Camera selection for real-time in vivo radiation treatment verification systems using Cherenkov imaging

PURPOSE

To identify achievable camera performance and hardware needs in a clinical Cherenkov imaging system for real-time, in vivo monitoring of the surface beam profile on patients, as novel visual information, documentation, and possible treatment verification for clinicians.

METHODS

Complementary metal-oxide-semiconductor (CMOS), charge-coupled device (CCD), intensified charge-coupled device (ICCD), and electron multiplying-intensified charge coupled device (EM-ICCD) cameras were investigated to determine Cherenkov imaging performance in a clinical radiotherapy setting, with one emphasis on the maximum supportable frame rate. Where possible, the image intensifier was synchronized using a pulse signal from the Linac in order to image with room lighting conditions comparable to patient treatment scenarios. A solid water phantom irradiated with a 6 MV photon beam was imaged by the cameras to evaluate the maximum frame rate for adequate Cherenkov detection. Adequate detection was defined as an average electron count in the background-subtracted Cherenkov image region of interest in excess of 0.5% (327 counts) of the 16-bit maximum electron count value. Additionally, an ICCD and an EM-ICCD were each used clinically to image two patients undergoing whole-breast radiotherapy to compare clinical advantages and limitations of each system.

RESULTS

Intensifier-coupled cameras were required for imaging Cherenkov emission on the phantom surface with ambient room lighting; standalone CMOS and CCD cameras were not viable. The EM-ICCD was able to collect images from a single Linac pulse delivering less than 0.05 cGy of dose at 30 frames/s (fps) and pixel resolution of 512 × 512, compared to an ICCD which was limited to 4.7 fps at 1024 × 1024 resolution. An intensifier with higher quantum efficiency at the entrance photocathode in the red wavelengths [30% quantum efficiency (QE) vs previous 19%] promises at least 8.6 fps at a resolution of 1024 × 1024 and lower monetary cost than the EM-ICCD.

CONCLUSIONS

The ICCD with an intensifier better optimized for red wavelengths was found to provide the best potential for real-time display (at least 8.6 fps) of radiation dose on the skin during treatment at a resolution of 1024 × 1024.

Is integrated transit planar portal dosimetry able to detect geometric changes in lung cancer patients treated with volumetric modulated arc therapy?

BACKGROUND

Geometric changes are frequent during the course of treatment of lung cancer patients. This may potentially result in deviations between the planned and actual delivered dose. Electronic portal imaging device (EPID)-based integrated transit planar portal dosimetry (ITPD) is a fast method for absolute in-treatment dose verification. The aim of this study was to investigate if ITPD could detect geometric changes in lung cancer patients.

MATERIALS AND METHODS

A total of 460 patients treated with volumetric modulated arc therapy (VMAT) following daily cone beam computed tomography (CT)-based setup were visually inspected for geometrical changes on a daily basis. Forty-six patients were subject to changes and had a re-CT and an adaptive treatment plan. The reasons for adaptation were: change in atelectasis (n = 18), tumor regression (n = 9), change in pleural effusion (n = 8) or other causes (n = 11). The ITPDs were calculated on both the initial planning CT and the re-CT and compared with a global gamma (γ) evaluation (criteria: 3%\3mm). A treatment fraction failed when the percentage of pixels failing in the radiation fields exceeded 10%. Dose-volume histograms (DVHs) were compared between the initial plan versus the plan re-calculated on the re-CT.

RESULTS

The ITPD threshold method detected 76% of the changes in atelectasis, while only 50% of the tumor regression cases and 42% of the pleural effusion cases were detected. Only 10% of the cases adapted for other reasons were detected with ITPD. The method has a 17% false-positive rate. No significant correlations were found between changes in DVH metrics and γ fail-rates.

CONCLUSIONS

This study showed that most cases with geometric changes caused by atelectasis could be captured by ITPD, however for other causes ITPD is not sensitive enough to detect the clinically relevant changes and no predictive power of ITPD was found.

A comprehensive study of the mechanical performance of gantry, EPID and the MLC assembly in Elekta linacs during gantry rotation

OBJECTIVE

In radiotherapy treatments, it is crucial to monitor the performance of linear accelerator (linac) components, including gantry, collimation system and electronic portal imaging device (EPID) during arc deliveries. In this study, a simple EPID-based measurement method is suggested in conjunction with an algorithm to investigate the stability of these systems at various gantry angles with the aim of evaluating machine-related errors in treatments.

METHODS

The EPID sag, gantry sag, changes in source-to-detector distance (SDD), EPID and collimator skewness, EPID tilt and the sag in leaf bank assembly owing to linac rotation were separately investigated by acquisition of 37 EPID images of a simple phantom with 5 ball bearings at various gantry angles. A fast and robust software package was developed for automated analysis of the image data. Nine Elekta AB (Stockholm, Sweden) linacs of different models and number of years in service were investigated.

RESULTS

The average EPID sag was within 2 mm for all tested linacs. Some machines showed >1-mm gantry sag. Changes in the SDD values were within 1.3 cm. EPID skewness and tilt values were <1° in all machines. The maximum sag in multileaf collimator leaf bank assemblies was around 1 mm. A meaningful correlation was found between the age of the linacs and their mechanical performance. Conclusions and Advances in knowledge: The method and software developed in this study provide a simple tool for effective investigation of the behaviour of Elekta linac components with gantry rotation. Such a comprehensive study has been performed for the first time on Elekta machines.

Comparison between an in-house 1D profile correction method and a 2D correction provided in Varian's PDPC Package for improving the accuracy of portal dosimetry images

While commissioning Varian's Portal Dose Image Prediction (PDIP) algorithm for portal dosimetry, an asymmetric radial response in the portal imager due to backscatter from the support arm was observed. This asymmetric response led to differences on the order of 2%–3% for simple square fields (<20×20 cm2) when comparing the measured to predicted portal fluences. A separate problem was that discrepancies of up to 10% were seen in measured to predicted portal fluences at increasing off‐axis distance (>10 cm). We have modified suggested methods from the literature to provide a 1D correction for the off‐axis response problem which adjusts the diagonal profile used in the portal imager calibration. This inherently cannot fix the 2D problem since the PDIP algorithm assumes a radially symmetric response and will lead to some uncertainty in portal dosimetry results. Varian has recently released generic “2D correction” files with their Portal Dosimetry Pre‐configuration (PDPC) package, but no independent testing has been published. We present the comparison between QA results using the Varian correction method to results using our 1D profile correction method using the gamma passing rates with a 3%, 3 mm criterion. The average, minimum, and maximum gamma pass rates for nine fixed‐field IMRT fields at gantry 0° using our profile correction method were 98.1%, 93.7%, and 99.8%, respectively, while the results using the PDPC correction method were 98.4%, 93.1%, and 99.8%. For four RapidArc fields, the average, minimum, and maximum gamma pass rates using our correction method were 99.6%, 99.4%, and 99.9%, respectively, while the results using the PDPC correction method were 99.8%, 99.5%, and 99.9%. The average gamma pass rates for both correction methods are quite similar, but both show improvement over the uncorrected results.

PACS numbers: 87.55.Qr, 87.55.N‐

Sensitivity of an Elekta iView GT a-Si EPID model to delivery errors for pre-treatment verification of IMRT fields

A Monte Carlo model of an Elekta iViewGT amorphous silicon electronic portal imaging device (a-Si EPID) has been validated for pre-treatment verification of clinical IMRT treatment plans. The simulations involved the use of the BEAMnrc and DOSXYZnrc Monte Carlo codes to predict the response of the iViewGT a-Si EPID model. The predicted EPID images were compared to the measured images obtained from the experiment. The measured EPID images were obtained by delivering a photon beam from an Elekta Synergy linac to the Elekta iViewGT a-Si EPID. The a-Si EPID was used with no additional build-up material. Frame averaged EPID images were acquired and processed using in-house software. The agreement between the predicted and measured images was analyzed using the gamma analysis technique with acceptance criteria of 3 %/3 mm. The results show that the predicted EPID images for four clinical IMRT treatment plans have a good agreement with the measured EPID signal. Three prostate IMRT plans were found to have an average gamma pass rate of more than 95.0 % and a spinal IMRT plan has the average gamma pass rate of 94.3 %. During the period of performing this work a routine MLC calibration was performed and one of the IMRT treatments re-measured with the EPID. A change in the gamma pass rate for one field was observed. This was the motivation for a series of experiments to investigate the sensitivity of the method by introducing delivery errors, MLC position and dosimetric overshoot, into the simulated EPID images. The method was found to be sensitive to 1 mm leaf position errors and 10 % overshoot errors.