Real-time measurement and correction of both B0 changes and subject motion in diffusion tensor imaging using a double volumetric navigated (DvNav) sequence

Prospective motion correction for R 2 * and susceptibility mapping using spherical navigators

PURPOSE

To perform prospective motion correction (PMC) for improvedR 2 * $$ {R}_2^{\ast } $$ and susceptibility mapping using a purely navigator-based approach.

METHODS

Spherical navigators (SNAVs) were combined with an additional FID readout for simultaneous measurement of motion and zeroth-order field shifts. The resulting FIDSNAVs were interleaved for PMC of a multi-echo gradient echo sequence with retrospectiveB 0 $$ {B}_0 $$ correction. Experiments were performed on a 3T scanner with a 32-channel head coil. Performance was assessed in five volunteers with motion prompts derived from real unintentional motion trajectories.

RESULTS

At short TEs, PMC alone was sufficient to achieve good image quality; at longer TEs, retrospectiveB 0 $$ {B}_0 $$ correction was often just as important for artifact reduction as motion correction. Both PMC and retrospectiveB 0 $$ {B}_0 $$ correction reduced error inR 2 * $$ {R}_2^{\ast } $$ and susceptibility maps for all participants. Residual artifacts were observed in the most severe motion case.

CONCLUSION

Combining SNAVs with an additional FID readout enables simultaneous motion and field correction with no additional hardware requirements, improving the fidelity of quantitative mapping in the presence of motion.

Chemical shift-encoded multishot EPI for navigator-free prostate DWI

PURPOSE

DWI is an important contrast for prostate MRI to enable early and accurate detection of cancer. This study introduces a Dixon 3-shot-EPI protocol with structured low-rank reconstruction for navigator-free DWI. The aim is to overcome the limitations of single-shot EPI (ssh-EPI), such as geometric distortions and fat signal interference, while addressing the motion-induced phase variations of multishot EPI and simultaneously allowing water/fat separation.

METHODS

DWI data were acquired from 7 healthy volunteers using both Dixon 3-shot EPI and standard fat-suppressed ssh-EPI with similar scan times for comparison. Two readers evaluated image quality using a 5-point Likert scale regarding different aspects. The ADC values were quantitatively compared between protocols. To show feasibility in a clinical setting, the protocol was applied to two patients.

RESULTS

From the reader scores, Dixon 3-shot EPI significantly reduced geometric distortion compared with ssh-EPI (p < 0.01), with no significant differences in edge definition, SNR, or overall image quality. There was no significant difference in ADC values between the two protocols. However, the Dixon multishot-EPI protocol offered advantages such as self-referenced B0 map-driven distortion correction, greater flexibility in imaging parameters, and superior fat suppression. In the patient data, the lesion could be clearly identified in both protocols and on the associated ADC maps.

CONCLUSION

The proposed Dixon 3-shot-EPI protocol shows promise as an alternative to ssh-EPI for prostate DWI, providing reduced geometric distortions and improved fat suppression. It addresses common DWI issues based on EPI and enhances scanning flexibility, indicating potential for optimized imaging.

Advancements in MR hardware systems and magnetic field control: B0 shimming, RF coils, and gradient techniques for enhancing magnetic resonance imaging and spectroscopy

High magnetic field homogeneity is critical for magnetic resonance imaging (MRI), functional MRI, and magnetic resonance spectroscopy (MRS) applications. B0 inhomogeneity during MR scans is a long-standing problem resulting from magnet imperfections and site conditions, with the main issue being the inhomogeneity across the human body caused by differences in magnetic susceptibilities between tissues, resulting in signal loss, image distortion, and poor spectral resolution. Through a combination of passive and active shim techniques, as well as technological advances employing multi-coil techniques, optimal coil design, motion tracking, and real-time modifications, improved field homogeneity and image quality have been achieved in MRI/MRS. The integration of RF and shim coils brings a high shim efficiency due to the proximity of participants. This technique will potentially be applied to high-density RF coils with a high-density shim array for improved B0 homogeneity. Simultaneous shimming and image encoding can be achieved using multi-coil array, which also enables the development of novel encoding methods using advanced magnetic field control. Field monitoring enables the capture and real-time compensation for dynamic field perturbance beyond the static background inhomogeneity. These advancements have the potential to better use the scanner performance to enhance diagnostic capabilities and broaden applications of MRI/MRS in a variety of clinical and research settings. The purpose of this paper is to provide an overview of the latest advances in B0 magnetic field shimming and magnetic field control techniques as well as MR hardware, and to emphasize their significance and potential impact on improving the data quality of MRI/MRS.

Dynamic field mapping and distortion correction using single-shot blip-rewound EPI and joint multi-echo reconstruction

PURPOSE

To develop a method for dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping and distortion correction.

METHODS

A blip-rewound EPI trajectory was developed to acquire multiple 2D EPI images in a single readout with an interleaved order, which allows a short TE difference. A joint multi-echo reconstruction was utilized to exploit the shared information between EPI images. The reconstructed images from each readout are combined to produce a final magnitude image. A∆ B 0 $$ \Delta {B}_0 $$ map is calculated from the phase of these images for distortion correction. The efficacy of the proposed method is assessed with phantom and in vivo experiments. The performance of the proposed method in the presence of subject motion is also investigated.

RESULTS

Compared to conventional multi-echo EPI, the proposed method allows dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping at matched resolution with a much shorter TR. Phantom and in vivo results show that the proposed method can provide a comparable magnitude image as conventional single-shot EPI. The∆ B 0 $$ \Delta {B}_0 $$ maps calculated from the proposed method are consistent with conventional multi-echo EPI in the phantom experiment. For in vivo experiments, the proposed method provides a more accurate estimation of∆ B 0 $$ \Delta {B}_0 $$ than conventional multi-echo EPI, which is prone to phase wrapping problems due to the long TE difference. In-vivo scan with subject motion shows the proposed dynamic field mapping method can improve the temporal stability of EPI time series compared to gradient echo (GRE) based static field mapping.

CONCLUSION

