Routine testing of magnetic field homogeneity on clinical MRI systems

MRI scarcity in low- and middle-income countries

Since the introduction of MRI as a sustainable diagnostic modality, global accessibility to its services has revealed a wide discrepancy between populations-leaving most of the population in LMICs without access to this important imaging modality. Several factors lead to the scarcity of MRI in LMICs; for example, inadequate infrastructure and the absence of a dedicated workforce are key factors in the scarcity observed. RAD-AID has contributed to the advancement of radiology globally by collaborating with our partners to make radiology more accessible for medically underserved communities. However, progress is slow and further investment is needed to ensure improved global access to MRI.

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.

Quantification of magnetic susceptibility fingerprint of a 3D linearity medical device

PURPOSE

The study investigates the numerical modelling as well as experimental validation of magnetic susceptibility effects with respect to a 3D linearity phantom used for the quantification of MR image distortions.

METHODS

Magnetic field numerical simulations based on finite difference methods were conducted to generate the susceptibility (χ) model of the MRID^3D phantom. Experimental data was acquired and analyzed for eight different MR scanners to include a wide range of scanning parameters. Distortion vector fields were generated by applying a harmonic analysis based on finite elements methods. Phantom scans for the same setup but with opposite polarities of the frequency encoding gradient were processed in conjunction with the susceptibility modelling to separately quantify three field components due to gradient non-linearities (GNL), B₀ inhomogeneities and χ perturbations.

RESULTS

The numerical modelling showed a significant range of χ value of up to 8.23 ppm, with a mean value of 2.9 ppm. The χ perturbations were found to be mostly present at the end plates of the cylindrical phantom design. The simulations also showed that setup rotations of up to 10° introduced only negligible variations in the χ model of less than 0.1 ppm. This allows for a straightforward practical implementation of the modelling as a single lookup table. After correcting for the χ perturbations, the B₀ inhomogeneities were derived and found to be in good agreement with either the MR system manufacturer specifications or experimental data available in the literature.

CONCLUSIONS

It is possible to accurately model the magnetic susceptibility signature of a 3D linearity device and remove it as a post-processing correction step. This is important as the procedure unlocks the ability of determining both the GNL field and B₀ map of the scanner without the need of extra acquisitions or phantoms.

Evaluation of the influence of susceptibility-induced magnetic field distortions on the precision of contouring intracranial organs at risk for stereotactic radiosurgery

•

45 data sets (18 on a 1.5 T MR and 27 on a 3 T MR) were evaluated for susceptibility induced distortions.

•

Maximum distortions of up to 1.7 mm were found for organs at risk in standard diagnostic settings.

•

Median distortions ranged between 0.1 and 0.2 mm for all organs at risk.

•

Active shimming was estimated to reduce distortions by a factor of 2.3 to 2.9.

•

A safety margin of 1 mm would have encompassed 99.8% of the distortions.

Background and purpose

Magnetic resonance imaging (MRI) is a crucial factor in optimal treatment planning for stereotactic radiosurgery. To further the awareness of possible errors in MRI, this work aimed to investigate the magnitude of susceptibility induced MRI distortions for intracranial organs at risk (OARs) and test the effectiveness of actively shimming these distortions.

Materials and methods

Distortion maps for 45 exams of 42 patients (18 on a 1.5 T MRI scanner, 27 on a 3 T MRI scanner) were calculated based on a high-bandwidth double-echo gradient echo sequence. The investigated OARs were brainstem, chiasm, eyes, and optic nerves. The influence of active shimming was investigated by comparing unshimmed 1.5 T data with shimmed 3 T data and comparing the results to a model based prediction.

Results

The median distortion for the different OARs was found to be between 0.13 and 0.18 mm for 1.5 T and between 0.11 and 0.13 mm for 3 T. The maximum distortion was found to be between 1.3 and 1.7 mm for 1.5 T and between 1.1 and 1.4 mm for 3 T. The variation of values was much higher for 1.5 T than for 3 T across all investigated OARs. Active shimming was found to reduce distortions by a factor of 2.3 to 2.9 compared to the expected values.

Conclusions

Using a safety margin for OARs of 1 mm would have encompassed 99.8% of the distortions. Since distortions are inversely proportional to the readout bandwidth, they can be further reduced by increasing the bandwidth. Additional error sources like gradient nonlinearities need to be addressed separately.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Characterization of hardware-related spatial distortions for IR-PETRA pulse sequence using a brain specific phantom

OBJECTIVE

Inversion recovery-pointwise encoding time reduction with radial acquisition (IR-PETRA) is an effective magnetic resonance (MR) pulse sequence in generating pseudo-CTs. The hardware-related spatial-distortion (HRSD) in MR images potentially deteriorates the accuracy of pseudo-CTs. Thus, we aimed at characterizing HRSD for IR-PETRA.

