Use of spherical harmonic deconvolution methods to compensate for nonlinear gradient effects on MRI images

Improving fat saturation robustness in outer extremity MRI with a local shim coil insert

PURPOSE

High-quality fat suppression is essential for various MRI applications. In musculoskeletal imaging, poor fat suppression caused by severe B0 inhomogeneity can obscure important lesions, potentially leading to inaccurate diagnoses. This problem is particularly exacerbated in off-isocenter imaging, where conventional shimming using second-order spherical harmonic shim coils often proves inadequate due to elevated B0 inhomogeneity. To address this challenge, we configured a simple local shim insert to provide additional localized B0 shimming for off-isocenter regions, offering a practical hardware solution.

METHODS

We designed and constructed a seven-channel shim coil and evaluated its performance in comparison to conventional second-order spherical harmonic shimming within a targeted volume near the scanner bore. The coil was tested with both phantom and in vivo studies using a clinical 3 T MRI scanner.

RESULTS

The improved B0 homogeneity achieved with the local shim coil significantly enhanced fat saturation (fat-sat) uniformity across the imaged volumes. This improvement was particularly beneficial for areas far from the scanner isocenter, where B0 inhomogeneity is most severe. Our results indicated a 40% reduction in RMS error of the B0 field for elbow imaging and 35% for hand imaging, highlighting substantial improvements in B0 field homogeneity. Additionally, the image quality score increased by 1 point for both hand and elbow images, reflecting enhanced fat-sat quality.

CONCLUSION

The simple local shim insert we configured improves fat-sat capability in both hand and elbow imaging. It offers the potential for improving off-isocenter musculoskeletal MRI.

Characterization and inter-scanner reproducibility of geometric distortion on a small footprint, high-performance, head-specific 0.5 T scanner

BACKGROUND

Magnetic resonance imaging (MRI) offers superior soft tissue contrast and essential imaging capabilities for modern medicine. MRI is increasingly being used in applications that require a high degree of spatial fidelity; however, distortions are a well-known limitation of the modality. The mid-field (0.3 T ≤ B0 < 1 T) has advantages in this respect due to being less susceptible to patient-induced distortions.

PURPOSE

The purpose of this work was to characterize the geometric fidelity of a short-bore, head-specific, 0.5T MRI system.

METHODS

Assessment of spatial fidelity was performed using a 3D gradient recalled echo (GRE) acquisition on a commercial distortion phantom using the validated distortion analysis software provided. B0-induced distortions were measured using a 3D field map. Inter-scanner reproducibility was assessed across four distinct systems of identical make and model, while intra-scanner repeatability was assessed at one site over six repeat measurements.

RESULTS

Inter-scanner reproducibility measured an average 95th percentile distortion over 100 and 180 mm DSV of 0.15 ± 0.03  and 0.33 ± 0.05 mm. Average 95th percentile distortions due to B0 field inhomogeneities over 100 and 180 mm DSV were 0.02 ± 0.01  and 0.07 ± 0.02 mm. Intra-scanner repeatability measured the uncertainty in distortion values to be 0.020 ± 0.005 mm.

CONCLUSION

The total residual distortions measured in this phantom study were less than half the recommended value required for radiosurgery and significantly better than data published from other MR systems. This demonstrates that in addition to the compact footprint of the Synaptive 0.5T scanner, it exceeds current standards for geometric accuracy.

MRI distortion correction is associated with improved local control in stereotactic radiotherapy for brain metastases

Distortions in brain MRI caused by gradient nonlinearities may reach several millimeters, thus distortion correction is strongly recommended for radiotherapy treatment planning. However, the significance of MRI distortion correction on actual clinical outcomes has not been described yet. Therefore, we investigated the impact of planning MRI distortion correction on subsequent local control in a historic series of 419 brain metastases in 189 patients treated with stereotactic radiotherapy between 01/2003 and 04/2015. Local control was evaluated using a volumetric extension of the RANO-BM criteria. The predictive significance of distortion correction was assessed using competing risk analysis. In this cohort, 2D distortion-corrected MRIs had been used for treatment planning in 52.5% (220/419) of lesions, while uncorrected MRIs had been employed in 47.5% (199/419) of metastases. 2D distortion correction was associated with improved local control (Cumulative incidence of local progression at 12 months: 14.3% vs. 21.2% and at 24 months: 18.7% vs. 28.6%, p = 0.038). In multivariate analysis, adjusting for histology, baseline tumor volume, interval between MRI and treatment delivery, year of planning MRI, biologically effective dose and adjuvant Whole-brain radiotherapy, use of distortion correction remained significantly associated with improved local control (HR 0.55, p = 0.020). This is the first study to clinically evaluate the impact of MRI gradient nonlinearity distortion correction on local control in stereotactic radiotherapy for brain metastases. In this historic series, we found significantly higher local control when using 2D corrected vs. uncorrected MRI studies for treatment planning. These results stress the importance of assuring that MR images used for radiotherapy treatment planning are properly distortion-corrected.

Three-Dimensional Variable Slab-Selective Projection Acquisition Imaging

Three-dimensional (3D) projection acquisition (PA) imaging has recently gained attention because of its advantages, such as achievability of very short echo time, less sensitivity to motion, and undersampled acquisition of projections without sacrificing spatial resolution. However, larger subjects require a stronger Nyquist criterion and are more likely to be affected by outer-volume signals outside the field of view (FOV), which significantly degrades the image quality. Here, we proposed a variable slab-selective projection acquisition (VSS-PA) method to mitigate the Nyquist criterion and effectively suppress aliasing streak artifacts in 3D PA imaging. The proposed method involves maintaining the vertical orientation of the slab-selective gradient for frequency-selective spin excitation and the readout gradient for data acquisition. As VSS-PA can selectively excite spins only in the width of the desired FOV in the projection direction during data acquisition, the effective size of the scanned object that determines the Nyquist criterion can be reduced. Additionally, unwanted signals originating from outside the FOV (e.g., aliasing streak artifacts) can be effectively avoided. The mitigation of the Nyquist criterion owing to VSS-PA was theoretically described and confirmed through numerical simulations and phantom and human lung experiments. These experiments further showed that the aliasing streak artifacts were nearly suppressed.

Improving MR axon radius estimation in human white matter using spiral acquisition and field monitoring

PURPOSE

To compare MR axon radius estimation in human white matter using a multiband spiral sequence combined with field monitoring to the current state-of-the-art echo-planar imaging (EPI)-based approach.

METHODS

A custom multiband spiral sequence was used for diffusion-weighted imaging at ultra-highb $$ b $$ -values. Field monitoring and higher order image reconstruction were employed to greatly reduce artifacts in spiral images. Diffusion weighting parameters were chosen to match a state-of-the art EPI-based axon radius mapping protocol. The spiral approach was compared to the EPI approach by comparing the image signal-to-noise ratio (SNR) and performing a test-retest study to assess the respective variability and repeatability of axon radius mapping. Effective axon radius estimates were compared over white matter voxels and along the left corticospinal tract.

RESULTS

Increased SNR and reduced artifacts in spiral images led to reduced variability in resulting axon radius maps, especially in low-SNR regions. Test-retest variability was reduced by a factor of approximately 1.5 using the spiral approach. Reduced repeatability due to significant bias was found for some subjects in both spiral and EPI approaches, and attributed to scanner instability, pointing to a previously unknown limitation of the state-of-the-art approach.

CONCLUSION

Combining spiral readouts with field monitoring improved mapping of the effective axon radius compared to the conventional EPI approach.

Open-source hardware and software for the measurement, characterization, reporting, and correction of geometric distortion in MRI

BACKGROUND

Geometric distortion is a serious problem in MRI, particularly in MRI guided therapy. A lack of affordable and adaptable tools in this area limits research progress and harmonized quality assurance.

PURPOSE

To develop and test a suite of open-source hardware and software tools for the measurement, characterization, reporting, and correction of geometric distortion in MRI.

