Real-time motion-including dose estimation of simulated multi-leaf collimator-tracked magnetic resonance-guided radiotherapy.

BACKGROUND: Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking. PURPOSE: We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions. METHODS: A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation. RESULTS: On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively. CONCLUSIONS: Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose.

Experimental comparison of linear regression and LSTM motion prediction models for MLC-tracking on an MRI-linac.

BACKGROUND: Magnetic resonance imaging (MRI)-guided radiotherapy with multileaf collimator (MLC)-tracking is a promising technique for intra-fractional motion management, achieving high dose conformality without prolonging treatment times. To improve beam-target alignment, the geometric error due to system latency should be reduced by using temporal prediction. PURPOSE: To experimentally compare linear regression (LR) and long-short-term memory (LSTM) motion prediction models for MLC-tracking on an MRI-linac using multiple patient-derived traces with different complexities. METHODS: Experiments were performed on a prototype 1.0 T MRI-linac capable of MLC-tracking. A motion phantom was programmed to move a target in superior-inferior (SI) direction according to eight lung cancer patient respiratory motion traces. Target centroid positions were localized from sagittal 2D cine MRIs acquired at 4 Hz using a template matching algorithm. The centroid positions were input to one of four motion prediction models. We used (1) a LSTM network which had been optimized in a previous study on patient data from another cohort (offline LSTM). We also used (2) the same LSTM model as a starting point for continuous re-optimization of its weights during the experiment based on recent motion (offline+online LSTM). Furthermore, we implemented (3) a continuously updated LR model, which was solely based on recent motion (online LR). Finally, we used (4) the last available target centroid without any changes as a baseline (no-predictor). The predictions of the models were used to shift the MLC aperture in real-time. An electronic portal imaging device (EPID) was used to visualize the target and MLC aperture during the experiments. Based on the EPID frames, the root-mean-square error (RMSE) between the target and the MLC aperture positions was used to assess the performance of the different motion predictors. Each combination of motion trace and prediction model was repeated twice to test stability, for a total of 64 experiments. RESULTS: The end-to-end latency of the system was measured to be (389 ± 15) ms and was successfully mitigated by both LR and LSTM models. The offline+online LSTM was found to outperform the other models for all investigated motion traces. It obtained a median RMSE over all traces of (2.8 ± 1.3) mm, compared to the (3.2 ± 1.9) mm of the offline LSTM, the (3.3 ± 1.4) mm of the online LR and the (4.4 ± 2.4) mm when using the no-predictor. According to statistical tests, differences were significant (p-value <0.05) among all models in a pair-wise comparison, but for the offline LSTM and online LR pair. The offline+online LSTM was found to be more reproducible than the offline LSTM and the online LR with a maximum deviation in RMSE between two measurements of 10%. CONCLUSIONS: This study represents the first experimental comparison of different prediction models for MRI-guided MLC-tracking using several patient-derived respiratory motion traces. We have shown that among the investigated models, continuously re-optimized LSTM networks are the most promising to account for the end-to-end system latency in MRI-guided radiotherapy with MLC-tracking.

Evaluation of real-time tumor contour prediction using LSTM networks for MR-guided radiotherapy.

BACKGROUND AND PURPOSE: Magnetic resonance imaging guided radiotherapy (MRgRT) with deformable multileaf collimator (MLC) tracking would allow to tackle both rigid displacement and tumor deformation without prolonging treatment. However, the system latency must be accounted for by predicting future tumor contours in real-time. We compared the performance of three artificial intelligence (AI) algorithms based on long short-term memory (LSTM) modules for the prediction of 2D-contours 500ms into the future. MATERIALS AND METHODS: Models were trained (52 patients, 3.1h of motion), validated (18 patients, 0.6h) and tested (18 patients, 1.1h) with cine MRs from patients treated at one institution. Additionally, we used three patients (2.9h) treated at another institution as second testing set. We implemented 1) a classical LSTM network (LSTM-shift) predicting tumor centroid positions in superior-inferior and anterior-posterior direction which are used to shift the last observed tumor contour. The LSTM-shift model was optimized both in an offline and online fashion. We also implemented 2) a convolutional LSTM model (ConvLSTM) to directly predict future tumor contours and 3) a convolutional LSTM combined with spatial transformer layers (ConvLSTM-STL) to predict displacement fields used to warp the last tumor contour. RESULTS: The online LSTM-shift model was found to perform slightly better than the offline LSTM-shift and significantly better than the ConvLSTM and ConvLSTM-STL. It achieved a 50% Hausdorff distance of 1.2mm and 1.0mm for the two testing sets, respectively. Larger motion ranges were found to lead to more substantial performance differences across the models. CONCLUSION: LSTM networks predicting future centroids and shifting the last tumor contour are the most suitable for tumor contour prediction. The obtained accuracy would allow to reduce residual tracking errors during MRgRT with deformable MLC-tracking.

