A field-monitoring-based approach for correcting eddy-current-induced artifacts of up to the 2nd spatial order in human-connectome-project-style multiband diffusion MRI experiment at 7T: A pilot study

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.

Maximizing SNR per unit time in diffusion MRI with multiband T-Hex spirals

PURPOSE

The characterization of tissue microstructure using diffusion MRI (dMRI) signals is rapidly evolving, with increasing sophistication of signal representations and microstructure models. However, this progress often requires signals to be acquired with very high b-values (e.g., b > 30 ms/μm2 ), along many directions, and using multiple b-values, leading to long scan times and extremely low SNR in dMRI images. The purpose of this work is to boost the SNR efficiency of dMRI by combining three particularly efficient spatial encoding techniques and utilizing a high-performance gradient system (Gmax  ≤ 300 mT/m) for efficient diffusion encoding.

METHODS

Spiral readouts, multiband imaging, and sampling on tilted hexagonal grids (T-Hex) are combined and implemented on a 3T MRI system with ultra-strong gradients. Image reconstruction is performed through an iterative cg-SENSE algorithm incorporating static off-resonance distributions and field dynamics as measured with an NMR field camera. Additionally, T-Hex multiband is combined with a more conventional EPI-readout and compared with state-of-the-art blipped-CAIPIRINHA sampling. The advantage of the proposed approach is furthermore investigated for clinically available gradient performance and diffusion kurtosis imaging.

RESULTS

High fidelity in vivo images with b-values up to 40 ms/μm2 are obtained. The approach provides superior SNR efficiency over other state-of-the-art multiband diffusion readout schemes.

CONCLUSION

The demonstrated gains hold promise for the widespread dissemination of advanced microstructural scans, especially in clinical populations.

Eddy current-induced artifact correction in high b-value ex vivo human brain diffusion MRI with dynamic field monitoring

PURPOSE

To investigate whether spatiotemporal magnetic field monitoring can correct pronounced eddy current-induced artifacts incurred by strong diffusion-sensitizing gradients up to 300 mT/m used in high b-value diffusion-weighted (DW) EPI.

METHODS

A dynamic field camera equipped with 16 1 H NMR field probes was first used to characterize field perturbations caused by residual eddy currents from diffusion gradients waveforms in a 3D multi-shot EPI sequence on a 3T Connectom scanner for different gradient strengths (up to 300 mT/m), diffusion directions, and shots. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-gradient strength, submillimeter resolution whole-brain ex vivo diffusion MRI. A 3D multi-shot image reconstruction framework was developed that incorporated the nonlinear phase evolution measured with the dynamic field camera.

RESULTS

Phase perturbations in the readout induced by residual eddy currents from strong diffusion gradients are highly nonlinear in space and time, vary among diffusion directions, and interfere significantly with the image encoding gradients, changing the k-space trajectory. During the readout, phase modulations between odd and even EPI echoes become non-static and diffusion encoding direction-dependent. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting reduction approaches such as navigator- and structured low-rank-based methods or MUSE followed by image-based distortion correction with the FSL tool "eddy."

CONCLUSION

Strong eddy current artifacts characteristic of high-gradient strength DW-EPI can be well corrected with dynamic field monitoring-based image reconstruction.

High-resolution diffusion-weighted imaging at 7 Tesla: single-shot readout trajectories and their impact on signal-to-noise ratio, spatial resolution and accuracy