The proposed method allows accurate dynamic∆ B 0 $$ \Delta {B}_0 $$ mapping for robust distortion correction without compromising spatial or temporal resolution.

Servo navigators: Linear regression and feedback control for rigid-body motion correction

PURPOSE

Navigator-based correction of rigid-body motion reconciling high precision with minimal acquisition, minimal calibration and simple, fast processing.

METHODS

A short orbital navigator (2.3 ms) is inserted in a three-dimensional (3D) gradient echo sequence for human head imaging. Head rotation and translation are determined by linear regression based on a complex-valued model built either from three reference navigators or in a reference-less fashion, from the first actual navigator. Optionally, the model is expanded by global phase and field offset. Run-time scan correction on this basis establishes servo control that maintains validity of the linear picture by keeping its expansion point stable in the head frame of reference. The technique is assessed in a phantom and demonstrated by motion-corrected imaging in vivo.

RESULTS

The proposed approach is found to establish stable motion control both with and without reference acquisition. In a phantom, it is shown to accurately detect motion mimicked by rotation of scan geometry as well as change in global B0 . It is demonstrated to converge to accurate motion estimates after perturbation well beyond the linear signal range. In vivo, servo navigation achieved motion detection with precision in the single-digit range of micrometers and millidegrees. Involuntary and intentional motion in the range of several millimeters were successfully corrected, achieving excellent image quality.

CONCLUSION

The combination of linear regression and feedback control enables prospective motion correction for head imaging with high precision and accuracy, short navigator readouts, fast run-time computation, and minimal demand for reference data.

Rapid and accurate navigators for motion and B0 tracking using QUEEN: Quantitatively enhanced parameter estimation from navigators

PURPOSE

To develop a framework that jointly estimates rigid motion and polarizing magnetic field (B0 ) perturbations (δ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ ) for brain MRI using a single navigator of a few milliseconds in duration, and to additionally allow for navigator acquisition at arbitrary timings within any type of sequence to obtain high-temporal resolution estimates.

THEORY AND METHODS

Methods exist that match navigator data to a low-resolution single-contrast image (scout) to estimate either motion orδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . In this work, called QUEEN (QUantitatively Enhanced parameter Estimation from Navigators), we propose combined motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimation from a fast, tailored trajectory with arbitrary-contrast navigator data. To this end, the concept of a quantitative scout (Q-Scout) acquisition is proposed from which contrast-matched scout data is predicted for each navigator. Finally, navigator trajectories, contrast-matched scout, andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ are integrated into a motion-informed parallel-imaging framework.

RESULTS

Simulations and in vivo experiments show the need to modelδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ to obtain accurate motion parameters estimated in the presence of strongδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Simulations confirm that tailored navigator trajectories are needed to robustly estimate both motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ . Furthermore, experiments show that a contrast-matched scout is needed for parameter estimation from multicontrast navigator data. A retrospective, in vivo reconstruction experiment shows improved image quality when using the proposed Q-Scout and QUEEN estimation.

CONCLUSIONS

We developed a framework to jointly estimate rigid motion parameters andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ from navigators. Combing a contrast-matched scout with the proposed trajectory allows for navigator deployment in almost any sequence and/or timing, which allows for higher temporal-resolution motion andδ B 0 $$ \delta {\mathbf{B}}_{\mathbf{0}} $$ estimates.

Prediction of motion induced magnetic fields for human brain MRI at 3 T

OBJECTIVE

Maps of B0 field inhomogeneities are often used to improve MRI image quality, even in a retrospective fashion. These field inhomogeneities depend on the exact head position within the static field but acquiring field maps (FM) at every position is time consuming. Here we propose a forward simulation strategy to obtain B0 predictions at different head-positions.

METHODS

FM were predicted by combining (1) a multi-class tissue model for estimation of tissue-induced fields, (2) a linear k-space model for capturing gradient imperfections, (3) a dipole estimation for quantifying lower-body perturbing fields (4) and a position-dependent tissue mask to model FM alterations caused by large motion effects. The performance of the combined simulation strategy was compared with an approach based on a rigid body transformation of the FM measured in the reference position to the new position.

RESULTS

The transformed FM provided inconsistent results for large head movements (> 5° rotation, approximately), while the simulation strategy had a superior prediction accuracy for all positions. The simulated FM was used to optimize B0 shims with up to 22.2% improvement with respect to the transformed FM approach.

CONCLUSION

The proposed simulation strategy is able to predict movement-induced B0 field inhomogeneities yielding more precise estimates of the ground truth field homogeneity than the transformed FM.

Water/fat separation for self-navigated diffusion-weighted multishot echo-planar imaging

The purpose of this study was to develop a self-navigation strategy to improve scan efficiency and image quality of water/fat-separated, diffusion-weighted multishot echo-planar imaging (ms-EPI). This is accomplished by acquiring chemical shift-encoded diffusion-weighted data and using an appropriate water-fat and diffusion-encoded signal model to enable reconstruction directly from k-space data. Multishot EPI provides reduced geometric distortion and improved signal-to-noise ratio in diffusion-weighted imaging compared with single-shot approaches. Multishot acquisitions require corrections for physiological motion-induced shot-to-shot phase errors using either extra navigators or self-navigation principles. In addition, proper fat suppression is important, especially in regions with large B₀ inhomogeneity. This makes the use of chemical shift encoding attractive. However, when combined with ms-EPI, shot-to-shot phase navigation can be challenging because of the spatial displacement of fat signals along the phase-encoding direction. In this work, a new model-based, self-navigated water/fat separation reconstruction algorithm is proposed. Experiments in legs and in the head-neck region of 10 subjects were performed to validate the algorithm. The results are compared with an image-based, two-dimensional (2D) navigated water/fat separation approach for ms-EPI and with a conventional fat saturation approach. Compared with the 2D navigated method, the use of self-navigation reduced the shot duration time by 30%-35%. The proposed algorithm provided improved diffusion-weighted water images in both leg and head-neck regions compared with the 2D navigator-based approach. The proposed algorithm also produced better fat suppression compared with the conventional fat saturation technique in the B₀ inhomogeneous regions. In conclusion, the proposed self-navigated reconstruction algorithm can produce superior water-only diffusion-weighted EPI images with less artefacts compared with the existing methods.