Rapid distortion correction enables accurate magnetic resonance imaging-guided real-time adaptive radiotherapy.

Background and purpose: Magnetic resonance imaging (MRI)-Linac systems combine simultaneous MRI with radiation delivery, allowing treatments to be guided by anatomically detailed, real-time images. However, MRI can be degraded by geometric distortions that cause uncertainty between imaged and actual anatomy. In this work, we develop and integrate a real-time distortion correction method that enables accurate real-time adaptive radiotherapy. Materials and methods: The method was based on the pre-treatment calculation of distortion and the rapid correction of intrafraction images. A motion phantom was set up in an MRI-Linac at isocentre (P0 ), the edge (P 1) and just outside (P 2) the imaging volume. The target was irradiated and tracked during real-time adaptive radiotherapy with and without the distortion correction. The geometric tracking error and latency were derived from the measurements of the beam and target positions in the EPID images. Results: Without distortion correction, the mean geometric tracking error was 1.3 mm at P 1 and 3.1 mm at P 2. When distortion correction was applied, the error was reduced to 1.0 mm at P 1 and 1.1 mm at P 2. The corrected error was similar to an error of 0.9 mm at P0 where the target was unaffected by distortion indicating that this method has accurately accounted for distortion during tracking. The latency was 319 ± 12 ms without distortion correction and 335 ± 34 ms with distortion correction. Conclusions: We have demonstrated a real-time distortion correction method that maintains accurate radiation delivery to the target, even at treatment locations with large distortion.

Optimising multi-target multileaf collimator tracking using real-time dose for locally advanced prostate cancer patients.

Objective. The accuracy of radiotherapy for patients with locally advanced cancer is compromised by independent motion of multiple targets. To date, MLC tracking approaches have used 2D geometric optimisation where the MLC aperture shape is simply translated to correspond to the target's motion, which results in sub-optimal delivered dose. To address this limitation, a dose-optimised multi-target MLC tracking method was developed and evaluated through simulated locally advanced prostate cancer treatments.Approach. A dose-optimised multi-target tracking algorithm that adapts the MLC aperture to minimise 3D dosimetric error was developed for moving prostate and static lymph node targets. A fast dose calculation algorithm accumulated the planned dose to the prostate and lymph node volumes during treatment in real time, and the MLC apertures were recalculated to minimise the difference between the delivered and planned dose with the included motion. Dose-optimised tracking was evaluated by simulating five locally advanced prostate plans and three prostate motion traces with a relative interfraction displacement. The same simulations were performed using geometric-optimised tracking and no tracking. The dose-optimised, geometric-optimised, and no tracking results were compared with the planned doses using a 2%/2 mmγcriterion.Main results. The mean dosimetric error was lowest for dose-optimised MLC tracking, withγ-failure rates of 12% ± 8.5% for the prostate and 2.2% ± 3.2% for the nodes. Theγ-failure rates for geometric-optimised MLC tracking were 23% ± 12% for the prostate and 3.6% ± 2.5% for the nodes. When no tracking was used, theγ-failure rates were 37% ± 28% for the prostate and 24% ± 3.2% for the nodes.Significance. This study developed a dose-optimised multi-target MLC tracking method that minimises the difference between the planned and delivered doses in the presence of intrafraction motion. When applied to locally advanced prostate cancer, dose-optimised tracking showed smaller errors than geometric-optimised tracking and no tracking for both the prostate and nodes.

First experimental demonstration of VMAT combined with MLC tracking for single and multi fraction lung SBRT on an MR-linac.

BACKGROUND AND PURPOSE: VMAT is not currently available on MR-linacs but could maximize plan conformality. To mitigate respiration without compromising delivery efficiency, MRI-guided MLC tumour tracking was recently developed for the 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden) in combination with IMRT. Here, we provide a first experimental demonstration of VMAT + MLC tracking for several lung SBRT indications. MATERIALS AND METHODS: We created central patient and phantom VMAT plans (8×7.5 Gy, 2 arcs) and we created peripheral phantom plans (3×18 & 1×34 Gy, 4 arcs). A motion phantom mimicked subject-recorded respiratory motion (A‾=11 mm, f‾=0.33 Hz, drift‾=0.3 mm/min). This was monitored using 2D-cine MRI at 4 Hz to continuously realign the beam with the target. VMAT + MLC tracking performance was evaluated using 2D film dosimetry and a novel motion-encoded and time-resolved pseudo-3D dosimetry approach. RESULTS: We found an MLC leaf and jaw end-to-end latency of 328.05(±3.78) ms and 317.33(±4.64) ms, which was mitigated by a predictor. The VMAT plans required maximum MLC speeds of 12.1 cm/s and MLC tracking superimposed an additional 1.48 cm/s. A local 2%/1 mm gamma analysis with a static measurement as reference, revealed pass-rates of 28-46% without MLC tracking and 88-100% with MLC tracking for the 2D film analysis. Similarly, the pseudo-3D gamma passing-rates increased from 22-77% to 92-100%. The dose area histograms showed that MLC tracking increased the GTV D98% by 5-20% and the PTV D95% by 7-24%, giving similar target coverage as their respective static reference. CONCLUSION: MRI-guided VMAT + MLC tracking is technically feasible on the MR-linac and results in highly conformal dose distribution.

