A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Three-dimensional assessment of image distortion induced by active cardiac implants in 3.0T CMR

CMR at 3.0T in the presence of active cardiac implants remains a challenge due to susceptibility artifacts. Beyond a signal void that cancels image information, magnetic field inhomogeneities may cause distorted appearances of anatomical structures. Understanding influencing factors and the extent of distortion are a first step towards optimizing the image quality of CMR with active implants at 3.0T. All measurements were obtained at a clinical 3.0T scanner. An in-house designed phantom with a 3D cartesian grid of water filled spheres was used to analyze the distortion caused by four representative active cardiac devices (cardiac loop recorder, pacemaker, 2 ICDs). For imaging a gradient echo (3D-TFE) sequence and a turbo spin echo (2D-TSE) sequence were used. The work defines metrics to quantify the different features of distortion such as changes in size, location and signal intensity. It introduces a specialized segmentation technique based on a reaction-diffusion-equation. The distortion features are dependent on the amount of magnetic material in the active implants and showed a significant increase when measured with the 3D TFE compared to the 2D TSE. This work presents a quantitative approach for the evaluation of image distortion at 3.0T caused by active cardiac implants and serves as foundation for both further optimization of sequences and devices but also for planning of imaging procedures.

Phantom assessment of three-dimensional geometric distortion of a dedicated wide-bore MR-simulator for radiotherapy

This study evaluated the machine-dependent three-dimensional geometric distortion images acquired from a 1.5T 700mm-wide bore MR-simulator based on a large geometric accuracy phantom. With the consideration of radiation therapy (RT) application requirements, every sequence was examined in various combinations of acquisition-orientations and receiver-bandwidths with console-integrated distortion correction enabled. Distortion was repeatedly measured over a six-month period. The distortion measured from the images acquired at the beginning of this period was employed to retrospectively correct the distortion in the subsequent acquisitions. Geometric distortion was analyzed within the largest field-of-view allowed. Six sequences were examined for comprehensive distortion analysis - VIBE, SPACE, TSE, FLASH, BLADE and PETRA. Based on optimal acquisition parameters, their diameter-sphere-volumes (DSVs) of CT-comparable geometric fidelity (where 1mm distortion was allowed) were 333.6mm, 315.1mm, 316.0mm, 318.9mm, 306.2mm and 314.5mm respectively. This was a significant increase from 254.0mm, 245.5mm, 228.9mm, 256.6mm, 230.8mm and 254.2mm DSVs respectively, when images were acquired using un-optimized parameters. The longitudinal stability of geometric distortion and the efficacy of retrospective correction of console-corrected images, based on prior distortion measurements, were inspected using VIBE and SPACE. The retrospectively corrected images achieved over 500mm DSVs with 1mm distortion allowed. The median distortion was below 1mm after retrospective correction, proving that obtaining prior distortion map for subsequent retrospective distortion correction is beneficial. The systematic evaluation of distortion using various combinations of sequence-type, acquisition-orientation and receiver-bandwidth in a six-month time span would be a valuable guideline for optimizing sequence for various RT applications.

Dosimetric impact of geometric distortions in an MRI-only proton therapy workflow for lung, liver and pancreas

In a radiation therapy workflow based on Magnetic Resonance Imaging (MRI), dosimetric errors may arise due to geometric distortions introduced by MRI. The aim of this study was to quantify the dosimetric effect of system-dependent geometric distortions in an MRI-only workflow for proton therapy applied at extra-cranial sites. An approach was developed, in which computed tomography (CT) images were distorted using an MRI displacement map, which represented the MR distortions in a spoiled gradient-echo sequence due to gradient nonlinearities and static magnetic field inhomogeneities. A retrospective study was conducted on 4DCT/MRI digital phantoms and 18 4DCT clinical datasets of the thoraco-abdominal site. The treatment plans were designed and separately optimized for each beam in a beam specific Planning Target Volume on the distorted CT, and the final dose distribution was obtained as the average. The dose was then recalculated in undistorted CT using the same beam geometry and beam weights. The analysis was performed in terms of Dose Volume Histogram (DVH) parameters. No clinically relevant dosimetric impact was observed on organs at risk, whereas in the target structure, geometric distortions caused statistically significant variations in the planned dose DVH parameters and dose homogeneity index (DHI). The dosimetric variations in the target structure were smaller in abdominal cases (ΔD_2%, ΔD_98%, and ΔD_mean all below 0.1% and ΔDHI below 0.003) compared to the lung cases. Indeed, lung patients with tumors isolated inside lung parenchyma exhibited higher dosimetric variations (ΔD_2%≥0.3%, ΔD_98%≥15.9%, ΔD_mean≥3.3% and ΔDHI≥0.102) than lung patients with tumor close to soft tissue (ΔD_2%≤0.4%, ΔD_98%≤5.6%, ΔD_mean≤0.9% and ΔDHI≤0.027) potentially due to higher density variations along the beam path. Results suggest the potential applicability of MRI-only proton therapy, provided that specific analysis is applied for isolated lung tumors.

