Assessment of Geometric Distortion in Six Clinical Scanners Using a 3D-Printed Grid Phantom

Characterization of spatial integrity with active and passive implants in a low-field magnetic resonance linear accelerator scanner

Background and Purpose

Standard imaging protocols can guarantee the spatial integrity of magnetic resonance (MR) images utilized in radiotherapy. However, the presence of metallic implants can significantly compromise this integrity. Our proposed method aims at characterizing the geometric distortions induced by both passive and active implants commonly encountered in planning images obtained from a low-field 0.35 T MR-linear accelerator (LINAC).

Materials and Methods

We designed a spatial integrity phantom defining 1276 control points and covering a field of view of 20x20x20 cm3. This phantom was scanned in a water tank with and without different implants used in hip and shoulder arthroplasty procedures as well as with active cardiac stimulators. The images were acquired with the clinical planning sequence (balanced steady-state free-precession, resolution 1.5x1.5x1.5 mm3). Spatial integrity was assessed by the Euclidian distance between the control point detected on the image and their theoretical locations. A first plane free of artefact (FPFA) was defined to evaluate the spatial integrity beyond the larger banding artefact.

Results

In the region extending up to 20 mm from the largest banding artefacts, the tested passive and active implants could cause distortions up to 2 mm and 3 mm, respectively. Beyond this region the spatial integrity was recovered and the image could be considered as unaffected by the implants.

Conclusions

We characterized the impact of common implants on a low field MR-LINAC planning sequence. These measurements could support the creation of extra margin while contouring organs at risk and target volumes in the vicinity of implants.

Apparent Diffusion Coefficient Reproducibility Across 3 T Scanners in a Breast Diffusion Phantom

BACKGROUND

To date, the accuracy and variability of diffusion-weighted MRI (DW-MRI) metrics have been reported in a limited number of scanner/protocol/coil combinations.

PURPOSE

To evaluate the reproducibility of DW-MRI estimates across multiple scanners and DW-MRI protocols and to assess the effects of using an 8-channel vs. 16-channel breast coil in a breast phantom.

STUDY TYPE

Prospective.

PHANTOM

Breast phantom containing tubes of water and differing polyvinylpyrrolidone (PVP) concentrations with apparent diffusion coefficients (ADCs) matching breast tissue.

FIELD STRENGTH/SEQUENCE

3 T (three standard and one wide bore), three DW-MRI single-shot echo planar imaging protocols of varying acquired spatial resolution.

ASSESSMENT

Accuracy of estimated ADCs was assessed using percent differences (PD) relative to nuclear magnetic resonance spectroscopy-derived reference values. Coefficients of variation (CV) were calculated to determine variation across scanners. CVs based on the median standard deviation (CV_M ) were used to evaluate tube-specific dispersion using 8- or 16-channel coils at a single scanner. ADCs of PVP-containing tubes were additionally normalized by the median water tube ADC to account for temperature effects.

STATISTICAL TESTS

Two-way repeated measures analysis of variance and post hoc tests were used to evaluate differences in ADC, CV, and CV_M across scanners and protocols (α = 0.05).

RESULTS

ADCs were within 11% (interquartile range [IQR] 7%) of reference values and significantly improved to 2% (IQR 7%) after normalization to an internal water reference. Normalization significantly reduced interscanner variability of ADC estimates from 7% to 4%. DW-MRI protocol did not affect ADC accuracy; however, the clinical and higher-resolution clinical protocols resulted in the greatest (9%) and least (6%) interscanner variability, respectively. The 8- and 16-channel receive coils yielded similar accuracy (PD: 12% vs. 16%) and precision (CV_M : 2.7% vs. 2.9%).

DATA CONCLUSION

Normalization by an internal reference improved interscanner ADC reproducibility. High-resolution protocols yielded comparably accurate and significantly less variable ADCs compared to a clinical standard protocol.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

MRI Distortion Correction and Robot-to-MRI Scanner Registration for an MRI-Guided Robotic System