A hybrid 2D/4D-MRI methodology using simultaneous multislice imaging for radiotherapy guidance.

PURPOSE: Respiratory motion management is important in abdominothoracic radiotherapy. Fast imaging of the tumor can facilitate multileaf collimator (MLC) tracking that allows for smaller treatment margins, while repeatedly imaging the full field-of-view is necessary for 4D dose accumulation. This study introduces a hybrid 2D/4D-MRI methodology that can be used for simultaneous MLC tracking and dose accumulation on a 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden). METHODS: We developed a hybrid 2D/4D-MRI methodology that uses a simultaneous multislice (SMS) accelerated MRI sequence, which acquires two coronal slices simultaneously and repeatedly cycles through slice positions over the image volume. As a result, the fast 2D imaging can be used prospectively for MLC tracking and the SMS slices can be sorted retrospectively into respiratory-correlated 4D-MRIs for dose accumulation. Data were acquired in five healthy volunteers with an SMS-bTFE and SMS-TSE MRI sequence. For each sequence, a prebeam dataset and a beam-on dataset were acquired simulating the two phases of MR-linac treatments. Prebeam data were used to generate a 4D-based motion model and a reference mid-position volume, while beam-on data were used for real-time motion extraction and reconstruction of beam-on 4D-MRIs. In addition, an in-silico computational phantom was used for validation of the hybrid 2D/4D-MRI methodology. MLC tracking experiments were performed with the developed methodology, for which real-time SMS data reconstruction was enabled on the scanner. A 15-beam 8× 7.5 Gy intensity-modulated radiotherapy plan for lung stereotactic body radiotherapy with isotropic 3 mm GTV-to-PTV margins was created. Dosimetry experiments were performed using a 4D motion phantom. The latency between target motion and updating the radiation beam was determined and compensated. Local gamma analyses were performed to quantify dose differences compared to a static reference delivery, and dose area histograms (DAHs) were used to quantify the GTV and PTV coverage. RESULTS: In-vivo data acquisition and MLC tracking experiments were successfully performed with the developed hybrid 2D/4D-MRI methodology. Real-time liver-lung interface motion estimation had a Pearson's correlation of 0.996 (in-vivo) and 0.998 (in-silico). A median (5th-95th percentile) error of 0.0 (-0.9 to 0.7) mm and 0.0 (-0.2 to 0.2) mm was found for real-time motion estimation for in-vivo and in-silico, respectively. Target motion prediction beyond the liver-lung interface had a median root mean square error of 1.6 mm (in-vivo) and 0.5 mm (in-silico). Beam-on 4D MRI reconstruction required a median amount of data equal to an acquisition time of 2:21-3:17 min, which was 20% less data compared to the prebeam-derived 4D-MRI. System latency was reduced from 501 ± 12 ms to -1 ± 3 ms (SMS-TSE) and from 398 ± 10 ms to -10 ± 4 ms (SMS-bTFE) by a linear regression prediction filter. The local gamma analysis agreed within - 3.8 % $-3.8\%$ to 3.3% (SMS-bTFE) and - 5.3 % $-5.3\%$ to 10% (SMS-TSE) with a reference MRI sequence. The DAHs revealed a relative D 98 % $D_{98\%}$ GTV coverage between 97% and 100% (SMS-bTFE) and 100% and 101% (SMS-TSE) compared to the static reference. CONCLUSIONS: The presented 2D/4D-MRI methodology demonstrated the potential for accurately extracting real-time motion for MLC tracking in abdominothoracic radiotherapy, while simultaneously reconstructing contiguous respiratory-correlated 4D-MRIs for dose accumulation.

Experimental investigation of dynamic real-time rotation-including dose reconstruction during prostate tracking radiotherapy.