MATERIALS AND METHODS

gross-HRSD_overall (Euclidean-sum of gross-HRSD_i (i = x, y, z)) for IR-PETRA was assessed using a brain-specific phantom for two MR scanners (1.5 T-Aera and 3.0 T-Prisma). Moreover, hardware imperfections were analyzed by determining gradient-nonlinearity spatial-distortion (GNSD) and B₀-inhomogeneity spatial-distortion (B₀ISD) for magnetization-prepared rapid acquisition gradient-echo (MP-RAGE) which has well-known distortion characteristics.

RESULTS

In 3.0 T, maximum of gross-GNSD_overall (Euclidean-sum of gross-GNSD_i) and gross-B₀ISD for MP-RAGE was 2.77 mm and 0.57 mm, respectively. For this scanner, the mean and maximum of gross-HRSD_overall for IR-PETRA were 0.63 ± 0.38 mm and 1.91 mm, respectively. In 1.5 T, maximum of gross-GNSD_overall and gross-B₀ISD for MP-RAGE was 3.41 mm and 0.78 mm, respectively. The mean and maximum of gross-HRSD_overall for IR-PETRA were 1.02 ± 0.50 mm and 3.12 mm, respectively.

DISCUSSION

The spatial accuracy of MR images, besides being impacted by hardware performance, scanner capabilities, and imaging parameters, is mainly affected by its imaging strategy and data acquisition scheme. In 3.0 T, even without applying vendor correction algorithms, spatial accuracy of IR-PETRA image is sufficient for generating pseudo-CTs. In 1.5 T, distortion-correction is required to provide this accuracy.

An Automatic Image Processing Workflow for Daily Magnetic Resonance Imaging Quality Assurance

The performance of magnetic resonance imaging (MRI) equipment is typically monitored with a quality assurance (QA) program. The QA program includes various tests performed at regular intervals. Users may execute specific tests, e.g., daily, weekly, or monthly. The exact interval of these measurements varies according to the department policies, machine setup and usage, manufacturer's recommendations, and available resources. In our experience, a single image acquired before the first patient of the day offers a low effort and effective system check. When this daily QA check is repeated with identical imaging parameters and phantom setup, the data can be used to derive various time series of the scanner performance. However, daily QA with manual processing can quickly become laborious in a multi-scanner environment. Fully automated image analysis and results output can positively impact the QA process by decreasing reaction time, improving repeatability, and by offering novel performance evaluation methods. In this study, we have developed a daily MRI QA workflow that can measure multiple scanner performance parameters with minimal manual labor required. The daily QA system is built around a phantom image taken by the radiographers at the beginning of day. The image is acquired with a consistent phantom setup and standardized imaging parameters. Recorded parameters are processed into graphs available to everyone involved in the MRI QA process via a web-based interface. The presented automatic MRI QA system provides an efficient tool for following the short- and long-term stability of MRI scanners.

Continuous table acquisition MRI for radiotherapy treatment planning: distortion assessment with a new extended 3D volumetric phantom

PURPOSE

Accurate geometry is required for radiotherapy treatment planning (RTP). When considering the use of magnetic resonance imaging (MRI) for RTP, geometric distortions observed in the acquired images should be considered. While scanner technology and vendor supplied correction algorithms provide some correction, large distortions are still present in images, even when considering considerably smaller scan lengths than those typically acquired with CT in conventional RTP. This study investigates MRI acquisition with a moving table compared with static scans for potential geometric benefits for RTP.

METHODS

A full field of view (FOV) phantom (diameter 500 mm; length 513 mm) was developed for measuring geometric distortions in MR images over volumes pertinent to RTP. The phantom consisted of layers of refined plastic within which vitamin E capsules were inserted. The phantom was scanned on CT to provide the geometric gold standard and on MRI, with differences in capsule location determining the distortion. MRI images were acquired with two techniques. For the first method, standard static table acquisitions were considered. Both 2D and 3D acquisition techniques were investigated. With the second technique, images were acquired with a moving table. The same sequence was acquired with a static table and then with table speeds of 1.1 mm/s and 2 mm/s. All of the MR images acquired were registered to the CT dataset using a deformable B-spline registration with the resulting deformation fields providing the distortion information for each acquisition.

RESULTS

MR images acquired with the moving table enabled imaging of the whole phantom length while images acquired with a static table were only able to image 50%-70% of the phantom length of 513 mm. Maximum distortion values were reduced across a larger volume when imaging with a moving table. Increased table speed resulted in a larger contribution of distortion from gradient nonlinearities in the through-plane direction and an increased blurring of capsule images, resulting in an apparent capsule volume increase by up to 170% in extreme axial FOV regions. Blurring increased with table speed and in the central regions of the phantom, geometric distortion was less for static table acquisitions compared to a table speed of 2 mm/s over the same volume. Overall, the best geometric accuracy was achieved with a table speed of 1.1 mm/s.

