Technical Note: Characterization and correction of gradient nonlinearity induced distortion on a 1.0 T open bore MR-SIM

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.

Patient specific distortion detection and mitigation in MR images used for stereotactic radiosurgery

Objective. In MRI-based radiation therapy planning, mitigating patient-specific distortion with standard high bandwidth scans can result in unnecessary sacrifices of signal to noise ratio. This study investigates a technique for distortion detection and mitigation on a patient specific basis. Approach. Fast B0 mapping was performed using a previously developed technique for high-resolution, large dynamic range field mapping without the need for phase unwrapping algorithms. A phantom study was performed to validate the method. Distortion mitigation was validated by reducing geometric distortion with increased acquisition bandwidth and confirmed by both the B0 mapping technique and manual measurements. Images and contours from 25 brain stereotactic radiosurgery patients and 95 targets were analyzed to estimate the range of geometric distortions expected in the brain and to estimate bandwidth required to keep all treatment targets within the ±0.5 mm iso-distortion contour. Main Results. The phantom study showed, at 3 T, the technique can measure distortions with a mean absolute error of 0.12 mm (0.18 ppm), and a maximum error of 0.37 mm (0.6 ppm). For image acquisition at 3 T and 1.0 mm resolution, mean absolute distortion under 0.5 mm in patients required bandwidths from 109 to 200 Hz px−1 for patients with the least and most distortion, respectively. Maximum absolute distortion under 0.5 mm required bandwidths from 120 to 390 Hz px−1. Significance. The method for B0 mapping was shown to be valid and may be applied to assess distortion clinically. Future work will adapt the readout bandwidth to prospectively mitigate distortion with the goal to improve radiosurgery treatment outcomes by reducing healthy tissue exposure.

Task group 284 report: magnetic resonance imaging simulation in radiotherapy: considerations for clinical implementation, optimization, and quality assurance

The use of dedicated magnetic resonance simulation (MR-SIM) platforms in Radiation Oncology has expanded rapidly, introducing new equipment and functionality with the overall goal of improving the accuracy of radiation treatment planning. However, this emerging technology presents a new set of challenges that need to be addressed for safe and effective MR-SIM implementation. The major objectives of this report are to provide recommendations for commercially available MR simulators, including initial equipment selection, siting, acceptance testing, quality assurance, optimization of dedicated radiation therapy specific MR-SIM workflows, patient-specific considerations, safety, and staffing. Major contributions include guidance on motion and distortion management as well as MRI coil configurations to accommodate patients immobilized in the treatment position. Examples of optimized protocols and checklists for QA programs are provided. While the recommendations provided here are minimum requirements, emerging areas and unmet needs are also highlighted for future development.

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.

ReUINet: A fast GNL distortion correction approach on a 1.0 T MRI-Linac scanner

PURPOSE

The hybrid system combining a magnetic resonance imaging (MRI) scanner with a linear accelerator (Linac) has become increasingly desirable for tumor treatment because of excellent soft tissue contrast and nonionizing radiation. However, image distortions caused by gradient nonlinearity (GNL) can have detrimental impacts on real-time radiotherapy using MRI-Linac systems, where accurate geometric information of tumors is essential.

METHODS

In this work, we proposed a deep convolutional neural network-based method to efficiently recover undistorted images (ReUINet) for real-time image guidance. The ReUINet, based on the encoder-decoder structure, was created to learn the relationship between the undistorted images and distorted images. The ReUINet was pretrained and tested on a publically available brain MR image dataset acquired from 23 volunteers. Then, transfer learning was adopted to implement the pretrained model (i.e., network with optimal weights) on the experimental three-dimensional (3D) grid phantom and in-vivo pelvis image datasets acquired from the 1.0 T Australian MRI-Linac system.

RESULTS

Evaluations on the phantom (768 slices) and pelvis data (88 slices) showed that the ReUINet achieved improvement over 15 times and 45 times on computational efficiency in comparison with standard interpolation and GNL-encoding methods, respectively. Moreover, qualitative and quantitative results demonstrated that the ReUINet provided better correction results than the standard interpolation method, and comparable performance compared to the GNL-encoding approach.

CONCLUSIONS

Validated by simulation and experimental results, the proposed ReUINet showed promise in obtaining accurate MR images for the implementation of real-time MRI-guided radiotherapy.

A phantom based laser marking workflow to visually assess geometric image distortion in magnetic resonance guided radiotherapy

