MRI quality control for low-field MR-IGRT systems: Lessons learned

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.

Receive coil quality assurance procedure and automated analysis for ViewRay MRIdian MR-Linac

PURPOSE

Regular receiving coil quality assurance (QA) is required to ensure image quality of an MRIdian Linac system. The manufacturer provides a spherical phantom and positioning tube for single-slice signal-to-noise ratio (SNR) and uniformity assessments. We aimed to improve imaging setup and coverage and eliminate inter-scan variability by employing multi-slice imaging of a stable phantom. Additionally, we strived to expedite analysis by developing objective, automated analysis software.

METHODS

A 5300 mL cylindrical plastic bottle placed in plastic bins was scanned at isocenter using a spin-echo sequence with NEMA-recommended parameters and 18 axial slices, avoiding phantom repositioning. Acquisition was repeated with and without prescan normalization filtering and by saving uncombined element images. Obtained data were analyzed using custom open-source MATLAB code. Signal and noise images were automatically assigned, and ROIs for SNR and uniformity calculations were defined using image thresholding. SNR and uniformity pass/fail decisions were made using baseline comparisons.

RESULTS

The proposed method was successfully implemented as monthly coil QA for 3.5 years. Setup and scanning took 41 min on average for a coil set. Automated image analysis was completed in a few minutes. Signal intensity peaked around +90 or -90 mm for Torso or Head/Neck coil unfiltered images. Noise peaked and minimized SNR inside ±30 mm from isocenter, while maximizing it around ±130 mm. Prescan normalization smoothed signal response, reduced SNR and increased uniformity. Individual coil element image analysis identified their position, signal or noise response and SNR. SNR and uniformity pass/fail thresholds were set for already tested and new coils. Conspicuous and subtle Torso coil malfunctions were detected considering baseline deviations of combined and individual element results.

CONCLUSIONS

Our QA method eliminated observer bias and provided insights into coil function, image filtering performance and coil element location. It provided SNR and uniformity thresholds and identified faulty coil elements.

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.

Comprehensive MR imaging QA of 0.35 T MR-Linac using a multi-purpose large FOV phantom: A single-institution experience

PURPOSE

Magnetic resonance-guided radiotherapy (MRgRT) is desired for the treatment of diseases in the abdominothoracic region, which has a broad imaging area and continuous motion. To ensure accurate treatment delivery, an effective image quality assurance (QA) program, with a phantom that covers the field of view (FOV) similar to a human torso, is required. However, routine image QA for a large FOV is not readily available at many MRgRT centers. In this work, we present the clinical experience of the large FOV MRgRT Insight phantom for periodic daily and monthly comprehensive magnetic resonance imaging (MRI)-QA and its feasibility compared to the existing institutional routine MRI-QA procedures in 0.35 T MRgRT.

METHODS

Three phantoms; ViewRay cylindrical water phantom, Fluke 76-907 uniformity and linearity phantom, and Modus QA large FOV MRgRT Insight phantom, were imaged on the 0.35 T MR-Linac. The measurements were made in MRI mode with the true fast imaging with steady-state free precession (TRUFI) sequence. The ViewRay cylindrical water phantom was imaged in a single-position setup whereas the Fluke phantom and Insight phantom were imaged in three different orientations: axial, sagittal, and coronal. Additionally, the phased array coil QA was performed using the horizontal base plate of the Insight phantom by placing the desired coil around the base section which was compared to an in-house built Polyurethane foam phantom for reference.

RESULT

The Insight phantom captured image artifacts across the entire planar field of view, up to 400 mm, in a single image acquisition, which is beyond the FOV of the conventional phantoms. The geometric distortion test showed a similar distortion of 0.45 ± 0.01  and 0.41 ± 0.01 mm near the isocenter, that is, within 300 mm lengths for Fluke and Insight phantoms, respectively, but showed higher geometric distortion of 0.8 ± 0.4 mm in the peripheral region between 300 and 400 mm of the imaging slice for the Insight phantom. The Insight phantom with multiple image quality features and its accompanying software utilized the modulation transform function (MTF) to evaluate the image spatial resolution. The average MTF values were 0.35 ± 0.01, 0.35 ± 0.01, and 0.34 ± 0.03 for axial, coronal, and sagittal images, respectively. The plane alignment and spatial accuracy of the ViewRay water phantom were measured manually. The phased array coil test for both the Insight phantom and the Polyurethane foam phantoms ensured the proper functionality of each coil element.

CONCLUSION

The multifunctional large FOV Insight phantom helps in tracking MR imaging quality of the system to a larger extent compared to the routine daily and monthly QA phantoms currently used in our institute. Also, the Insight phantom is found to be more feasible for routine QA with easy setup.

