First clinical experience of correcting phantom-based image distortion related to gantry position on a 0.35T MR-Linac

Characterization of real-time cine MR imaging distortion on 0.35 T MRgRT with concentric cine imaging QA phantom

Objective. Real-time MRgRT uses 2D-cine imaging for target tracking and motion evaluation. Rotation of gantry inducedB0off-resonance, resulting in image artifacts and imaging isocenter-shift precluding MR-guided arc therapy. Standard MRI phantoms designed for higher resolution images face challenges when low-resolution cine imaging is needed to achieve high frame rates. This work aimed to examine the spatial accuracy including geometric distortion and isocenter shift in real-time during gantry rotation on a 0.35 T MR-Linac using the concentric Cine imaging quality assurance (QA) phantom and its associated image analysis software.Approach. The Cine imaging QA phantom consists of two concentric shells of low-T1mineral oil and a central alignment structure. The phantom was scanned on three different MRI systems; 0.55 T Siemens Free.Max, 1.5 T Philips Ingenia, and 0.35 T ViewRay MRIdian MR-Linac using 2D balanced steady-state free precession (bSSFP) imaging sequence. In addition, bSSFP cine MRI with the banding artifact correction was tested on 0.35 T ViewRay MR-Linac. Images from the MR-Linac were acquired with the Linac gantry stationary and rotating from gantry 300°→ 0° and vice versa. Three orthogonal image planes were scanned excluding the 1.5 T Philips Ingenia, where only the axial plane was scanned. The image analysis software calculated the distortion values as well as the isocenter position for each cine frame.Main results. The geometric distortion of cine imaging on MRIs and MR-Linac at gantry stationary are within 1 mm while the substantial geometric distortion of 2 and 2.2 mm were observed on 0.35 T MR-Linac while rotating the gantry clockwise (300°→ 0°) and counterclockwise 0°→ 300° respectively. The average imaging isocenter shift was 0.1 mm for both MRIs and the static gantry and imaging isocenter shift of ≤1.5 mm was observed during the gantry rotation. The imaging isocenter shift decreased by 1 ± 0.2 mm clockwise and counterclockwise withB0compensation.Significance. The concentric Cine imaging QA phantom and its associated software effectively demonstrate the image distortion on real-time cine imaging on regular MRIs and 0.35 T MR-Linac. The results of significant geometric distortion with a rotating gantry in the MR-Linac system require further investigation to alleviate the extent of the image distortion.

MRgRT Quality Assurance for a Low-Field MR-Linac

The introduction of MR-guided treatment machines into the radiation oncology clinic has provided unique challenges for the radiotherapy QA program. These MR-linac systems require that existing QA procedures be adapted to verify linac performance within the magnetic field environment and that new procedures be added to ensure acceptable image quality for the MR system. While both high and low-field MR-linac options exist, this chapter is intended to provide a structure for implementing a QA program within the low-field MR environment. This review is divided into three sections. The first section focuses on machine QA tasks including mechanical and dosimetric verification. The second section is concentrated on the procedures implemented for imaging QA. Finally, the last section covers patient specific QA tasks including special considerations related to the performance of patient specific QA within the framework of online adaptive radiotherapy.

Characterisation of a split gradient coil design induced systemic imaging artefact on 0.35 T MR-linac systems

Objective. The aim of this work was to highlight and characterize a systemic 'star-like' artefact inherent to the low field 0.35 T MRIdian MR-linac system, a magnetic resonance guided radiotherapy device. This artefact is induced by the original split gradients coils design. This design causes a surjection of the intensity gradient inZ(or head-feet) direction. This artefact appears on every sequence with phase encoding in the head-feet direction.Approach. Basic gradient echo sequence and clinical mandatory bSSFP sequence were used. Three setups using manufacturer provided QA phantoms were designed: two including the linearity control grid used for the characterisation and a third including two homogeneity control spheres dedicated to the artefact management in a more clinical like situation. The presence of the artefact was checked in four different MRidian sites. The tested parameters based on the literature were: phase encoding orientation, slab selectivity, excitation bandwidth (BW_RF), acceleration factor (R) and phase/slab oversampling (PO/SO).Main results. The position of this artefact is constant and reproducible over the tested MRIdian sites. The typical singularity saturated dot or star is visible even with the 3D slab-selection enabled. A management is proposed by decreasing the BW_RF, theRin head-feet direction and increasing the PO/SO. The oversampling can be optimized using a formula to anticipate the location of artefact in the field of view.Significance. The star-like artefact has been well characterised. A manageable solution comes at the cost of acquisition time. Observed in clinical cases, the artefact may degrade the images used for the RT planning and repositioning during the treatment unless corrected.

System-dependent image distortion related to gantry positions of a 0.35 T MRgRT: Characterization and the corresponding correction

PURPOSE