Magnetic resonance (MR)-only workflows require quality assurance due to potential dosimetric impacts of using geometry distorted MR images in radiotherapy planning. MR-visible silicone-based fiducials were arranged in regular 3D structures to cover extended imaging volumes. The scanner’s patient marking workflow with a 2-axes movable laser bridge allowed to visually check geometric distortions of each MR reconstructed fiducial against its true position in 3D space. A measurement resolution and uncertainty of the order of 0.5 mm in sagittal and coronal, and 1 mm in transversal direction was found. The proposed workflow required 1 min of evaluation time per fiducial position, and a 9 min 3D MR volume acquisition.

Audit feasibility for geometric distortion in magnetic resonance imaging for radiotherapy

BACKGROUND AND PURPOSE

Magnetic Resonance Imaging (MRI) is increasingly being used in radiotherapy (RT). However, geometric distortions are a known challenge of using MRI in RT. The aim of this study was to demonstrate feasibility of a national audit of MRI geometric distortions. This was achieved by assessing large field of view (FOV) MRI distortions on a number of scanners used clinically for RT.

MATERIALS AND METHODS

MRI scans of a large FOV MRI geometric distortion phantom were acquired on 11 MRI scanners that are used clinically for RT in the UK. The mean and maximum distortions and variance between scanners were reported at different distances from the isocentre.

RESULTS

For a small FOV representing a brain (100-150 mm from isocentre) all distortions were < 2 mm except for the maximum distortion of one scanner. For a large FOV representing a head and neck/pelvis (200-250 mm from isocentre) mean distortions were < 2 mm except for one scanner, maximum distortions were > 10 mm in some cases. The variance between scanners was low and was found to increase with distance from isocentre.

CONCLUSIONS

This study demonstrated feasibility of the technique to be repeated in a country wide geometric distortion audit of all MRI scanners used clinically for RT. Recommendations were made for performing such an audit and how to derive acceptable limits of distortion in such an audit.

A Systematic Review of the Clinical Implementation of Pelvic Magnetic Resonance Imaging-Only Planning for External Beam Radiation Therapy

The use of magnetic resonance (MR) imaging scans alone for radiation therapy treatment planning (MR-only planning) has been highlighted as one method of improving patient outcomes. Recent technologic advances have meant that introducing MR-only planning to the clinic is becoming a reality, with several specialist radiation therapy clinics using this technique for treatment. As such, substantial efforts are being made to introduce this technique into wide-spread clinical implementation. A systematic review of publications investigating the clinical implementation of pelvic MR-only radiation therapy treatment planning was undertaken following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. The Medline, Embase, Scopus, Science Direct, Cumulative Index to Nursing and Allied Health Literature, and Web of Science databases were searched (timespan: all years to January 2, 2019). Twenty-six articles met the inclusion criteria. The studies were grouped into the following categories: (1) MR acquisition and synthetic computed tomography generation verification, (2) MR distortion quantification and phantom development, (3) clinical validation of patient treatment positioning in an MR-only workflow, and (4) MR-only commissioning processes. Key conclusions from this review are (1) MR-only planning has been implemented clinically for prostate cancer treatments; (2) a substantial amount of work remains to translate MR-only planning into widespread clinical implementation for all pelvic sites; (3) MR scanner distortions are no longer a barrier to MR-only planning, but they must be managed appropriately; (4) MR-only-based patient positioning verification shows promise, but limited evidence is reported in the literature and further investigation is required; and (5) a number of MR-only commissioning processes have been reported, which can aid centers as they undertake local commissioning; however, this needs to be formalized in guidance from national bodies.

Modern Radiation Therapy Planning and Delivery

Modern radiation therapy treatment planning and delivery is a complex process that relies on advanced imaging and computing technology as well as expertise from the medical team. The process begins with simulation imaging, in which three-dimensional computed tomography images (or magnetic resonance images in some cases) are used to characterize the patient anatomy. From there, the radiation oncologist delineates the relevant target/tumor volumes and normal tissue and communicates the goals for treatment planning. The planning process attempts to generate a radiation therapy treatment plan that will deliver a therapeutic dose of radiation to the tumor while sparing nearby normal tissue.

MRI-based treatment planning for liver stereotactic body radiotherapy: validation of a deep learning-based synthetic CT generation method

OBJECTIVE

The purpose of this work is to develop and validate a learning-based method to derive electron density from routine anatomical MRI for potential MRI-based SBRT treatment planning.

METHODS