A deep learning approach to generate synthetic CT in low field MR-guided radiotherapy for lung cases

INTRODUCTION

This study aims to apply a conditional Generative Adversarial Network (cGAN) to generate synthetic Computed Tomography (sCT) from 0.35 Tesla Magnetic Resonance (MR) images of the thorax.

METHODS

Sixty patients treated for lung lesions were enrolled and divided into training (32), validation (8), internal (10,T_A) and external (10,T_B) test set. Image accuracy of generated sCT was evaluated computing the mean absolute (MAE) and mean error (ME) with respect the original CT. Three treatment plans were calculated for each patient considering MRI as reference image: original CT, sCT (pure sCT) and sCT with GTV density override (hybrid sCT) were used as Electron Density (ED) map. Dose accuracy was evaluated comparing treatment plans in terms of gamma analysis and Dose Volume Histogram (DVH) parameters.

RESULTS

No significant difference was observed between the test sets for image and dose accuracy parameters. Considering the whole test cohort, a MAE of 54.9 ± 10.5 HU and a ME of 4.4 ± 7.4 HU was obtained. Mean gamma passing rates for 2%/2mm, and 3%/3mm tolerance criteria were 95.5 ± 5.9% and 98.2 ± 4.1% for pure sCT, 96.1 ± 5.1% and 98.5 ± 3.9% for hybrid sCT: the difference between the two approaches was significant (p = 0.01). As regards DVH analysis, differences in target parameters estimation were found to be within 5% using hybrid approach and 20% using pure sCT.

CONCLUSION

The DL algorithm here presented can generate sCT images in the thorax with good image and dose accuracy, especially when the hybrid approach is used. The algorithm does not suffer from inter-scanner variability, making feasible the implementation of MR-only workflows for palliative treatments.

Practical guidelines of online MR-guided adaptive radiotherapy

The first magnetic resonance (MR)-guided radiotherapy system in Japan was installed in May 2017. Implementation of online MR-guided adaptive radiotherapy (MRgART) began in February 2018. Online MRgART offers greater treatment accuracy owing to the high soft-tissue contrast in MR-images (MRI), compared to that in X-ray imaging. The Japanese Society for Magnetic Resonance in Medicine (JSMRM), Japan Society of Medical Physics (JSMP), Japan Radiological Society (JRS), Japanese Society of Radiological Technology (JSRT), and Japanese Society for Radiation Oncology (JASTRO) jointly established the comprehensive practical guidelines for online MRgART. These guidelines propose the essential requirements for clinical implementation of online MRgART with respect to equipment, personnel, institutional environment, practice guidance, and quality assurance/quality control (QA/QC). The minimum requirements for related equipment and QA/QC tools, recommendations for safe operation of MRI system, and the implementation system are described. The accuracy of monitor chamber and detector in dose measurements should be confirmed because of the presence of magnetic field. The ionization chamber should be MR-compatible. Non-MR-compatible devices should be used in an area that is not affected by the static magnetic field (outside the five Gauss line), and their operation should be checked to ensure that they do not affect the MR image quality. Dose verification should be performed using an independent dose verification system that has been confirmed to be reliable through commissioning. This guideline proposes the checklists to ensure the safety of online MRgART. Successful clinical implementation of online MRgART requires close collaboration between physician, radiological technologist, nurse, and medical physicist.

Real-time B0 compensation during gantry rotation in a 0.35 T MRI-Linac

BACKGROUND

Rotation of the ferromagnetic gantry of a low magnetic field MRI-Linac was previously demonstrated to cause large center frequency offsets of ±400 Hz. The B₀ off-resonances cause image artifacts and imaging isocenter shifts that would preclude MRI-guided arc therapy.

PURPOSE

The purpose of this study was to measure and compensate for center frequency offsets in real-time during gantry rotation on a 0.35 T MRI-Linac using a free induction decay (FID) navigator.

METHODS

