Tumor-site specific geometric distortions in high field integrated magnetic resonance linear accelerator radiotherapy

Geometric distortions in clinical MRI sequences for radiotherapy: insights gained from a multicenter investigation

BACKGROUND

As magnetic resonance imaging (MRI) becomes increasingly integrated into radiotherapy (RT) for enhanced treatment planning and adaptation, the inherent geometric distortion in acquired MR images pose a potential challenge to treatment accuracy. This study aimed to evaluate the geometric distortion levels in the clinical MRI protocols used across Danish RT centers and discuss influence of specific sequence parameters. Based on the variety in geometric performance across centers, we assess if harmonization of MRI sequences is a relevant measure.

MATERIALS AND METHODS

Nine centers participated with 12 MRI scanners and MRI-Linacs (MRL). Using a travelling phantom approach, a reference MRI sequence was used to assess variation in baseline distortion level between scanners. The phantom was also scanned with local clinical MRI sequences for brain, head/neck (H/N), abdomen, and pelvis. The influence of echo time, receiver bandwidth, image weighting, and 2D/3D acquisition was investigated.

RESULTS

We found a large variation in geometric accuracy across 93 clinical sequences examined, exceeding the baseline variation found between MRI scanners (σ = 0.22 mm), except for abdominal sequences where the variation was lower. Brain and abdominal sequences showed lowest distortion levels ([0.22, 2.26] mm), and a large variation in performance was found for H/N and pelvic sequences ([0.19, 4.07] mm). Post hoc analyses revealed that distortion levels decreased with increasing bandwidth and a less clear increase in distortion levels with increasing echo time. 3D MRI sequences had lower distortion levels than 2D (median of 1.10 and 2.10 mm, respectively), and in DWI sequences, the echo-planar imaging read-out resulted in highest distortion levels.

CONCLUSION

There is a large variation in the geometric distortion levels of clinical MRI sequences across Danish RT centers, and between anatomical sites. The large variation observed makes harmonization of MRI sequences across institutions and adoption of practices from well-performing anatomical sites, a relevant measure within RT.

Multi-parametric MRI for radiotherapy simulation

Magnetic resonance imaging (MRI) has become an important imaging modality in the field of radiotherapy (RT) in the past decade, especially with the development of various novel MRI and image-guidance techniques. In this review article, we will describe recent developments and discuss the applications of multi-parametric MRI (mpMRI) in RT simulation. In this review, mpMRI refers to a general and loose definition which includes various multi-contrast MRI techniques. Specifically, we will focus on the implementation, challenges, and future directions of mpMRI techniques for RT simulation.

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 (P₀ ), the edge (P ₁) and just outside (P ₂) 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 ₁ and 3.1 mm at P ₂. When distortion correction was applied, the error was reduced to 1.0 mm at P ₁ and 1.1 mm at P ₂. The corrected error was similar to an error of 0.9 mm at P₀ 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.

Toward MR-integrated proton therapy: modeling the potential benefits for liver tumors

Magnetic resonance imaging (MRI)-integrated proton therapy (MRiPT) is envisioned to improve treatment quality for many cancer patients. However, given the availability of alternative image-guided strategies, its clinical need is yet to be justified. This study aims to compare the expected clinical outcomes of MRiPT with standard of practice cone-beam CT (CBCT)-guided PT, and other MR-guided methods, i.e. offline MR-guided PT and MR-linac, for treatment of liver tumors. Clinical outcomes were assessed by quantifying the dosimetric and biological impact of target margin reduction enabled by each image-guided approach. Planning target volume (PTV) margins were calculated using random and systematic setup, delineation and motion uncertainties, which were quantified by analyzing longitudinal MRI data for 10 patients with liver tumors. Proton treatment plans were created using appropriate PTV margins for each image-guided PT method. Photon plans with margins equivalent to MRiPT were generated to represent MR-linac. Normal tissue complication probabilities (NTCP) of the uninvolved liver were compared. We found that PTV margin can be reduced by 20% and 40% for offline MR-guided PT and MRiPT, respectively, compared with CBCT-guided PT. Furthermore, clinical target volume expansion could be largely alleviated when delineating on MRI rather than CT. Dosimetric implications included decreased equivalent mean dose of the uninvolved liver, i.e. up to 24.4 Gy and 27.3 Gy for offline MR-guided PT and MRiPT compared to CBCT-guided PT, respectively. Considering Child-Pugh score increase as endpoint, NTCP of the uninvolved liver was significantly decreased for MRiPT compared to CBCT-guided PT (up to 48.4%,p < 0.01), offline MR-guided PT (up to 12.9%,p < 0.01) and MR-linac (up to 30.8%,p < 0.05). Target underdose was possible in the absence of MRI-guidance (D90 reduction up to 4.2 Gy in 20% of cases). In conclusion, MRiPT has the potential to significantly reduce healthy liver toxicities in patients with liver tumors. It is superior to other image-guided techniques currently available.