We proposed to integrate dense block into cycle generative adversarial network (GAN) to effectively capture the relationship between the CT and MRI for CT synthesis. A cohort of 21 patients with co-registered CT and MR pairs were used to evaluate our proposed method by the leave-one-out cross-validation. Mean absolute error, peak signal-to-noise ratio and normalized cross-correlation were used to quantify the imaging differences between the synthetic CT (sCT) and CT. The accuracy of Hounsfield unit (HU) values in sCT for dose calculation was evaluated by comparing the dose distribution in sCT-based and CT-based treatment planning. Clinically relevant dose-volume histogram metrics were then extracted from the sCT-based and CT-based plans for quantitative comparison.

RESULTS

The mean absolute error, peak signal-to-noise ratio and normalized cross-correlation of the sCT were 72.87 ± 18.16 HU, 22.65 ± 3.63 dB and 0.92 ± 0.04, respectively. No significant differences were observed in the majority of the planning target volume and organ at risk dose-volume histogram metrics ( p  >  0.05). The average pass rate of γ analysis was over 99% with 1%/1 mm acceptance criteria on the coronal plane that intersects with isocenter.

CONCLUSION

The image similarity and dosimetric agreement between sCT and original CT warrant further development of an MRI-only workflow for liver stereotactic body radiation therapy.

ADVANCES IN KNOWLEDGE

This work is the first deep-learning-based approach to generating abdominal sCT through dense-cycle-GAN. This method can successfully generate the small bony structures such as the rib bones and is able to predict the HU values for dose calculation with comparable accuracy to reference CT images.

Geometric and dosimetric impact of anatomical changes for MR-only radiation therapy for the prostate

PURPOSE

With the move towards magnetic resonance imaging (MRI) as a primary treatment planning modality option for men with prostate cancer, it becomes critical to quantify the potential uncertainties introduced for MR-only planning. This work characterized geometric and dosimetric intra-fractional changes between the prostate, seminal vesicles (SVs), and organs at risk (OARs) in response to bladder filling conditions.

MATERIALS AND METHODS

T2-weighted and mDixon sequences (3-4 time points/subject, at 1, 1.5 and 3.0 T with totally 34 evaluable time points) were acquired in nine subjects using a fixed bladder filling protocol (bladder void, 20 oz water consumed pre-imaging, 10 oz mid-session). Using mDixon images, Magnetic Resonance for Calculating Attenuation (MR-CAT) synthetic computed tomography (CT) images were generated by classifying voxels as muscle, adipose, spongy, and compact bone and by assignment of bulk Hounsfield Unit values. Organs including the prostate, SVs, bladder, and rectum were delineated on the T2 images at each time point by one physician. The displacement of the prostate and SVs was assessed based on the shift of the center of mass of the delineated organs from the reference state (fullest bladder). Changes in dose plans at different bladder states were assessed based on volumetric modulated arc radiotherapy (VMAT) plans generated for the reference state.

RESULTS

Bladder volume reduction of 70 ± 14% from the final to initial time point (relative to the final volume) was observed in the subject population. In the empty bladder condition, the dose delivered to 95% of the planning target volume (PTV) (D95%) reduced significantly for all cases (11.53 ± 6.00%) likely due to anterior shifts of prostate/SVs relative to full bladder conditions. D15% to the bladder increased consistently in all subjects (42.27 ± 40.52%). Changes in D15% to the rectum were patient-specific, ranging from -23.93% to 22.28% (-0.76 ± 15.30%).

CONCLUSIONS

Variations in the bladder and rectal volume can significantly dislocate the prostate and OARs, which can negatively impact the dose delivered to these organs. This warrants proper preparation of patients during treatment and imaging sessions, especially when imaging required longer scan times such as MR protocols.

Large field of view distortion assessment in a low-field MR-linac

PURPOSE

MR-guided radiation therapy (RT) offers unparalleled soft tissue contrast for localization and target tracking. However, MRI distortions may be detrimental to high precision RT. This work characterizes the gradient nonlinearity (GNL) and total distortions over the first year of clinical operation of a 0.35T MR-linac.

METHODS

For GNL characterization, an in-house large field of view (FOV) phantom (60 × 42.5 × 55 cm³ , >6000 spherical landmarks) was configured and scanned at four timepoints with forward/reverse read polarities (Gradient Echo sequence, FA/TR/TE = 28°/30 ms/6 ms). GNL was measured in Anterior-Posterior (AP), Left-Right (LR), and Superior-Inferior (SI) frequency-encoding directions based on deviation of the auto-segmented landmark centroids between rigidly registered MR and CT images and assessed based on radial distance from magnet isocenter. Total distortion was assessed using a 30 × 30 cm² grid phantom oriented along the cardinal axes over >1 year of operation.