METHODS

An open-source python library was developed, comprising modules for parametric phantom design, data processing, spherical harmonics, distortion correction, and interactive reporting. The code was used to design and manufacture a distortion phantom consisting of 618 oil filled markers covering a sphere of radius 150 mm. This phantom was imaged on a CT scanner and a novel split-bore 1.0 T MRI magnet. The CT images provide distortion-free dataset. These data were used to test all modules of the open-source software.

RESULTS

All markers were successfully extracted from all images. The distorted MRI markers were mapped to undistorted CT data using an iterative search approach. Spherical harmonics reconstructed the fitted gradient data to 1.0 ± 0.6% of the input data. High resolution data were reconstructed via spherical harmonics and used to generate an interactive report. Finally, distortion correction on an independent data set reduced distortion inside the DSV from 5.5 ± 3.1 to 1.6 ± 0.8 mm.

CONCLUSION

Open-source hardware and software for the measurement, characterization, reporting, and correction of geometric distortion in MRI have been developed. The utility of these tools has been demonstrated via their application on a novel 1.0 T split bore magnet.

Nonlinear Gradient Field Mapping Using a Spherical Grid Phantom for 3 and 7 Tesla MR Imaging Systems Equipped with High-performance Gradient Coils

Recent high-performance gradient coils are fabricated mainly at the expense of spatial linearity. In this study, we measured the spatial nonlinearity of the magnetic field generated by the gradient coils of two MRI systems with high-performance gradient coils. The nonlinearity of the gradient fields was measured using 3D gradient echo sequences and a spherical phantom with a built-in lattice structure. The spatial variation of the gradient field was approximated to the 3rd order polynomials. The coefficients of the polynomials were calculated using the steepest descent method. The geometric distortion of the acquired 3D MR images was corrected using the polynomials and compared with the 3D images corrected using the harmonic functions provided by the MRI venders. As a result, it was found that the nonlinearity correction formulae provided by the vendors were insufficient and needed to be verified or corrected using a geometric phantom such as used in this study.

Quantitative characterization of breast lesions and normal fibroglandular tissue using compartmentalized diffusion-weighted model: comparison of intravoxel incoherent motion and restriction spectrum imaging

BACKGROUND

To compare the compartmentalized diffusion-weighted models, intravoxel incoherent motion (IVIM) and restriction spectrum imaging (RSI), in characterizing breast lesions and normal fibroglandular tissue.

METHODS

This prospective study enrolled 152 patients with 157 histopathologically verified breast lesions (41 benign and 116 malignant). All patients underwent a full-protocol preoperative breast MRI, including a multi-b-value DWI sequence. The diffusion parameters derived from the mono-exponential model (ADC), IVIM model (Dt, Dp, f), and RSI model (C1, C2, C3, C1C2, F1, F2, F3, F1F2) were quantitatively measured and then compared among malignant lesions, benign lesions and normal fibroglandular tissues using Kruskal-Wallis test. The Mann-Whitney U-test was used for the pairwise comparisons. Diagnostic models were built by logistic regression analysis. The ROC analysis was performed using five-fold cross-validation and the mean AUC values were calculated and compared to evaluate the discriminative ability of each parameter or model.

RESULTS

Almost all quantitative diffusion parameters showed significant differences in distinguishing malignant breast lesions from both benign lesions (other than C2) and normal fibroglandular tissue (all parameters) (all P < 0.0167). In terms of the comparisons of benign lesions and normal fibroglandular tissues, the parameters derived from IVIM (Dp, f) and RSI (C1, C2, C1C2, F1, F2, F3) showed significant differences (all P < 0.005). When using individual parameters, RSI-derived parameters-F1, C1C2, and C2 values yielded the highest AUCs for the comparisons of malignant vs. benign, malignant vs. normal tissue and benign vs. normal tissue (AUCs = 0.871, 0.982, and 0.863, respectively). Furthermore, the combined diagnostic model (IVIM + RSI) exhibited the highest diagnostic efficacy for the pairwise discriminations (AUCs = 0.893, 0.991, and 0.928, respectively).

CONCLUSIONS

Quantitative parameters derived from the three-compartment RSI model have great promise as imaging indicators for the differential diagnosis of breast lesions compared with the bi-exponential IVIM model. Additionally, the combined model of IVIM and RSI achieves superior diagnostic performance in characterizing breast lesions.

Quality requirements for MRI simulation in cranial stereotactic radiotherapy: a guideline from the German Taskforce "Imaging in Stereotactic Radiotherapy"

Accurate Magnetic Resonance Imaging (MRI) simulation is fundamental for high-precision stereotactic radiosurgery and fractionated stereotactic radiotherapy, collectively referred to as stereotactic radiotherapy (SRT), to deliver doses of high biological effectiveness to well-defined cranial targets. Multiple MRI hardware related factors as well as scanner configuration and sequence protocol parameters can affect the imaging accuracy and need to be optimized for the special purpose of radiotherapy treatment planning. MRI simulation for SRT is possible for different organizational environments including patient referral for imaging as well as dedicated MRI simulation in the radiotherapy department but require radiotherapy-optimized MRI protocols and defined quality standards to ensure geometrically accurate images that form an impeccable foundation for treatment planning. For this guideline, an interdisciplinary panel including experts from the working group for radiosurgery and stereotactic radiotherapy of the German Society for Radiation Oncology (DEGRO), the working group for physics and technology in stereotactic radiotherapy of the German Society for Medical Physics (DGMP), the German Society of Neurosurgery (DGNC), the German Society of Neuroradiology (DGNR) and the German Chapter of the International Society for Magnetic Resonance in Medicine (DS-ISMRM) have defined minimum MRI quality requirements as well as advanced MRI simulation options for cranial SRT.

A correction algorithm for improved magnetic field monitoring with distal field probes

PURPOSE

A significant source of artifacts in MRI are field fluctuations. Field monitoring is a new technology that allows measurement of field dynamics during a scan via "field probes," which can be used to improve image reconstruction. Ideally, probes are located within the volume where gradients produce nominally linear field patterns. However, in some situations probes must be located far from isocenter where rapid field variation can arise, leading to erroneous field-monitoring characterizations and images. This work aimed to develop an algorithm that improves the robustness of field dynamics in these situations.

METHODS

The algorithm is split into three components. Component 1 calculates field dynamics one spatial order at a time, whereas the second implements a weighted least squares solution based on probe distance. Component 3 then calculates phase residuals and removes the residual phase for distant probes before recalculation. Two volunteers and a phantom were scanned on a 7T MRI using diffusion-weighted sequences, and field monitoring was performed. Image reconstructions were informed with field dynamics calculated conventionally, and with the correction algorithm, after which in vivo images were compared qualitatively and phantom image error was quantitatively assessed.

RESULTS

The algorithm was able to correct corrupted field dynamics, resulting in image-quality improvements. Significant artifact reduction was observed when correcting higher-order fits. Stepwise fitting provided the most correction benefit, which was marginally improved when adding the other correction strategies.

CONCLUSION

The proposed algorithm can mitigate effects of phase errors in field monitoring, providing improved characterization of field dynamics.

Relationship between imaging parameters and distortion in magnetic resonance images for gamma knife treatment planning

In magnetic resonance imaging (MRI), it is necessary to reduce image distortion as much as possible because it suppresses the increase in the planning target volume. This study investigated the relationship between imaging parameters and image distortion when using G-frames. The images were obtained using a 1.5-T MRI system with a 09-101 Pro-MRI phantom. Image distortion was measured by changing the RF pulse mode, gradient mode, asymmetric echo, and bandwidth (BW). The image distortion was increased in the high RF mode than in the Normal mode. The image distortion increased in the following order: Whisper ≦ Normal < Fast in the different gradient modes. The image distortion increased in the following order: Without ≦ Weak < Strong in the different asymmetric echo modes. The image distortion increased in the following order: 300 Hz/pixel > 670 Hz/pixel ≧ REF (150 Hz/pixel) in the different Bw. The relationship between parameters and image distortion was clarified in this study when G-frames were used for gamma knife therapy. There is had relationship between the parameters causing variation in the gradient magnetic field and image distortion. Therefore, these parameters should be adjusted to minimize distortion.