Real-time shimming with FID navigators

PURPOSE

To implement a method for real-time field control using rapid FID navigator (FIDnav) measurements and evaluate the efficacy of the proposed approach for mitigating dynamic field perturbations and improvingT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality.

METHODS

FIDnavs were embedded in a gradient echo sequence and a subject-specific linear calibration model was generated on the scanner to facilitate rapid shim updates in response to measured FIDnav signals. To confirm the accuracy of FID-navigated field updates, phantom and volunteer scans were performed with online updates of the scanner B0 shim settings. To evaluate improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality with real-time shimming, 10 volunteers were scanned at 3T while performing deep-breathing and nose-touching tasks designed to modulate the B0 field. Quantitative image quality metrics were compared with and without FID-navigated field control. An additional volunteer was scanned at 7T to evaluate performance at ultra-high field.

RESULTS

Applying measured FIDnav shim updates successfully compensated for applied global and linear field offsets in phantoms and across all volunteers. FID-navigated real-time shimming led to a substantial reduction in field fluctuations and a consequent improvement inT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality in volunteers performing deep-breathing and nose-touching tasks, with 7.57% ± 6.01% and 8.21% ± 10.90% improvement in peak SNR and structural similarity, respectively.

CONCLUSION

FIDnavs facilitate rapid measurement and application of field coefficients for slice-wise B0 shimming. The proposed approach can successfully counteract spatiotemporal field perturbations and substantially improvesT 2 * $$ {\mathrm{T}}_2^{\ast } $$ -weighted image quality, which is important for a variety of clinical and research applications, particularly at ultra-high field.

Data‐driven motion‐corrected brain MRI incorporating pose‐dependent B 0 fields

Purpose

To develop a fully data‐driven retrospective intrascan motion‐correction framework for volumetric brain MRI at ultrahigh field (7 Tesla) that includes modeling of pose‐dependent changes in polarizing magnetic (B₀) fields.

Theory and Methods

Tissue susceptibility induces spatially varying B₀ distributions in the head, which change with pose. A physics‐inspired B₀ model has been deployed to model the B₀ variations in the head and was validated in vivo. This model is integrated into a forward parallel imaging model for imaging in the presence of motion. Our proposal minimizes the number of added parameters, enabling the developed framework to estimate dynamic B₀ variations from appropriately acquired data without requiring navigators. The effect on data‐driven motion correction is validated in simulations and in vivo.

Results

The applicability of the physics‐inspired B₀ model was confirmed in vivo. Simulations show the need to include the pose‐dependent B₀ fields in the reconstruction to improve motion‐correction performance and the feasibility of estimating B₀ evolution from the acquired data. The proposed motion and B₀ correction showed improved image quality for strongly corrupted data at 7 Tesla in simulations and in vivo.

Conclusion

We have developed a motion‐correction framework that accounts for and estimates pose‐dependent B₀ fields. The method improves current state‐of‐the‐art data‐driven motion‐correction techniques when B₀ dependencies cannot be neglected. The use of a compact physics‐inspired B₀ model together with leveraging the parallel imaging encoding redundancy and previously proposed optimized sampling patterns enables a purely data‐driven approach.

Self-navigated prospective motion correction for 3D-EPI acquisition

PURPOSE

Although 3D EPI is more susceptible to motion artifacts than 2D EPI, it presents some benefits for functional MRI, including the absence of spin-history artifacts, greater potential for parallel imaging acceleration, and better functional sensitivity in high-resolution imaging. Here we present a self-navigated 3D-EPI sequence suitable for prospective motion-corrected functional MRI without additional hardware or pulses.

METHODS

For each volume acquisition, the first 24 of the 52 partitions being acquired are accumulated to a new feedback block that was added to the image reconstruction pipeline. After zero-filling the remaining partitions, the feedback block constructs a volumetric self-navigator (vSNav), co-registers it to the reference vSNav acquired during the first volume acquisition, and sends motion estimates to the sequence. The sequence then updates its FOV and acquires subsequent partitions with the adjusted FOV, until the next update is received. The sequence was validated without and with intentional motion in phantom and in vivo on a 3T Skyra.

RESULTS

For phantom scans, the FOV was updated 0.704 s after acquisition of the vSNav partitions, and for in vivo scans after 0.768 s. Both phantom and in vivo data demonstrated stable motion estimates in the absence of motion. For in vivo acquisitions, prospective head-pose estimates using the vSNav's and retrospective estimates with FLIRT (FMRIB's Linear Image Registration Tool) agreed to within 0.23 mm (< 10% of the slice thickness) and 0.14° in all directions.

CONCLUSION

Depending when motion occurs during a volume acquisition, the proposed method fully corrects the FOV and recovers image quality within one volume acquisition.

Model-based dynamic off-resonance correction for improved accelerated fMRI in awake behaving nonhuman primates

Purpose

To estimate dynamic off‐resonance due to vigorous body motion in accelerated fMRI of awake behaving nonhuman primates (NHPs) using the echo‐planar imaging reference navigator, in order to attenuate the effects of time‐varying off‐resonance on the reconstruction.

Methods

In NHP fMRI, the animal’s head is usually head‐posted, and the dynamic off‐resonance is mainly caused by motion in body parts that are distant from the brain and have low spatial frequency. Hence, off‐resonance at each frame can be approximated as a spatially linear perturbation of the off‐resonance at a reference frame, and is manifested as a relative linear shift in k‐space. Using GRAPPA operators, we estimated these shifts by comparing the navigator at each time frame with that at the reference frame. Estimated shifts were then used to correct the data at each frame. The proposed method was evaluated in phantom scans, simulations, and in vivo data.

Results

The proposed method is shown to successfully estimate spatially low‐order dynamic off‐resonance perturbations, including induced linear off‐resonance perturbations in phantoms, and is able to correct retrospectively corrupted data in simulations. Finally, it is shown to reduce ghosting artifacts and geometric distortions by up to 20% in simultaneous multislice in vivo acquisitions in awake‐behaving NHPs.