RESULTS

The scanner's spatial integrity within the first ~10 months was stable (maximum total distortion variation = 10/6/8%, maximum distortion = 1.41/0.99/1.56 mm in Axial/Coronal/Sagittal planes, respectively). GNL distortions measured during this time period <10 cm from isocenter were (-0.74, 0.45), (-0.67, 0.53), and (-0.86, 0.70) mm in AP/LR/SI directions. In the 10-20 cm range, <1.5% of the distortions exceeded 2 mm in the AP and LR axes while <4% of the distortions exceeded 2 mm for SI. After major repairs and magnet re-shim, detectable changes were observed in total and GNL distortions (20% reduction in AP and 36% increase in SI direction in the 20-25 cm range). Across all timepoints and axes, 38-53% of landmarks in the 20-25 cm range were displaced by >1 mm.

CONCLUSIONS

GNL distortions were negligible within a 10 cm radius from isocenter. However, in the periphery, non-negligible distortions of up to ~7 mm were observed, which may necessitate GNL corrections for MR-IGRT for treatment sites distant from magnet isocenter.

FMEA of MR-Only Treatment Planning in the Pelvis

Purpose

To evaluate the implementation of a magnetic resonance (MR)-only workflow (ie, implementing MR simulation as the primary planning modality) using failure mode and effects analysis (FMEA) in comparison with a conventional multimodality (MR simulation in conjunction with computed tomography simulation) workflow for pelvis external beam planning.

Methods and Materials

To perform the FMEA, a multidisciplinary 9-member team was assembled and developed process maps, identified potential failure modes (FMs), and assigned numerical values to the severity (S), frequency of occurrence (O), and detectability (D) of those FMs. Risk priority numbers (RPNs) were calculated via the product of S, O, and D as a metric for evaluating relative patient risk. An alternative 3-digit composite number (SOD) was computed to emphasize high-severity FMs. Fault tree analysis identified the causality chain leading to the highest-severity FM.

Results

Seven processes were identified, 3 of which were shared between workflows. Image fusion and target delineation subprocesses using the conventional workflow added 9 and 10 FMs, respectively, with 6 RPNs >100. By contrast, synthetic computed tomography generation introduced 3 major subprocesses and propagated 46 unique FMs, 15 with RPNs >100. For the conventional workflow, the largest RPN scores were introduced by image fusion (RPN range, 120-192). For the MR-only workflow, the highest RPN scores were from inaccuracies in target delineation resulting from misinterpretation of MR images (RPN = 240) and insufficient management of patient- and system-level distortions (RPN = 210 and 168, respectively). Underestimation (RPN = 140) or overestimation (RPN = 192) of bone volume produced higher RPN scores. The highest SODs for both workflows were related to changes in target location because of internal anatomy changes (conventional = 961, MR-only = 822).

Conclusions

FMEA identified areas for mitigating risk in MR-only pelvis RTP, and SODs identified high-severity process modes. Efforts to develop a quality management program to mitigate high FMs are underway.

Using synthetic CT for partial brain radiation therapy: Impact on image guidance

PURPOSE

Recent advancements in synthetic computed tomography (synCT) from magnetic resonance (MR) imaging data have made MRI-only treatment planning feasible in the brain, although synCT performance for image guided radiation therapy (IGRT) is not well understood. This work compares geometric equivalence of digitally reconstructed radiographs (DRRs) from CTs and synCTs for brain cancer patients and quantifies performance for partial brain IGRT.

METHODS AND MATERIALS

Ten brain cancer patients (12 lesions, 7 postsurgical) underwent MR-SIM and CT-SIM. SynCTs were generated by combining ultra-short echo time, T1, T2, and fluid attenuation inversion recovery datasets using voxel-based weighted summation. SynCT and CT DRRs were compared using patient-specific thresholding and assessed via overlap index, Dice similarity coefficient, and Jaccard index. Planar IGRT images for 22 fractions were evaluated to quantify differences between CT-generated DRRs and synCT-generated DRRs in 6 quadrants. Previously validated software was implemented to perform 2-dimensional (2D)-2D rigid registrations using normalized mutual information. Absolute (planar image/DRR registration) and relative (differences between synCT and CT DRR registrations) shifts were calculated for each axis and 3-dimensional vector difference. A total of 1490 rigid registrations were assessed.

