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.

Commissioning and quality assurance of a novel solution for respiratory-gated PBS proton therapy based on optical tracking of surface markers

We present the commissioning and quality assurance of our clinical protocol for respiratory gating in pencil beam scanning proton therapy for cancer patients with moving targets. In a novel approach, optical tracking has been integrated in the therapy workflow and used to monitor respiratory motion from multiple surrogates, applied on the patients' chest. The gating system was tested under a variety of experimental conditions, specific to proton therapy, to evaluate reaction time and reproducibility of dose delivery control. The system proved to be precise in the application of beam gating and allowed the mitigation of dose distortions even for large (1.4cm) motion amplitudes, provided that adequate treatment windows were selected. The total delivered dose was not affected by the use of gating, with measured integral error within 0.15cGy. Analysing high-resolution images of proton transmission, we observed negligible discrepancies in the geometric location of the dose as a function of the treatment window, with gamma pass rate greater than 95% (2%/2mm) compared to stationary conditions. Similarly, pass rate for the latter metric at the 3%/3mm level was observed above 97% for clinical treatment fields, limiting residual movement to 3mm at end-exhale. These results were confirmed in realistic clinical conditions using an anthropomorphic breathing phantom, reporting a similarly high 3%/3mm pass rate, above 98% and 94%, for regular and irregular breathing, respectively. Finally, early results from periodic QA tests of the optical tracker have shown a reliable system, with small variance observed in static and dynamic measurements.

Anthropomorphic phantom for deformable lung and liver CT and MR imaging for radiotherapy

In this study, a functioning and ventilated anthropomorphic phantom was further enhanced for the purpose of CT and MR imaging of the lung and liver. A deformable lung, including respiratory tract was 3D printed. Within the lung's inner structures is a solid region shaped from a patient's lung tumour and six nitro-glycerine capsules as reference landmarks. The full internal mesh was coated, and the tumour filled, with polyorganosiloxane based gel. A moulded liver was created with an external casing of silicon filled with polyorganosiloxane gel and flexible plastic internal structures. The liver, fitted to the inferior portion of the right lung, moves along with the lung's ventilation. In the contralateral side, a cavity is designed to host a dosimeter, whose motion is correlated to the lung pressure. A 4DCT of the phantom was performed along with static and 4D T1 weighted MR images. The CT Hounsfield units (HU) for the flexible 3D printed material were -600-100 HU (lung and liver structures), for the polyorganosiloxane gel 30-120 HU (lung coating and liver filling) and for the silicon 650-800 HU (liver casing). The MR image intensity units were 0-40, 210-280 and 80-130, respectively. The maximum range of motion in the 4D imaging for the superior lung was 1-3.5 mm and 3.5-8 mm in the inferior portion. The liver motion was 5.5-8.0 mm at the tip and 5.7-10.0 mm at the dome. No measurable drift in motion was observed over a 2 h session and motion was reproducible over three different sessions for sin²(t), sin⁴(t) and a patient-like breathing curve with the interquartile range of amplitudes for all breathing cycles within 0.5 mm. The addition of features within the lung and of a deformable liver will allow the phantom to be used for imaging studies such as validation of 4DMRI and pseudo CT methods.

Technical Note: In silico and experimental evaluation of two leaf-fitting algorithms for MLC tracking based on exposure error and plan complexity

PURPOSE

Multileaf collimator (MLC) tracking is being clinically pioneered to continuously compensate for thoracic and pelvic motion during radiotherapy. The purpose of this work was to characterize the performance of two MLC leaf-fitting algorithms, direct optimization and piecewise optimization, for real-time motion compensation with different plan complexity and tumor trajectories.

METHODS

To test the algorithms, both in silico and phantom experiments were performed. The phantom experiments were performed on a Trilogy Varian linac and a HexaMotion programmable motion platform. High and low modulation VMAT plans for lung and prostate cancer cases were used along with eight patient-measured organ-specific trajectories. For both MLC leaf-fitting algorithms, the plans were run with their corresponding patient trajectories. To compare algorithms, the average exposure errors, i.e., the difference in shape between ideal and fitted MLC leaves by the algorithm, plan complexity and system latency of each experiment were calculated.

RESULTS