Conclusion

A method is proposed that does not need sequence modification or extra acquisitions and makes accelerated awake behaving NHP imaging more robust and reliable, reducing the gap between what is possible with NHP protocols and state‐of‐the‐art human imaging.

Real-time motion and retrospective coil sensitivity correction for CEST using volumetric navigators (vNavs) at 7T

PURPOSE

To explore the impact of temporal motion-induced coil sensitivity changes on CEST-MRI at 7T and its correction using interleaved volumetric EPI navigators, which are applied for real-time motion correction.

METHODS

Five healthy volunteers were scanned via CEST. A 4-fold correction pipeline allowed the mitigation of (1) motion, (2) motion-induced coil sensitivity variations, ΔB1- , (3) motion-induced static magnetic field inhomogeneities, ΔB₀ , and (4) spatially varying transmit RF field fluctuations, ΔB1+ . Four CEST measurements were performed per session. For the first 2, motion correction was turned OFF and then ON in absence of voluntary motion, whereas in the other 2 controlled head rotations were performed. During post-processing ΔB1- was removed additionally for the motion-corrected cases, resulting in a total of 6 scenarios to be compared. In all cases, retrospective ∆B₀ and - ΔB1+ corrections were performed to compute artifact-free magnetization transfer ratio maps with asymmetric analysis (MTR_asym ).

RESULTS

Dynamic ΔB1- correction successfully mitigated signal deviations caused by head motion. In 2 frontal lobe regions of volunteer 4, induced relative signal errors of 10.9% and 3.9% were reduced to 1.1% and 1.0% after correction. In the right frontal lobe, the motion-corrected MTR_asym contrast deviated 0.92%, 1.21%, and 2.97% relative to the static case for Δω = 1, 2, 3 ± 0.25 ppm. The additional application of ΔB1- correction reduced these deviations to 0.10%, 0.14%, and 0.42%. The fully corrected MTR_asym values were highly consistent between measurements with and without intended head rotations.

CONCLUSION

Temporal ΔB1- cause significant CEST quantification bias. The presented correction pipeline including the proposed retrospective ΔB1- correction significantly reduced motion-related artifacts on CEST-MRI.

Characterization of B0-field fluctuations in prostate MRI

Multi-parametric MRI is increasingly used for prostate cancer detection. Improving information from current sequences, such as T2-weighted and diffusion-weighted (DW) imaging, and additional sequences, such as magnetic resonance spectroscopy (MRS) and chemical exchange saturation transfer (CEST), may enhance the performance of multi-parametric MRI. The majority of these techniques are sensitive to B0-field variations and may result in image distortions including signal pile-up and stretching (echo planar imaging (EPI) based DW-MRI) or unwanted shifts in the frequency spectrum (CEST and MRS). Our aim is to temporally and spatially characterize B0-field changes in the prostate. Ten male patients are imaged using dual-echo gradient echo sequences with varying repetitions on a 3 T scanner to evaluate the temporal B0-field changes within the prostate. A phantom is also imaged to consider no physiological motion. The spatial B0-field variations in the prostate are reported as B0-field values (Hz), their spatial gradients (Hz/mm) and the resultant distortions in EPI based DW-MRI images (b-value = 0 s/mm2 and two oppositely phase encoded directions). Over a period of minutes, temporal changes in B0-field values were ≤19 Hz for minimal bowel motion and ≥30 Hz for large motion. Spatially across the prostate, the B0-field values had an interquartile range of ≤18 Hz (minimal motion) and ≤44 Hz (large motion). The B0-field gradients were between -2 and 5 Hz/mm (minimal motion) and 2 and 12 Hz/mm (large motion). Overall, B0-field variations can affect DW, MRS and CEST imaging of the prostate.

Simultaneous feedback control for joint field and motion correction in brain MRI

T2*-weighted gradient-echo sequences count among the most widely used techniques in neuroimaging and offer rich magnitude and phase contrast. The susceptibility effects underlying this contrast scale with B₀, making T2*-weighted imaging particularly interesting at high field. High field also benefits baseline sensitivity and thus facilitates high-resolution studies. However, enhanced susceptibility effects and high target resolution come with inherent challenges. Relying on long echo times, T2*-weighted imaging not only benefits from enhanced local susceptibility effects but also suffers from increased field fluctuations due to moving body parts and breathing. High resolution, in turn, renders neuroimaging particularly vulnerable to motion of the head. This work reports the implementation and characterization of a system that aims to jointly address these issues. It is based on the simultaneous operation of two control loops, one for field stabilization and one for motion correction. The key challenge with this approach is that the two loops both operate on the magnetic field in the imaging volume and are thus prone to mutual interference and potential instability. This issue is addressed at the levels of sensing, timing, and control parameters. Performance assessment shows the resulting system to be stable and exhibit adequate loop decoupling, precision, and bandwidth. Simultaneous field and motion control is then demonstrated in examples of T2*-weighted in vivo imaging at 7T.

Dynamic distortion correction for functional MRI using FID navigators

PURPOSE

To develop a method for slice-wise dynamic distortion correction for EPI using rapid spatiotemporal B₀ field measurements from FID navigators (FIDnavs) and to evaluate the efficacy of this new approach relative to an established data-driven technique.

METHODS

A low-resolution reference image was used to create a forward model of FIDnav signal changes to enable estimation of spatiotemporal B₀ inhomogeneity variations up to second order from measured FIDnavs. Five volunteers were scanned at 3 T using a 64-channel coil with FID-navigated EPI. The accuracy of voxel shift measurements and geometric distortion correction was assessed for experimentally induced magnetic field perturbations. The temporal SNR was evaluated in EPI time-series acquired at rest and with a continuous nose-touching action, before and after image realignment.

RESULTS

Field inhomogeneity coefficients and voxel shift maps measured using FIDnavs were in excellent agreement with multi-echo EPI measurements. The FID-navigated distortion correction accurately corrected image geometry in the presence of induced magnetic field perturbations, outperforming the data-driven approach in regions with large field offsets. In functional MRI scans with nose touching, FIDnav-based correction yielded temporal SNR gains of 30% in gray matter. Following image realignment, which accounted for global image shifts, temporal SNR gains of 3% were achieved.