Efficient approximate signal reconstruction for correction of gradient nonlinearities in diffusion-weighted imaging

In diffusion weighted MRI (DW-MRI), hardware nonlinearities lead to spatial variations in the orientation and magnitude of diffusion weighting. While the correction of these spatial distortions has been well established for analyses of DW-MRI, the existing voxel-wise empirical correction for gradient nonlinearities requires reimplementation of existing models, as the resultant gradients vary by voxel. Herein, we propose a two-step signal approximation after voxel-wise correction of gradient nonlinearity effects in DW-MRI. The proposed technique (1) scales the diffusion signal and (2) resamples the gradient orientations. This results in uniform gradients across the corrected image and provides the key advantage of seamless integration into current diffusion workflows. We investigated the validity of our technique by fitting a multi-compartment neurite orientation dispersion and density imaging (NODDI) model to the empirical correction and proposed approximation in five subjects from the MASiVar pediatric dataset. We evaluated intra-cellular volume fraction (iVF), CSF volume fraction (cVF), and orientation dispersion index (ODI) from NODDI. The Cohen's d of iVF, cVF and ODI between the techniques was <0.2 indicating the proposed technique does not exhibit significant differences from the voxel-wise correction technique. Our two-step signal approximation is an efficient representation of the voxel-wise gradient table correction. Using this approximation, correction of gradient nonlinearities can be easily incorporated into existing diffusion preprocessing pipelines and is implemented in "PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images".

Distortion-corrected image reconstruction with deep learning on an MRI-Linac

PURPOSE

MRI is increasingly utilized for image-guided radiotherapy due to its outstanding soft-tissue contrast and lack of ionizing radiation. However, geometric distortions caused by gradient nonlinearities (GNLs) limit anatomical accuracy, potentially compromising the quality of tumor treatments. In addition, slow MR acquisition and reconstruction limit the potential for effective image guidance. Here, we demonstrate a deep learning-based method that rapidly reconstructs distortion-corrected images from raw k-space data for MR-guided radiotherapy applications.

METHODS

We leverage recent advances in interpretable unrolling networks to develop a Distortion-Corrected Reconstruction Network (DCReconNet) that applies convolutional neural networks (CNNs) to learn effective regularizations and nonuniform fast Fourier transforms for GNL-encoding. DCReconNet was trained on a public MR brain dataset from 11 healthy volunteers for fully sampled and accelerated techniques, including parallel imaging (PI) and compressed sensing (CS). The performance of DCReconNet was tested on phantom, brain, pelvis, and lung images acquired on a 1.0T MRI-Linac. The DCReconNet, CS-, PI-and UNet-based reconstructed image quality was measured by structural similarity (SSIM) and RMS error (RMSE) for numerical comparisons. The computation time and residual distortion for each method were also reported.

RESULTS

Imaging results demonstrated that DCReconNet better preserves image structures compared to CS- and PI-based reconstruction methods. DCReconNet resulted in the highest SSIM (0.95 median value) and lowest RMSE (<0.04) on simulated brain images with four times acceleration. DCReconNet is over 10-times faster than iterative, regularized reconstruction methods.

CONCLUSIONS

DCReconNet provides fast and geometrically accurate image reconstruction and has the potential for MRI-guided radiotherapy applications.

Quality assurance and temporal stability of a 1.5 T MRI scanner for MR-guided Photon and Particle Therapy

PURPOSE

To describe performance measurements, adaptations and time stability over 20 months of a diagnostic MR scanner for integration into MR-guided photon and particle radiotherapy.

MATERIAL AND METHODS

For realization of MR-guided photon and particle therapy (MRgRT/MRgPT), a 1.5 T MR scanner was installed at the Heidelberg Ion Beam Therapy Center. To integrate MRI into the treatment process, a flat tabletop and dedicated coil holders for flex coils were used, which prevent deformation of the patient external contour and allow for the use of immobilization tools for reproducible positioning. The signal-to-noise ratio (SNR) was compared for the diagnostic and therapy-specific setup using the flat couch top and flexible coils for the a) head & neck and b) abdominal region as well as for different bandwidths and clinical pulse sequences. Additionally, a quality assurance (QA) protocol with monthly measurements of the ACR phantom and measurement of geometric distortions for a large field-of-view (FOV) was implemented to assess the imaging quality parameters of the device over the course of 20 months.

RESULTS

The SNR measurements showed a decreased SNR for the RT-specific as compared to the diagnostic setup of (a) 26% to 34% and (b) 11% to 33%. No significant bandwidth dependency for this ratio was found. The longitudinal assessment of the image quality parameters with the ACR and distortion phantom confirmed the long-term stability of the MRI device.

CONCLUSION

A diagnostic MRI was commissioned for use in MR-guided particle therapy. Using a radiotherapy specific setup, a high geometric accuracy and signal homogeneity was obtained after some adaptions and the measured parameters were shown to be stable over a period of 20 months.

The Value of Various Post-Processing Modalities of Diffusion Weighted Imaging in the Detection of Multiple Sclerosis

Diffusion tensor imaging (DTI) showed its adequacy in evaluating the normal-appearing white matter (NAWM) and lesions in the brain that are difficult to evaluate with routine clinical magnetic resonance imaging (MRI) in multiple sclerosis (MS). Recently, MRI systems have been developed with regard to software and hardware, leading to different proposed diffusion analysis methods such as diffusion tensor imaging, q-space imaging, diffusional kurtosis imaging, neurite orientation dispersion and density imaging, and axonal diameter measurement. These methods have the ability to better detect in vivo microstructural changes in the brain than DTI. These different analysis modalities could provide supplementary inputs for MS disease characterization and help in monitoring the disease’s progression as well as treatment efficacy. This paper reviews some of the recent diffusion MRI methods used for the assessment of MS in vivo.

Eddy current-induced artifacts correction in high gradient strength diffusion MRI with dynamic field monitoring: demonstration in ex vivo human brain imaging

Purpose

To demonstrate the advantages of spatiotemporal magnetic field monitoring to correct eddy current-induced artifacts (ghosting and geometric distortions) in high gradient strength diffusion MRI (dMRI).

Methods

A dynamic field camera with 16 NMR field probes was used to characterize eddy current fields induced from diffusion gradients for different gradients strengths (up to 300 mT/m), diffusion directions, and shots in a 3D multi-shot EPI sequence on a 3T Connectom scanner. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-resolution whole brain ex vivo dMRI. A 3D multi-shot image reconstruction framework was informed with the actual nonlinear phase evolution measured with the dynamic field camera, thereby accounting for high-order eddy currents fields on top of the image encoding gradients in the image formation model.

Results

Eddy current fields from diffusion gradients at high gradient strength in a 3T Connectom scanner are highly nonlinear in space and time, inducing high-order spatial phase modulations between odd/even echoes and shots that are not static during the readout. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting approaches such as navigator- and structured low-rank-based methods or MUSE, followed by image-based distortion correction with eddy. Improved dMRI analysis is demonstrated with diffusion tensor imaging and high-angular resolution diffusion imaging.

Conclusion

Strong eddy current artifacts characteristic of high gradient strength dMRI can be well corrected with dynamic field monitoring-based image reconstruction, unlike the two-step approach consisting of ghosting correction followed by geometric distortion reduction with eddy.

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 (P0 ), the edge (P 1) and just outside (P 2) 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 1 and 3.1 mm at P 2. When distortion correction was applied, the error was reduced to 1.0 mm at P 1 and 1.1 mm at P 2. The corrected error was similar to an error of 0.9 mm at P0 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.

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.