RESULTS

DRR agreements in anteroposterior and lateral views for overlap index, Dice similarity coefficient, and Jaccard index were >0.95. Normalized mutual information results were equivalent in 75% of quadrants. Rotational registration results were negligible (<0.07°). Statistically significant differences between CT and synCT registrations were observed in 9/18 matched quadrants/axes (P < .05). The population average absolute shifts were 0.77 ± 0.58 and 0.76 ± 0.59 mm for CT and synCT, respectively, for all axes/quadrants. Three-dimensional vectors were <2 mm in 77.7 ± 10.8% and 76.5 ± 7.2% of CT and synCT registrations, respectively. SynCT DRRs were sensitive in postsurgical cases (vector displacements >2 mm in affected quadrants).

CONCLUSIONS

DRR synCT geometry was robust. Although statistically significant differences were observed between CT and synCT registrations, results were not clinically significant. Future work will address synCT generation in postsurgical settings.

Per-organ assessment of subject-induced susceptibility distortion for MR-only male pelvis treatment planning

Background

Patient-specific distortions, particularly near tissue/air interfaces, require assessment for magnetic resonance (MR) only radiation treatment planning (RTP). However, patients are dynamic due to changes in physiological status during imaging sessions. This work investigated changes in subject-induced susceptibility distortions to pelvic organs at different bladder states to support pelvis MR-only RTP.

Methods

Pelvises of 9 healthy male volunteers were imaged at 1.0 Tesla (T), 1.5 T, and 3.0 T. Subject-induced susceptibility distortion field maps were generated using a dual-echo gradient-recalled echo (GRE) sequence with B₀ field maps obtained from the phase difference between the two echoes acquired at several bladder volume states (3–4/subject, 32 overall). T2 turbo spin echo images were also acquired at each bladder state for organ delineation. Magnet central frequency was tracked over time. Distortion map differences and boxplots were computed to characterize changes within the clinical target volume (CTV), bladder, seminal vesicles, and prostate volumes.

Results

The time between the initial and final B0 maps was 42.6 ± 13.9 (range: 13.2–62.1) minutes with minimal change in magnet central frequency (0.02 ± 0.05 mm (range: − 0.06 – 0.12 mm)). Subject-induced susceptibility distortion across all bladder states, field strengths, and subjects was relatively small (1.4–1.9% of all voxels in the prostate and seminal vesicles were distorted > 0.5 mm). In the bladder, no voxels exhibited distortions > 1 mm. An extreme case acquired at 3.0 T with a large volume of rectal air yielded 27.4–34.6% of voxels within the CTVs had susceptibility-induced distortions > 0.5 mm across all time points.

Conclusions

Our work suggests that subject-induced susceptibility distortions caused by bladder/rectal conditions are generally small and subject-dependent. Local changes may be non-negligible within the CTV, thus proper management of filling status is warranted. Future work evaluating the impact of multiple models to accommodate for extreme status changes may be advantageous.

MRI-only treatment planning: benefits and challenges

Over the past decade, the application of magnetic resonance imaging (MRI) has increased, and there is growing evidence to suggest that improvements in the accuracy of target delineation in MRI-guided radiation therapy may improve clinical outcomes in a variety of cancer types. However, some considerations should be recognized including patient motion during image acquisition and geometric accuracy of images. Moreover, MR-compatible immobilization devices need to be used when acquiring images in the treatment position while minimizing patient motion during the scan time. Finally, synthetic CT images (i.e. electron density maps) and digitally reconstructed radiograph images should be generated from MRI images for dose calculation and image guidance prior to treatment. A short review of the concepts and techniques that have been developed for implementation of MRI-only workflows in radiation therapy is provided in this document.

Magnetic Resonance Imaging-Guided Adaptive Radiation Therapy: A "Game Changer" for Prostate Treatment?

Radiation therapy to the prostate involves increasingly sophisticated delivery techniques and changing fractionation schedules. With a low estimated α/β ratio, a larger dose per fraction would be beneficial, with moderate fractionation schedules rapidly becoming a standard of care. The integration of a magnetic resonance imaging (MRI) scanner and linear accelerator allows for accurate soft tissue tracking with the capacity to replan for the anatomy of the day. Extreme hypofractionation schedules become a possibility using the potentially automated steps of autosegmentation, MRI-only workflow, and real-time adaptive planning. The present report reviews the steps involved in hypofractionated adaptive MRI-guided prostate radiation therapy and addresses the challenges for implementation.