BACKGROUND: Hypofractionation in prostate radiotherapy is of increasing interest. Steep dose gradients and a large weight on each individual fraction emphasize the need for motion management. Real-time motion management techniques such as multileaf collimator (MLC) tracking or couch tracking typically adjust for translational motion while rotations remain uncompensated with unknown dosimetric impact. PURPOSE: The purpose of this study is to demonstrate and validate dynamic real-time rotation-including dose reconstruction during radiotherapy experiments with and without MLC and couch tracking. METHODS: Real-time dose reconstruction was performed using the in-house developed software DoseTracker. DoseTracker receives streamed target positions and accelerator parameters during treatment delivery and uses a pencil beam algorithm with water density assumption to reconstruct the dose in a moving target. DoseTracker's ability to reconstruct motion-induced dose errors in a dynamically rotating and translating target was investigated during three different scenarios: (1) no motion compensation and translational motion correction with (2) MLC tracking and (3) couch tracking. In each scenario, dose reconstruction was performed online and in real time during delivery of two dual-arc volumetric-modulated arc therapy prostate plans with a prescribed fraction dose of 7 Gy to the prostate and simultaneous intraprostatic lesion boosts with doses of at least 8 Gy, but up to 10 Gy as long as the organs at risk dose constraints were fulfilled. The plans were delivered to a pelvis phantom that replicated three patient-measured motion traces using a rotational insert with 21 layers of EBT3 film spaced 2.5 mm apart. DoseTracker repeatedly calculated the actual motion-including dose increment and the planned static dose increment since the last calculation in 84 500 points in the film stack. The experiments were performed with a TrueBeam accelerator with MLC and couch tracking based on electromagnetic transponders embedded in the film stack. The motion-induced dose error was quantified as the difference between the final cumulative dose with motion and without motion using the 2D 2%/2 mm γ-failure rate and the difference in dose to 95% of the clinical target volume (CTV ΔD95% ) and the gross target volume (GTV ΔD95% ) as well as the difference in dose to 0.1 cm3 of the urethra, bladder, and rectum (ΔD0.1CC ). The motion-induced errors were compared between dose reconstructions and film measurements. RESULTS: The dose was reconstructed in all calculation points at a mean frequency of 4.7 Hz. The root-mean-square difference between real-time reconstructed and film-measured motion-induced errors was 3.1%-points (γ-failure rate), 0.13 Gy (CTV ΔD95% ), 0.23 Gy (GTV ΔD95% ), 0.19 Gy (urethra ΔD0.1CC ), 0.09 Gy (bladder ΔD0.1CC ), and 0.07 Gy (rectum ΔD0.1CC ). CONCLUSIONS: In a series of phantom experiments, online real-time rotation-including dose reconstruction was performed for the first time. The calculated motion-induced errors agreed well with film measurements. The dose reconstruction provides a valuable tool for monitoring dose delivery and investigating the efficacy of advanced motion-compensation techniques in the presence of translational and rotational motion.

First experimental exploration of real-time cardiorespiratory motion management for future stereotactic arrhythmia radioablation treatments on the MR-linac.

Objective.Stereotactic arrhythmia radioablation (STAR) is a novel, non-invasive treatment for refractory ventricular tachycardia (VT). The VT isthmus is subject to both respiratory and cardiac motion. Rapid cardiac motion presents a unique challenge. In this study, we provide first experimental evidence for real-time cardiorespiratory motion-mitigated MRI-guided STAR on the 1.5 T Unity MR-linac (Elekta AB, Stockholm, Sweden) aimed at simultaneously compensating cardiac and respiratory motions.Approach.A real-time cardiorespiratory motion-mitigated radiotherapy workflow was developed on the Unity MR-linac in research mode. A 15-beam intensity-modulated radiation therapy treatment plan (1 × 25 Gy) was created in Monaco v.5.40.01 (Elekta AB) for the Quasar MRI4Dphantom (ModusQA, London, ON). A film dosimetry insert was moved by combining either artificial (cos4, 70 bpm, 10 mm peak-to-peak) or subject-derived (59 average bpm, 15.3 mm peak-to-peak) cardiac motion with respiratory (sin, 12 bpm, 20 mm peak-to-peak) motion. A balanced 2D cine MRI sequence (13 Hz, field-of-view = 400 × 207 mm2, resolution = 3 × 3 × 15 mm3) was developed to estimate cardiorespiratory motion. Cardiorespiratory motion was estimated by rigid registration and then deconvoluted into cardiac and respiratory components. For beam gating, the cardiac component was used, whereas the respiratory component was used for MLC-tracking. In-silico dose accumulation experiments were performed on three patient data sets to simulate the dosimetric effect of cardiac motion on VT targets.Main results.Experimentally, a duty cycle of 57% was achieved when simultaneously applying respiratory MLC-tracking and cardiac gating. Using film, excellent agreement was observed compared to a static reference delivery, resulting in a 1%/1 mm gamma pass rate of 99%. The end-to-end gating latency was 126 ms on the Unity MR-linac. Simulations showed that cardiac motion decreased the target's D98% dose between 0.1 and 1.3 Gy, with gating providing effective mitigation.Significance.Real-time MRI-guided cardiorespiratory motion management greatly reduces motion-induced dosimetric uncertainty and warrants further research and development for potential future use in STAR.

First experimental evaluation of multi-target multileaf collimator tracking during volumetric modulated arc therapy for locally advanced prostate cancer.