Cancellation of streak artifacts in radial abdominal imaging using interference null space projection

PURPOSE

In radial abdominal imaging, it has been commonly observed that signal from the arms cause streaks due to system imperfections. We previously introduced a streak removal technique (B-STAR), which is inherently spatially variant and limited to work in image space. In this work, we propose a spatially invariant streak cancellation technique (CACTUS), which can be applied in either image space or k-space and is compatible with iterative reconstructions.

THEORY AND METHODS

Streak sources are typically spatially localized and can be represented using a low-dimensional subspace. CACTUS identifies the streak subspace by leveraging the spatial redundancy of receiver coils and projects the data onto the streak null space to eliminate the streaks. When applied in k-space, CACTUS can be combined with iterative reconstructions. CACTUS was tested in phantoms and in vivo abdominal imaging using a radial turbo spin-echo pulse sequence.

RESULTS

In phantoms, CACTUS improved T2 estimation in comparison to previous de-streaking methods. In vivo experiments showed that CACTUS reduced streaks and yielded T2 estimation, in regions affected by streaks, closer to a streak-free reference. Evaluation using a clinical abdominal dataset (n = 20) showed that CACTUS is comparable to B-STAR and yields significantly better signal preservation and streak cancellation than coil removal and suppression methods.

CONCLUSION

CACTUS provides superior signal preservation and streak reduction performance compared to coil removal and suppression methods. As a clear advantage over B-STAR, CACTUS can be integrated with iterative reconstruction methods. In abdominal T2 mapping, CACTUS improves the accuracy of parameter estimation in areas affected by streaks.

Evaluation of the ability of the Brainlab Elements Cranial Distortion Correction algorithm to correct clinically relevant MRI distortions for cranial SRT

PURPOSE

Cranial stereotactic radiotherapy (SRT) requires highly accurate lesion delineation. However, MRI can have significant inherent geometric distortions. We investigated how well the Elements Cranial Distortion Correction algorithm of Brainlab (Munich, Germany) corrects the distortions in MR image-sets of a phantom and patients.

METHODS

A non-distorted reference computed tomography image-set of a CIRS Model 603-GS (CIRS, Norfolk, VA, USA) phantom was acquired. Three-dimensional T1-weighted images were acquired with five MRI scanners and reconstructed with vendor-derived distortion correction. Some were reconstructed without correction to generate heavily distorted image-sets. All MR image-sets were corrected with the Brainlab algorithm relative to the computed tomography acquisition. CIRS Distortion Check software measured the distortion in each image-set. For all uncorrected and corrected image-sets, the control points that exceeded the 0.5-mm clinically relevant distortion threshold and the distortion maximum, mean, and standard deviation were recorded. Empirical cumulative distribution functions (eCDF) were plotted. Intraclass correlation coefficient (ICC) was calculated. The algorithm was evaluated with 10 brain metastases using Dice similarity coefficients (DSC).

RESULTS

The algorithm significantly reduced mean and standard deviation distortion in all image-sets. It reduced the maximum distortion in the heavily distorted image-sets from 2.072 to 1.059 mm and the control points with > 0.5-mm distortion fell from 50.2% to 4.0%. Before and especially after correction, the eCDFs of the four repeats were visually similar. ICC was 0.812 (excellent-good agreement). The algorithm increased the DSCs for all patients and image-sets.

CONCLUSION

The Brainlab algorithm significantly and reproducibly ameliorated MRI distortion, even with heavily distorted images. Thus, it increases the accuracy of cranial SRT lesion delineation. After further testing, this tool may be suitable for SRT of small lesions.

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.

MaxGIRF: Image reconstruction incorporating concomitant field and gradient impulse response function effects

Purpose

To develop and evaluate an improved strategy for compensating concomitant field effects in non‐Cartesian MRI at the time of image reconstruction.

Theory

We present a higher‐order reconstruction method, denoted as MaxGIRF, for non‐Cartesian imaging that simultaneously corrects off‐resonance, concomitant fields, and trajectory errors without requiring specialized hardware. Gradient impulse response functions are used to predict actual gradient waveforms, which are in turn used to estimate the spatiotemporally varying concomitant fields based on analytic expressions. The result, in combination with a reference field map, is an encoding matrix that incorporates a correction for all three effects.

Methods

The MaxGIRF reconstruction is applied to noiseless phantom simulations, spiral gradient‐echo imaging of an International Society for Magnetic Resonance in Medicine/National Institute of Standards and Technology phantom, and axial and sagittal multislice spiral spin‐echo imaging of a healthy volunteer at 0.55 T. The MaxGIRF reconstruction was compared against previously established concomitant field‐compensation and image‐correction methods. Reconstructed images are evaluated qualitatively and quantitatively using normalized RMS error. Finally, a low‐rank approximation of MaxGIRF is used to reduce computational burden. The accuracy of the low‐rank approximation is studied as a function of minimum rank.

Results

The MaxGIRF reconstruction successfully mitigated blurring artifacts both in phantoms and in vivo and was effective in regions where concomitant fields counteract static off‐resonance, superior to the comparator method. A minimum rank of 8 and 30 for axial and sagittal scans, respectively, gave less than 2% error compared with the full‐rank reconstruction.

Conclusions

The MaxGIRF reconstruction simultaneously corrects off‐resonance, trajectory errors, and concomitant field effects. The impact of this method is greatest when imaging with longer readouts and/or at lower field strength.

Click here for author‐reader discussions

Combining navigator and optical prospective motion correction for high-quality 500 μm resolution quantitative multi-parameter mapping at 7T

PURPOSE

High-resolution quantitative multi-parameter mapping shows promise for non-invasively characterizing human brain microstructure but is limited by physiological artifacts. We implemented corrections for rigid head movement and respiration-related B0-fluctuations and evaluated them in healthy volunteers and dementia patients.

METHODS

Camera-based optical prospective motion correction (PMC) and FID navigator correction were implemented in a gradient and RF-spoiled multi-echo 3D gradient echo sequence for mapping proton density (PD), longitudinal relaxation rate (R1) and effective transverse relaxation rate (R2*). We studied their effectiveness separately and in concert in young volunteers and then evaluated the navigator correction (NAVcor) with PMC in a group of elderly volunteers and dementia patients. We used spatial homogeneity within white matter (WM) and gray matter (GM) and scan-rescan measures as quality metrics.

RESULTS

NAVcor and PMC reduced artifacts and improved the homogeneity and reproducibility of parameter maps. In elderly participants, NAVcor improved scan-rescan reproducibility of parameter maps (coefficient of variation decreased by 14.7% and 11.9% within WM and GM respectively). Spurious inhomogeneities within WM were reduced more in the elderly than in the young cohort (by 9% vs. 2%). PMC increased regional GM/WM contrast and was especially important in the elderly cohort, which moved twice as much as the young cohort. We did not find a significant interaction between the two corrections.

CONCLUSION

Navigator correction and PMC significantly improved the quality of PD, R1, and R2* maps, particularly in less compliant elderly volunteers and dementia patients.

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Long-Term Stability of Gradient Characteristics Warrants Model-Based Correction of Diffusion Weighting Bias

