Estimation of geometrically undistorted B(0) inhomogeneity maps

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.

Quantitative Assessment of the Impact of Geometric Distortions and Their Correction on fMRI Data Analyses

Functional magnetic resonance imaging (fMRI) data is typically collected with gradient-echo echo-planar imaging (GE-EPI) sequences, which are particularly prone to the susceptibility artifact as a result of B₀ field inhomogeneity. The component derived from in-plane spin dephasing induces pixel intensity variations and, more critically, geometric distortions. Despite the physical mechanisms underlying the susceptibility artifact being well established, a systematic investigation on the impact of the associated geometric distortions, and the direct comparison of different approaches to tackle them, on fMRI data analyses is missing. Here, we compared two different distortion correction approaches, by acquiring additional: (1) EPI data with reversed phase encoding direction (TOPUP), and (2) standard (and undistorted) GE data at two different echo times (GRE). We first characterized the geometric distortions and the correction approaches based on the estimated ΔB₀ field offset and voxel shift maps, and then conducted three types of analyses on the distorted and corrected fMRI data: (1) registration into structural data, (2) identification of resting-state networks (RSNs), and (3) mapping of task-related brain regions of interest. GRE estimated the largest voxel shifts and more positively impacted the quality of the analyses, in terms of the (significantly lower) cost function of the registration, the (higher) spatial overlap between the RSNs and appropriate templates, and the (significantly higher) sensitivity of the task-related mapping based on the Z-score values of the associated activation maps, although also evident when considering TOPUP. fMRI data should thus be corrected for geometric distortions, with the choice of the approach having a modest, albeit positive, impact on the fMRI analyses.

Optimizing MRI sequences and images for MRI-based stereotactic radiosurgery treatment planning

AIM

Development of MRI sequences and processing methods for the production of images appropriate for direct use in stereotactic radiosurgery (SRS) treatment planning.

BACKGROUND

MRI is useful in SRS treatment planning, especially for patients with brain lesions or anatomical targets that are poorly distinguished by CT, but its use requires further refinement. This methodology seeks to optimize MRI sequences to generate distortion-free and clinically relevant MR images for MRI-only SRS treatment planning.

MATERIALS AND METHODS

We used commercially available SRS MRI-guided radiotherapy phantoms and eight patients to optimize sequences for patient imaging. Workflow involved the choice of correct MRI sequence(s), optimization of the sequence parameters, evaluation of image quality (artifact free and clinically relevant), measurement of geometrical distortion, and evaluation of the accuracy of our offline correction algorithm.

RESULTS

CT images showed a maximum deviation of 1.3 mm and minimum deviation of 0.4 mm from true fiducial position for SRS coordinate definition. Interestingly, uncorrected MR images showed maximum deviation of 1.2 mm and minimum of 0.4 mm, comparable to CT images used for SRS coordinate definition. After geometrical correction, we observed a maximum deviation of 1.1 mm and minimum deviation of only 0.3 mm.

CONCLUSION

Our optimized MRI pulse sequences and image correction technique show promising results; MR images produced under these conditions are appropriate for direct use in SRS treatment planning.

A Robust Method for Estimating B0 Inhomogeneity Field in the Liver by Mitigating Fat Signals and Phase-Wrapping

We developed an optimized and robust method to estimate liver B0 field inhomogeneity for monitoring and correcting susceptibility-induced geometric distortion in magnetic resonance images for precision therapy. A triple-gradient-echo acquisition was optimized for the whole liver B0 field estimation within a single-exhale breath-hold scan on a 3 T scanner. To eliminate chemical-shift artifacts, fat signals were chosen in-phase between 2 echoes with an echo time difference (ΔTE) of 2.3 milliseconds. To avoid phase-wrapping, other 2 echoes provided a large field dynamic range (1/ΔTE) to cover the B0 field inhomogeneity. In addition, using high parallel imaging factor of 4 and a readout-bandwidth of 1955 Hz/pixel, an ∼18-second acquisition time for breath-held scans was achieved. A 2-step, 1-dimensional regularized method for the ΔB0 field map estimation was developed, tested and validated in phantom and patient studies. Our method was validated on a water phantom with fat components and air pockets; it yielded ΔB0-field maps that had no chemical-shift and phase-wrapping artifacts, and it had a <0.5 mm of geometric distortion near the air pockets. The ΔB0-field maps of the patients' abdominal regions were also free from phase-wrapping and chemical-shift artifacts. The maximum field inhomogeneity was found near the lung–liver interface, up to ∼300 Hz, resulting in ∼2 mm of distortions in anatomical images with a readout-bandwidth of 440 Hz/pixel. The field mapping method in the abdominal region is robust; it can be easily integrated in clinical workflow for patient-based quality control of magnetic resonance imaging geometric integrity.