Diffusion MRI (dMRI) is a valuable imaging technique to study the connectivity and microstructure of the brain in vivo. However, the resolution of dMRI is limited by the low signal-to-noise ratio (SNR) of this technique. Various multi-shot acquisition strategies have been developed to achieve sub-millimeter resolution, but they require long scan times which can be restricting for patient scans. Alternatively, the SNR of single-shot acquisitions can be increased by using a spiral readout trajectory to minimize the sequence echo time. Imaging at ultra-high fields (UHF) could further increase the SNR of single-shot dMRI; however, the shorter T2* of brain tissue and the greater field non-uniformities at UHFs will degrade image quality, causing image blurring, distortions, and signal loss. In this study, we investigated the trade-off between the SNR and resolution of different k-space trajectories, including echo planar imaging (EPI), partial Fourier EPI, and spiral trajectories, over a range of dMRI resolutions at 7T. The effective resolution, spatial specificity and sharpening effect were measured from the point spread function (PSF) of the simulated diffusion sequences for a nominal resolution range of 0.6-1.8 mm. In-vivo partial brain scans at a nominal resolution of 1.5 mm isotropic were acquired using the three readout trajectories to validate the simulation results. Field probes were used to measure dynamic magnetic fields offline up to the 3^rd order of spherical harmonics. Image reconstruction was performed using static ΔB₀ field maps and the measured trajectories to correct image distortions and artifacts, leaving T2* effects as the primary source of blurring. The effective resolution was examined in fractional anisotropy (FA) maps calculated from a multi-shell dataset with b-values of 300, 1000, and 2000 s/mm² in 5, 16, and 48 directions, respectively. In-vivo scans at nominal resolutions of 1, 1.2, and 1.5 mm were acquired and the SNR of the different trajectories calculated using the multiple replica method to investigate the SNR. Finally, in-vivo whole brain scans with an effective resolution of 1.5 mm isotropic were acquired to explore the SNR and efficiency of different trajectories at a matching effective resolution. FA and intra-cellular volume fraction (ICVF) maps calculated using neurite orientation dispersion and density imaging (NODDI) were used for the comparison. The simulations and in vivo imaging results showed that for matching nominal resolutions, EPI trajectories had the highest specificity and effective resolution with maximum image sharpening effect. However, spirals have a significantly higher SNR, in particular at higher resolutions and even when the effective image resolutions are matched. Overall, this work shows that the higher SNR of single-shot spiral trajectories at 7T allows us to achieve higher effective resolutions compared to EPI and PF-EPI to map the microstructure and connectivity of small brain structures.

Reconstruction for 7T high-resolution whole-brain diffusion MRI using two-stage N/2 ghost correction and L1-SPIRiT without single-band reference

PURPOSE

To combine a new two-stage N/2 ghost correction and an adapted L1-SPIRiT method for reconstruction of 7T highly accelerated whole-brain diffusion MRI (dMRI) using only autocalibration scans (ACS) without the need of additional single-band reference (SBref) scans.

METHODS

The proposed ghost correction consisted of a 3-line reference approach in stage 1 and the reference-free entropy method in stage 2. The adapted L1-SPIRiT method was formulated within the 3D k-space framework. Its efficacy was examined by acquiring two dMRI data sets at 1.05-mm isotropic resolutions with a total acceleration of 6 or 9 (i.e., 2-fold or 3-fold slice and 3-fold in-plane acceleration). Diffusion analysis was performed to derive DTI metrics and estimate fiber orientation distribution functions (fODFs). The results were compared with those of 3D k-space GRAPPA using only ACS, all in reference to 3D k-space GRAPPA using both ACS and SBref (serving as a reference).

RESULTS

The proposed ghost correction eliminated artifacts more robustly than conventional approaches. Our adapted L1-SPIRiT method outperformed 3D k-space GRAPPA when using only ACS, improving image quality to what was achievable with 3D k-space GRAPPA using both ACS and SBref scans. The improvement in image quality further resulted in an improvement in estimation performances for DTI and fODFs.

CONCLUSION

The combination of our new ghost correction and adapted L1-SPIRiT method can reliably reconstruct 7T highly accelerated whole-brain dMRI without the need of SBref scans, increasing acquisition efficiency and reducing motion sensitivity.

Fast measurement of the gradient system transfer function at 7 T

PURPOSE

In this work, a new method to determine the gradient system transfer function (GSTF) with high frequency resolution and high SNR is presented, using fast and simple phantom measurements. The GSTF is an effective instrument for hardware characterization and calibration, which can be used to correct for gradient distortions, or enhance gradient fidelity.

METHODS

The thin-slice approach for phantom-based measurements of the GSTF is expanded by adding excitations that are shifted after the application of the probing gradient, to capture long-lasting field fluctuations with high SNR. A physics-informed regularization procedure is implemented to derive high-quality transfer functions from a small number of measurements. The resulting GSTFs are evaluated by means of gradient time-course estimation and pre-emphasis of a trapezoidal test gradient on a 7T scanner.

RESULTS

The GSTFs determined with the proposed method capture sharp mechanical resonances with a high level of detail. The measured trapezoidal gradient progressions are authentically reproduced by the GSTF estimations on all three axes. The GSTF-based pre-emphasis considerably improves the gradient fidelity in the plateau phase of the test gradient and almost completely eliminates lingering field oscillations.

CONCLUSION

The presented approach allows fast and simple characterization of gradient field fluctuations caused by long-living eddy current and vibration effects, which become more pronounced at ultrahigh field strengths.

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.

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.

Integration of an RF coil and commercial field camera for ultrahigh-field MRI

PURPOSE