The study aims to test the long-term stability of gradient characteristics for model-based correction of diffusion weighting (DW) bias in an apparent diffusion coefficient (ADC) for multisite imaging trials. Single spin echo (SSE) DWI of a long-tube ice-water phantom was acquired quarterly on six MR scanners over two years for individual diffusion gradient channels, along with B0 mapping, as a function of right-left (RL) and superior-inferior (SI) offsets from the isocenter. Additional double spin-echo (DSE) DWI was performed on two systems. The offset dependences of derived ADC were fit to 4th-order polynomials. Chronic shim gradients were measured from spatial derivatives of B0 maps along the tube direction. Gradient nonlinearity (GNL) was modeled using vendor-provided gradient field descriptions. Deviations were quantified by root-mean-square differences (RMSD), normalized to reference ice-water ADC, between the model and reference (RMSDREF), measurement and model (RMSDEXP), and temporal measurement variations (RMSDTMP). Average RMSDREF was 4.9 ± 3.2 (%RL) and –14.8 ± 3.8 (%SI), and threefold larger than RMSDEXP. RMSDTMP was close to measurement errors (~3%). GNL-induced bias across gradient systems varied up to 20%, while deviation from the model accounted at most for 6.5%, and temporal variation for less than 3% of ADC reproducibility error. Higher SSE RMSDEXP = 7.5–11% was reduced to 2.5–4.8% by DSE, consistent with the eddy current origin. Measured chronic shim gradients below 0.1 mT/m had a minor contribution to ADC bias. The demonstrated long-term stability of spatial ADC profiles and consistency with system GNL models justifies retrospective and prospective DW bias correction based on system gradient design models. Residual errors due to eddy currents and shim gradients should be corrected independent of GNL.

B1 -gradient-based MRI using frequency-modulated Rabi-encoded echoes

PURPOSE

Reduce expense and increase accessibility of MRI by eliminating pulsed field (B0 ) gradient hardware.

METHODS

A radiofrequency imaging method is described that enables spatial encoding without B0 gradients. This method, herein referred to as frequency-modulated Rabi-encoded echoes (FREE), utilizes adiabatic full passage pulses and a gradient in the RF field (B1 ) to produce spatially dependent phase modulation, equivalent to conventional phase encoding. In this work, Cartesian phase encoding was accomplished using FREE in a multi-shot double spin-echo sequence. Theoretical analysis and computer simulations investigated the influence of resonance offset and B1 -gradient steepness and magnitude on reconstruction quality, which limit other radiofrequency imaging methodologies. Experimentally, FREE was compared to conventional phase-encoded MRI on human visual cortex using a simple surface transceiver coil.

RESULTS

Image distortions occurred in FREE when using nonlinear B1 fields where the phase dependence becomes nonlinear, but with minimal change in signal intensity. Resonance offset effects were minimal for Larmor frequencies within the adiabatic full-passage pulse bandwidth.

CONCLUSION

For the first time, FREE enabled slice-selective 2D imaging of the human brain without a B0 gradient in the y-direction. FREE achieved high resolution in regions where the B1 gradient was steepest, whereas images were distorted in regions where nonlinearity in the B1 gradient was significant. Given that FREE experiences no significant signal loss due to B1 nonlinearities and resonance offset, image distortions shown in this work might be corrected in the future based on B1 and B0 maps.

Gradient nonlinearity correction in liver DWI using motion-compensated diffusion encoding waveforms

OBJECTIVE : To experimentally characterize the effectiveness of a gradient nonlinearity correction method in removing ADC bias for different motion-compensated diffusion encoding waveforms.

METHODS

The diffusion encoding waveforms used were the standard monopolar Stejskal-Tanner pulsed gradient spin echo (pgse) waveform, the symmetric bipolar velocity-compensated waveform (sym-vc), the asymmetric bipolar velocity-compensated waveform (asym-vc) and the asymmetric bipolar partial velocity-compensated waveform (asym-pvc). The effectiveness of the gradient nonlinearity correction method using the spherical harmonic expansion of the gradient coil field was tested with the aforementioned waveforms in a phantom and in four healthy subjects.

RESULTS

The gradient nonlinearity correction method reduced the ADC bias in the phantom experiments for all used waveforms. The range of the ADC values over a distance of ± 67.2 mm from isocenter reduced from 1.29 × 10^-4 to 0.32 × 10^-4 mm²/s for pgse, 1.04 × 10^-4 to 0.22 × 10^-4 mm²/s for sym-vc, 1.22 × 10^-4 to 0.24 × 10^-4 mm²/s for asym-vc and 1.07 × 10^-4 to 0.11 × 10^-4 mm²/s for asym-pvc. The in vivo results showed that ADC overestimation due to motion or bright vessels can be increased even further by the gradient nonlinearity correction.

CONCLUSION

The investigated gradient nonlinearity correction method can be used effectively with various motion-compensated diffusion encoding waveforms. In coronal liver DWI, ADC errors caused by motion and residual vessel signal can be increased even further by the gradient nonlinearity correction.

Technical note: Unexpected external markers artifact in 3D k-space based parallel imaging turbo spin-echo magnetic resonance imaging

PURPOSE

MRI for radiotherapy planning requires spatial referencing using immobilization devices and markers. Clinical images of a difficult-to-interpret artifact are presented, resembling a metastasis, which occurs when combining CAIPIRINHA k-space-based parallel imaging (PI), 3D distortion correction, and external markers.

METHODS

A 3D variable flip angle Turbo Spin Echo sequence was used on a 1.5 T and 3 T MRI using flexible and head and neck coils. Two types of markers were tested: Liquimark LM1 and Spee-D-Mark. A silicone oil phantom was used that represents low signal intensity, such as gray matter. 3D Fourier transforms were also used to show the issue's origin.

RESULTS

The markers can appear in an unexpected region of a patient, not in the same original or reconstructed slice nor in a rectilinear direction in a slice, especially when using CAIPIRINHA acceleration with 3D distortion correction. The probability of occurrence was respectively 13% and 80% for distances of <=2 mm and >2 mm between marker and patient, for example when using thermoplastic masks. Clinical cases are shown where this semi-randomly occurring artifact appears post contrast only, and thus can be interpreted as metastases. The artifact did not appear when using compressed sensing acceleration.

CONCLUSION

Markers used for radiotherapy MRI application can introduce additional artifacts that can be interpreted as metastases. However, other high signal intensity structures on the surface of a patient, such as the ear, can lead to an equivalent error.

Recent Advances and Future Directions in Brain MR Imaging Studies in Schizophrenia: Toward Elucidating Brain Pathology and Developing Clinical Tools

Schizophrenia is a common severe psychiatric disorder that affects approximately 1% of general population through the life course. Historically, in Kraepelin's time, schizophrenia was a disease unit conceptualized as dementia praecox; however, since then, the disease concept has changed. Recent MRI studies had shown that the neuropathology of the brain in this disorder was characterized by mild progression before and after the onset of the disease, and that the brain alterations were relatively smaller than assumed. Although genetic factors contribute to the brain alterations in schizophrenia, which are thought to be trait differences, other changes include factors that are common in psychiatric diseases. Furthermore, it has been shown that the brain differences specific to schizophrenia were relatively small compared to other changes, such as those caused by brain development, aging, and gender. In addition, compared to the disease and participant factors, machine and imaging protocol differences could affect MRI signals, which should be addressed in multi-site studies. Recent advances in MRI modalities, such as multi-shell diffusion-weighted imaging, magnetic resonance spectroscopy, and multimodal brain imaging analysis, may be candidates to sharpen the characterization of schizophrenia-specific factors and provide new insights. The Brain/MINDS Beyond Human Brain MRI (BMB-HBM) project has been launched considering the differences and noises irrespective of the disease pathologies and includes the future perspectives of MRI studies for various psychiatric and neurological disorders. The sites use restricted MRI machines and harmonized multi-modal protocols, standardized image preprocessing, and traveling subject harmonization. Data sharing to the public will be planned in FY 2024. In the future, we believe that combining a high-quality human MRI dataset with genetic data, randomized controlled trials, and MRI for non-human primates and animal models will enable us to understand schizophrenia, elucidate its neural bases and therapeutic targets, and provide tools for clinical application at bedside.

Comparison of traveling-subject and ComBat harmonization methods for assessing structural brain characteristics