Quantifying the Effect of 3T Magnetic Resonance Imaging Residual System Distortions and Patient-Induced Susceptibility Distortions on Radiation Therapy Treatment Planning for Prostate Cancer

PURPOSE

To investigate the effect of magnetic resonance system- and patient-induced susceptibility distortions from a 3T scanner on dose distributions for prostate cancers.

METHODS AND MATERIALS

Combined displacement fields from the residual system and patient-induced susceptibility distortions were used to distort 17 prostate patient CT images. VMAT dose plans were initially optimized on distorted CT images and the plan parameters transferred to the original patient CT images to calculate a new dose distribution.

RESULTS

Maximum residual mean distortions of 3.19 mm at a radial distance of 25 cm and maximum mean patient-induced susceptibility shifts of 5.8 mm were found using the lowest bandwidth of 122 Hz per pixel. There was a dose difference of <0.5% between distorted and undistorted treatment plans. The 90% confidence intervals of the mean difference between the dCT and CT treatment plans were all within an equivalence interval of (-0.5, 0.5) for all investigated plan quality measures.

CONCLUSIONS

Patient-induced susceptibility distortions at high field strengths in closed bore magnetic resonance scanners are larger than residual system distortions after using vendor-supplied 3-dimensional correction for the delineated regions studied. However, errors in dose due to disturbed patient outline and shifts caused by patient-induced susceptibility effects are below 0.5%.

Prospective analysis of in vivo landmark point-based MRI geometric distortion in head and neck cancer patients scanned in immobilized radiation treatment position: Results of a prospective quality assurance protocol

Purpose

Uncertainties related to geometric distortion are a major obstacle for effectively utilizing MRI in radiation oncology. We aim to quantify the geometric distortion in patient images by comparing their in-treatment position MRIs with the corresponding planning CTs, using CT as the non-distorted gold standard.

Methods

Twenty-one head and neck cancer patients were imaged with MRI as part of a prospective Institutional Review Board approved study. MR images were acquired with a T2 SE sequence (0.5 × 0.5 × 2.5 mm voxel size) in the same immobilization position as in the CTs. MRI to CT rigid registration was then done and geometric distortion comparison was assessed by measuring the corresponding anatomical landmarks on both the MRI and the CT images. Several landmark measurements were obtained including; skin to skin (STS), bone to bone, and soft tissue to soft tissue at specific levels in horizontal and vertical planes of both scans. Inter-observer variability was assessed and interclass correlation (ICC) was calculated.

Results

A total of 430 landmark measurements were obtained. The median distortion for all landmarks in all scans was 1.06 mm (IQR 0.6–1.98). For each patient 48% of the measurements were done in the right-left direction and 52% were done in the anteroposterior direction. The measured geometric distortion was not statistically different in the right-left direction compared to the anteroposterior direction (1.5 ± 1.6 vs. 1.6 ± 1.7 mm, respectively, p = 0.4). The magnitude of distortion was higher in the STS peripheral landmarks compared to the more central landmarks (2.0 ± 1.9 vs. 1.2 ± 1.3 mm, p < 0.0001). The mean distortion measured by observer one was not significantly different compared to observer 2, 3, and 4 (1.05, 1.23, 1.06 and 1.05 mm, respectively, p = 0.4) with ICC = 0.84.

Conclusion

MRI geometric distortions were quantified in radiotherapy planning applications with a clinically insignificant error of less than 2 mm compared to the gold standard CT.

Quantitative assessment of the susceptibility artefact and its interaction with motion in diffusion MRI