PURPOSE: Locally advanced and oligometastatic cancer patients require radiotherapy treatment to multiple independently moving targets. There is no existing commercial solution that can simultaneously track and treat multiple targets. This study experimentally implemented and evaluated a real-time multi-target tracking system for locally advanced prostate cancer. METHODS: Real-time multi-target MLC tracking was integrated with 3D x-ray image guidance on a standard linac. Three locally advanced prostate cancer treatment plans were delivered to a static lymph node phantom and dynamic prostate phantom that reproduced three prostate trajectories. Treatments were delivered using multi-target MLC tracking, single-target MLC tracking, and no tracking. Doses were measured using Gafchromic film placed in the dynamic and static phantoms. Dosimetric error was quantified by the 2%/2 mm gamma failure rate. Geometric error was evaluated as the misalignment between target and aperture positions. The multi-target tracking system latency was measured. RESULTS: The mean (range) gamma failure rates for the prostate and lymph nodes, were 18.6% (5.2%, 28.5%) and 7.5% (1.1%, 13.7%) with multi-target tracking, 7.9% (0.7%, 15.4%) and 37.8% (18.0%, 57.9%) with single-target tracking, and 38.1% (0.6%, 75.3%) and 37.2% (29%, 45.3%) without tracking. Multi-target tracking had the lowest geometric error with means and standard deviations within 0.2 ± 1.5 for the prostate and 0.0 ± 0.3 mm for the lymph nodes. The latency was 730 ± 20 ms. CONCLUSION: This study presented the first experimental implementation of multi-target tracking to independently track prostate and lymph node displacement during VMAT. Multi-target tracking reduced dosimetric and geometric errors compared to single-target tracking and no tracking.

Dosimetric evaluation of MRI-guided multi-leaf collimator tracking and trailing for lung stereotactic body radiation therapy.

PURPOSE: The treatment margins for lung stereotactic body radiotherapy (SBRT) are often large to cover the tumor excursions resulting from respiration, such that underdosage of the tumor can be avoided. Magnetic resonance imaging (MRI)-guided multi-leaf collimator (MLC) tracking can potentially reduce the influence of respiration to allow for smaller treatment margins. However, tracking is accompanied by system latency that may induce residual tracking errors. Alternatively, a simpler mid-position delivery combined with trailing can be used. Trailing reduces influences of respiration by compensating for baseline motion, to potentially improve target coverage. In this study, we aim to show the feasibility of MRI-guided tracking and trailing to reduce influences of respiration during lung SBRT. METHODS: We implemented MRI-guided tracking on the MR-linac using an Elekta research tracking interface to track tumor motion during intensity modulated radiotherapy (IMRT). A Quasar MRI 4 D phantom was used to generate Lujan motion ( cos 4 , 4 s period, 20 mm peak-to-peak amplitude) with and without 1.0 mm/min cranial drift. Phantom tumor positions were estimated from sagittal 2D cine-MRI (4 or 8 Hz) using cross-correlation-based template matching. To compensate the anticipated system latency, a linear ridge regression predictor was optimized for online MRI by comparing two predictor training approaches: training on multiple traces and training on a single trace. We created 15-beam clinical-grade lung SBRT plans for central targets (8 × 7.5 Gy) and peripheral targets (3 × 18 Gy) with different PTV margins for mid-position based motion management (3-5 mm) and for MLC tracking (3 mm). We used a film insert with a 3 cm spherical target to measure the spatial distribution and quantity of the delivered dose. A 1%/1 mm local gamma-analysis quantified dose differences between motion management strategies and reference cases. Additionally, a dose area histogram (DAH) revealed the target coverage relative to the reference scenario. RESULTS: The prediction filter gain was on average 25% when trained on multiple traces and 44% when trained on a single trace. The filter reduced system latency from 313 ± 2 ms to 0 ± 5 ms for 4 Hz imaging and from 215 ± 3 ms to 3 ± 3 ms for 8 Hz. The local gamma analysis for the central delivery showed that tracking improved the gamma pass-rate from 23% to 96% for periodic motion and from 14% to 93% when baseline drift was applied. For the peripheral delivery during periodic motion, delivery pass-rates improved from 22% to 93%. Comparing mid-position delivery to trailing for periodic+drift motion increased the local gamma pass rate from 15% to 98% for a central delivery and from 8% to 98% for a peripheral delivery. Furthermore, the DAHs revealed a relative D 98 % GTV coverage of 101% and 97% compared to the reference scenario for, respectively, central and peripheral tracking of periodic+drift motion. For trailing, a relative D 98 % of 99% for central and 98% for peripheral trailing was found. CONCLUSIONS: We provided a first experimental demonstration of the technical feasibility of MRI-guided MLC tracking and trailing for central and peripheral lung SBRT. Tracking maximizes the sparing of healthy tissue, while trailing is highly effective in mitigating baseline motion.