CONCLUSIONS

Our proposed application of FIDnavs enables slice-wise dynamic distortion correction with high temporal efficiency. We achieved improved signal stability by leveraging the encoding information from multichannel coils. This approach can be easily adapted to other EPI-based sequences to improve temporal SNR for a variety of clinical and research applications.

Improved detection of molecularly targeted iron oxide particles in mouse brain using B0 field stabilised high resolution MRI

PURPOSE

High resolution multi-gradient echo (MGE) scanning is typically used for detection of molecularly targeted iron oxide particles. The images of individual echoes are often combined to generate a composite image with improved SNR from the early echoes and boosted contrast from later echoes. In 3D implementations prolonged scanning at high gradient duty cycles induces a B0 shift that predominantly affects image alignment in the slow phase encoding dimension of 3D MGE images. The effect corrupts the composite echo image and limits the image resolution that is realised. A real-time adaptive B0 stabilisation during respiration gated 3D MGE scanning is shown to reduce image misalignment and improve detection of molecularly targeted iron oxide particles in composite images of the mouse brain.

METHODS

An optional B0 measurement block consisting of a 16 μs hard pulse with FA 1°, an acquisition delay of 3.2 ms, followed by gradient spoiling in all three axes was added to a respiration gated 3D MGE scan. During the acquisition delay of each B0 measurement block the NMR signal was routed to a custom built B0 stabilisation unit which mixed the signal to an audio frequency nominally centred around 1000 Hz to enable an Arduino based single channel receiver to measure frequency shifts. The frequency shift was used to effect correction to the main magnetic field via the B0 coil. The efficacy of B0 stabilisation and respiration gating was validated in vivo and used to improve detection of molecularly targeted microparticles of iron oxide (MPIO) in a mouse model of acute neuroinflammation.

RESULTS

Without B0 stabilisation 3D MGE image data exhibit varying mixtures of translation, scaling and blurring, which compromise the fidelity of the composite image. The real-time adaptive B0 stabilisation minimises corruption of the composite image as the images from the different echoes are properly aligned. The improved detection of molecularly targeted MPIO easily compensates for the scan time penalty of 14% incurred by the B0 stabilisation method employed. Respiration gating of the B0 measurement and the MRI scan was required to preserve high resolution detail, especially towards the back of the brain.

CONCLUSIONS

High resolution imaging for the detection of molecularly targeted iron oxide particles in the mouse brain requires good stabilisation of the main B0 field, and can benefit from a respiration gated image acquisition strategy.

Intra-session and inter-subject variability of 3D-FID-MRSI using single-echo volumetric EPI navigators at 3T

Purpose

In this study, we demonstrate the first combination of 3D FID proton MRSI and spatial encoding via concentric‐ring trajectories (CRTs) at 3T. FID‐MRSI has many benefits including high detection sensitivity, in particular for J‐coupled metabolites (e.g., glutamate/glutamine). This makes it highly attractive, not only for clinical, but also for, potentially, functional MRSI. However, this requires excellent reliability and temporal stability. We have, therefore, augmented this 3D‐FID‐MRSI sequence with single‐echo, imaging‐based volumetric navigators (se‐vNavs) for real‐time motion/shim‐correction (SHMOCO), which is 2× quicker than the original double‐echo navigators (de‐vNavs), hence allowing more efficient integration also in short‐TR sequences.

Methods

The tracking accuracy (position and B₀‐field) of our proposed se‐vNavs was compared to the original de‐vNavs in phantoms (rest and translation) and in vivo (voluntary head rotation). Finally, the intra‐session stability of a 5:40 min 3D‐FID‐MRSI scan was evaluated with SHMOCO and no correction (NOCO) in 5 resting subjects. Intra/inter‐subject coefficients of variation (CV) and intra‐class correlations (ICC) over the whole 3D volume and in selected regions of interest ROI were assessed.

Results

Phantom and in vivo scans showed highly consistent tracking performance for se‐vNavs compared to the original de‐vNavs, but lower frequency drift. Up to ~30% better intra‐subject CVs were obtained for SHMOCO (P < 0.05), with values of 9.3/6.9/6.5/7.8% over the full VOI for Glx/tNAA/tCho/m‐Ins ratios to tCr. ICCs were good‐to‐high (91% for Glx/tCr in motor cortex), whereas the inter‐subject variability was ~11–19%.

Conclusion

Real‐time motion/shim corrected 3D‐FID‐MRSI with time‐efficient CRT‐sampling at 3T allows reliable, high‐resolution metabolic imaging that is fast enough for clinical use and even, potentially, for functional MRSI.

Repeatability and reproducibility of prospective motion- and shim corrected 2D glycoCEST MRI

Background

Repeated glycoCEST MRI measurements on the same subject should produce similar results under the same environmental and experimental conditions. However, fluctuations in the static B₀ field, which may occur between and within measurements due to heating of the shim iron or subject motion, may alter results and affect reproducibility. Here we investigate the repeatability and reproducibility of glycoCEST measurements and examine the effectiveness of a real-time shim- and motion navigated chemical exchange saturation transfer (CEST) sequence to improve reproducibility.

Methods

In nine subjects, double volumetric navigated (DvNav)-CEST acquisitions in the calf muscle were repeated five times in each of two sessions-the first without correction, and the second with real-time shim- and motion correction applied. In both sessions a dynamically changing field was introduced by running a 5-minute gradient intensive diffusion sequence. We evaluated the effect of the introduced B₀ inhomogeneity on the reproducibility of glycoCEST, where the small chemical shift difference between the hydroxyl and bulk water protons at 3 T makes CEST quantification extremely sensitive to magnetic field inhomogeneities.

Results

With real-time shim- and motion correction, glycoCEST results were relatively consistent with mean coefficient of variation (CoV) 2.7%±1.4% across all subjects, whereas without correction the results were less consistent with CoV 84%±71%.

Conclusions

Our results demonstrate that real-time shim- and motion correction can mitigate effects of B0 field fluctuations and improve reproducibility of glycoCEST data. This is important when conducting longitudinal studies or when using glycoCEST MRI to assess treatment or physiological responses over time.