In this paper we evaluate the three main methods for correcting the susceptibility-induced artefact in diffusion-weighted magnetic-resonance (DW-MR) data, and assess how correction is affected by the susceptibility field's interaction with motion. The susceptibility artefact adversely impacts analysis performed on the data and is typically corrected in post-processing. Correction strategies involve either registration to a structural image, the application of an acquired field-map or the use of additional images acquired with different phase-encoding. Unfortunately, the choice of which method to use is made difficult by the absence of any systematic comparisons of them. In this work we quantitatively evaluate these methods, by extending and employing a recently proposed framework that allows for the simulation of realistic DW-MR datasets with artefacts. Our analysis separately evaluates the ability for methods to correct for geometric distortions and to recover lost information in regions of signal compression. In terms of geometric distortions, we find that registration-based methods offer the poorest correction. Field-mapping techniques are better, but are influenced by noise and partial volume effects, whilst multiple phase-encode methods performed best. We use our simulations to validate a popular surrogate metric of correction quality, the comparison of corrected data acquired with AP and LR phase-encoding, and apply this surrogate to real datasets. Furthermore, we demonstrate that failing to account for the interaction of the susceptibility field with head movement leads to increased errors when analysing DW-MR data. None of the commonly used post-processing methods account for this interaction, and we suggest this may be a valuable area for future methods development.

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.

MR-guided radiation therapy: transformative technology and its role in the central nervous system

This review article describes advancement of magnetic resonance imaging technologies in radiation therapy planning, guidance, and adaptation of brain tumors. The potential for MR-guided radiation therapy to improve outcomes and the challenges in its adoption are discussed.

Multiatlas approach with local registration goodness weighting for MRI-based electron density mapping of head and neck anatomy

PURPOSE

The growing use of magnetic resonance imaging (MRI) as a substitute for computed tomography-based treatment planning requires the development of effective algorithms to generate electron density maps for treatment planning and patient setup verification. The purpose of this work was to develop a method to synthesize computerized tomography (CT) for MR-only radiotherapy of head and neck cancer patients.

METHODS

The algorithm is based on registration of multiple patient datasets containing both MRI and CT images (a "multiatlas" algorithm). Twelve matched pairs of good quality CT and MRI scans (those without apparent motion and blurring artifacts) were selected from a pool of head and neck cancer patients to form the atlas. All atlas MRI scans were preprocessed to reduce scanner- and patient-induced intensity inhomogeneities and to standardize their intensity histograms. Atlas CT and MRIs were coregistered using a novel bone-to-air replacement technique applied to the CT scans that improves the similarity between CTs and MRIs and facilitates the registration process. For each new patient, all atlas MRIs are deformed initially onto the new patients' MRI. We introduce a generalized registration error (GRE) metric that automatically measures the goodness of local registration between MRI pairs. The final synthetic CT value at each point is a nonlinear GRE-weighted average of the atlas CTs. For evaluation, the leave-one-out technique was used for synthetic CT generation and the mean absolute error (MAE) between the original and synthetic CT was computed over the entire CT image. The impact of our proposed CT-MR registration scheme on the accuracy of the final synthetic CT was also studied. The original treatment plans were also recomputed on the new synthetic CTs and dose-volume histogram metrics were compared. In addition, the two-dimensional (2D) gamma analysis at 1%/1 mm and 2%/2 mm dose difference/distance to agreement was also performed to study the dose distribution at the isocenter.

RESULTS

MAE error (± standard deviation) between the original and the synthetic CTs was 64 ± 10, 113 ± 12, and 130 ± 28 Hounsfield Unit (HU) for the entire image, air, and bone regions respectively. Our results showed that our proposed bone-suppression based CT-MR fusion and GRE-weighted strategy could lower the overall MAE error between the original and synthetic CTs by ~69% and ~34% respectively. Dose recalculation comparison showed highly consistent results between plans based on the synthetic vs. the original CTs. The 2D gamma analysis revealed the pass rate of 95.44 ± 2.5 and 99.36 ± 0.71 for 1%/1 mm and 2%/2 mm criteria respectively. Due to local registration weighting, the method is robust with respect to MRI imaging artifacts.

CONCLUSION

We developed a novel image analysis technique to synthesize CT for head and neck anatomy. Novel methods were introduced to accurately register atlas CTs and MRIs as well as to weight the final electron density maps using local registration goodness estimates. The resulting accuracy is clinically acceptable, at least for these atlas patients.

Stereotactic radiosurgery planning of vestibular schwannomas: Is MRI at 3 Tesla geometrically accurate?

Purpose