Comparison of exposure errors for the in silico and phantom experiments showed minor differences between the two algorithms. The average exposure errors for in silico experiments with low/high plan complexity were 0.66/0.88 cm² for direct optimization and 0.66/0.88 cm² for piecewise optimization, respectively. The average exposure errors for the phantom experiments with low/high plan complexity were 0.73/1.02 cm² for direct and 0.73/1.02 cm² for piecewise optimization, respectively. The measured latency for the direct optimization was 226 ± 10 ms and for the piecewise algorithm was 228 ± 10 ms. In silico and phantom exposure errors quantified for each treatment plan demonstrated that the exposure errors from the high plan complexity (0.96 cm² mean, 2.88 cm² 95% percentile) were all significantly different from the low plan complexity (0.70 cm² mean, 2.18 cm² 95% percentile) (P < 0.001, two-tailed, Mann-Whitney statistical test).

CONCLUSIONS

The comparison between the two leaf-fitting algorithms demonstrated no significant differences in exposure errors, neither in silico nor with phantom experiments. This study revealed that plan complexity impacts the overall exposure errors significantly more than the difference between the algorithms.

Validation of fast motion-including dose reconstruction for proton scanning therapy in the liver

This study validates a method of fast motion-including dose reconstruction for proton pencil beam scanning in the liver. The method utilizes a commercial treatment planning system (TPS) and calculates the delivered dose for any translational 3D target motion. Data from ten liver patients previously treated with photon radiotherapy with intrafraction tumour motion monitoring were used. The dose reconstruction method utilises an in-house developed program to incorporate beam's-eye-view tumour motion by shifting each spot in the opposite direction of the tumour and in-depth motion as beam energy changes for each spot. The doses are then calculated on a single CT phase in the TPS. Two aspects of the dose reconstruction were assessed: (1) The accuracy of reconstruction, by comparing dose reconstructions created using 4DCT motion with ground truth doses obtained by calculating phase specific doses in all 4DCT phases and summing up these partial doses. (2) The error caused by assuming 4DCT motion, by comparing reconstructions with 4DCT motion and actual tumour motion. The CTV homogeneity index (HI) and the root-mean-square (rms) dose error for all dose points receiving  >70%, >80% and  >90% of the prescribed dose were calculated. The dose reconstruction resulted in mean (range) absolute CTV HI errors of 1.0% (0.0-3.0)% and rms dose errors of 2.5% (1.0%-5.3%), 2.1% (0.9%-4.5%), and 1.8% (0.7%-3.7%) for  >70%, >80% and  >90% doses, respectively, when compared with the ground truth. The assumption of 4DCT motion resulted in mean (range) absolute CTV HI errors of 5.9% (0.0-15.0)% and rms dose errors of 6.3% (3.9%-12.6%), 5.9% (3.4%-12.5%), and 5.4% (2.6%-12.1%) for  >70%, >80% and  >90% doses, respectively. The investigated method allows tumour dose reconstruction with the actual tumour motion and results in significantly smaller dose errors than those caused by assuming that motion at treatment is identical to the 4DCT motion.

A dosimetric comparison of real-time adaptive and non-adaptive radiotherapy: A multi-institutional study encompassing robotic, gimbaled, multileaf collimator and couch tracking

PURPOSE

A study of real-time adaptive radiotherapy systems was performed to test the hypothesis that, across delivery systems and institutions, the dosimetric accuracy is improved with adaptive treatments over non-adaptive radiotherapy in the presence of patient-measured tumor motion.

METHODS AND MATERIALS

Ten institutions with robotic(2), gimbaled(2), MLC(4) or couch tracking(2) used common materials including CT and structure sets, motion traces and planning protocols to create a lung and a prostate plan. For each motion trace, the plan was delivered twice to a moving dosimeter; with and without real-time adaptation. Each measurement was compared to a static measurement and the percentage of failed points for γ-tests recorded.

RESULTS

For all lung traces all measurement sets show improved dose accuracy with a mean 2%/2mm γ-fail rate of 1.6% with adaptation and 15.2% without adaptation (p<0.001). For all prostate the mean 2%/2mm γ-fail rate was 1.4% with adaptation and 17.3% without adaptation (p<0.001). The difference between the four systems was small with an average 2%/2mm γ-fail rate of <3% for all systems with adaptation for lung and prostate.

CONCLUSIONS

The investigated systems all accounted for realistic tumor motion accurately and performed to a similar high standard, with real-time adaptation significantly outperforming non-adaptive delivery methods.

Prostate motion during radiotherapy of prostate cancer patients with and without application of a hydrogel spacer: a comparative study

Background and purpose

The use of a tissue expander (hydrogel) for sparing of the rectum from increased irradiation during prostate radiotherapy is becoming popular. The goal of this study is to investigate the effect of a tissue expander (hydrogel) on the intrafraction prostate motion during radiotherapy.