A modular phantom and software to characterize 3D geometric distortion in MRI

Magnetic resonance imaging (MRI) offers outstanding soft tissue contrast that may reduce uncertainties in target and organ-at-risk delineation and enable online adaptive image-guided treatment. Spatial distortions resulting from non-linearities in the gradient fields and non-uniformity in the main magnetic field must be accounted for across the imaging field-of-view to prevent systematic errors during treatment delivery. This work presents a modular phantom and software application to characterize geometric distortion (GD) within the large field-of-view MRI images required for radiation therapy simulation. The modular phantom is assembled from a series of rectangular foam blocks containing high-contrast fiducial markers in a known configuration. The modular phantom design facilitates transportation of the phantom between different MR scanners and MR-guided linear accelerators and allows the phantom to be adapted to fit different sized bores or coils. The phantom was evaluated using a 1.5 T MR-guided linear accelerator (MR-Linac) and 1.5 T and 3.0 T diagnostic scanners. Performance was assessed by varying acquisition parameters to induce image distortions in a known manner. Imaging was performed using T1 and T2 weighted pulse sequences with 2D and 3D distortion correction algorithms and the receiver bandwidth (BW) varied as 250-815 Hz pixel-1. Phantom set-up reproducibility was evaluated across independent set-ups. The software was validated by comparison with a non-modular phantom. Average geometric distortion was 0.94 ± 0.58 mm for the MR-Linac, 0.90 ± 0.53 mm for the 1.5 T scanner, and 1.15 ± 0.62 mm for the 3.0 T scanner, for a 400 mm diameter volume-of-interest. GD increased, as expected, with decreasing BW, and with the 2D versus 3D correction algorithm. Differences in GD attributed to phantom set-up were 0.13 mm or less. Differences in GD for the two software applications were less than 0.07 mm. A novel modular phantom was developed to evaluate distortions in MR images for radiation therapy applications.

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.

Effect of gradient field nonlinearity distortions in MRI-based attenuation maps for PET reconstruction

PURPOSE

Attenuation correction is a requirement for quantification of the activity distribution in PET. The need to base attenuation correction on MRI instead of CT has arisen with the introduction of integrated PET/MRI systems. The aim was to describe the effect of residual gradient field nonlinearity distortions on PET attenuation correction.

METHODS

MRI distortions caused by gradient field nonlinearity were simulated in CT images used for attenuation correction in PET reconstructions. The simulations yielded radial distortion of up to ±2.3mm at 15cm from the scanner isocentre for distortion corrected images. The mean radial distortion of uncorrected images were 6.3mm at the same distance. Reconstructions of PET data were performed using the distortion corrected images as well as the images where no correction had been applied.

RESULTS

The mean relative difference in reconstructed PET uptake intensity due to incomplete distortion correction was less than ±5%. The magnitude of this difference varied between patients and the size of the distortions remaining after distortion correction.

CONCLUSIONS

Radial distortions of 2mm at 15cm radius from the scanner isocentre lead to PET attenuation correction errors smaller than 5%. Keeping the gradient field nonlinearity distortions below this limit can be a reasonable goal for MRI systems used for attenuation correction in PET for quantification purposes. A higher geometrical accuracy may, however, be warranted for quantification of peripheral lesions. These distortions can, e.g., be controlled at acceptance testing and subsequent quality assurance intervals.

Characterization of 3D geometric distortion of magnetic resonance imaging scanners commissioned for radiation therapy planning

OBJECTIVE

To develop a method for the assessment and characterization of 3D geometric distortion as part of routine quality assurance for MRI scanners commissioned for Radiation Therapy planning.

MATERIALS AND METHODS

In this study, the in-plane and through-plane geometric distortions on a 1.5T GE MRI-SIM unit are characterized and the 2D and 3D correction algorithms provided by the vendor are evaluated. We used a phantom developed by GE Healthcare that covers a large field of view of 500mm, and consists of layers of foam embedded with a matrix of ellipsoidal markers. An in-house Java-based software module was developed to automatically assess the geometric distortion by calculating the center of each marker using the center of mass method, correcting of gross rotation errors and comparing the corrected positions with a CT gold standard data set. Spatial accuracy of typical pulse sequences used in RT planning was assessed (2D T1/T2 FSE, 3D CUBE, T1 SPGR) using the software. The accuracy of vendor specific geometric distortion correction (GDC) algorithms was quantified by measuring distortions before and after the application of the 2D and 3D correction algorithms.

RESULTS

Our algorithm was able to accurately calculate geometric distortion with sub-pixel precision. For all typical MR sequences used in Radiotherapy, the vendor's GDC was able to substantially reduce the distortions. Our results showed also that the impact of the acquisition produced a maximum variation of 0.2mm over a radial distance of 200mm. It has been shown that while the 2D correction algorithm remarkably reduces the in-plane geometric distortion, 3D geometric distortion further reduced the geometric distortion by correcting both in-plane and through-plane distortions in all acquisitions.

CONCLUSION