MRI-Related Geometric Distortions in Stereotactic Radiotherapy Treatment Planning: Evaluation and Dosimetric Impact

In view of their superior soft tissue contrast compared to computed tomography, magnetic resonance images are commonly involved in stereotactic radiosurgery/radiotherapy applications for target delineation purposes. It is known, however, that magnetic resonance images are geometrically distorted, thus deteriorating dose delivery accuracy. The present work focuses on the assessment of geometric distortion inherent in magnetic resonance images used in stereotactic radiosurgery/radiotherapy treatment planning and attempts to quantitively evaluate the consequent impact on dose delivery. The geometric distortions for 3 clinical magnetic resonance protocols (at both 1.5 and 3.0 T) used for stereotactic radiosurgery/radiotherapy treatment planning were evaluated using a recently proposed phantom and methodology. Areas of increased distortion were identified at the edges of the imaged volume which was comparable to a brain scan. Although mean absolute distortion did not exceed 0.5 mm on any spatial axis, maximum detected control point disposition reached 2 mm. In an effort to establish what could be considered as acceptable geometric uncertainty, highly conformal plans were utilized to irradiate targets of different diameters (5-50 mm). The targets were mispositioned by 0.5 up to 3 mm, and dose–volume histograms and plan quality indices clinically used for plan evaluation and acceptance were derived and used to investigate the effect of geometrical uncertainty (distortion) on dose delivery accuracy and plan quality. The latter was found to be strongly dependent on target size. For targets less than 20 mm in diameter, a spatial disposition of the order of 1 mm could significantly affect (>5%) plan acceptance/quality indices. For targets with diameter greater than 2 cm, the corresponding disposition was found greater than 1.5 mm. Overall results of this work suggest that efficacy of stereotactic radiosurgery/radiotherapy applications could be compromised in case of very small targets lying distant from the scanner’s isocenter (eg, the periphery of the brain).

Dosimetric evaluation of synthetic CT for magnetic resonance-only based radiotherapy planning of lung cancer

Background

Interest in MR-only treatment planning for radiation therapy is growing rapidly with the emergence of integrated MRI/linear accelerator technology. The purpose of this study was to evaluate the feasibility of using synthetic CT images generated from conventional Dixon-based MRI scans for radiation treatment planning of lung cancer.

Methods

Eleven patients who underwent whole-body PET/MR imaging following a PET/CT exam were randomly selected from an ongoing prospective IRB-approved study. Attenuation maps derived from the Dixon MR Images and atlas-based method was used to create CT data (synCT). Treatment planning for radiation treatment of lung cancer was optimized on the synCT and subsequently copied to the registered CT (planCT) for dose calculation. Planning target volumes (PTVs) with three sizes and four different locations in the lung were planned for irradiation. The dose-volume metrics comparison and 3D gamma analysis were performed to assess agreement between the synCT and CT calculated dose distributions.

Results

Mean differences between PTV doses on synCT and CT across all the plans were −0.1% ± 0.4%, 0.1% ± 0.5%, and 0.4% ± 0.5% for D95, D98 and D100, respectively. Difference in dose between the two datasets for organs at risk (OARs) had average differences of −0.14 ± 0.07 Gy, 0.0% ± 0.1%, and −0.1% ± 0.2% for maximum spinal cord, lung V20, and heart V40 respectively. In patient groups based on tumor size and location, no significant differences were observed in the PTV and OARs dose-volume metrics (p > 0.05), except for the maximum spinal-cord dose when the target volumes were located at the lung apex (p = 0.001). Gamma analysis revealed a pass rate of 99.3% ± 1.1% for 2%/2 mm (dose difference/distance to agreement) acceptance criteria in every plan.

Conclusions

The synCT generated from Dixon-based MRI allows for dose calculation of comparable accuracy to the standard CT for lung cancer treatment planning. The dosimetric agreement between synCT and CT calculated doses warrants further development of a MR-only workflow for radiotherapy of lung cancer.

Optimization of a novel large field of view distortion phantom for MR-only treatment planning

PURPOSE

MR-only treatment planning requires images of high geometric fidelity, particularly for large fields of view (FOV). However, the availability of large FOV distortion phantoms with analysis software is currently limited. This work sought to optimize a modular distortion phantom to accommodate multiple bore configurations and implement distortion characterization in a widely implementable solution.