To develop an RF coil with an integrated commercial field camera for ultrahigh field (7T) neuroimaging. The RF coil would operate within a head-only gradient coil and be subject to the corresponding design constraints. The RF coil can thereafter be used for subject-specific correction of k-space trajectories-notably in gradient-sensitive sequences such as single-shot spiral imaging.

METHODS

The transmit and receive performance was evaluated before and after the integration of field probes, whereas field probes were evaluated when in an optimal configuration external to the coil and after their integration. Diffusion-weighted EPI and single-shot spiral acquisitions were employed to evaluate the efficacy of correcting higher order field perturbations and the consequent effect on image quality.

RESULTS

Field probes had a negligible effect on RF-coil performance, including the transmit efficiency, transmit uniformity, and mean SNR over the brain. Modest reductions in field-probe signal lifetimes were observed, caused primarily by nonidealities in the gradient and shim fields of the head-only gradient coil at the probe positions. The field-monitoring system could correct up to second-order field perturbations in single-shot spiral imaging.

CONCLUSION

The integrated RF coil and field camera was capable of concurrent-field monitoring within a 7T head-only scanner and facilitated the subsequent correction of k-space trajectories during spiral imaging.

Feasibility of spiral fMRI based on an LTI gradient model

Spiral imaging is very well suited for functional MRI, however its use has been limited by the fact that artifacts caused by gradient imperfections and B₀ inhomogeneity are more difficult to correct compared to EPI. Effective correction requires accurate knowledge of the traversed k-space trajectory. With the goal of making spiral fMRI more accessible, we have evaluated image reconstruction using trajectories predicted by the gradient impulse response function (GIRF), which can be determined in a one-time calibration step. GIRF-predicted reconstruction was tested for high-resolution (0.8 mm) fMRI at 7T. Image quality and functional results of the reconstructions using GIRF-prediction were compared to reconstructions using the nominal trajectory and concurrent field monitoring. The reconstructions using nominal spiral trajectories contain substantial artifacts and the activation maps contain misplaced activation. Image artifacts are substantially reduced when using the GIRF-predicted reconstruction, and the activation maps for the GIRF-predicted and monitored reconstructions largely overlap. The GIRF reconstruction provides a large increase in the spatial specificity of the activation compared to the nominal reconstruction. The GIRF-reconstruction generates image quality and fMRI results similar to using a concurrently monitored trajectory. The presented approach does not prolong or complicate the fMRI acquisition. Using GIRF-predicted trajectories has the potential to enable high-quality spiral fMRI in situations where concurrent trajectory monitoring is not available.

Diffusion-PEPTIDE: Distortion- and blurring-free diffusion imaging with self-navigated motion-correction and relaxometry capabilities

PURPOSE

To implement the time-resolved relaxometry PEPTIDE technique into a diffusion acquisition to provide self-navigated, distortion- and blurring-free diffusion imaging that is robust to motion, while simultaneously providing T₂ and T 2 ∗ mapping.

THEORY AND METHODS

The PEPTIDE readout was implemented into a spin-echo diffusion acquisition, enabling reconstruction of a time-series of T₂ - and T 2 ∗ -weighted images, free from conventional echo planar imaging (EPI) distortion and blurring, for each diffusion-encoding. Robustness of PEPTIDE to motion and shot-to-shot phase variation was examined through a deliberate motion-corrupted diffusion experiment. Two diffusion-relaxometry in vivo brain protocols were also examined: (1)1 × 1 × 3 mm³ across 32 diffusion directions in 20 min, (2)1.5 × 1.5 × 3.0 mm³ across 6 diffusion-weighted images in 3.4 min. T₂ , T 2 ∗ , and diffusion parameter maps were calculated from these data. As initial exploration of the rich diffusion-relaxometry data content for use in multi-compartment modeling, PEPTIDE data were acquired of a gadolinium-doped asparagus phantom. These datasets contained two compartments with different relaxation parameters and different diffusion orientation properties, and T₂ relaxation variations across these diffusion directions were explored.

RESULTS

Diffusion-PEPTIDE showed the capability to provide high quality diffusion images and T₂ and T 2 ∗ maps from both protocols. The reconstructions were distortion-free, avoided potential resolution losses exceeding 100% in equivalent EPI acquisitions, and showed tolerance to nearly 30° of rotational motion. Expected variation in T₂ values as a function of diffusion direction was observed in the two-compartment asparagus phantom (P < .01), demonstrating potential to explore diffusion-PEPTIDE data for multi-compartment modeling.

CONCLUSIONS

Diffusion-PEPTIDE provides highly robust diffusion and relaxometry data and offers potential for future applications in diffusion-relaxometry multi-compartment modeling.