MR-guided radiotherapy with high accuracy treatment planning requires addressing MR imaging artifacts that originate from system imperfections. This work presents the characterization and corresponding correction of gantry-related imaging distortions including geometric distortion and isocenter shift in a 0.35 T magnetic resonance imaging (MRI)-guided radiotherapy (MRgRT) system using distortion vector fields (DVFs).

METHODS

Two phantoms, the magnetic resonance imaging distortion in 3D (MRID^3D ) phantom and the Fluke phantom, along with a human volunteer were imaged at different gantry angles on a 0.35 T MR-Linac. The geometric distortion and isocenter shift were characterized for both phantom images. DVFs with a field of view extended beyond the physical boundary of the MRID^3D phantom were extracted from images taken at 30° gantry angle increments, with vendor-provided distortion correction turned on and off (DstOff). These extended DVFs were then applied to the relevant phantom images to correct their geometric distortions and isocenter shift at the respective gantry angles. The extended DVFs produced from the MRID^3D phantom were also applied to Fluke phantom and human MR images at their respective gantry angles. The resampled images were evaluated using structural similarity index measure (SSIM) comparison with the vendor corrected images from the MRgRT system.

RESULTS

Geometric distortion with "mean (± SD) distortion" of 3.2 ± 0.02, 2.9 ± 0.02, and 1.8 ± 0.01 mm and isocenter shift (±SD) of 0.49 ± 0.3, 0.05 ± 0.2, and 0.01 ± 0.03 mm were present in the DstOff MRID^3D phantom images in right-left (RL), anterior-posterior (AP), and superior-inferior (SI) directions, respectively. After resampling the originally acquired images by applying extended DVFs, the distortion was corrected to 0.18 ± 0.02, 0.09 ± 0.01, 0.15 ± 0.01 mm, and isocenter shift was corrected to 0.14 ± 0.05, -0.02 ± 0.04, and -0.07 ± 0.05 mm in RL, AP, and SI directions, respectively. The Fluke phantom average geometric distortion with "mean (± SD) distortion" of 2.7 ± 0.1 mm was corrected to 0.2 ± 0. 1 mm and the average isocenter shift (± SD) of 0.51 ± 0.2 mm, and 0.05 ± 0.03 was corrected to -0.08 ± 0.03 mm, and -0.05 ± 0.01 in RL and AP directions, respectively. SSIM (mean ± SD) of the original images to resampled images was increased from 0.49 ± 0.02 to 0.78 ± 0.01, 0.45 ± 0.02 to 0.75 ± 0.01, and 0.86 ± 0.25 to 0.98 ± 0.08 for MRID^3D phantom, Fluke phantom, and human MR images, respectively, for all the gantry angles compared to the vendor corrected images.

CONCLUSION

The gantry-related MR imaging distortion including geometric distortion and isocenter shift was characterized and a corresponding correction was demonstrated using extended DVFs on 0.35 T MRgRT system. The characterized gantry-related isocenter shift can be combined with geometric distortion correction to provide a technique for the correction of the full system-dependent distortion in an MRgRT system.

First clinical experience of correcting phantom-based image distortion related to gantry position on a 0.35T MR-Linac

MR‐guided radiotherapy requires strong imaging spatial integrity to deliver high quality plans and provide accurate dose calculation. The MRI system, however, can be compromised by the integrated linear accelerator (Linac), resulting in inaccurate imaging isocenter position and geometric distortion. Dependence on gantry position further complicates the correction of distortions. This work presents a new clinical application of a commercial phantom and software system that quantifies isocenter alignment and geometric distortion, as well as providing a deformation vector field (DVF). A large distortion phantom and a smaller grid phantom were imaged at multiple gantry angles from 0 to 330° on a 0.35 T integrated MR‐Linac. The software package was used to assess geometric distortion and generate DVFs to correct distortions within the phantom volume. The DVFs were applied to the grid phantom with resampling software then evaluated using structural similarity index measure (SSIM). Scans were also performed with a ferromagnetic clip near the phantom to investigate the correction of more severe artifacts. The mean magnitude isocenter shift was 0.67 mm, ranging from 0.25 to 1.04 mm across all angles. The DVF had a mean component value of 0.27 ± 0.02, 0.24 ± 0.01, and 0.19 ± 0.01 mm in the right‐left (RL), anterior‐posterior (AP), and superior‐inferior (SI) directions. The ferromagnetic clip increased isocenter position error from 1.98 mm to 2.20 mm and increased mean DVF component values in the RL and AP directions. The resampled grid phantom had an increased SSIM for all gantry angles compared to original images, increasing from 0.26 ± 0.001 to 0.70 ± 0.004. Through this clinical assessment, we were able to correct geometric distortion and isocenter shift related to gantry position on a 0.35 T MR‐Linac using the distortion phantom and software package. This provides encouragement that it could be used for quality assurance and clinically to correct systematic distortion caused by imaging at different gantry angles.