A nonselective FID navigator was added before each 2D balanced steady-state free precession (bSSFP) cine image acquisition on a 0.35 T MRI-Linac. Images were acquired at 7.3 frames per second. Phase data from the initial FID navigator (while the gantry was stationary) was used as a reference. The phase data from each subsequent FID navigator was used to calculate the real-time B₀ off-resonance. The transmitter/receiver phase and the phase accrual over the adjacent image acquisition were adjusted to correct for the center frequency offset. Measurements were performed using a MRI-Linac dynamic phantom prior to and while the gantry rotated clockwise and counterclockwise. Image quality and signal-to-noise ratio were compared between uncorrected and B₀ corrected MRIs using a reference image acquired while the gantry was stationary. Four targets in the phantom were manually contoured on the first image frame and an active contouring algorithm was used retrospectively on each subsequent frame to assess image variations and calculate Dice coefficients. Additionally, three healthy volunteers were imaged using the same pulse sequences with and without real-time B₀ compensation during gantry rotation. Normalized root mean square errors (nRMSEs) were calculated for the phantom and in vivo to assess the efficacy of the B₀ compensation on image quality. The measured center frequency offsets from the volunteer and MRI dynamic phantom navigator data were also compared. The sinusoidal behavior of the center frequency offsets was modeled based on the gantry layout and long time constant eddy currents resulting from gantry rotation.

RESULTS

The duration of the FID navigator and processing was 4.5 ms. The FID navigator resulted in a ≤11% drop in signal-to-noise ratio (SNR) in the phantom and in vivo (liver). Dice coefficients from the MR-IGRT phantom contour measurements remained above 0.8 with B₀ compensation. Without B₀ compensation, the Dice coefficients dropped below 0.8 for up to 21% of the time depending on the contour. Real-time B₀ compensation resulted in mean reductions in nRMSE of 51% and 16% for the MR-IGRT phantom and in vivo, respectively. Peak-to-peak center frequency offsets ranged from 757 to 773 Hz in the phantom and 670 to 871 Hz in vivo.

CONCLUSION

Dynamic real-time B₀ compensation significantly improved image quality and reduced artifacts during gantry rotation in the phantom and in vivo. However, the FID navigator resulted in a small drop in the imaging duty cycle and SNR. This article is protected by copyright. All rights reserved.

Technical Note: Effects of rotating gantry on magnetic field and eddy currents in 0.35 T MRI-guided radiotherapy (MR-IGRT) system

PURPOSE

The purpose of this study was to identify the cause of severe image artifacts that occurred during gantry rotation in a 0.35 T MRI-Linac by comparing measurements of eddy currents, center frequency, and field inhomogeneities made with the gantry in motion and stationary.

METHODS

Gradient and B₀ eddy currents were calculated from the free induction decays (FIDs) resulting from selective excitation at a temporal resolution of 200 ms/measurement. B₀ eddy currents were also calculated from FIDs acquired with nonselective excitation at a temporal resolution of 100 ms/measurement. Center frequencies and B₀ inhomogeneities were measured by acquiring FIDs with a repetition time (TR) of 290 ms. Cartesian and radial 2D true fast imaging with steady-state precession (TrueFISP) pulse sequences used in real-time MRI-guided radiation therapy (MR-IGRT) were acquired. To assess artifact severity, the normalized root mean square error (nRMSE) was calculated between a reference MRI (static gantry) and MRIs acquired during gantry rotation for each serial acquisition. Image artifacts were qualitatively graded as nominal, minor, or severe. Measurements were conducted while the gantry was rotated through its entire range for both clockwise and counterclockwise. Measurements during gantry rotation were compared to measurements with a stationary gantry (every 30°).

RESULTS

Severe image artifacts were observed 22-35% of the time while the gantry was rotating. Short time constant eddy currents were not affected by gantry rotation. The peak to peak center frequency and FWHM rose by factors of 13.2-14.5 and 1.1-1.6, respectively, for the rotating versus stationary gantry. The magnitude of the center frequency offset and field inhomogeneities depended on the direction of the gantry rotation.

CONCLUSIONS

Image artifacts during gantry rotation were primarily caused by center frequency variations and field inhomogeneities. Therefore, dynamic B₀ compensation techniques should be able to reduce artifacts during gantry rotation.

Clinical assessment of geometric distortion for a 0.35T MR-guided radiotherapy system

Purpose

To estimate the overall spatial distortion on clinical patient images for a 0.35 T MR‐guided radiotherapy system.

Methods

Ten patients with head‐and‐neck cancer underwent CT and MR simulations with identical immobilization. The MR images underwent the standard systematic distortion correction post‐processing. The images were rigidly registered and landmark‐based analysis was performed by an anatomical expert. Distortion was quantified using Euclidean distance between each landmark pair and tagged by tissue interface: bone‐tissue, soft tissue, or air‐tissue. For baseline comparisons, an anthropomorphic phantom was imaged and analyzed.

Results

The average spatial discrepancy between CT and MR landmarks was 1.15 ± 1.14 mm for the phantom and 1.46 ± 1.78 mm for patients. The error histogram peaked at 0–1 mm. 66% of the discrepancies were <2 mm and 51% <1 mm. In the patient data, statistically significant differences (p‐values < 0.0001) were found between the different tissue interfaces with averages of 0.88 ± 1.24 mm, 2.01 ± 2.20 mm, and 1.41 ± 1.56 mm for the air/tissue, bone/tissue, and soft tissue, respectively. The distortion generally correlated with the in‐plane radial distance from the image center along the longitudinal axis of the MR.