Multisite magnetic resonance imaging (MRI) is increasingly used in clinical research and development. Measurement biases—caused by site differences in scanner/image‐acquisition protocols—negatively influence the reliability and reproducibility of image‐analysis methods. Harmonization can reduce bias and improve the reproducibility of multisite datasets. Herein, a traveling‐subject (TS) dataset including 56 T1‐weighted MRI scans of 20 healthy participants in three different MRI procedures—20, 19, and 17 subjects in Procedures 1, 2, and 3, respectively—was considered to compare the reproducibility of TS‐GLM, ComBat, and TS‐ComBat harmonization methods. The minimum participant count required for harmonization was determined, and the Cohen's d between different MRI procedures was evaluated as a measurement‐bias indicator. The measurement‐bias reduction realized with different methods was evaluated by comparing test–retest scans for 20 healthy participants. Moreover, the minimum subject count for harmonization was determined by comparing test–retest datasets. The results revealed that TS‐GLM and TS‐ComBat reduced measurement bias by up to 85 and 81.3%, respectively. Meanwhile, ComBat showed a reduction of only 59.0%. At least 6 TSs were required to harmonize data obtained from different MRI scanners, complying with the imaging protocol predetermined for multisite investigations and operated with similar scan parameters. The results indicate that TS‐based harmonization outperforms ComBat for measurement‐bias reduction and is optimal for MRI data in well‐prepared multisite investigations. One drawback is the small sample size used, potentially limiting the applicability of ComBat. Investigation on the number of subjects needed for a large‐scale study is an interesting future problem.

RAVE-T2/T1 - Feasibility of a new hybrid MR-sequence for free-breathing abdominal MRI in children and adolescents

BACKGROUND

The new radial volumetric encoding RAVE-T₂/T₁ hybrid sequence is a modern three-dimensional sequence with multiparametric approach, which includes T₂- and T₁-weighted contrasts obtained in identical slice position during one measurement. However, the RAVE-T₂/T₁ hybrid sequence is not yet being used in clinical routine.

PURPOSE

The aim of this study was to evaluate the RAVE-T₂/T₁ hybrid sequence in a pediatric population with a clinical indication for an abdominal MRI examination to demonstrate that the hybrid imaging may be less challenging to perform on children.

MATERIALS AND METHODS

Our retrospective observational study included pediatric patients of all age groups and required for an abdominal MRI examination. Non-contrast standard axial T₁ DIXON and non-contrast RAVE-T₂/T₁ hybrid sequence were obtained at 3 T. MRI studies were analyzed independently by two pediatric radiologists using a 5-point Likert-type scale in five different categories. T₁- and T₂-weighted sequences were each compared with the RAVE-T₂/T₁-sequence using a Wilcoxon signed-rank test.

RESULTS

The analysis included 15 children (mean age, 11 years and 4 months, 7 girls and 8 boys). The Cohens Kappa of interrater agreement measured 0.62. The T₂ weighted part of the RAVE-T₂/T₁ sequence was significantly better than the standard T₂ HASTE sequence in four of five image quality categories: overall image quality (2.2 ± 0.7 vs 1.8 ± 0,7, p = 0.03), respiratory motion artefacts (3.8 ± 0.4 vs 2.0 ± 0.7, p <= 0.01), portal vein clarity (3.3 ± 0.8 vs 2.2 ± 0.7, p <= 0.01), hepatic margin sharpness (2.4 ± 1,0 vs 1.8 ± 0.7, p <= 0.01). The T₁ weighted part of the RAVE-T₂/T₁ sequence was significantly better than the standard T₁ DIXON weighted sequence in three of five image quality categories: respiratory motion artefacts (4.0 ± 0.2 vs 3.6 ± 0.8, p = 0.01), portal vein clarity (2.7 ± 0.9 vs 2.1 ± 0.7, p <= 0.01), hepatic margin sharpness (3.2 ± 0.7 vs 2.6 ± 0.9, p <= 0.01).

CONCLUSIONS

The RAVE-T₂/T₁ hybrid sequence is feasible and equal compared to standard T₁- and T₂-weighted sequences in the assessment of abdominal organs in a pediatric population. Due to non-inferiority to the current standard sequences for abdominal imaging, the RAVE-T₂/T₁ hybrid sequence is a good alternative for children who cannot be examined in breath-hold technique.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

Empirical validation of gradient field models for an accurate ADC measured on clinical 3T MR systems in body oncologic applications

PURPOSE

To empirically corroborate vendor-provided gradient nonlinearity (GNL) characteristics and demonstrate efficient GNL bias correction for human brain apparent diffusion coefficient (ADC) across 3T MR systems and spatial locations.

METHODS

Spatial distortion vector fields (DVF) were mapped in 3D using a surface fiducial array phantom for individual gradient channels on three 3T MR platforms from different vendors. Measured DVF were converted into empirical 3D GNL tensors and compared with their theoretical counterparts derived from vendor-provided spherical harmonic (SPH) coefficients. To illustrate spatial impact of GNL on ADC, diffusion weighted imaging using three orthogonal gradient directions was performed on a volunteer brain positioned at isocenter (as a reference) and offset superiorly by 10-17 cm (>10% predicted GNL bias). The SPH tensor-based GNL correction was applied to individual DWI gradient directions, and derived ADC was compared with low-bias reference for human brain white matter (WM) ROIs.

RESULTS

Empiric and predicted GNL errors were comparable for all three studied 3T MR systems, with <1.0% differences in the median and width of spatial histograms for individual GNL tensor elements. Median (±width) of ADC (10-3mm2/s) histograms measured at isocenter in WM reference ROIs from three MR systems were: 0.73 ± 0.11, 0.71 ± 0.14, 0.74 ± 0.17, and at off-isocenters (before versus after GNL correction) were respectively 0.63 ± 0.14 versus 0.72 ± 0.11, 0.53 ± 0.16 versus 0.74 ± 0.18, and 0.65 ± 0.16 versus 0.76 ± 0.18.

CONCLUSION

The phantom-based spatial distortion measurements validated vendor-provided gradient fields, and accurate WM ADC was recovered regardless of spatial locations and clinical MR platforms using system-specific tensor-based GNL correction for routine DWI.

Simultaneous multi-slice accelerated 4D-MRI for radiotherapy guidance

4D-MRI is becoming increasingly important for daily guidance of thoracic and abdominal radiotherapy. This study exploits the simultaneous multi-slice (SMS) technique to accelerate the acquisition of a balanced turbo field echo (bTFE) and a turbo spin echo (TSE) coronal 4D-MRI sequence performed on 1.5 T MRI scanners. SMS single-shot bTFE and TSE sequences were developed to acquire a stack of 52 coronal 2D images over 30 dynamics. Simultaneously excited slices were separated by half the field of view. Slices intersecting with the liver-lung interface were used as navigator slices. For each navigator slice location, an end-exhale dynamic was automatically identified, and used to derive the self-sorting signal by rigidly registering the remaining dynamics. Navigator slices were sorted into 10 amplitude bins, and the temporal relationship of simultaneously excited slices was used to generate sorted 4D-MRIs for 12 healthy volunteers. The self-sorting signal was validated using anin vivopeak-to-peak motion analysis. The smoothness of the liver-lung interface was quantified by comparing to sagittal cine images acquired directly after the SMS-4D-MRI sequence. To ensure compatibility with the MR-linac radiotherapy workflow, the 4D-MRIs were transformed into 3D mid-position (MidP) images using deformable image registration. Consistency of the deformable vector fields was quantified in terms of the distance discordance metric (DDM) in the body. The SMS-4D-TSE sequence was additionally acquired for 3 lung cancer patients to investigate tumor visibility. SMS-4D-MRI acquisition and processing took approximately 7 min. 4D-MRI reconstruction was possible for 26 out of 27 acquired datasets. Missing data in the sorted 4D-MRIs varied from 4%-26% for the volunteers and varied from 8%-24% for the patients. Peak-to-peak (SD) amplitudes analysis agreed within 1.8 (1.1) mm and 0.9 (0.4) mm between the sorted 4D-MRIs and the self-sorting signals of the volunteers and patients, respectively. Liver-lung interface smoothness was found to be in the range of 0.6-3.1 mm for volunteers. The percentage of DDM values smaller than 2 mm was in the range of 85%-89% and 86%-92% for the volunteers and patients, respectively. Lung tumors were clearly visibility in the SMS-4D-TSE images and MidP images. Two fast SMS-accelerated 4D-MRI sequences were developed resulting in T₂/T₁or T₂weighted contrast. The SMS-4D-MRIs and derived 3D MidP-MRIs yielded anatomically plausible images and good tumor visibility. SMS-4D-MRI is therefore a strong candidate to be used for treatment simulation and daily guidance of thoracic and abdominal MR-guided radiotherapy.

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.