METHOD AND MATERIALS

To determine candidate materials, 1.0 T MR and CT images were acquired of twelve urethane foam samples of various densities and strengths. Samples were precision-machined to accommodate 6 mm diameter paintballs used as landmarks. Final material candidates were selected by balancing strength, machinability, weight, and cost. Bore sizes and minimum aperture width resulting from couch position were tabulated from the literature (14 systems, 5 vendors). Bore geometry and couch position were simulated using MATLAB to generate machine-specific models to optimize the phantom build. Previously developed software for distortion characterization was modified for several magnet geometries (1.0 T, 1.5 T, 3.0 T), compared against previously published 1.0 T results, and integrated into the 3D Slicer application platform.

RESULTS

All foam samples provided sufficient MR image contrast with paintball landmarks. Urethane foam (compressive strength ∼1000 psi, density ~20 lb/ft³ ) was selected for its accurate machinability and weight characteristics. For smaller bores, a phantom version with the following parameters was used: 15 foam plates, 55 × 55 × 37.5 cm³ (L×W×H), 5,082 landmarks, and weight ~30 kg. To accommodate > 70 cm wide bores, an extended build used 20 plates spanning 55 × 55 × 50 cm³ with 7,497 landmarks and weight ~44 kg. Distortion characterization software was implemented as an external module into 3D Slicer's plugin framework and results agreed with the literature.

CONCLUSION

The design and implementation of a modular, extendable distortion phantom was optimized for several bore configurations. The phantom and analysis software will be available for multi-institutional collaborations and cross-validation trials to support MR-only planning.

Image-based gradient non-linearity characterization to determine higher-order spherical harmonic coefficients for improved spatial position accuracy in magnetic resonance imaging

PURPOSE

Spatial position accuracy in magnetic resonance imaging (MRI) is an important concern for a variety of applications, including radiation therapy planning, surgical planning, and longitudinal studies of morphologic changes to study neurodegenerative diseases. Spatial accuracy is strongly influenced by gradient linearity. This work presents a method for characterizing the gradient non-linearity fields on a per-system basis, and using this information to provide improved and higher-order (9th vs. 5th) spherical harmonic coefficients for better spatial accuracy in MRI.

METHODS

A large fiducial phantom containing 5229 water-filled spheres in a grid pattern is scanned with the MR system, and the positions all the fiducials are measured and compared to the corresponding ground truth fiducial positions as reported from a computed tomography (CT) scan of the object. Systematic errors from off-resonance (i.e., B0) effects are minimized with the use of increased receiver bandwidth (±125kHz) and two acquisitions with reversed readout gradient polarity. The spherical harmonic coefficients are estimated using an iterative process, and can be subsequently used to correct for gradient non-linearity. Test-retest stability was assessed with five repeated measurements on a single scanner, and cross-scanner variation on four different, identically-configured 3T wide-bore systems.

RESULTS

A decrease in the root-mean-square error (RMSE) over a 50cm diameter spherical volume from 1.80mm to 0.77mm is reported here in the case of replacing the vendor's standard 5th order spherical harmonic coefficients with custom fitted 9th order coefficients, and from 1.5mm to 1mm by extending custom fitted 5th order correction to the 9th order. Minimum RMSE varied between scanners, but was stable with repeated measurements in the same scanner.

CONCLUSIONS

The results suggest that the proposed methods may be used on a per-system basis to more accurately calibrate MR gradient non-linearity coefficients when compared to vendor standard corrections.

Distortion inherent to magnetic resonance imaging can lead to geometric miss in radiosurgery planning

PURPOSE

Anatomic distortion is present in all magnetic resonance imaging (MRI) data because of nonlinearity of gradient fields; it measures up to several millimeters. We evaluated the potential for uncorrected MRI to lead to geometric miss of the target volume in stereotactic radiosurgery (SRS).

METHODS AND MATERIALS

Twenty-eight SRS cases were studied retrospectively. MRI scans were corrected for gradient nonlinearity distortion in 3 dimensions, and gross tumor volumes (GTVs) were contoured. The manufacturer-specified distortion field was then reapplied to GTV masks to allow measurement of GTV displacement in uncorrected images. The uncorrected GTV was used for SRS planning, and the dose received by the true (corrected) GTV was measured.

RESULTS

Median displacement of the GTV resulting from gradient distortion was 1.2 mm (interquartile range, 0.1-2.3 mm), with a minimum of 0 mm and a maximum of 3.9 mm. Eight of the 28 cases met a priori criteria for "geometric miss."