Rapid measurement and correction of spatiotemporal B0 field changes using FID navigators and a multi-channel reference image

PURPOSE

To measure spatiotemporal B₀ field changes in real time using FID navigators (FIDnavs) and to demonstrate the efficacy of retrospectively correcting high-resolution T 2 * -weighted images using a novel FIDnav framework.

METHODS

A forward model of the complex FIDnav signals was generated by simulating the effect of changes in the underlying B₀ inhomogeneity coefficients, with spatial encoding provided by a multi-channel reference image. Experiments were performed at 3T to assess the accuracy of B₀ field estimates from FIDnavs acquired from a 64-channel head coil under different shim settings and in 5 volunteers performing deep-breathing and nose-touching tasks designed to modulate the B₀ field. Second-order, in-plane spherical harmonic (SH) inhomogeneity coefficients estimated from FIDnavs were incorporated into an iterative reconstruction to retrospectively correct 2D gradient-echo images acquired in both axial and sagittal planes.

RESULTS

Spatiotemporal B₀ field changes measured from rapidly acquired FIDnavs were in good agreement with the results of second-order SH fitting to the measured field maps. FIDnav field estimates accounted for a significant proportion of the ΔB₀ variance induced by deep breathing (64 ± 21%) and nose touching (67 ± 34%) across all volunteers. Ghosting, blurring, and intensity modulation artifacts in T 2 * -weighted images, induced by spatiotemporal field changes, were visibly reduced following retrospective correction with FIDnav inhomogeneity coefficients.

CONCLUSIONS

Spatially resolved B₀ inhomogeneity changes up to second order can be characterized in real time using the proposed approach. Retrospective FIDnav correction substantially improves T 2 * -weighted image quality in the presence of strong B₀ field modulations, with potential for real-time shimming.

A comparison of static and dynamic ∆B0 mapping methods for correction of CEST MRI in the presence of temporal B0 field variations

PURPOSE

To assess the performance, in the presence of scanner instabilities, of three dynamic correction methods which integrate ∆B₀ mapping into the chemical exchange saturation transfer (CEST) measurement and three established static ∆B₀ -correction approaches.

METHODS

A homogeneous phantom and five healthy volunteers were scanned with a CEST sequence at 7 T. The in vivo measurements were performed twice: first with unaltered system frequency and again applying frequency shifts during the CEST acquisition. In all cases, retrospective voxel-wise ∆B₀ -correction was performed using one intrinsic and two extrinsic [prescans with dual-echo gradient-echo and water saturation shift referencing (WASSR)] static approaches. These were compared with two intrinsic [using phase data directly generated by single-echo or double-echo GRE (gradient-echo) CEST readout (CEST-GRE-2TE)] and one extrinsic [phase from interleaved dual-echo EPI (echo planar imaging) navigator (NAV-EPI-2TE)] dynamic ∆B₀ -correction approaches [allowing correction of each Z-spectral point before magnetization transfer ratio asymmetry (MTR_asym) analysis].

RESULTS

All three dynamic methods successfully mapped the induced drift. The intrinsic approaches were affected by the CEST labeling near water (∆ω < |0.3| ppm). The MTR_asym contrast was distorted by the frequency drift in the brain by up to 0.21%/Hz when static ∆B₀ -corrections were applied, whereas the dynamic ∆B₀ corrections reduced this to <0.01%/Hz without the need of external scans. The CEST-GRE-2TE and NAV-EPI-2TE resulted in highly consistent MTR_asym values with/without drift for all subjects.

CONCLUSION

Reliable correction of scanner instabilities is essential to establish clinical CEST MRI. The three dynamic approaches presented improved the ∆B₀ -correction performance significantly in the presence of frequency drift compared to established static methods. Among them, the self-corrected CEST-GRE-2TE was the most accurate and straightforward to implement.

Whole-slice mapping of GABA and GABA+ at 7T via adiabatic MEGA-editing, real-time instability correction, and concentric circle readout

An adiabatic MEscher-GArwood (MEGA)-editing scheme, using asymmetric hyperbolic secant editing pulses, was developed and implemented in a B1+-insensitive, 1D-semiLASER (Localization by Adiabatic SElective Refocusing) MR spectroscopic imaging (MRSI) sequence for the non-invasive mapping of γ-aminobutyric acid (GABA) over a whole brain slice. Our approach exploits the advantages of edited-MRSI at 7T while tackling challenges that arise with ultra-high-field-scans. Spatial-spectral encoding, using density-weighted, concentric circle echo planar trajectory readout, enabled substantial MRSI acceleration and an improved point-spread-function, thereby reducing extracranial lipid signals. Subject motion and scanner instabilities were corrected in real-time using volumetric navigators optimized for 7T, in combination with selective reacquisition of corrupted data to ensure robust subtraction-based MEGA-editing. Simulations and phantom measurements of the adiabatic MEGA-editing scheme demonstrated stable editing efficiency even in the presence of ±0.15 ppm editing frequency offsets and B1+ variations of up to ±30% (as typically encountered in vivo at 7T), in contrast to conventional Gaussian editing pulses. Volunteer measurements were performed with and without global inversion recovery (IR) to study regional GABA levels and their underlying, co-edited, macromolecular (MM) signals at 2.99 ppm. High-quality in vivo spectra allowed mapping of pure GABA and MM-contaminated GABA+ (GABA + MM) along with Glx (Glu + Gln), with high-resolution (eff. voxel size: 1.4 cm3) and whole-slice coverage in 24 min scan time. Metabolic ratio maps of GABA/tNAA, GABA+/tNAA, and Glx/tNAA were correlated linearly with the gray matter fraction of each voxel. A 2.15-fold increase in gray matter to white matter contrast was observed for GABA when enabling IR, which we attribute to the higher abundance of macromolecules at 2.99 ppm in the white matter than in the gray matter. In conclusion, adiabatic MEGA-editing with 1D-semiLASER selection is as a promising approach for edited-MRSI at 7T. Our sequence capitalizes on the benefits of ultra-high-field MRSI while successfully mitigating the challenges related to B0/B1+ inhomogeneities, prolonged scan times, and motion/scanner instability artifacts. Robust and accurate 2D mapping has been shown for the neurotransmitters GABA and Glx.

