Respiratory motion introduces substantial uncertainties in abdominal radiotherapy for which traditionally large margins are used. The MR-Linac will open up the opportunity to acquire high resolution ...MR images just prior to radiation and during treatment. However, volumetric MRI time series are not able to characterize 3D tumor and organ-at-risk motion with sufficient temporal resolution. In this study we propose a method to estimate 3D deformation vector fields (DVFs) with high spatial and temporal resolution based on fast 2D imaging and a subject-specific motion model based on respiratory correlated MRI. In a pre-beam phase, a retrospectively sorted 4D-MRI is acquired, from which the motion is parameterized using a principal component analysis. This motion model is used in combination with fast 2D cine-MR images, which are acquired during radiation, to generate full field-of-view 3D DVFs with a temporal resolution of 476 ms. The geometrical accuracies of the input data (4D-MRI and 2D multi-slice acquisitions) and the fitting procedure were determined using an MR-compatible motion phantom and found to be 1.0-1.5 mm on average. The framework was tested on seven healthy volunteers for both the pancreas and the kidney. The calculated motion was independently validated using one of the 2D slices, with an average error of 1.45 mm. The calculated 3D DVFs can be used retrospectively for treatment simulations, plan evaluations, or to determine the accumulated dose for both the tumor and organs-at-risk on a subject-specific basis in MR-guided radiotherapy.
•The Elekta Unity MR-linac adaptive radiotherapy concept is explained.•The adapt to shape and adapt to position workflows are compared.•Different methods for dose re-calculation and optimization are ...discussed.•Full online re-planning is the most robust adaptive planning method for the Unity.•Faster methods are available, but should be dosimetrically explored per use case.
The promise of the MR-linac is that one can visualize all anatomical changes during the course of radiotherapy and hence adapt the treatment plan in order to always have the optimal treatment. Yet, there is a trade-off to be made between the time spent for adapting the treatment plan against the dosimetric gain. In this work, the various daily plan adaptation methods will be presented and applied on a variety of tumour sites. The aim is to provide an insight in the behavior of the state-of-the-art 1.5 T MRI guided on-line adaptive radiotherapy methods.
To explore the different available plan adaptation workflows and methods, we have simulated online plan adaptation for five cases with varying levels of inter-fraction motion, regions of interest and target sizes: prostate, rectum, esophagus and lymph node oligometastases (single and multiple target). The plans were evaluated based on the clinical dose constraints and the optimization time was measured.
The time needed for plan adaptation ranged between 17 and 485 s. More advanced plan adaptation methods generally resulted in more plans that met the clinical dose criteria. Violations were often caused by insufficient PTV coverage or, for the multiple lymph node case, a too high dose to OAR in the vicinity of the PTV. With full online replanning it was possible to create plans that met all clinical dose constraints for all cases.
Daily full online replanning is the most robust adaptive planning method for Unity. It is feasible for specific sites in clinically acceptable times. Faster methods are available, but before applying these, the specific use cases should be explored dosimetrically.
•A comprehensive commissioning protocol is designed to assess the MRI performance of integrated MR-linac systems.•The imaging performance of four recently installed 1.5T MR-linac systems is assessed ...in a multi-institutional study.•A benchmark data set is provided by making the results as well as the corresponding measurement protocols publically available.
Magnetic Resonance linear accelerator (MR-linac) systems represent a new type of technology that allows for online MR-guidance for high precision radiotherapy (RT). Currently, the first MR-linac installations are being introduced clinically. Since the imaging performance of these integrated MR-linac systems is critical for their application, a thorough commissioning of the MRI performance is essential. However, guidelines on the commissioning of MR-guided RT systems are not yet defined and data on the performance of MR-linacs are not yet available.
Here we describe a comprehensive commissioning protocol, which contains standard MRI performance measurements as well as dedicated hybrid tests that specifically assess the interactions between the Linac and the MRI system. The commissioning results of four MR-linac systems are presented in a multi-center study.
Although the four systems showed similar performance in all the standard MRI performance tests, some differences were observed relating to the hybrid character of the systems. Field homogeneity measurements identified differences in the gantry shim configuration, which was later confirmed by the vendor.
Our results highlight the importance of dedicated hybrid commissioning tests and the ability to compare the machines between institutes at this very early stage of clinical introduction. Until formal guidelines and tolerances are defined the tests described in this study may be used as a practical guideline. Moreover, the multi-center results provide initial bench mark data for future MR-linac installations.
•5 patients with pelvic lymph node metastases received SBRT using a 1.5 T MR-linac.•Session time was <60 min for all 25 treatment fractions.•All quality assurance tests were passed (dose calculations ...& film measurements).
Online adaptive radiotherapy using the 1.5 Tesla MR-linac is feasible for SBRT (5 × 7 Gy) of pelvic lymph node oligometastases. The workflow allows full online planning based on daily anatomy. Session duration is less than 60 min. Quality assurance tests, including independent 3D dose calculations and film measurements were passed.
To enable magnetic resonance imaging (MRI)-guided radiotherapy with real-time adaptation, motion must be quickly estimated with low latency. The motion estimate is used to adapt the radiation beam to ...the current anatomy, yielding a more conformal dose distribution. As the MR acquisition is the largest component of latency, deep learning (DL) may reduce the total latency by enabling much higher undersampling factors compared to conventional reconstruction and motion estimation methods. The benefit of DL on image reconstruction and motion estimation was investigated for obtaining accurate deformation vector fields (DVFs) with high temporal resolution and minimal latency. 2D cine MRI acquired at 1.5 T from 135 abdominal cancer patients were retrospectively included in this study. Undersampled radial golden angle acquisitions were retrospectively simulated. DVFs were computed using different combinations of conventional- and DL-based methods for image reconstruction and motion estimation, allowing a comparison of four approaches to achieve real-time motion estimation. The four approaches were evaluated based on the end-point-error and root-mean-square error compared to a ground-truth optical flow estimate on fully-sampled images, the structural similarity (SSIM) after registration and time necessary to acquire k-space, reconstruct an image and estimate motion. The lowest DVF error and highest SSIM were obtained using conventional methods up to R≤10. For undersampling factors R>10, the lowest DVF error and highest SSIM were obtained using conventional image reconstruction and DL-based motion estimation. We have found that, with this combination, accurate DVFs can be obtained up to R=25 with an average root-mean-square error up to 1 millimeter and an SSIM greater than 0.8 after registration, taking 60 milliseconds. High-quality 2D DVFs from highly undersampled k-space can be obtained with a high temporal resolution with conventional image reconstruction and a deep learning-based motion estimation approach for real-time adaptive MRI-guided radiotherapy.
This work aimed to quantify any principal magnetic field (B
) inhomogeneity and changes in MR image geometric distortion with continuous linac gantry rotation on an Elekta Unity MR-linac. This ...situation occurs for around a second between treatment beams during current image guided radiotherapy treatment and would occur frequently in foreseeable real-time adaptive radiotherapy treatment. Pixel by pixel maps of B
inhomogeneity were obtained via repeated high temporal resolution pulse sequences with the linac gantry static at 36 gantry angles spaced ten degrees apart, and in continuous rotation at both 1 and 2 rpm. Individual B
maps were subtracted from average maps across all data and the residual peak to peak inhomogeneity was calculated for each. The bulk geometric shift and change in physical extent of a 10 cm diameter spherical flood phantom during continuous linac gantry rotation at 1 and 2 rpm was compared to the static gantry case for two pulse sequences: the real-time clinical monitoring bFFE sequence and a non-clinical EPI sequence, chosen for its susceptibility to geometric distortion. The peak to peak inhomogeneity in the deviation-from-average ppm maps, plotted against gantry angle with the gantry in continuous rotation at 1 and 2 rpm were negligibly different from equivalent data obtained with the gantry static. The real-time clinical monitoring pulse sequence was shown to give negligible geometric distortion during continuous gantry motion, whilst a non-clinical EPI sequence showed bulk shifts of the order of one pixel and gantry angle dependent changes in extent, demonstrating the sensitivity of the chosen method. MR imaging on the Elekta Unity MR-Linac with the gantry in continuous motion is negligibly different from the static gantry case with current clinical pulse sequences. Real-time tracking and treatment plan adaptation using MR images obtained with the linac gantry in motion is possible.
Abstract Background and purpose To characterize pancreatic tumor motion and to develop a gating scheme for radiotherapy in pancreatic cancer. Materials and methods Two cine MRIs of 60 s each were ...performed in fifteen pancreatic cancer patients, one in sagittal direction and one in coronal direction. A Minimum Output Sum of Squared Error (MOSSE) adaptive correlation filter was used to quantify tumor motion in craniocaudal, lateral and anteroposterior directions. To develop a gating scheme, stability of the breathing phases was examined and a gating window assessment was created, incorporating tumor motion, treatment time and motion margins. Results The largest tumor motion was found in craniocaudal direction, with an average peak-to-peak amplitude of 15 mm (range 6–34 mm). Amplitude of the tumor in the anteroposterior direction was on average 5 mm (range 1–13 mm). The least motion was seen in lateral direction (average 3 mm, range 2–5 mm). The end exhale position was the most stable position in the breathing cycle and tumors spent more time closer to the end exhale position than to the end inhale position. On average, a margin of 25% of the maximum craniocaudal breathing amplitude was needed to achieve full target coverage with a duty cycle of 50%. When reducing the duty cycle to 50%, a margin of 5 mm was sufficient to cover the target in 11 out of 15 patients. Conclusion Gated delivery for radiotherapy of pancreatic cancer is best performed around the end exhale position as this is the most stable position in the breathing cycle. Considerable margin reduction can be established at moderate duty cycles, yielding acceptable treatment efficiency. However, motion patterns and amplitude do substantially differ between individual patients. Therefore, individual treatment strategies should be considered for radiotherapy in pancreatic cancer.
To determine the optimum sampling strategy for retrospective reconstruction of 4-dimensional (4D) MR data for nonrigid motion characterization of tumor and organs at risk for radiation therapy ...purposes.
For optimization, we compared 2 surrogate signals (external respiratory bellows and internal MRI navigators) and 2 MR sampling strategies (Cartesian and radial) in terms of image quality and robustness. Using the optimized protocol, 6 pancreatic cancer patients were scanned to calculate the 4D motion. Region of interest analysis was performed to characterize the respiratory-induced motion of the tumor and organs at risk simultaneously.
The MRI navigator was found to be a more reliable surrogate for pancreatic motion than the respiratory bellows signal. Radial sampling is most benign for undersampling artifacts and intraview motion. Motion characterization revealed interorgan and interpatient variation, as well as heterogeneity within the tumor.
A robust 4D-MRI method, based on clinically available protocols, is presented and successfully applied to characterize the abdominal motion in a small number of pancreatic cancer patients.
•Machine characteristics must be considered when selecting b-values for measuring the ADC.•Spatial variation of the ADC requires positioning of a ROI near the iso-center.•The ADC can be measured on ...the Unity system during gantry rotation and irradiation.•Includes example acquisition parameters to facilitate multi-center research.
Diffusion-weighted imaging (DWI) for treatment response monitoring is feasible on hybrid magnetic resonance linear accelerator (MR-linac) systems. The MRI scanner of the Elekta Unity system has an adjusted design compared to diagnostic scanners. We investigated its impact on measuring the DWI-derived apparent diffusion coefficient (ADC) regarding three aspects: the choice of b-values, the spatial variation of the ADC, and scanning during radiation treatment. The aim of this study is to give recommendations for accurate ADC measurements on Unity systems.
Signal-to-noise ratio (SNR) measurements with increasing b-values were done to determine the highest bvalue that can be measured reliably. The spatial variation of the ADC was assessed on six Unity systems with a cylindrical phantom of 40 cm diameter. The influence of gantry rotation and irradiation was investigated by acquiring DWI images before and during treatment of 11 prostate cancer patients.
On the Unity system, a maximum b-value of 500 s/mm2 should be used for ADC quantification, as a trade-off between SNR and diffusion weighting. Accurate ADC values were obtained within 7 cm from the iso-center, while outside this region ADC values deviated more than 5%. The ADC was not influenced by the rotating linac or irradiation during treatment.
We provide Unity system specific recommendations for measuring the ADC. This will increase the consistency of ADC values acquired in different centers on the Unity system, enabling large cohort studies for biomarker discovery and treatment response monitoring.