CONCLUSIONS

The phantom designed enables full FOV imaging for distortion assessment for the purposes of RTP. MRI acquisition with a moving table extends the imaging volume in the z direction with reduced distortions which could be useful particularly if considering MR-only planning. If utilizing MR images to provide additional soft tissue information to the planning CT, standard acquisition sequences over a smaller volume would avoid introducing additional blurring or distortions from the through-plane table movement.

Patient-induced susceptibility effect on geometric distortion of clinical brain MRI for radiation treatment planning on a 3T scanner

Concerns about the geometric accuracy of MRI in radiation therapy (RT) have been present since its invention. Although modern scanners typically have system levels of geometric accuracy that meet requirements of RT, subject-specific distortion is variable, and methods to in vivo assess and control patient-induced geometric distortion are not yet resolved. This study investigated the nature and magnitude of the subject-induced susceptibility effect on geometric distortions in clinical brain MRI, and tested the feasibility of in vivo quality control using field inhomogeneity mapping. For 19 consecutive patients scanned on a dedicated 3T MR scanner, B0 field inhomogeneity maps were acquired and analyzed to determine subject-induced distortions. For 3D T1 weighted images frequency-encoded with a bandwidth of 180 Hz/pixel, 86.9% of the estimated displacements were <0.5 mm, 97.4% <1 mm, and only 0.1% of displacements > 2 mm. The maximum displacement was <4 mm. The greatest distortions were observed at the interfaces with air at the sinuses. Displacements decayed to less than 1 mm over a distance of 8 mm. Metal surgical wires generated smaller distortions, with an averaged maximum displacement of 0.76 mm. Repeat acquisition of the field maps in 17 patients revealed a within-subject standard deviation of 0.25 ppm, equivalent to 0.22 mm displacement in the frequency-encoding direction in the 3D T1 weighted images. Susceptibility-induced voxel displacements in the brain are generally small, but should be monitored for precision RT. These effects are manageable at 3T and lower fields, and the methods applied can be used to monitor for potential local errors in individual patients, as well as to correct for local distortions as needed.

Characterization of a 3D MEMS fabricated micro-solenoid at 9.4 T

We present for the first time a complete characterization of a micro-solenoid for high resolution MR imaging of mass- and volume-limited samples based on three-dimensional B(0), B(1) per unit current (B(1)(unit)) and SNR maps. The micro-solenoids are fabricated using a fully micro-electromechanical systems (MEMS) compatible process in conjunction with an automatic wire-bonder. We present 15 μm isotropic resolution 3D B(0) maps performed using the phase difference method. The resulting B(0) variation in the range of [-0.07 ppm to -0.157 ppm] around the coil center, compares favorably with the 0.5 ppm limit accepted for MR microscopy. 3D B(1)(unit) maps of 40 μm isotropic voxel size were acquired according to the extended multi flip angle (ExMFA) method. The results demonstrate that the characterized microcoil provides a high and uniform sensitivity distribution around its center (B(1)(unit) = 3.4 mT/A ± 3.86%) which is in agreement with the corresponding 1D theoretical data computed along the coil axis. The 3D SNR maps reveal a rather uniform signal distribution around the coil center with a mean value of 53.69 ± 19%, in good agreement with the analytical 1D data along coil axis in the axial slice. Finally, we prove the microcoil capabilities for MR microscopy by imaging Eremosphaera viridis cells with 18 μm isotropic resolution.

Feasibility of using real-time MRI to measure joint kinematics in 1.5T and open-bore 0.5T systems

PURPOSE

To test the feasibility and accuracy of measuring joint motion with real-time MRI in a 1.5T scanner and in a 0.5T open-bore scanner and to assess the dependence of measurement accuracy on movement speed.

MATERIALS AND METHODS

We developed an MRI-compatible motion phantom to evaluate the accuracy of tracking bone positions with real-time MRI for varying movement speeds. The measurement error was determined by comparing phantom positions estimated from real-time MRI to those measured using optical motion capture techniques. To assess the feasibility of measuring in vivo joint motion, we calculated 2D knee joint kinematics during knee extension in six subjects and compared them to previously reported measurements.

RESULTS

Measurement accuracy decreased as the phantom's movement speed increased. The measurement accuracy was within 2 mm for velocities up to 217 mm/s in the 1.5T scanner and 38 mm/s in the 0.5T scanner. We measured knee joint kinematics with small intraobserver variation (variance of 0.8 degrees for rotation and 3.6% of patellar width for translation).

CONCLUSION

Our results suggest that real-time MRI can be used to measure joint kinematics when 2 mm accuracy is sufficient. They can also be used to prescribe the speed of joint motion necessary to achieve certain measurement accuracy.