Aliasing artifact reduction in spiral real-time MRI

PURPOSE

To mitigate a common artifact in spiral real-time MRI, caused by aliasing of signal outside the desired FOV. This artifact frequently occurs in midsagittal speech real-time MRI.

METHODS

Simulations were performed to determine the likely origin of the artifact. Two methods to mitigate the artifact are proposed. The first approach, denoted as "large FOV" (LF), keeps an FOV that is large enough to include the artifact signal source during reconstruction. The second approach, denoted as "estimation-subtraction" (ES), estimates the artifact signal source before subtracting a synthetic signal representing that source in multicoil k-space raw data. Twenty-five midsagittal speech-production real-time MRI data sets were used to evaluate both of the proposed methods. Reconstructions without and with corrections were evaluated by two expert readers using a 5-level Likert scale assessing artifact severity. Reconstruction time was also compared.

RESULTS

The origin of the artifact was found to be a combination of gradient nonlinearity and imperfect anti-aliasing in spiral sampling. The LF and ES methods were both able to substantially reduce the artifact, with an averaged qualitative score improvement of 1.25 and 1.35 Likert levels for LF correction and ES correction, respectively. Average reconstruction time without correction, with LF correction, and with ES correction were 160.69 ± 1.56, 526.43 ± 5.17, and 171.47 ± 1.71 ms/frame.

CONCLUSION

Both proposed methods were able to reduce the spiral aliasing artifacts, with the ES-reduction method being more effective and more time efficient.

Porcine lung phantom-based validation of estimated 4D-MRI using orthogonal cine imaging for low-field MR-Linacs

Real-time motion monitoring of lung tumors with low-field magnetic resonance imaging-guided linear accelerators (MR-Linacs) is currently limited to sagittal 2D cine magnetic resonance imaging (MRI). To provide input data for improved intrafractional and interfractional adaptive radiotherapy, the 4D anatomy has to be inferred from data with lower dimensionality. The purpose of this study was to experimentally validate a previously proposed propagation method that provides continuous time-resolved estimated 4D-MRI based on orthogonal cine MRI for a low-field MR-Linac. Ex vivo porcine lungs were injected with artificial nodules and mounted in a dedicated phantom that allows for the simulation of periodic and reproducible breathing motion. The phantom was scanned with a research version of a commercial 0.35 T MR-Linac. Respiratory-correlated 4D-MRI were reconstructed and served as ground truth images. Series of interleaved orthogonal slices in sagittal and coronal orientation, intersecting the injected targets, were acquired at 7.3 Hz. Estimated 4D-MRI at 3.65 Hz were created in post-processing using the propagation method and compared to the ground truth 4D-MRI. Eight datasets at different breathing frequencies and motion amplitudes were acquired for three porcine lungs. The overall median (95[Formula: see text] percentile) deviation between ground truth and estimated deformation vector fields was 2.3 mm (5.7 mm), corresponding to 0.7 (1.6) times the in-plane imaging resolution (3.5 × 3.5 mm²). Median (95[Formula: see text] percentile) estimated nodule position errors were 1.5 mm (3.8 mm) for nodules intersected by orthogonal slices and 2.1 mm (7.1 mm) for nodules located more than 2 cm away from either of the orthogonal slices. The estimation error depended on the breathing phase, the motion amplitude and the location of the estimated position with respect to the orthogonal slices. By using the propagation method, the 4D motion within the porcine lung phantom could be accurately and robustly estimated. The method could provide valuable information for treatment planning, real-time motion monitoring, treatment adaptation, and post-treatment evaluation of MR-guided radiotherapy treatments.

Magnetic modeling of actively shielded rotating MRI magnets in the presence of environmental steel

Rotating MRI systems could enable novel integrated medical devices such as MRI-Linacs, MRI-xray-angiography systems, and MRI-proton therapy systems. This work aimed to investigate the feasibility of rotating actively shielded superconducting MRI magnets in the presence of environmental steel-in particular, construction steel in the floor of the installation site. Two magnets were investigated: a 1.0 T split bore magnet, and a 1.5 T closed bore magnet. Each magnet was scaled to emulate field strengths of 0.5, 1.0, and 1.5 T. Finite Element Modeling was used to simulate these magnets in the presence of a 3 × 4 m steel plate located 1250 mm or 1400 mm below the isocenter. There are two possible rotation directions: around the longitudinal (z) axis or around the transverse (x) axis. Each model was solved for rotation angles between 0 and 360° in 30° intervals around each of these axes. For each simulation, a 300 mm DSV was extracted and decomposed into spherical harmonics. For the closed-bore magnet, total induced perturbation for the zero degree rotation angle was 223, 432, and 562 μT peak-to-peak (pk-pk) for the 0.5, 1.0, and 1.5 T models respectively (steel at 1250 mm). For the split-bore magnet, the same numbers were 1477, 16747, and 1766 μT. The substantially higher perturbation for the split-bore magnet can be traced to its larger fringe field. For rotation around the z-axis, total perturbation does not change as a function of angle but is exchanged between different harmonics. For rotation around the x-axis, total perturbation is different at each rotation angle. For the closed bore magnet, maximum perturbations occurred for a 90° rotation around the transverse axis. For the split-bore magnet, the opposite was observed, with the same 90° rotation yielding total perturbation lower than the conventional position. In all cases, at least 95% of the total perturbation was composed of 1st and 2nd order harmonics. The presence of environmental steel poses a major challenge to the realization of an actively shielded rotating superconducting MRI system, requiring some novel form of shimming. Possible shimming strategies are discussed at length.

Empirical field mapping for gradient nonlinearity correction of multi-site diffusion weighted MRI

BACKGROUND

Achieving inter-site / inter-scanner reproducibility of diffusion weighted magnetic resonance imaging (DW-MRI) metrics has been challenging given differences in acquisition protocols, analysis models, and hardware factors.

PURPOSE

Magnetic field gradients impart scanner-dependent spatial variations in the applied diffusion weighting that can be corrected if the gradient nonlinearities are known. However, retrieving manufacturer nonlinearity specifications is not well supported and may introduce errors in interpretation of units or coordinate systems. We propose an empirical approach to mapping the gradient nonlinearities with sequences that are supported across the major scanner vendors.

STUDY TYPE

Prospective observational study.

SUBJECTS

A spherical isotropic diffusion phantom, and a single human control volunteer.

FIELD STRENGTH/SEQUENCE

3 T (two scanners). Stejskal-Tanner spin echo sequence with b-values of 1000, 2000 s/mm² with 12, 32, and 384 diffusion gradient directions per shell.

ASSESSMENT

We compare the proposed correction with the prior approach using manufacturer specifications against typical diffusion pre-processing pipelines (i.e., ignoring spatial gradient nonlinearities). In phantom data, we evaluate metrics against the ground truth. In human and phantom data, we evaluate reproducibility across scans, sessions, and hardware.

STATISTICAL TESTS

Wilcoxon rank-sum test between uncorrected and corrected data.

RESULTS

In phantom data, our correction method reduces variation in mean diffusivity across sessions over uncorrected data (p < 0.05). In human data, we show that this method can also reduce variation in mean diffusivity across scanners (p < 0.05).

CONCLUSION