Magnetic resonance imaging (MRI) guided robotic procedures require safe robotic instrument navigation and precise target localization. This depends on reliable tracking of the instrument from MR images, which requires accurate registration of the robot to the scanner. A novel differential image based robot-to-MRI scanner registration approach is proposed that utilizes a set of active fiducial coils, where background subtraction method is employed for coil detection. In order to use the presented preoperative registration approach jointly with the real-time high speed MRI image acquisition and reconstruction methods in real-time interventional procedures, the effects of the geometric MRI distortion in robot to scanner registration is analyzed using a custom distortion mapping algorithm. The proposed approach is validated by a set of target coils placed within the workspace, employing multi-planar capabilities of the scanner. Registration and validation errors are respectively 2.05 mm and 2.63 mm after the distortion correction showing an improvement of respectively 1.08 mm and 0.14 mm compared to the results without distortion correction.

Increased brain volumetric measurement precision from multi-site 3D T1-weighted 3 T magnetic resonance imaging by correcting geometric distortions

PURPOSE

Magnetic resonance imaging (MRI) scanner-specific geometric distortions may contribute to scanner induced variability and decrease volumetric measurement precision for multi-site studies. The purpose of this study was to determine whether geometric distortion correction increases the precision of brain volumetric measurements in a multi-site multi-scanner study.

METHODS

Geometric distortion variation was quantified over a one-year period at 10 sites using the distortion fields estimated from monthly 3D T1-weighted MRI geometrical phantom scans. The variability of volume and distance measurements were quantified using synthetic volumes and a standard quantitative MRI (qMRI) phantom. The effects of geometric distortion corrections on MRI derived volumetric measurements of the human brain were assessed in two subjects scanned on each of the 10 MRI scanners and in 150 subjects with cerebrovascaular disease (CVD) acquired across imaging sites.

RESULTS

Geometric distortions were found to vary substantially between different MRI scanners but were relatively stable on each scanner over a one-year interval. Geometric distortions varied spatially, increasing in severity with distance from the magnet isocenter. In measurements made with the qMRI phantom, the geometric distortion correction decreased the standard deviation of volumetric assessments by 35% and distance measurements by 42%. The average coefficient of variance decreased by 16% in gray matter and white matter volume estimates in the two subjects scanned on the 10 MRI scanners.

CONCLUSION

Geometric distortion correction using an up-to-date correction field is recommended to increase precision in volumetric measurements made from MRI images.

Characterizing geometric distortions of 3D sequences in clinical head MRI

Objective

Phantoms are often used to estimate the geometric accuracy in magnetic resonance imaging (MRI). However, the distortions may differ between anatomical and phantom images. This study aimed to investigate the applicability of a phantom-based and a test-subject-based method in evaluating geometric distortion present in clinical head-imaging sequences.

Materials and methods

We imaged a 3D-printed phantom and test subjects with two MRI scanners using two clinical head-imaging 3D sequences with varying patient-table positions and receiver bandwidths. The geometric distortions were evaluated through nonrigid registrations: the displaced acquisitions were compared against the ideal isocenter positioning, and the varied bandwidth volumes against the volume with the highest bandwidth. The phantom acquisitions were also registered to a computed tomography scan.

Results

Geometric distortion magnitudes increased with larger table displacements and were in good agreement between the phantom and test-subject acquisitions. The effect of increased distortions with decreasing receiver bandwidth was more prominent for test-subject acquisitions.

Conclusion

Presented results emphasize the sensitivity of the geometric accuracy to positioning and imaging parameters. Phantom limitations may become an issue with some sequence types, encouraging the use of anatomical images for evaluating the geometric accuracy.

Differences between 3D isovoxel fat suppression VIBE MRI and CT models of proximal femur osseous anatomy: A preliminary study for bone tumor resection planning

Purpose

To evaluate the osseous anatomy of the proximal femur extracted from a 3D-MRI volumetric interpolated breath-hold (VIBE) sequence using either a Dixon or water excitation (WE) fat suppression method, and to measure the overall difference using CT as a reference standard.

Material and methods

This retrospective study reviewed imaging of adult patients with hip pain who underwent 3D hip MRI and CT. A semi-automatically segmented CT model served as the reference standard, and MRI segmentation was performed manually for each unilateral hip joint. The differences between Dixon-VIBE-3D-MRI vs. CT, and WE-VIBE-3D-MRI vs. CT, were measured. Equivalence tests between Dixon-VIBE and WE-VIBE models were performed with a threshold of 0.1 mm. Bland–Altman plots and Lin’s concordance-correlation coefficient were used to analyze the agreement between WE and Dixon sequences. Subgroup analyses were performed for the femoral head/neck, intertrochanteric, and femoral shaft areas.