Adapting to the motion of multiple independent targets using multileaf collimator tracking for locally advanced prostate cancer: Proof of principle simulation study.

PURPOSE: For patients with locally advanced cancer, multiple targets are treated simultaneously with radiotherapy. Differential motion between targets can compromise the treatment accuracy, yet there are currently no methods able to adapt to independent target motion. This study developed a multileaf collimator (MLC) tracking algorithm for differential motion adaptation and evaluated it in simulated treatments of locally advanced prostate cancer. METHODS: A multi-target MLC tracking algorithm was developed that consisted of three steps: (a) dividing the MLC aperture into two possibly overlapping sections assigned to the prostate and lymph nodes, (b) calculating the ideally shaped MLC aperture as a union of the individually translated sections, and (c) fitting the MLC positions to the ideal aperture shape within the physical constraints of the MLC leaves. The multi-target tracking method was evaluated and compared with two existing motion management methods: single-target tracking and no tracking. Treatment simulations of six locally advanced prostate cancer patients with three prostate motion traces were performed for all three motion adaptation methods. The geometric error for each motion adaptation method was calculated using the area of overexposure and underexposure of each field. The dosimetric error was estimated by calculating the dose delivered to the prostate, lymph nodes, bladder, rectum, and small bowel using a motion-encoded dose reconstruction method. RESULTS: Multi-target MLC tracking showed an average improvement in geometric error of 84% compared to single-target tracking, and 83% compared to no tracking. Multi-target tracking maintained dose coverage to the prostate clinical target volume (CTV) D98% and planning target volume (PTV) D95% to within 4.8% and 3.9% of the planned values, compared to 1.4% and 0.7% with single-target tracking, and 20.4% and 31.8% with no tracking. With multi-target tracking, the node CTV D95%, PTV D90%, and gross tumor volume (GTV) D95% were within 0.3%, 0.6%, and 0.3% of the planned values, compared to 9.1%, 11.2%, and 21.1% for single-target tracking, and 0.8%, 2.0%, and 3.2% with no tracking. The small bowel V57% was maintained within 0.2% to the plan using multi-target tracking, compared to 8% and 3.5% for single-target tracking and no tracking, respectively. Meanwhile, the bladder and rectum V50% increased by up to 13.6% and 5.2%, respectively, using multi-target tracking, compared to 2.7% and 1.9% for single-target tracking, and 11.2% and 11.5% for no tracking. CONCLUSIONS: A multi-target tracking algorithm was developed and tracked the prostate and lymph nodes independently during simulated treatments. As the algorithm optimizes for target coverage, tracking both targets simultaneously may increase the dose delivered to the organs at risk.

Dosimetric impact of intrafraction prostate rotation and accuracy of gating, multi-leaf collimator tracking and couch tracking to manage rotation: An end-to-end validation using volumetric film measurements.

BACKGROUND AND PURPOSE: Both gating and tracking can mitigate the deteriorating dosimetric impact of intrafraction translation during prostate stereotactic body radiotherapy (SBRT). However, their ability to manage intrafraction rotation has not yet been thoroughly investigated. The dosimetric accuracy of gating, MLC tracking and couch tracking to manage intrafraction prostate rotation was investigated. MATERIALS AND METHODS: Treatment plans for end-to-end tests of prostate SBRT with focal boosting were generated for a dynamic anthropomorphic pelvis phantom. The phantom applied internal lateral rotation (up to 25°) and coupled vertical and longitudinal translation of a radiochromic film stack that was used for dose measurements. Dose was delivered for each plan while the phantom applied motion according to three typical prostate motion traces without compensation (i), with gating (ii), with MLC tracking (iii) or with couch tracking (iv). Measured doses for the four motion compensation strategies were compared with the planned dose in terms of γ-index analysis, target coverage and organs at risk (OAR) sparing. RESULTS: Intrafraction rotation reduced the 3%(global)/2mm γ-index passing rate (γPR) for the prostate target volume by median (range) -33.2% (-68.6%, -4.1%) when no motion compensation was applied. The use of motion compensation improved the γPR by 13.2% (-0.4%, 32.9%) for gating, by 6.0% (-0.8%, 27.7%) for MLC tracking and by 11.1% (1.2%, 22.9%) for couch tracking. The three compensation techniques improved the target coverage in most cases. Gating showed better OAR sparing than MLC tracking or couch tracking. CONCLUSIONS: Compensation of intrafraction prostate rotation with gating, MLC tracking and couch tracking was investigated experimentally for the first time. All three techniques improved the dosimetric accuracy, but residual motion-related dose errors remained due to the lack of rotation correction.

MLC tracking for lung SABR is feasible, efficient and delivers high-precision target dose and lower normal tissue dose.