The presented methods represent a valuable tool for routine quality assurance of MR applications that require stringent spatial accuracy assessment such as radiotherapy. The phantom used in this study provides three dimensional arrays of control points. These tools and the detailed results can be also used for developing new geometric distortion correction algorithms or improving the existing ones.

Phantom-based characterization of distortion on a magnetic resonance imaging simulator for radiation oncology

One of the major issues potentially limiting treatment planning with solely MR images is the possibility of geometric distortion inherent in MR images. We designed a large distortion phantom containing a 3D array of spheres and proposed a three-dimensional (3D) approach to determine the distortion of MR image volume. The approach to overcome partially filled spheres is also presented. The phantom was assembled with a 3D array of spheres filled with contrast and was scanned with a 3T MRI simulator. A 3D whole-sphere or half-sphere template is used to match the image pattern. The half-sphere template is used when the normalized cross-correlation value for the whole-sphere template is below a predetermined threshold. Procrustes method was applied to remove the shift induced by rotation and translation of the phantom. Then the distortion map was generated. Accuracy of the method was verified using CT images of a small phantom of the same design. The analysis of the small phantom showed that the method is accurate with an average offset of estimated sphere center 0.12 ± 0.04 mm. The Procrustes analysis estimated the rotation angle to be 1.95° and 0.01°, respectively, when the phantom was placed at 2° and 0° from the ceiling laser. The analysis showed that on the central plane through the magnet center, the average displacement is less than 1 mm for all radii. At distal planes, when the radius is less than 18 cm, the average displacement is less than 1 mm. However, the average displacement is over 1 mm but still less than 1.5 mm for larger radii. A large distortion phantom was assembled and analysis software was developed to characterize distortions in MRI scans. The use of two templates helps reduce the potential impact of residual air bubbles in some of the spheres.

Fast concomitant gradient field and field inhomogeneity correction for spiral cardiac imaging

Non-Cartesian imaging provides many advantages in terms of flexibility, functionality, and speed. However, a major drawback to these imaging methods is off-resonance distortion artifacts. These artifacts manifest as blurring in spiral imaging. Common techniques that remove the off-resonance field inhomogeneity distortion effects are not sufficient, because the high order concomitant gradient fields are nontrivial for common imaging conditions, such as imaging 5 cm off isocenter in an 1.5 T scanner. Previous correction algorithms are either slow or do not take into account the known effects of concomitant gradient fields along with the field inhomogeneities. To ease the correction, the distortion effects are modeled as a non-stationary convolution problem. In this work, two fast and accurate postgridding algorithms are presented and analyzed. These methods account for both the concomitant field effects and the field inhomogeneities. One algorithm operates in the frequency domain and the other in the spatial domain. To take advantage of their speed and accuracy, the algorithms are applied to a real-time cardiac study and a high-resolution cardiac study. Both of the presented algorithms provide for a practical solution to the off-resonance problem in spiral imaging.

Geometric distortion in clinical MRI systems Part I: evaluation using a 3D phantom

Recently, a 3D phantom that can provide a comprehensive and accurate measurement of the geometric distortion in MRI has been developed. Using this phantom, a full assessment of the geometric distortion in a number of clinical MRI systems (GE and Siemens) has been carried out and detailed results are presented in this paper. As expected, the main source of geometric distortion in modern superconducting MRI systems arises from the gradient field nonlinearity. Significantly large distortions with maximum absolute geometric errors ranged between 10 and 25 mm within a volume of 240 x 240 x 240 mm(3) were observed when imaging with the new generation of gradient systems that employs shorter coils. By comparison, the geometric distortion was much less in the older-generation gradient systems. With the vendor's correction method, the geometric distortion measured was significantly reduced but only within the plane in which these 2D correction methods were applied. Distortion along the axis normal to the plane was, as expected, virtually unchanged. Two-dimensional correction methods are a convenient approach and in principle they are the only methods that can be applied to correct geometric distortion in a single slice or in multiple noncontiguous slices. However, these methods only provide an incomplete solution to the problem and their value can be significantly reduced if the distortion along the normal of the correction plane is not small.

Method for a detailed measurement of image intensity nonuniformity in magnetic resonance imaging

In magnetic resonance imaging (MRI), the MR signal intensity can vary spatially and this spatial variation is usually referred to as MR intensity nonuniformity. Although the main source of intensity nonuniformity arises from B1 inhomogeneity of the coil acting as a receiver and/or transmitter, geometric distortion also alters the MR signal intensity. It is useful on some occasions to have these two different sources be separately measured and analyzed. In this paper, we present a practical method for a detailed measurement of the MR intensity nonuniformity. This method is based on the same three-dimensional geometric phantom that was recently developed for a complete measurement of the geometric distortion in MR systems. In this paper, the contribution to the intensity nonuniformity from the geometric distortion can be estimated and thus, it provides a mechanism for estimation of the intensity nonuniformity that reflects solely the spatial characteristics arising from B1. Additionally, a comprehensive scheme for characterization of the intensity nonuniformity based on the new measurement method is proposed. To demonstrate the method, the intensity nonuniformity in a 1.5 T Sonata MR system was measured and is used to illustrate the main features of the method.