Results

The mean and maximum differences between Dixon-VIBE-3D-MRI vs. CT were 0.2917 and 3.4908 mm, respectively, whereas for WE-VIBE-3D-MRI vs. CT they were 0.3162 and 3.1599 mm. The mean differences of the WE and Dixon methods were equivalent (P = 0.0292). However, the maximum difference was not equivalent between the two methods and it was higher in WE method. Lin’s concordance-correlation coefficient showed poor agreement between Dixon and WE methods. The mean differences between the CT and 3D-MRI models were significantly higher in the femoral shaft area (P = 0.0004 for WE and P = 0.0015 for Dixon) than in the other areas. The maximum difference was greatest in the intertrochanteric area for both techniques.

Conclusion

The difference between 3D-MR and CT models were acceptable with a maximal difference below 3.5mm. WE and Dixon fat suppression methods were equivalent. The mean difference was highest at the femoral shaft area, which was off-center from the magnetization field.

A 3D printed phantom to assess MRI geometric distortion

Magnetic Resonance has become a standard imaging modality for target volume delineation and treatment planning in radiation oncology. Geometric distortions, however, have the potential to detrimentally affect both tumour definition and the dose delivered to the target volume. We report the design, fabrication and imaging of a 3D printed unibody MR distortion phantom along with quantitative image analysis.

METHODS

The internal cavity of the phantom is an orthogonal three-dimensional planar lattice, composed of 3mm diameter rods spaced equidistantly at a 20mm centre-centre offset repeating along the X, Y and Z axes. The phantom featured an overall length of 308.5 mm, a width of 246 mm and a height of 264 mm with lines on the external surface for phantom positioning matched to external lasers. The MR phantom was 3D printed in Nylon-12 using an advancement on traditional selective laser sintering (SLS) (HP Jet Fusion 3D - 4200 machine). The phantom was scanned on a Toshiba Aquilion CT scanner to check the integrity of the 3D print and to correct for any resultant issues. The phantom was then filled with NiSO4 solution and scanned on a 3T PET-MR Siemens scanner for selected T1 and T2 sequences, from which distortion vectors were generated and analysed using in-house software written in Python.

RESULTS

All deviations were less than 1 mm, with an average displacement of 0.228 mm. The majority of the deviations are smaller than the 0.692 mm pixel size for this dataset.

CONCLUSION

A cost-effective, 3D printed MRI-phantom was successfully printed and tested for assessing geometric distortion on MRI scanners. The custom phantom with markings for phantom alignment may be considered for radiotherapy departments looking to add MR scanners for simulation and image guidance.

Impact of 3-Dimensional Versus 2-Dimensional Image Distortion Correction on Stereotactic Neurosurgical Navigation Image Fusion Reliability for Images Acquired With Intraoperative Magnetic Resonance Imaging

BACKGROUND

Fusion of preoperative and intraoperative magnetic resonance imaging (iMRI) studies during stereotactic navigation may be very useful for procedures such as tumor resections but can be subject to error because of image distortion.

OBJECTIVE

To assess the impact of 3-dimensional (3D) vs 2-dimensional (2D) image distortion correction on the accuracy of auto-merge image fusion for stereotactic neurosurgical images acquired with iMRI using a head phantom in different surgical positions.

METHODS

T1-weighted intraoperative images of the head phantom were obtained using 1.5T iMRI. Images were postprocessed with 2D and 3D image distortion correction. These studies were fused to T1-weighted preoperative MRI studies performed on a 1.5T diagnostic MRI. The reliability of the auto-merge fusion of these images for 2D and 3D correction techniques was assessed both manually using the stereotactic navigation system and via image analysis software.

RESULTS

Eight surgical positions of the head phantom were imaged with iMRI. Greater image distortion occurred with increased distance from isocenter in all 3 axes, reducing accuracy of image fusion to preoperative images. Visually reliable image fusions were accomplished in 2/8 surgical positions using 2D distortion correction and 5/8 using 3D correction. Three-dimensional correction yielded superior image registration quality as defined by higher maximum mutual information values, with improvements ranging between 2.3% and 14.3% over 2D correction.