Conclusion

Spatial distortion remains in the MR images after systematic distortion corrections. Although the average errors were relatively small, large distortions observed at bone/tissue interfaces emphasize the need for quantitative methods for assessing and correcting patient‐specific spatial distortions.

VMAT-like plans for magnetic resonance guided radiotherapy: Addressing unmet needs

PURPOSE

VMAT delivery technique is currently not applicable to Magnetic Resonance-guided radiotherapy (MRgRT) hybrid systems. Aim of this study is to evaluate an innovative VMAT-like (VML) delivery technique.

MATERIAL AND METHODS

First, planning and dosimetric evaluation of the MRgRT VML treatment have been performed on 10 different disease sites and the results have been compared with the corresponding IMRT plans. Then, in the second phase, 10 of the most dosimetrically challenging locally advanced pancreas treatment plans have been retrospectively re-planned using the VML approach to explore the potentiality of this new delivery technique. Finally, VML robustness was evaluated and compared with the IMRT plans, considering a lateral positioning error of ± 5 mm.

RESULTS

In phase one, all VML plans were within constraint for all OARs. When PTV coverage is considered, in the 50% of the cases VML PTV coverage is equal or higher than in IMRT plan. In the remaining 50%, the highest target under coverage difference in comparison with IMRT plan is -1.71%. The mean and maximum treatment time differences (VML-IMRT) is 0.2 min and 3.1 min respectively. In phase two, the treatment time variation (VML-IMRT), shows a mean, maximum and minimum variations of 1.3, 4.6 and -0.6 min respectively. All VML plans have a better target coverage if compared with IMRT plans, keeping in any case the OARs constraints within tolerance. VML doesn't increase plan robustness.

CONCLUSION

VMAT-like treatment approach appeared to be an efficient planning solution and it was decided to clinically implement it in daily practice, especially in the frame of hypo fractionated treatments.

Effects of B0 eddy currents on imaging isocenter shifts in 0.35-T MRI-guided radiotherapy (MR-IGRT) system

PURPOSE

The purpose of this study was to measure gantry angle-related eddy currents in a 0.35-T MRI-Linac and determine if B₀ (zeroth order) eddy currents are the primary cause of gantry angle-dependent imaging isocenter shifts vs other potential causes like B₀ inhomogeneities and gradient (first order) eddy currents. For conventional Cartesian acquisitions, B₀ eddy currents can cause imaging isocenter shifts along both phase encode and readout directions. Gradient eddy currents can cause spatial distortion along both the phase encode and readout directions. Center frequency offsets can cause imaging isocenter shifts along the readout direction that vary with readout gradient polarity.

METHODS

MRI-related eddy currents and imaging isocenter shifts were measured on a 0.35-T MRI-Linac at gantry angles from 0° to 330° in increments of 30^° . All measurements were made after gradient shimming and center frequency tuning at each planned gantry angle. Eddy current and field homogeneity measurements were conducted using a 24-cm diameter spherical phantom. Gradient and B₀ eddy currents were calculated from the free induction decays (FIDs) resulting from selective excitation of slices located ±5 cm from isocenter. B₀ eddy currents were also calculated from FIDs acquired with nonselective excitation and compared with B₀ eddy current values derived using selective excitation. B₀ inhomogeneities and center frequency offsets were measured by acquiring FIDs with nonselective excitation. Imaging isocenter shifts were measured using a 33x33x10.5 cm³ uniformity linearity (grid) phantom and a 3D true fast imaging with steady-state precession (TrueFISP) sequence used in MRI-guided radiation therapy. Eddy currents were compared to vendor specifications and correlated with the imaging isocenter shifts. Measurements were conducted before and after the MRI-Linac's waveguide was replaced with an updated design to reduce eddy currents.

RESULTS

B₀ eddy currents were highly correlated (r = 0.986, P << 0.001) for measurements made with vs without selective excitation. Transverse (X and Y) axis B₀ eddy currents before and after the waveguide upgrade were out of specification (specification: ≤0.1 μT m/mT for delays < 10 ms) for most of the measured gantry angles. Gradient eddy currents before and after the upgrade were within specifications for the measured gantry angles (≤0.1% for delays < 10 ms). B₀ eddy currents and imaging isocenter shifts were highly correlated (r = 0.965, P << 0.001). After the Linac waveguide upgrade, root mean square (RMS) peak B₀ and gradient eddy currents dropped 45% and 11%, respectively, for delays <10 ms, while imaging isocenter shifts dropped 53%. Isocenter shifts were observed in both phase encode and readout directions. Center frequency offsets were <26 Hz while B₀ inhomogeneities were <33 Hz full width at half maximum (FWHM).