CONCLUSIONS

Although MRI distortion is often subtle on visual inspection, there is a significant clinical impact of this distortion on SRS planning. Distortion-corrected MRI should uniformly be used for intracranial radiosurgery planning because uncorrected MRI can lead to potential geometric miss.

Image Guided Radiation Therapy Using Synthetic Computed Tomography Images in Brain Cancer

PURPOSE

The development of synthetic computed tomography (CT) (synCT) derived from magnetic resonance (MR) images supports MR-only treatment planning. We evaluated the accuracy of synCT and synCT-generated digitally reconstructed radiographs (DRRs) relative to CT and determined their performance for image guided radiation therapy (IGRT).

METHODS AND MATERIALS

Magnetic resonance simulation (MR-SIM) and CT simulation (CT-SIM) images were acquired of an anthropomorphic skull phantom and 12 patient brain cancer cases. SynCTs were generated using fluid attenuation inversion recovery, ultrashort echo time, and Dixon data sets through a voxel-based weighted summation of 5 tissue classifications. The DRRs were generated from the phantom synCT, and geometric fidelity was assessed relative to CT-generated DRRs through bounding box and landmark analysis. An offline retrospective analysis was conducted to register cone beam CTs (n=34) to synCTs and CTs using automated rigid registration in the treatment planning system. Planar MV and KV images (n=37) were rigidly registered to synCT and CT DRRs using an in-house script. Planar and volumetric registration reproducibility was assessed and margin differences were characterized by the van Herk formalism.

RESULTS

Bounding box and landmark analysis of phantom synCT DRRs were within 1 mm of CT DRRs. Absolute planar registration shift differences ranged from 0.0 to 0.7 mm for phantom DRRs on all treatment platforms and from 0.0 to 0.4 mm for volumetric registrations. For patient planar registrations, the mean shift differences were 0.4 ± 0.5 mm (range, -0.6 to 1.6 mm), 0.0 ± 0.5 mm (range, -0.9 to 1.2 mm), and 0.1 ± 0.3 mm (range, -0.7 to 0.6 mm) for the superior-inferior (S-I), left-right (L-R), and anterior-posterior (A-P) axes, respectively. The mean shift differences in volumetric registrations were 0.6 ± 0.4 mm (range, -0.2 to 1.6 mm), 0.2 ± 0.4 mm (range, -0.3 to 1.2 mm), and 0.2 ± 0.3 mm (range, -0.2 to 1.2 mm) for the S-I, L-R, and A-P axes, respectively. The CT-SIM and synCT derived margins were <0.3 mm different.

CONCLUSION

DRRs generated by synCT were in close agreement with CT-SIM. Planar and volumetric image registrations to synCT-derived targets were comparable with CT for phantom and patients. This validation is the next step toward MR-only planning for the brain.

Spatial Precision in Magnetic Resonance Imaging-Guided Radiation Therapy: The Role of Geometric Distortion

Because magnetic resonance imaging-guided radiation therapy (MRIgRT) offers exquisite soft tissue contrast and the ability to image tissues in arbitrary planes, the interest in this technology has increased dramatically in recent years. However, intrinsic geometric distortion stemming from both the system hardware and the magnetic properties of the patient affects MR images and compromises the spatial integrity of MRI-based radiation treatment planning, given that for real-time MRIgRT, precision within 2 mm is desired. In this article, we discuss the causes of geometric distortion, describe some well-known distortion correction algorithms, and review geometric distortion measurements from 12 studies, while taking into account relevant imaging parameters. Eleven of the studies reported phantom measurements quantifying system-dependent geometric distortion, while 2 studies reported simulation data quantifying magnetic susceptibility-induced geometric distortion. Of the 11 studies investigating system-dependent geometric distortion, 5 reported maximum measurements less than 2 mm. The simulation studies demonstrated that magnetic susceptibility-induced distortion is typically smaller than system-dependent distortion but still nonnegligible, with maximum distortion ranging from 2.1 to 2.6 mm at a field strength of 1.5 T. As expected, anatomic landmarks containing interfaces between air and soft tissue had the largest distortions. The evidence indicates that geometric distortion reduces the spatial integrity of MRI-based radiation treatment planning and likely diminishes the efficacy of MRIgRT. Better phantom measurement techniques and more effective distortion correction algorithms are needed to achieve the desired spatial precision.