Methods and material

Real time prostate motion was analysed for 26 patients and 742 fractions; 12 patients with and 14 patients without hydrogel (SpaceOAR™). The intra-fraction motion was quantified and compared between the two groups.

Results

The average (±standard deviation) of the mean motion during the treatment for patients with and without hydrogel was 1.5 (±0.8 mm) and 1.1 (±0.9 mm) respectively (p < 0.05). The average time of motion >3 mm for patients with and without hydrogel was 7.7 % (±1.1 %) and 4.5 % (±0.9 %) respectively (p > 0.05). The hydrogel age, fraction number and treatment time were found to have no effect (R² < 0.05) on the prostate motion.

Conclusions

Differences in intrafraction motion in patients with hydrogel and without hydrogel were within measurement uncertainty (<1 mm). This result confirms that the addition of a spacer does not negate the need for intrafraction motion management if clinically indicated.

Quality assurance for the clinical implementation of kilovoltage intrafraction monitoring for prostate cancer VMAT

PURPOSE

Kilovoltage intrafraction monitoring (KIM) is a real-time 3D tumor monitoring system for cancer radiotherapy. KIM uses the commonly available gantry-mounted x-ray imager as input, making this method potentially more widely available than dedicated real-time 3D tumor monitoring systems. KIM is being piloted in a clinical trial for prostate cancer patients treated with VMAT (NCT01742403). The purpose of this work was to develop clinical process and quality assurance (QA) practices for the clinical implementation of KIM.

METHODS

Informed by and adapting existing guideline documents from other real-time monitoring systems, KIM-specific QA practices were developed. The following five KIM-specific QA tests were included: (1) static localization accuracy, (2) dynamic localization accuracy, (3) treatment interruption accuracy, (4) latency measurement, and (5) clinical conditions accuracy. Tests (1)-(4) were performed using KIM to measure static and representative patient-derived prostate motion trajectories using a 3D programmable motion stage supporting an anthropomorphic phantom with implanted gold markers to represent the clinical treatment scenario. The threshold for system tolerable latency is <1 s. The tolerances for all other tests are that both the mean and standard deviation of the difference between the programmed trajectory and the measured data are <1 mm. The (5) clinical conditions accuracy test compared the KIM measured positions with those measured by kV/megavoltage (MV) triangulation from five treatment fractions acquired in a previous pilot study.

RESULTS

For the (1) static localization, (2) dynamic localization, and (3) treatment interruption accuracy tests, the mean and standard deviation of the difference are <1.0 mm. (4) The measured latency is 350 ms. (5) For the tests with previously acquired patient data, the mean and standard deviation of the difference between KIM and kV/MV triangulation are <1.0 mm.

CONCLUSIONS

Clinical process and QA practices for the safe clinical implementation of KIM, a novel real-time monitoring system using commonly available equipment, have been developed and implemented for prostate cancer VMAT.

The first clinical implementation of electromagnetic transponder-guided MLC tracking

PURPOSE

We report on the clinical process, quality assurance, and geometric and dosimetric results of the first clinical implementation of electromagnetic transponder-guided MLC tracking which occurred on 28 November 2013 at the Northern Sydney Cancer Centre.

METHODS

An electromagnetic transponder-based positioning system (Calypso) was modified to send the target position output to in-house-developed MLC tracking code, which adjusts the leaf positions to optimally align the treatment beam with the real-time target position. Clinical process and quality assurance procedures were developed and performed. The first clinical implementation of electromagnetic transponder-guided MLC tracking was for a prostate cancer patient being treated with dual-arc VMAT (RapidArc). For the first fraction of the first patient treatment of electromagnetic transponder-guided MLC tracking we recorded the in-room time and transponder positions, and performed dose reconstruction to estimate the delivered dose and also the dose received had MLC tracking not been used.

RESULTS

The total in-room time was 21 min with 2 min of beam delivery. No additional time was needed for MLC tracking and there were no beam holds. The average prostate position from the initial setup was 1.2 mm, mostly an anterior shift. Dose reconstruction analysis of the delivered dose with MLC tracking showed similar isodose and target dose volume histograms to the planned treatment and a 4.6% increase in the fractional rectal V60. Dose reconstruction without motion compensation showed a 30% increase in the fractional rectal V60 from that planned, even for the small motion.

CONCLUSIONS

The real-time beam-target correction method, electromagnetic transponder-guided MLC tracking, has been translated to the clinic. This achievement represents a milestone in improving geometric and dosimetric accuracy, and by inference treatment outcomes, in cancer radiotherapy.