Our method is relatively simple, fast, and can be applied retroactively. We advocate incorporating voxel-specific b-value and b-vector maps should be incorporated in DW-MRI harmonization preprocessing pipelines to improve quantitative accuracy of measured diffusion parameters.

On the use of low-dimensional temporal subspace constraints to reduce reconstruction time and improve image quality of accelerated 4D-MRI

BACKGROUND AND PURPOSE

The purpose of this work is to investigate the use of low-dimensional temporal subspace constraints for 4D-MRI reconstruction from accelerated data in the context of MR-guided online adaptive radiation therapy (MRgOART).

MATERIALS AND METHODS

Subspace basis functions are derived directly from the accelerated golden angle radial stack-of-stars 4D-MRI data. The reconstruction times, image quality, and motion estimates are investigated as a function of the number of subspace coefficients and compared with a conventional frame-by-frame reconstruction. These experiments were performed in five patients with four 4D-MRI scans per patient on a 1.5T MR-Linac.

RESULTS

If two or three subspace coefficients are used, the iterative reconstruction time is reduced by 32% and 18%, respectively, compared to conventional parallel imaging with compressed sensing reconstructions. No significant difference was found between motion estimates made with the subspace-constrained reconstructions (p > 0.08). Qualitative improvements in image quality included reduction in apparent noise and reductions in streaking artifacts from the radial k-space coverage.

CONCLUSION

Incorporating subspace constraints for accelerated 4D-MRI reconstruction reduces noise and residual undersampling artifacts in the images while reducing computation time, making it a strong candidate for use in clinical MRgOART workflows.

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.

Quantitative Evaluations of Geometrical Distortion Corrections in Cortical Surface-Based Analysis of High-Resolution Functional MRI Data at 7T

Background

Although 7T functional MRI (fMRI) provides better signal‐to‐noise ratio and higher spatial resolution than 3T fMRI, geometric distortions become more challenging because fMRI is more susceptible to distortions than structural MRI. Accurate alignment of 7T fMRI to structural MRI data is critical for precise cortical surface‐based analysis.

Purpose

To quantify the effectiveness of distortion corrections of 7T fMRI data.

Study Type

Prospective.

Subjects

Fifteen healthy individuals aged 19–26 years (mean: 21.9 years).

Field Strength/Sequence

Multiband gradient‐echo echo‐planar imaging sequence at 7T; 3D T₁/T₂‐weighted sequences (magnetization prepared rapid acquisition with gradient echo [MPRAGE] and sampling perfection with application optimized contrast using different flip angle evolution [SPACE]) at 3T.

Assessment

fMRI data at 7T were registered to cortical surfaces reconstructed from 3T structural data acquired in the same subjects. Distortions induced by B₀ inhomogeneity and gradient nonlinearity (B₀ and gradient distortions) were evaluated as cortical fallout (misregistration of noncortical areas) and displacement (misregistration along gray matter).

Statistical Tests

Repeated measures analyses of variance with post‐hoc t‐tests with Bonferroni correction.

Results

The accuracy of fully corrected fMRI images based on the intensity distribution was 89.2%. Without any corrections, 9.7% of vertices in the whole surfaces were fallout and the average displacement was 0.96 mm for the rest of the vertices. B₀ and gradient distortion corrections significantly reduced the fallout (to 2.1% and 8.7%) and displacement (to 0.29 mm and 0.86 mm). These corrections were effective even around regions with moderate distortions (the somatosensory and visual cortices for B₀ distortion, and the anterior frontal, inferior temporal, and posterior occipital cortices for gradient distortion).

Data Conclusion

B₀ distortion correction is crucial for surface‐based analysis of fine‐resolution fMRI at 7T. Gradient distortion correction should be considered when regions of interest include regions distant from the isocenter of scanners.

Evidence Level

1

Technical Efficacy Stage

1

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.

Low-field magnetic resonance imaging using multiplicative regularization

In this paper we present a magnetic resonance imaging (MRI) technique that is based on multiplicative regularization. Instead of adding a regularizing objective function to a data fidelity term, we multiply by such a regularizing function. By following this approach, no regularization parameter needs to be determined for each new data set that is acquired. Reconstructions are obtained by iteratively updating the images using short-term conjugate gradient-type update formulas and Polak-Ribière update directions. We show that the algorithm can be used as an image reconstruction algorithm and as a denoising algorithm. We illustrate the performance of the algorithm on two-dimensional simulated low-field MR data that is corrupted by noise and on three-dimensional measured data obtained from a low-field MR scanner. Our reconstruction results show that the algorithm effectively suppresses noise and produces accurate reconstructions even for low-field MR signals with a low signal-to-noise ratio.

Test-retest of automated segmentation with different motion correction strategies: A comparison of prospective versus retrospective methods

Test-retest of automated image segmentation algorithms (FSL FAST, FSL FIRST, and FREESURFER) are computed on magnetic resonance images from 12 unsedated children aged 9.4±2.6 years ([min,max] ​= ​[6.5 years, 13.8 years]) using different approaches to motion correction (prospective versus retrospective). The prospective technique, PROMO MPRAGE, dynamically estimates motion using specially acquired navigator images and adjusts the remaining acquisition accordingly, whereas the retrospective technique, MPnRAGE, uses a self-navigation property to retrospectively estimate and account for motion during image reconstruction. To increase the likelihood and range of motions, participants heads were not stabilized with padding during repeated scans. When motion was negligible both techniques had similar performance. When motion was not negligible, the automated image segmentation and anatomical labeling software tools showed the most consistent performance with the retrospectively corrected MPnRAGE technique (≥80% volume overlaps for 15 of 16 regions for FIRST and FREESURFER, with greater than 90% volume overlaps for 12 regions with FIRST and 11 regions with FREESURFER). Prospectively corrected MPRAGE with linear view-ordering also demonstrated lower performance than MPnRAGE without retrospective motion correction.

Magnetic resonance imaging for brain stereotactic radiotherapy : A review of requirements and pitfalls

Due to its superior soft tissue contrast, magnetic resonance imaging (MRI) is essential for many radiotherapy treatment indications. This is especially true for treatment planning in intracranial tumors, where MRI has a long-standing history for target delineation in clinical practice. Despite its routine use, care has to be taken when selecting and acquiring MRI studies for the purpose of radiotherapy treatment planning. Requirements on MRI are particularly demanding for intracranial stereotactic radiotherapy, where accurate imaging has a critical role in treatment success. However, MR images acquired for routine radiological assessment are frequently unsuitable for high-precision stereotactic radiotherapy as the requirements for imaging are significantly different for radiotherapy planning and diagnostic radiology. To assure that optimal imaging is used for treatment planning, the radiation oncologist needs proper knowledge of the most important requirements concerning the use of MRI in brain stereotactic radiotherapy. In the present review, we summarize and discuss the most relevant issues when using MR images for target volume delineation in intracranial stereotactic radiotherapy.

[Evaluation of Long-term Fluctuation of Geometric Distortion in MRI for Radiation Therapy Planning by Using an Automatic Analysis Tool]

High tissue contrast in magnetic resonance imaging (MRI) allows better radiotherapy planning. However, geometric distortion in MRI induces inaccuracies affecting such planning, making it necessary to evaluate the characteristics of such geometric distortion. Although many studies have considered geometric distortion, most of these involved measurements performed only a few times. In this study, we evaluated MRI device-specific geometric distortion over long term and measured its variation by using an automatic analysis tool. The result showed that geometric distortion increased with distance from the center along both lateral and longitudinal directions. Specifically, the average distortion rate and average diameter error over the full measurement period increased by up to 1.02% and 1.96 mm, respectively, when using T1 weighted Image (WI) 3D fast spoiled gradient echo (FSPGR) at R15. In the case of T2 WI 2D fast spin echo (FSE) at R15, the standard deviation of the distortion rate and diameter error increased up to 0.38%, 0.72 mm, respectively. We conclude that periodic quality assurance of geometric distortion should be performed in order to maintain geometric distortion within allowable values.