CONCLUSIONS

Imaging isocenter shifts measured in a 0.35-T MRI-Linac were highly correlated with B₀ eddy currents. The eddy currents and imaging isocenter shifts decreased after the MRI-Linac's waveguide was replaced.

Machine QA for the Elekta Unity system: A Report from the Elekta MR-linac consortium

Over the last few years, magnetic resonance image‐guided radiotherapy systems have been introduced into the clinic, allowing for daily online plan adaption. While quality assurance (QA) is similar to conventional radiotherapy systems, there is a need to introduce or modify measurement techniques. As yet, there is no consensus guidance on the QA equipment and test requirements for such systems. Therefore, this report provides an overview of QA equipment and techniques for mechanical, dosimetric, and imaging performance of such systems and recommendation of the QA procedures, particularly for a 1.5T MR‐linac device. An overview of the system design and considerations for QA measurements, particularly the effect of the machine geometry and magnetic field on the radiation beam measurements is given. The effect of the magnetic field on measurement equipment and methods is reviewed to provide a foundation for interpreting measurement results and devising appropriate methods. And lastly, a consensus overview of recommended QA, appropriate methods, and tolerances is provided based on conventional QA protocols. The aim of this consensus work was to provide a foundation for QA protocols, comparative studies of system performance, and for future development of QA protocols and measurement methods.

MRI prostate contouring is not impaired by the use of a radiotherapy image acquisition set-up. An intra- and inter-observer paired comparative analysis with diagnostic set-up images

PURPOSE

The use of MRI for radiotherapy planning purposes is growing but image acquisition using radiotherapy set-ups has impaired image quality. Whether differences in image acquisition set-up could modify organ contouring has not been evaluated. Therefore, we aimed to evaluate differences in contouring between paired of image sets that were acquired in the same scanning session using different parameters.

MATERIAL AND METHODS

Ten patients underwent RT treatment planning with MRI co-registration. MRI was carried out using two different set-ups during the same session, MRI radiotherapy set-ups and MRI diagnostic set-ups. Prostates and rectums were retrospectively contoured in both image sets by 5 radiation oncologists and 4 radiologists. Intra-observer analysis was carried out comparing organ volumes, the Dice coefficient and hausdorff distance values between two contouring rounds. Inter-observer analysis was carried out by comparing individual contours to a generated STAPLE consensus contour, which is considered the gold standard reference.

RESULTS

No significant differences were observed between MRI acquisition set-ups. Significant differences were observed for the dice and hausdorff parameters, comparing individual contours to the STAPLE consensus contour, when analysing diagnostic images between rounds, although raw values were similar.

CONCLUSION

Prostate and rectum contours did not differ significantly when using diagnostic or radiotherapy MRI acquisition set-ups.

Integration of spatial distortion effects in a 4D computational phantom for simulation studies in extra-cranial MRI-guided radiation therapy: Initial results

PURPOSE

Spatial distortions in magnetic resonance imaging (MRI) are mainly caused by inhomogeneities of the static magnetic field, nonlinearities in the applied gradients, and tissue-specific magnetic susceptibility variations. These factors may significantly alter the geometrical accuracy of the reconstructed MR image, thus questioning the reliability of MRI for guidance in image-guided radiation therapy. In this work, we quantified MRI spatial distortions and created a quantitative model where different sources of distortions can be separated. The generated model was then integrated into a four-dimensional (4D) computational phantom for simulation studies in MRI-guided radiation therapy at extra-cranial sites.

METHODS

A geometrical spatial distortion phantom was designed in four modules embedding laser-cut PMMA grids, providing 3520 landmarks in a field of view of (345 × 260 × 480) mm³ . The construction accuracy of the phantom was verified experimentally. Two fast MRI sequences for extra-cranial imaging at 1.5 T were investigated, considering axial slices acquired with online distortion correction, in order to mimic practical use in MRI-guided radiotherapy. Distortions were separated into their sources by acquisition of images with gradient polarity reversal and dedicated susceptibility calculations. Such a separation yielded a quantitative spatial distortion model to be used for MR imaging simulations. Finally, the obtained spatial distortion model was embedded into an anthropomorphic 4D computational phantom, providing registered virtual CT/MR images where spatial distortions in MRI acquisition can be simulated.

RESULTS