Real-time simultaneous shim and motion measurement and correction in glycoCEST MRI using double volumetric navigators (DvNavs)

PURPOSE

CEST MRI allows for indirect detection of molecules with exchangeable protons, measured as a reduction in water signal because of continuous transfer of saturated protons. CEST requires saturation pulses on the order of a second, as well as repeated acquisitions at different offset frequencies. The resulting extended scan time makes CEST susceptible to subject motion, which introduces field inhomogeneity, shifting offset frequencies and causing distortions in CEST spectra that resemble true CEST effects. This is a particular problem for molecules that resonate close to water, such as hydroxyl group in glycogen. To address this, a technique for real-time measurement and correction of motion and field inhomogeneity is proposed.

METHODS

A CEST sequence was modified to include double volumetric navigators (DvNavs) for real-time simultaneous motion and shim correction. Phantom tests were conducted to investigate the effects of motion and shim changes on CEST quantification and to validate the accuracy of DvNav motion and shim estimates. To evaluate DvNav shim and motion correction in vivo, acquisitions including 5 experimental conditions were performed in the calf muscle of 2 volunteers.

RESULTS

Phantom data show that DvNav-CEST accurately estimates frequency and linear gradient changes because of motion and corrects resulting image distortions. In addition, DvNav-CEST improves CEST quantification in vivo in the presence of motion.

CONCLUSION

The proposed technique allows for real-time simultaneous motion and shim correction with no additional scanning time, enabling accurate CEST quantification even in the presence of motion and field inhomogeneity.

The impact of in-scanner head motion on structural connectivity derived from diffusion MRI

Multiple studies have shown that data quality is a critical confound in the construction of brain networks derived from functional MRI. This problem is particularly relevant for studies of human brain development where important variables (such as participant age) are correlated with data quality. Nevertheless, the impact of head motion on estimates of structural connectivity derived from diffusion tractography methods remains poorly characterized. Here, we evaluated the impact of in-scanner head motion on structural connectivity using a sample of 949 participants (ages 8-23 years old) who passed a rigorous quality assessment protocol for diffusion magnetic resonance imaging (dMRI) acquired as part of the Philadelphia Neurodevelopmental Cohort. Structural brain networks were constructed for each participant using both deterministic and probabilistic tractography. We hypothesized that subtle variation in head motion would systematically bias estimates of structural connectivity and confound developmental inference, as observed in previous studies of functional connectivity. Even following quality assurance and retrospective correction for head motion, eddy currents, and field distortions, in-scanner head motion significantly impacted the strength of structural connectivity in a consistency- and length-dependent manner. Specifically, increased head motion was associated with reduced estimates of structural connectivity for network edges with high inter-subject consistency, which included both short- and long-range connections. In contrast, motion inflated estimates of structural connectivity for low-consistency network edges that were primarily shorter-range. Finally, we demonstrate that age-related differences in head motion can both inflate and obscure developmental inferences on structural connectivity. Taken together, these data delineate the systematic impact of head motion on structural connectivity, and provide a critical context for identifying motion-related confounds in studies of structural brain network development.

Effect of head motion on MRI B0 field distribution

PURPOSE

To identify and characterize the sources of B₀ field changes due to head motion, to reduce motion sensitivity in human brain MRI.

METHODS

B₀ fields were measured in 5 healthy human volunteers at various head poses. After measurement of the total field, the field originating from the subject was calculated by subtracting the external field generated by the magnet and shims. A subject-specific susceptibility model was created to quantify the contribution of the head and torso. The spatial complexity of the field changes was analyzed using spherical harmonic expansion.

RESULTS

Minor head pose changes can cause substantial and spatially complex field changes in the brain. For rotations and translations of approximately 5 º and 5 mm, respectively, at 7 T, the field change that is associated with the subject's magnetization generates a standard deviation (SD) of about 10 Hz over the brain. The stationary torso contributes to this subject-associated field change significantly with a SD of about 5 Hz. The subject-associated change leads to image-corrupting phase errors in multi-shot <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> -weighted acquisitions.

CONCLUSION

The B₀ field changes arising from head motion are problematic for multishot <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow><mml:msubsup><mml:mi>T</mml:mi> <mml:mn>2</mml:mn> <mml:mo>*</mml:mo></mml:msubsup> </mml:mrow> </mml:math> -weighted imaging. Characterization of the underlying sources provides new insights into mitigation strategies, which may benefit from individualized predictive field models in addition to real-time field monitoring and correction strategies.

Susceptibility-induced distortion that varies due to motion: Correction in diffusion MR without acquiring additional data

Because of their low bandwidth in the phase-encode (PE) direction, the susceptibility-induced off-resonance field causes distortions in echo planar imaging (EPI) images. It is therefore crucial to correct for susceptibility-induced distortions when performing diffusion studies using EPI. The susceptibility-induced field is caused by the object (head) disrupting the field and it is typically assumed that it remains constant within a framework defined by the object, (i.e. it follows the object as it moves in the scanner).

However, this is only approximately true. When a non-spherical object rotates around an axis other than that parallel with the magnetic flux (the z-axis) it changes the way it disrupts the field, leading to different distortions. Hence, if using a single field to correct for distortions there will be residual distortions in the volumes where the object orientation is substantially different to that when the field was measured.

In this paper we present a post-processing method for estimating the field as it changes with motion during the course of an experiment. It only requires a single measured field and knowledge of the orientation of the subject when that field was acquired. The volume-to-volume changes of the field as a consequence of subject movement are estimated directly from the diffusion data without the need for any additional or special acquisitions. It uses a generative model that predicts how each volume would look predicated on field change and inverts that model to yield an estimate of the field changes. It has been validated on both simulations and experimental data. The results show that we are able to track the field with high accuracy and that we are able to correct the data for the adverse effects of the changing field.

Real-time motion analytics during brain MRI improve data quality and reduce costs