BACKGROUND AND PURPOSE: The purpose of this work is to present the clinical experience from the first-in-human trial of real-time tumor targeting via MLC tracking for stereotactic ablative body radiotherapy (SABR) of lung lesions. METHODS AND MATERIALS: Seventeen patients with stage 1 non-small cell lung cancer (NSCLC) or lung metastases were included in a study of electromagnetic transponder-guided MLC tracking for SABR (NCT02514512). Patients had electromagnetic transponders inserted near the tumor. An MLC tracking SABR plan was generated with planning target volume (PTV) expanded 5 mm from the end-exhale gross tumor volume (GTV). A clinically approved comparator plan was generated with PTV expanded 5 mm from a 4DCT-derived internal target volume (ITV). Treatment was delivered using a standard linear accelerator to continuously adapt the MLC based on transponder motion. Treated volumes and reconstructed delivered dose were compared between MLC tracking and comparator ITV-based treatment. RESULTS: All seventeen patients were successfully treated with MLC tracking (70 successful fractions). MLC tracking treatment delivery time averaged 8 minutes. The time from the start of CBCT to the end of treatment averaged 22 minutes. The MLC tracking PTV for 16/17 patients was smaller than the ITV-based PTV (range -1.6% to 44% reduction, or -0.6 to 18 cc). Reductions in mean lung dose (27 cGy) and V20Gy (50 cc) were statistically significant (p < 0.02). Reconstruction of treatment doses confirmed a statistically significant improvement in delivered GTV D98% (p < 0.05) from planned dose compared with the ITV-based plans. CONCLUSION: The first treatments with lung MLC tracking have been successfully performed in seventeen SABR patients. MLC tracking for lung SABR is feasible, efficient and delivers high-precision target dose and lower normal tissue dose.

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 (RMSError ). For simulated field size errors, the RMSError was 0.47 cm2 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 cm2 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.

First experimental investigation of simultaneously tracking two independently moving targets on an MRI-linac using real-time MRI and MLC tracking.

PURPOSE: High quality radiotherapy is challenging in cases where multiple targets with independent motion are simultaneously treated. A real-time tumor tracking system that can simultaneously account for the motion of two targets was developed and characterized. METHODS: The multitarget tracking system was implemented on a magnetic resonance imaging (MRI)-linac and utilized multi-leaf collimator (MLC) tracking to adapt the radiation beam to phantom targets reproducing motion with prostate and lung motion traces. Multitarget tracking consisted of three stages: (a) pretreatment aperture segmentation where the treatment aperture was divided into segments corresponding to each target, (b) MR imaging where the positions of the two targets were localized, and (c) MLC tracking where an updated treatment aperture was calculated. Electronic portal images (EPID) acquired during irradiation were analyzed to characterize geometric uncertainty and tracking latency. RESULTS: Multitarget MLC tracking effectively accounted for the motion of both targets during treatment. The root-mean-square error between the centers of the targets and the centers of the corresponding MLC leaves were reduced from 5.5 mm without tracking to 2.7 mm with tracking for lung motion traces and reduced from 4.2 to 1.4 mm for prostate motion traces. The end-to-end latency of tracking was measured to be 328 ± 44 ms. CONCLUSIONS: We have demonstrated the first experimental implementation of MLC tracking for multiple targets having independent motion. This technology takes advantage of the imaging capabilities of MRI-linacs and would allow treatment margins to be reduced in cases where multiple targets are simultaneously treated.

Couch and multileaf collimator tracking: A clinical feasibility study for pancreas and liver treatment.