The manufacturing accuracy of the geometrical distortion phantom was quantified to be within 0.2 mm in the grid planes and 0.5 mm in depth, including thickness variations and bending effects of individual grids. Residual spatial distortions after MRI distortion correction were strongly influenced by the applied correction mode, with larger effects in the trans-axial direction. In the axial plane, gradient nonlinearities caused the main distortions, with values up to 3 mm in a 1.5 T magnet, whereas static field and susceptibility effects were below 1 mm. The integration in the 4D anthropomorphic computational phantom highlighted that deformations can be severe in the region of the thoracic diaphragm, especially when using axial imaging with 2D distortion correction. Adaptation of the phantom based on patient-specific measurements was also verified, aiming at increased realism in the simulation.

CONCLUSIONS

The implemented framework provides an integrated approach for MRI spatial distortion modeling, where different sources of distortion can be quantified in time-dependent geometries. The computational phantom represents a valuable platform to study motion management strategies in extra-cranial MRI-guided radiotherapy, where the effects of spatial distortions can be modeled on synthetic images in a virtual environment.

Characterization of radiotherapy component impact on MR imaging quality for an MRgRT system

Radiotherapy components of an magnetic resonnace-guided radiotherapy (MRgRT) system can alter the magnetic fields, causing spatial distortion and image deformation, altering imaging and radiation isocenter coincidence and the accuracy of dose calculations. This work presents a characterization of radiotherapy component impact on MR imaging quality in terms of imaging isocenter variation and spatial integrity changes on a 0.35T MRgRT system, pre- and postupgrade of the system. The impact of gantry position, MLC field size, and treatment table power state on imaging isocenter and spatial integrity were investigated. A spatial integrity phantom was used for all tests. Images were acquired for gantry angles 0-330° at 30° increments to assess the impact of gantry position. For MLC and table power state tests all images were acquired at the home gantry position (330°). MLC field sizes ranged from 1.66 to 27.4 cm edge length square fields. Imaging isocenter shift caused by gantry position was reduced from 1.7 mm at gantry 150° preupgrade to 0.9 mm at gantry 120° postupgrade. Maximum spatial integrity errors were 0.5 mm or less pre- and postupgrade for all gantry angles, MLC field sizes, and treatment table power states. However, when the treatment table was powered on, there was significant reduction in SNR. This study showed that gantry position can impact imaging isocenter, but spatial integrity errors were not dependent on gantry position, MLC field size, or treatment table power state. Significant isocenter variation, while reduced postupgrade, is cause for further investigation.

Dosimetric feasibility of brain stereotactic radiosurgery with a 0.35 T MRI-guided linac and comparison vs a C-arm-mounted linac

PURPOSE

MRI is the gold-standard imaging modality for brain tumor diagnosis and delineation. The purpose of this work was to investigate the feasibility of performing brain stereotactic radiosurgery (SRS) with a 0.35 T MRI-guided linear accelerator (MRL) equipped with a double-focused multileaf collimator (MLC). Dosimetric comparisons were made vs a conventional C-arm-mounted linac with a high-definition MLC.

METHODS

The quality of MRL single-isocenter brain SRS treatment plans was evaluated as a function of target size for a series of spherical targets with diameters from 0.6 cm to 2.5 cm in an anthropomorphic head phantom and six brain metastases (max linear dimension = 0.7-1.9 cm) previously treated at our clinic on a conventional linac. Each target was prescribed 20 Gy to 99% of the target volume. Step-and-shoot IMRT plans were generated for the MRL using 11 static coplanar beams equally spaced over 360° about an isocenter placed at the center of the target. Couch and collimator angles are fixed for the MRL. Two MRL planning strategies (VR1 and VR2) were investigated. VR1 minimized the 12 Gy isodose volume while constraining the maximum point dose to be within ±1 Gy of 25 Gy which corresponded to normalization to an 80% isodose volume. VR2 minimized the 12 Gy isodose volume without the maximum dose constraint. For the conventional linac, the TB1 method followed the same strategy as VR1 while TB2 used five noncoplanar dynamic conformal arcs. Plan quality was evaluated in terms of conformity index (CI), conformity/gradient index (CGI), homogeneity index (HI), and volume of normal brain receiving ≥12 Gy (V_12Gy ). Quality assurance measurements were performed with Gafchromic EBT-XD film following an absolute dose calibration protocol.

RESULTS