MRI is a mandatory requirement to accurately plan Stereotactic Radiosurgery (SRS) for Vestibular Schwannomas. However, MRI may be distorted due not only to inhomogeneity of the static magnetic field and gradients but also due to susceptibility‐induced effects, which are more prominent at higher magnetic fields. We assess geometrical distortions around air spaces and consider MRI protocol requirements for SRS planning at 3 T.

Methods

Hardware‐related distortion and the effect of incorrect shimming were investigated with structured test objects. The magnetic field was mapped over the head on five volunteers to assess susceptibility‐related distortion in the naso‐oro‐pharyngeal cavities (NOPC) and around the internal ear canal (IAC).

Results

Hardware‐related geometric displacements were found to be less than 0.45 mm within the head volume, after distortion correction. Shimming errors can lead to displacements of up to 4 mm, but errors of this magnitude are unlikely to arise in practice. Susceptibility‐related field inhomogeneity was under 3.4 ppm, 2.8 ppm, and 2.7 ppm for the head, NOPC region and IAC region, respectively. For the SRS planning protocol (890 Hz/pixel, approximately 1 mm³ isotropic), susceptibility‐related displacements were less than 0.5 mm (head), and 0.4 mm (IAC and NOPC). Large displacements are possible in MRI examinations undertaken with lower receiver bandwidth values, commonly used in clinical MRI. Higher receiver bandwidth makes the protocol less vulnerable to sub‐optimal shimming. The shimming volume and the CT‐MR co‐registration must be considered jointly.

Conclusion

Geometric displacements can be kept under 1 mm in the vicinity of air spaces within the head at 3 T with appropriate setting of the receiver bandwidth, correct shimming and employing distortion correction.

A new gradient shimming method based on undistorted field map of B0 inhomogeneity

Most existing gradient shimming methods for NMR spectrometers estimate field maps that resolve B0 inhomogeneity spatially from dual gradient-echo (GRE) images acquired at different echo times. However, the distortions induced by B0 inhomogeneity that always exists in the GRE images can result in estimated field maps that are distorted in both geometry and intensity, leading to inaccurate shimming. This work proposes a new gradient shimming method based on undistorted field map of B0 inhomogeneity obtained by a more accurate field map estimation technique. Compared to the traditional field map estimation method, this new method exploits both the positive and negative polarities of the frequency encoded gradients to eliminate the distortions caused by B0 inhomogeneity in the field map. Next, the corresponding automatic post-data procedure is introduced to obtain undistorted B0 field map based on knowledge of the invariant characteristics of the B0 inhomogeneity and the variant polarity of the encoded gradient. The experimental results on both simulated and real gradient shimming tests demonstrate the high performance of this new method.

Assessing the Dosimetric Accuracy of Magnetic Resonance-Generated Synthetic CT Images for Focal Brain VMAT Radiation Therapy

PURPOSE

The purpose of this study was to assess the dosimetric accuracy of synthetic CT (MRCT) volumes generated from magnetic resonance imaging (MRI) data for focal brain radiation therapy.

METHODS AND MATERIALS

A study was conducted in 12 patients with gliomas who underwent both MR and CT imaging as part of their simulation for external beam treatment planning. MRCT volumes were generated from MR images. Patients' clinical treatment planning directives were used to create 12 individual volumetric modulated arc therapy (VMAT) plans, which were then optimized 10 times on each of their respective CT and MRCT-derived electron density maps. Dose metrics derived from optimization criteria, as well as monitor units and gamma analyses, were evaluated to quantify differences between the imaging modalities.

RESULTS

Mean differences between planning target volume (PTV) doses on MRCT and CT plans across all patients were 0.0% (range: -0.1 to 0.2%) for D(95%); 0.0% (-0.7 to 0.6%) for D(5%); and -0.2% (-1.0 to 0.2%) for D(max). MRCT plans showed no significant changes in monitor units (-0.4%) compared to CT plans. Organs at risk (OARs) had average D(max) differences of 0.0 Gy (-2.2 to 1.9 Gy) over 85 structures across all 12 patients, with no significant differences when calculated doses approached planning constraints.

CONCLUSIONS

Focal brain VMAT plans optimized on MRCT images show excellent dosimetric agreement with standard CT-optimized plans. PTVs show equivalent coverage, and OARs do not show any overdose. These results indicate that MRI-derived synthetic CT volumes can be used to support treatment planning of most patients treated for intracranial lesions.