PURPOSE: Real-time tumor tracking through active correction by the multileaf collimator or treatment couch offers a promising strategy to mitigate delivery uncertainty due to intrafractional tumor motion. This study evaluated the performance of MLC and couch tracking using the prototype iTools Tracking system in TrueBeam Developer Mode and the application for abdominal cancer treatments. METHODS: Experiments were carried out using a phantom with embedded Calypso transponders and a motion simulation platform. Geometric evaluations were performed using a circular conformal field with sinusoidal traces and pancreatic tumor motion traces. Geometric tracking accuracy was retrospectively calculated by comparing the compensational MLC or couch motion extracted from machine log files to the target motion reconstructed from real-time MV and kV images. Dosimetric tracking accuracy was measured with radiochromic films using clinical abdominal VMAT plans and pancreatic tumor traces. RESULTS: Geometrically, the root-mean-square errors for MLC tracking were 0.5 and 1.8 mm parallel and perpendicular to leaf travel direction, respectively. Couch tracking, in contrast, showed an average of 0.8 mm or less geometric error in all directions. Dosimetrically, both MLC and couch tracking reduced motion-induced local dose errors compared to no tracking. Evaluated with five pancreatic tumor motion traces, the average 2%/2 mm global gamma pass rate of eight clinical abdominal VMAT plans was 67.4% (range: 26.4%-92.7%) without tracking, which was improved to 86.0% (range: 67.9%-95.6%) with MLC tracking, and 98.1% (range: 94.9%-100.0%) with couch tracking. In 16 out of 40 deliveries with different plans and motion traces, MLC tracking did not achieve clinically acceptable dosimetric accuracy with 3%/3mm gamma pass rate below 95%. CONCLUSIONS: This study demonstrated the capability of MLC and couch tracking to reduce motion-induced dose errors in abdominal cases using a prototype tracking system. Clinically significant dose errors were observed with MLC tracking for certain plans which could be attributed to the inferior MLC tracking accuracy in the direction perpendicular to leaf travel, as well as the interplay between motion tracking and plan delivery for highly modulated plans. Couch tracking outperformed MLC tracking with consistently high dosimetric accuracy in all plans evaluated, indicating its clinical potential in the treatment of abdominal cancers.

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 cm2 for direct optimization and 0.66/0.88 cm2 for piecewise optimization, respectively. The average exposure errors for the phantom experiments with low/high plan complexity were 0.73/1.02 cm2 for direct and 0.73/1.02 cm2 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 cm2 mean, 2.88 cm2 95% percentile) were all significantly different from the low plan complexity (0.70 cm2 mean, 2.18 cm2 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.

Technical Note: Real-time image-guided adaptive radiotherapy of a rigid target for a prototype fixed beam radiotherapy system.

PURPOSE: Fixed beam radiotherapy systems utilize couch movement and rotation instead of gantry rotation in order to simplify linear accelerator design. We investigate the ability to deliver fixed beam treatments with the same level of clinical accuracy as conventional (rotating beam) treatments using real-time image guidance to maintain this accuracy in the presence of rigid target motion. METHODS: A prototype fixed beam radiotherapy system was built using a standard linac with the beam fixed in the vertical position and a computer controlled rotation stage that rotated a rigid phantom about the superior-inferior axis. Kilovoltage Intrafraction Monitoring (KIM) and real-time beam adaptation with MLC tracking was applied to a five-field IMRT treatment plan with motion introduced to the phantom. The same IMRT treatment was also delivered with real-time adaptation using the conventional rotating beam geometry. Film dosimetry was used to measure the dose delivered with a fixed beam compared to a rotating beam, as well as to compare treatments delivered with and without real-time adaptation. RESULTS: The dose distributions were found to be equivalent between the fixed beam and rotating beam geometry for real-time adaptive radiotherapy using KIM and MLC tracking beam adaptation. Gamma analysis on the films showed agreement >98% using a 2%/2 mm criteria with adaptation for static shifts and periodic motion. CONCLUSIONS: Fixed beam treatments with real-time beam adaptation are dosimetrically equivalent to conventional treatments with a rotating beam, even in the presence of rigid target motion. This suggests that, for a rigid target, the high clinical accuracy of real-time adaptive radiotherapy can be achieved with simpler beam geometry.

Real-time high spatial resolution dose verification in stereotactic motion adaptive arc radiotherapy.

PURPOSE: Radiation treatments delivered with real-time multileaf collimator (MLC) tracking currently lack fast pretreatment or real-time quality assurance. The purpose of this study is to test a 2D silicon detector, MagicPlate-512 (MP512), in a complex clinical environment involving real-time reconfiguration of the MLC leaves during target tracking. METHODS: MP512 was placed in the center of a solid water phantom and mounted on a motion platform used to simulate three different patient motions. Electromagnetic target tracking was implemented using the Calypso system (Varian Medical Systems, Palo Alto, CA, USA) and an MLC tracking software. A two-arc VMAT plan was delivered and 2D dose distributions were reconstructed by MP512, EBT3 film, and the Eclipse treatment planning system (TPS). Dose maps were compared using gamma analysis with 2%/2 mm and 3%/3 mm acceptance criteria. Dose profiles were generated in sup-inf and lateral directions to show the agreement of MP512 to EBT3 and to highlight the efficacy of the MLC tracking system in mitigating the effect of the simulated patient motion. RESULTS: Using a 3%/3 mm acceptance criterion for 2D gamma analysis, MP512 to EBT3 film agreement was 99% and MP512 to TPS agreement was 100%. For a 2%/2 mm criterion, the agreement was 95% and 98%, respectively. Full width at half maximum and 80%/20% penumbral width of the MP512 and EBT3 dose profiles agreed within 1 mm and 0.5 mm, respectively. Patient motion increased the measured dose profile penumbral width by nearly 2 mm (with respect to the no-motion case); however, the MLC tracking strategy was able to mitigate 80% of this effect. CONCLUSIONS: MP512 is capable of high spatial resolution 2D dose reconstruction during adaptive MLC tracking, including arc deliveries. It shows potential as an effective tool for 2D small field dosimetry and pretreatment quality assurance for MLC tracking modalities. These results provide confidence that detector-based pretreatment dosimetry is clinically feasible despite fast real-time MLC reconfigurations.