For the phantom study, the CI of MRL plans was not significantly different compared to a conventional linac (P > 0.05). The use of dynamic conformal arcs and noncoplanar beams with a conventional linac spared significantly more normal brain (P = 0.027) and maximized the CGI, as expected. The mean CGI was 95.9 ± 4.5 for TB2 vs 86.6 ± 3.7 (VR1), 88.2 ± 4.8 (VR2), and 88.5 ± 5.9 (TB1). Each method satisfied a normal brain V_12Gy  ≤ 10.0 cm³ planning goal for targets with diameter ≤2.25 cm. The mean V_12Gy was 3.1 cm³ for TB2 vs 5.5 cm³ , 5.0 cm³ and 4.3 cm³ , for VR1, VR2, and TB1, respectively. For a 2.5-cm diameter target, only TB2 met the V_12Gy planning objective. The MRL clinical brain plans were deemed acceptable for patient treatment. The normal brain V_12Gy was ≤6.0 cm³ for all clinical targets (maximum target volume = 3.51 cm³ ). CI and CGI ranged from 1.12-1.65 and 81.2-88.3, respectively. Gamma analysis pass rates (3%/1mm criteria) exceeded 97.6% for six clinical targets planned and delivered on the MRL. The mean measured vs computed absolute dose difference was -0.1%.

CONCLUSIONS

The MRL system can produce clinically acceptable brain SRS plans for spherical lesions with diameter ≤2.25 cm. Large lesions (>2.25 cm) should be treated with a linac capable of delivering noncoplanar beams.

B0 field homogeneity recommendations, specifications, and measurement units for MRI in radiation therapy

PURPOSE

The purpose is: (a) Relate magnetic resonance imaging (MRI) quality recommendations for radiation therapy (RT) to B₀ field homogeneity; (b) Evaluate manufacturer specifications of B₀ homogeneity for 34 commercial whole-body MRI systems based on the MRI quality recommendations and RT application; (c) Measure field homogeneity in five commercial MRI systems and one commercial MRI-Linac used in RT and compare the results with their B₀ homogeneity specifications.

METHODS

Magnetic resonance imaging quality recommendations for spatial integrity, image blurring, fat saturation, and null banding in RT were developed based on the literature. Guaranteed (maximum) and typical B₀ field homogeneity specifications for various diameter spherical volumes (DSVs) were provided by GE, Philips, Siemens, and Canon. For each system, the DSV that conforms to each MRI quality recommendation and anatomical RT application was estimated based on the manufacturer specifications. B₀ field homogeneity was measured on six MRI systems including Philips (1.5 T), Siemens (1.5 and 3 T), and ViewRay MRI (0.35 T) systems using 24 and 35 cm DSV spherical phantoms. Two measurement techniques were used: (a) MRI using phase contrast field mapping to measure peak-to-peak (pk-pk), volume root mean square (VRMS), and standard deviation (SD); and (b) Magnetic resonance (MR) spectroscopy by acquiring a volumetric free induction decay (FID) to measure full width at half maximum (FWHM). The measurements were used to assess: (a) conformance with the manufacturer specifications; and (b) the relationship between the various field homogeneity measurement units. Measurements were made with and without gradient shimming (gradshim) or second-order active shimming. Multiple comparisons, analysis of variance (ANOVA), and Pearson correlations were performed to assess the dependence of pk-pk, VRMS, SD, and FWHM measurements of field homogeneity on shim volume, level of shim, and MRI system.

RESULTS

For a 40 cm DSV, the B₀ homogeneity specifications ranged from 0.35 to 5 ppm (median = 0.75 ppm) VRMS for 1.5 T systems and 0.2 to 1.4 ppm (median = 0.5 ppm) VRMS for 3 T systems. The usable DSVs ranged from 16 to 49 cm (median = 35 cm) based on the image quality recommendations and the manufacturer specifications. There was general compliance between the six measured field homogeneities and manufacturer specifications although signal dephasing was observed in two systems at < 35 cm DSV. The relationships between pk-pk, VRMS, SD, and FWHM varied based on MRI system, shim volume, and quality of shim. However, VRMS and SD measurements were highly correlated.

CONCLUSIONS

The delineation of the diseased lesion from organs at risk is the main priority for RT. Therefore, field homogeneity performance for RT must minimize image blurring and image artifacts (null bands and signal dephasing) while optimizing spatial integrity and fat saturation. Based on the specifications and recommendations for field homogeneity, some MRI systems are not well suited to meet the strict demands of RT particularly for the large imaging volumes used in body MRI. VRMS and SD measurements of B₀ field homogeneity tend to be more stable and sensitive to field inhomogeneities in RT applications than pk-pk and FWHM.

An assessment of the ExacTrac intrafraction imaging capabilities for flattening filter free prostate stereotactic body radiotherapy