Head motion systematically distorts clinical and research MRI data. Motion artifacts have biased findings from many structural and functional brain MRI studies. An effective way to remove motion artifacts is to exclude MRI data frames affected by head motion. However, such post-hoc frame censoring can lead to data loss rates of 50% or more in our pediatric patient cohorts. Hence, many scanner operators collect additional 'buffer data', an expensive practice that, by itself, does not guarantee sufficient high-quality MRI data for a given participant. Therefore, we developed an easy-to-setup, easy-to-use Framewise Integrated Real-time MRI Monitoring (FIRMM) software suite that provides scanner operators with head motion analytics in real-time, allowing them to scan each subject until the desired amount of low-movement data has been collected. Our analyses show that using FIRMM to identify the ideal scan time for each person can reduce total brain MRI scan times and associated costs by 50% or more.

Quantitative assessment of the susceptibility artefact and its interaction with motion in diffusion MRI

In this paper we evaluate the three main methods for correcting the susceptibility-induced artefact in diffusion-weighted magnetic-resonance (DW-MR) data, and assess how correction is affected by the susceptibility field's interaction with motion. The susceptibility artefact adversely impacts analysis performed on the data and is typically corrected in post-processing. Correction strategies involve either registration to a structural image, the application of an acquired field-map or the use of additional images acquired with different phase-encoding. Unfortunately, the choice of which method to use is made difficult by the absence of any systematic comparisons of them. In this work we quantitatively evaluate these methods, by extending and employing a recently proposed framework that allows for the simulation of realistic DW-MR datasets with artefacts. Our analysis separately evaluates the ability for methods to correct for geometric distortions and to recover lost information in regions of signal compression. In terms of geometric distortions, we find that registration-based methods offer the poorest correction. Field-mapping techniques are better, but are influenced by noise and partial volume effects, whilst multiple phase-encode methods performed best. We use our simulations to validate a popular surrogate metric of correction quality, the comparison of corrected data acquired with AP and LR phase-encoding, and apply this surrogate to real datasets. Furthermore, we demonstrate that failing to account for the interaction of the susceptibility field with head movement leads to increased errors when analysing DW-MR data. None of the commonly used post-processing methods account for this interaction, and we suggest this may be a valuable area for future methods development.

Analysis and correction of field fluctuations in fMRI data using field monitoring

This work investigates the role of magnetic field fluctuations as a confound in fMRI. In standard fMRI experiments with single-shot EPI acquisition at 3 Tesla the uniform and gradient components of the magnetic field were recorded with NMR field sensors. By principal component analysis it is found that differences of field evolution between the EPI readouts are explainable by few components relating to slow and within-shot field dynamics of hardware and physiological origin. The impact of fluctuating field components is studied by selective data correction and assessment of its influence on image fluctuation and SFNR. Physiological field fluctuations, attributed to breathing, were found to be small relative to those of hardware origin. The dominant confounds were hardware-related and attributable to magnet drift and thermal changes. In raw image time series, field fluctuation caused significant SFNR loss, reflected by a 67% gain upon correction. Large part of this correction can be accomplished by traditional image realignment, which addresses slow and spatially uniform field changes. With realignment, explicit field correction increased the SFNR on the order of 6%. In conclusion, field fluctuations are a relevant confound in fMRI and can be addressed effectively by retrospective data correction. Based on the physics involved it is anticipated that the advantage of full field correction increases with field strength, with non-Cartesian readouts, and upon phase-sensitive BOLD analysis.

Retrospective 3D motion correction using spherical navigator echoes

PURPOSE

To develop and evaluate a rapid spherical navigator echo (SNAV) motion correction technique, then apply it for retrospective correction of brain images.

METHODS

The pre-rotated, template matching SNAV method (preRot-SNAV) was developed in combination with a novel hybrid baseline strategy, which includes acquired and interpolated templates. Specifically, the SNAV templates are only rotated around X- and Y-axis; for each rotated SNAV, simulated baseline templates that mimic object rotation about the Z-axis were interpolated. The new method was first evaluated with phantom experiments. Then, a customized SNAV-interleaved gradient echo sequence was used to image three volunteers performing directed head motion. The SNAV motion measurements were used to retrospectively correct the brain images. Experiments were performed using a 3.0T whole-body MRI scanner and both single and 8-channel head coils.

RESULTS

Phantom rotations and translations measured using the hybrid baselines agreed to within 0.9° and 1mm compared to those measured with the original preRot-SNAV method. Retrospective motion correction of in vivo images using the hybrid preRot-SNAV effectively corrected for head rotation up to 4° and 4mm.

CONCLUSIONS

The presented hybrid approach enables the acquisition of pre-rotated baseline templates in as little as 2.5s, and results in accurate measurement of rotations and translations. Retrospective 3D motion correction successfully reduced motion artifacts in vivo.

Assessing the performance of different DTI motion correction strategies in the presence of EPI distortion correction

Diffusion tensor imaging (DTI) is susceptible to several artifacts due to eddy currents, echo planar imaging (EPI) distortion and subject motion. While several techniques correct for individual distortion effects, no optimal combination of DTI acquisition and processing has been determined. Here, the effects of several motion correction techniques are investigated while also correcting for EPI distortion: prospective correction, using navigation; retrospective correction, using two different popular packages (FSL and TORTOISE); and the combination of both methods. Data from a pediatric group that exhibited incidental motion in varying degrees are analyzed. Comparisons are carried while implementing eddy current and EPI distortion correction. DTI parameter distributions, white matter (WM) maps and probabilistic tractography are examined. The importance of prospective correction during data acquisition is demonstrated. In contrast to some previous studies, results also show that the inclusion of retrospective processing also improved ellipsoid fits and both the sensitivity and specificity of group tractographic results, even for navigated data. Matches with anatomical WM maps are highest throughout the brain for data that have been both navigated and processed using TORTOISE. The inclusion of both prospective and retrospective motion correction with EPI distortion correction is important for DTI analysis, particularly when studying subject populations that are prone to motion. Hum Brain Mapp 37:4405-4424, 2016. © 2016 Wiley Periodicals, Inc.