Synthetic CT for single-fraction neoadjuvant partial breast irradiation on an MRI-linac

A synthetic computed tomography (sCT) is required for daily plan optimization on an MRI-linac. Yet, only limited information is available on the accuracy of dose calculations on sCT for breast radiotherapy. This work aimed to (1) evaluate dosimetric accuracy of treatment plans for single-fraction neoadjuvant partial breast irradiation (PBI) on a 1.5 T MRI-linac calculated on a) bulk-density sCT mimicking the current MRI-linac workflow and b) deep learning-generated sCT, and (2) investigate the number of bulk-density levels required. For ten breast cancer patients we created three bulk-density sCTs of increasing complexity from the planning-CT, using bulk-density for: (1) body, lungs, and GTV (sCT_BD1); (2) volumes for sCT_BD1plus chest wall and ipsilateral breast (sCT_BD2); (3) volumes for sCT_BD2plus ribs (sCT_BD3); and a deep learning-generated sCT (sCT_DL) from a 1.5 T MRI in supine position. Single-fraction neoadjuvant PBI treatment plans for a 1.5 T MRI-linac were optimized on each sCT and recalculated on the planning-CT. Image evaluation was performed by assessing mean absolute error (MAE) and mean error (ME) in Hounsfield Units (HU) between the sCTs and the planning-CT. Dosimetric evaluation was performed by assessing dose differences, gamma pass rates, and dose-volume histogram (DVH) differences. The following results were obtained (median across patients for sCT_BD1/sCT_BD2/sCT_BD3/sCT_DLrespectively): MAE inside the body contour was 106/104/104/75 HU and ME was 8/9/6/28 HU, mean dose difference in the PTV_GTVwas 0.15/0.00/0.00/-0.07 Gy, median gamma pass rate (2%/2 mm, 10% dose threshold) was 98.9/98.9/98.7/99.4%, and differences in DVH parameters were well below 2% for all structures except for the skin in the sCT_DL. Accurate dose calculations for single-fraction neoadjuvant PBI on an MRI-linac could be performed on both bulk-density and deep learning sCT, facilitating further implementation of MRI-guided radiotherapy for breast cancer. Balancing simplicity and accuracy, sCT_BD2showed the optimal number of bulk-density levels for a bulk-density approach.

Technical Challenges of Real-Time Adaptive MR-Guided Radiotherapy

In the past few years, radiotherapy (RT) has experienced a major technological innovation with the development of hybrid machines combining magnetic resonance (MR) imaging and linear accelerators. This new technology for MR-guided cancer treatment has the potential to revolutionize the field of adaptive RT due to the opportunity to provide high-resolution, real-time MR imaging before and during treatment application. However, from a technical point of view, several challenges remain which need to be tackled to ensure safe and robust real-time adaptive MR-guided RT delivery. In this manuscript, several technical challenges to MR-guided RT are discussed. Starting with magnetic field strength tradeoffs, the potential and limitations for purely MR-based RT workflows are discussed. Furthermore, the current status of real-time 3D MR imaging and its potential for real-time RT are summarized. Finally, the potential of quantitative MR imaging for future biological RT adaptation is highlighted.

End-to-end validation of the geometric dose delivery performance of MR linac adaptive radiotherapy

The clinical introduction of hybrid magnetic resonance (MR) guided radiotherapy (RT) delivery systems has led to the need to validate the end-to-end dose delivery performance on such machines. In the current study, an MR visible phantom was developed and used to test the spatial deviation between planned and delivered dose at two 1.5 T MR linear accelerator (MR linac) systems, including pre-treatment imaging, dose planning, online imaging, image registration, plan adaptation, and dose delivery. The phantom consisted of 3D printed plastic and MR visible silicone rubber. It was designed to minimise air gaps close to the radiochromic film used as a dosimeter. Furthermore, the phantom was designed to allow submillimetre, reproducible positioning of the film in the phantom. At both MR linac systems, 54 complete adaptive, MR guided RT workflow sessions were performed. To test the dose delivery performance of the MR linac systems in various adaptive RT (ART) scenarios, the sessions comprised a range of systematic positional shifts of the phantom and imaging or plan adaptation conditions. In each workflow session, the positional translation between the film and the adaptive planned dose was determined. The results showed that the accuracy of the MR linac systems was between 0.1 and 0.9 mm depending on direction. The highest mean deviance observed was in the posterior-anterior direction, and the direction of the error was consistent between centres. The precision of the systems was related to whether the workflow utilized the internal image registration algorithm of the MR linac. Workflows using the internal registration algorithm led to a worse precision (0.2-0.7 mm) compared to workflows where the algorithm was decoupled (0.2 mm). In summary, the spatial deviation between planned and delivered dose of MR-guided ART at the two MR linac systems was well below 1 mm and thus acceptable for clinical use.