The uncertainties associated with image matching using the ExacTrac® system (BrainLab, Munich, Germany) have been the subject of investigation in the literature for extra-cranial sites. However, the uncertainties involved in the use of ExacTrac in the presence of higher scatter conditions like that for intrafraction imaging of prostate stereotactic radiotherapy utilising unflattened beams is yet to be determined. A prostate phantom was created with 3 implanted gold fiducial markers. This phantom was shifted by 1 mm and 2 mm amounts in the translational planes and by 1° and 2° amounts in the rotational planes and subsequently imaged by ExacTrac during delivery of a clinical SBRT plan. ExacTrac auto-match results were compared to the known offsets with uncertainties calculated. Calculated shifts were shown to be accurate within one standard deviation of the known offsets. Uncertainties were found to vary considerably among the 6 dimensions with matching in the vertical and angle vertical directions having standard deviations of 0.7 mm and 1.3°, respectively. These results agreed with the literature cases for pre-treatment setup and lower scatter condition IMRT intrafraction delivery. Based on these values, probabilities of intrafraction inhibits were calculated based on patient movement and possible fusion tolerances. While the measured uncertainties are adequately defined in order to calculate appropriate target margins, their relatively large magnitudes made choice of intrafraction fusion tolerances problematic. A degree of compromise between the rate of false positives and false negatives is required when implementing ExacTrac into a SBRT prostate protocol.

A review of Image Guided Radiation Therapy in head and neck cancer from 2009-201 - Best Practice Recommendations for RTTs in the Clinic

Radiation therapy (RT) is beneficial in Head and Neck Cancer (HNC) in both the definitive and adjuvant setting. Highly complex and conformal planning techniques are becoming standard practice in delivering increased doses in HNC. A sharp falloff in dose outside the high dose area is characteristic of highly complex techniques and geometric uncertainties must be minimised to prevent under dosage of the target volume and possible over dosage of surrounding critical structures. CTV-PTV margins are employed to account for geometric uncertainties such as set up errors and both interfraction and intrafraction motion. Robust immobilisation and Image Guided Radiation Therapy (IGRT) is also essential in this group of patients to minimise discrepancies in patient position during the treatment course. IGRT has evolved with increased 2-Dimensional (2D) and 3-Dimensional (3D) IGRT modalities available for geometric verification. 2D and 3D IGRT modalities are both beneficial in geometric verification while 3D imaging is a valuable tool in assessing volumetric changes that may have dosimetric consequences for this group of patients. IGRT if executed effectively and efficiently provides clinicians with confidence to reduce CTV-PTV margins thus limiting treatment related toxicities in patients. Accumulated exposure dose from IGRT vary considerably and may be incorporated into the treatment plan to avoid excess dose. However, there are considerable variations in the application of IGRT in RT practice. This paper aims to summarise the advances in IGRT in HNC treatment and provide clinics with recommendations for an IGRT strategy for HNC in the clinic.

Characterizing MR Imaging isocenter variation in MRgRT

We characterized MRI isocenter variation at various gantry positions in two 0.35 T MRgRT systems using two independent methods. First, image center-based quantification was employed on 3D volumetric and 2D cine images of a 24 cm diameter spherical phantom at various gantry positions in the MRI QA mode. The center of the phantom images was identified to quantify the variation of the imaging center at each gantry position. Second, image registration-based quantification was used in radiotherapy mode. 3D volumetric MRIs of a cylindrical phantom were acquired and corresponding image registration from MRI to planning CT was performed. The shifts of the couch were identified to quantify the variation of the imaging center. For verification of noticeable MRI isocenter variation, star-shot pattern measurements with five beams were delivered on the radio-chromic film inserted into the phantom after the couch was shifted. The center of the star-shot pattern was identified to quantify the variation of the imaging center. The proposed methods for measuring MRI isocenter variation were demonstrated with MR-LINAC and MR-60Co systems. Both of the MRgRT systems had field inhomogeneities <5 ppm over a 24 cm diameter spherical volume (DSV) and spatial integrity distortion: <1 mm within 100 mm radius and <2 mm within 175 mm radius. The MRI isocenter of the MR-LINAC system showed noticeable 3D variation (max magnitude: 1.8 mm) compared to that of MR-60Co system (max magnitude: 0.9 mm) relative to the reference gantry positions. In addition, 2D variations (max magnitude) of the MRI isocenter from sagittal cine images were 0.9 mm for the MR-LINAC system and 0.5 mm for the MR-60Co system. Two proposed methods quantified the MRI isocenter variation for various gantry positions in two 0.35 T MRgRT systems. The results of significant isocenter variation in the MR-LINAC system requires further investigation to determine the cause.