CONCLUSION

Using 3D distortion correction enhanced the reliability of surgical navigation auto-merge fusion of phantom images acquired with iMRI across a wider range of head positions and may improve the accuracy of stereotactic navigation using iMRI images.

Systematic review of pre-clinical and clinical devices for magnetic resonance-guided radiofrequency hyperthermia

Clinical trials have demonstrated the therapeutic benefits of adding radiofrequency (RF) hyperthermia (HT) as an adjuvant to radio- and chemotherapy. However, maximum utilization of these benefits is hampered by the current inability to maintain the temperature within the desired range. RF HT treatment quality is usually monitored by invasive temperature sensors, which provide limited data sampling and are prone to infection risks. Magnetic resonance (MR) temperature imaging has been developed to overcome these hurdles by allowing noninvasive 3D temperature monitoring in the target and normal tissues. To exploit this feature, several approaches for inserting the RF heating devices into the MR scanner have been proposed over the years. In this review, we summarize the status quo in MR-guided RF HT devices and analyze trends in these hybrid hardware configurations. In addition, we discuss the various approaches, extract best practices and identify gaps regarding the experimental validation procedures for MR - RF HT, aimed at converging to a common standard in this process.

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.

Geometric inaccuracy and co-registration errors for CT, DynaCT and MRI images used in robotic stereotactic radiosurgery treatment planning

PURPOSE

To measure the combined errors due to geometric inaccuracy and image co-registration on secondary images (dynamic CT angiography (dCTA), 3D DynaCT angiography (DynaCTA), and magnetic resonance images (MRI)) that are routinely used to aid in target delineation and planning for stereotactic radiosurgery (SRS).

METHODS

Three phantoms (one commercial and two in-house built) and two different analysis approaches (commercial and MATLAB based) were used to quantify the magnitude of geometric image distortion and co-registration errors for different imaging modalities within CyberKnife's MultiPlan treatment planning software. For each phantom, the combined errors were reported as a mean target registration error (TRE). The mean TRE's for different intramodality imaging parameters (e.g., mAs, kVp, and phantom set-ups) and for dCTA, DynaCTA, and MRI systems were measured.

RESULTS

Only X-ray based imaging can be performed with the commercial phantom, and the mean TRE ± standard deviation values were large compared to the in-house analysis using MATLAB. With the 3D printed phantom, even drastic changes in treatment planning CT imaging protocols did not greatly influence the mean TRE (<0.5 mm for a 1 mm slice thickness CT). For all imaging modalities, the largest mean TRE was found on DynaCT, followed by T2-weighted MR images (albeit all <1 mm).

CONCLUSIONS

The user may overestimate the mean TRE if the commercial phantom and MultiPlan were used solely. The 3D printed phantom design is a sensitive and suitable quality assurance tool for measuring 3D geometric inaccuracy and co-registration errors across all imaging modalities.

Measuring geometric accuracy in magnetic resonance imaging with 3D-printed phantom and nonrigid image registration

Objective

We aimed to develop a vendor-neutral and interaction-free quality assurance protocol for measuring geometric accuracy of head and brain magnetic resonance (MR) images. We investigated the usability of nonrigid image registration in the analysis and looked for the optimal registration parameters.

Materials and methods

We constructed a 3D-printed phantom and imaged it with 12 MR scanners using clinical sequences. We registered a geometric-ground-truth computed tomography (CT) acquisition to the MR images using an open-source nonrigid-registration-toolbox with varying parameters. We applied the transforms to a set of control points in the CT image and compared their locations to the corresponding visually verified reference points in the MR images.

Results

With optimized registration parameters, the mean difference (and standard deviation) of control point locations when compared to the reference method was (0.17 ± 0.02) mm for the 12 studied scanners. The maximum displacements varied from 0.50 to 1.35 mm or 0.89 to 2.30 mm, with vendors’ distortion correction on or off, respectively.

Discussion

Using nonrigid CT–MR registration can provide a robust and relatively test-object-agnostic method for estimating the intra- and inter-scanner variations of the geometric distortions.