Motion Correction in fMRI by Mapping Slice-to-Volume with Concurrent Field-Inhomogeneity Correction

Head motion is the major source of error in measuring intensity changes related to given stimuli in fMRI. The effects of head motion are image shifts and field inhomogeneity variations which cause local changes in geometric distortions. The previously developed motion correction method, mapping slice-to-volume (MSV), retrospectively remaps slices that are shifted by head motion to their spatially correct locations in an anatomical reference. Images exhibiting spatially varying geometric distortions require non-linear mapping solutions. An accurate field map can be used for the correction of such spatial distortions. However, field-map changes with head motion and, in practice, only one field-map is available typically. This work evaluates the improved motion correction capability of MSV with concurrent iterative field-corrected reconstruction using only an initial field-map. The results from simulated motion data show effective convergence and accuracy in image registration for the correction of image artifacts complicated by the motion induced field effects.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85942/1/Fessler206.pd

Advances in Concurrent Motion and Field-Inhomogeneity Correction in Functional MRI.

Head motion and static magnetic field (B0) inhomogeneity are two important sources of image intensity variability in functional MRI (fMRI). Ideally, in MRI, any deviation in B0 homogeneity in an object occurs only by design. However, due to imperfections in the main magnet and gradient coils, and, magnetic susceptibility differences in the object, undesired B0 deviations may occur. This causes geometric distortion in Cartesian EPI images. In addition to spatial shifts and rotations of images, head motion during an fMRI experiment may induce time-varying field-inhomogeneity changes in the brain. As a result, correcting for motion and field-inhomogeneity effects independently of each other with a static field map may be insufficient, especially in the presence of large out-of-plane rotations. Our primary concern is the correction of the combined effects of motion and field-inhomogeneity induced geometric distortion in Cartesian EPI fMRI images. We formulate a concurrent field-inhomogeneity with map-slice-to-volume motion correction, and develop a motion-robust dual-echo bipolar gradient echo static field map estimation method. We also propose and evaluate a penalized weighted least squares approach to dynamic field map estimation using the susceptibility voxel convolution method. This technique accounts for field changes due to out-of-plane rotations, and estimates dynamic field maps from a high resolution static field map without requiring accurate image segmentation, or the use of literature susceptibility values. Experiments with simulated data suggest that the technique is promising, and the method will be applied to real data in future work. In a separate clinical fMRI project, which is independent of the above work, we also formulate a current density weighted index to quantify correspondence between electrocortical stimulation and fMRI maps for brain presurgical planning. The proposed index is formulated with the broader goal of defining safe limits for lesion resection, and is characterized extensively with simulated data. The index is also computed for real human datasets.Ph.D.Electrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/60787/1/tbyeo_1.pd

Advances in image acquisition and filtering for MRI neuroimaging at 7 tesla

Performing magnetic resonance imaging at high magnetic field strength promises many improvements over low fields that are of direct benefit in functional neuroimaging. This includes the possibility of improved signal-to-noise levels, and increased BOLD functional contrast and spatial specificity. However, human MRI at 7T and above suffers from unique engineering challenges that limit the achievable gains. In this thesis, three technological developments are introduced, all of which address separate issues associated with functional magnetic resonance neuroimaging at very high magnetic field strengths. First, the image homogeneity problem is addressed by investigating methods of RF shimming — modifying the excitation portion of the MRI experiment for use with multi-channel RF coils. It is demonstrated that in 2D MRI experiments, shimming on a slice-by slice basis allows utilization of an extra degree of freedom available from the slice dimension, resulting in significant gains in image homogeneity and reduced RF power requirements. After acceptable images are available, we move to address complications of high field imaging that manifest in the fMRI time series. In the second paper, the increased physiological noise present in BOLD time series at high field is addressed with a unique data-driven noise regressor scheme based upon information in the phase component of the MRI signal. It is demonstrated that this method identifies and removes a significant portion of physiological signals, and performs as good or better than other popular data driven methods that use only the magnitude signal information. Lastly, the BOLD phase signal is again leveraged to address the confounding role of veins in resting state BOLD fMRI experiments. The phase regressor technique (previously developed by Dr. Menon) is modified and applied to resting state fMRI to remove macro vascular contributions in the datasets, leading to changes in spatial extent and connectivity of common resting state networks on single subjects and at the group level

Mitigating susceptibility-induced distortions in high-resolution 3DEPI fMRI at 7T

Geometric distortion is a major limiting factor for spatial specificity in high-resolution fMRI using EPI readouts and is exacerbated at higher field strengths due to increased B0 field inhomogeneity. Prominent correction schemes are based on B0 field-mapping or acquiring reverse phase-encoded (reversed-PE) data. However, to date, comparisons of these techniques in the context of fMRI have only been performed on 2DEPI data, either at lower field or lower resolution. In this study, we investigate distortion compensation in the context of sub-millimetre 3DEPI data at 7T. B0 field-mapping and reversed-PE distortion correction techniques were applied to both partial coverage BOLD-weighted and whole brain MT-weighted 3DEPI data with matched distortion. Qualitative assessment showed overall improvement in cortical alignment for both correction techniques in both 3DEPI fMRI and whole-brain MT-3DEPI datasets. The distortion-corrected MT-3DEPI images were quantitatively evaluated by comparing cortical alignment with an anatomical reference using dice coefficient (DC) and correlation ratio (CR) measures. These showed that B0 field-mapping and reversed-PE methods both improved correspondence between the MT-3DEPI and anatomical data, with more substantial improvements consistently obtained using the reversed-PE approach. Regional analyses demonstrated that the largest benefit of distortion correction, and in particular of the reversed-PE approach, occurred in frontal and temporal regions where susceptibility-induced distortions are known to be greatest, but had not led to complete signal dropout. In conclusion, distortion correction based on reversed-PE data has shown the greater capacity for achieving faithful alignment with anatomical data in the context of high-resolution fMRI at 7T using 3DEPI

High resolution whole brain diffusion imaging at 7 T for the Human Connectome Project

Mapping structural connectivity in healthy adults for the Human Connectome Project (HCP) benefits from high quality, high resolution, multiband (MB)-accelerated whole brain diffusion MRI (dMRI). Acquiring such data at ultrahigh fields (7 T and above) can improve intrinsic signal-to-noise ratio (SNR), but suffers from shorter T2 and T2⁎ relaxation times, increased B1+ inhomogeneity (resulting in signal loss in cerebellar and temporal lobe regions), and increased power deposition (i.e. specific absorption rate (SAR)), thereby limiting our ability to reduce the repetition time (TR). Here, we present recent developments and optimizations in 7 T image acquisitions for the HCP that allow us to efficiently obtain high quality, high resolution whole brain in-vivo dMRI data at 7 T. These data show spatial details typically seen only in ex-vivo studies and complement already very high quality 3 T HCP data in the same subjects. The advances are the result of intensive pilot studies aimed at mitigating the limitations of dMRI at 7 T. The data quality and methods described here are representative of the datasets that will be made freely available to the community in 2015

Three-dimensional echo-shifted EPI with simultaneous blip-up and blip-down acquisitions for correcting geometric distortion

Purpose: Echo-planar imaging (EPI) with blip-up/down acquisition (BUDA) can provide high-quality images with minimal distortions by using two readout trains with opposing phase-encoding gradients. Because of the need for two separate acquisitions, BUDA doubles the scan time and degrades the temporal resolution when compared to single-shot EPI, presenting a major challenge for many applications, particularly functional MRI (fMRI). This study aims at overcoming this challenge by developing an echo-shifted EPI BUDA (esEPI-BUDA) technique to acquire both blip-up and blip-down datasets in a single shot. Methods: A three-dimensional (3D) esEPI-BUDA pulse sequence was designed by using an echo-shifting strategy to produce two EPI readout trains. These readout trains produced a pair of k-space datasets whose k-space trajectories were interleaved with opposite phase-encoding gradient directions. The two k-space datasets were separately reconstructed using a 3D SENSE algorithm, from which time-resolved B0-field maps were derived using TOPUP in FSL and then input into a forward model of joint parallel imaging reconstruction to correct for geometric distortion. In addition, Hankel structured low-rank constraint was incorporated into the reconstruction framework to improve image quality by mitigating the phase errors between the two interleaved k-space datasets. Results: The 3D esEPI-BUDA technique was demonstrated in a phantom and an fMRI study on healthy human subjects. Geometric distortions were effectively corrected in both phantom and human brain images. In the fMRI study, the visual activation volumes and their BOLD responses were comparable to those from conventional 3D echo-planar images. Conclusion: The improved imaging efficiency and dynamic distortion correction capability afforded by 3D esEPI-BUDA are expected to benefit many EPI applications.Comment: 8 figures, peer-reviewed journal pape

Motion Robust Magnetic Susceptibility and Field Inhomogeneity Estimation Using Regularized Image Restoration Techniques for fMRI

In functional MRI, head motion may cause dynamic nonlinear field-inhomogeneity changes, especially with large out-of-plane rotations. This may lead to dynamic geometric distortion or blurring in the time series, which may reduce activation detection accuracy. The use of image registration to estimate dynamic field inhomogeneity maps from a static field map is not sufficient in the presence of such rotations. This paper introduces a retrospective approach to estimate magnetic susceptibility induced field maps of an object in motion, given a static susceptibility induced field map and the associated object motion parameters. It estimates a susceptibility map from a static field map using regularized image restoration techniques, and applies rigid body motion to the former. The dynamic field map is then computed using susceptibility voxel convolution. The method addresses field map changes due to out-of-plane rotations during time series acquisition and does not involve real time field map acquisitions.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85944/1/Fessler233.pd

Cortical lamina-dependent blood volume changes in human brain at 7T

Cortical layer-dependent high (sub-millimeter) resolution functional magnetic resonance imaging (fMRI) in human or animal brain can be used to address questions regarding the functioning of cortical circuits, such as the effect of different afferent and efferent connectivities on activity in specific cortical layers. The sensitivity of gradient echo (GE) blood oxygenation level-dependent (BOLD) responses to large draining veins reduces its local specificity and can render the interpretation of the underlying laminar neural activity impossible. The application of the more spatially specific cerebral blood volume (CBV)-based fMRI in humans has been hindered by the low sensitivity of the noninvasive modalities available. Here, a vascular space occupancy (VASO) variant, adapted for use at high field, is further optimized to capture layer-dependent activity changes in human motor cortex at sub-millimeter resolution. Acquired activation maps and cortical profiles show that the VASO signal peaks in gray matter at 0.8–1.6 mm depth, and deeper compared to the superficial and vein-dominated GE-BOLD responses. Validation of the VASO signal change versus well-established iron-oxide contrast agent based fMRI methods in animals showed the same cortical profiles of CBV change, after normalization for lamina-dependent baseline CBV. In order to evaluate its potential of revealing small lamina-dependent signal differences due to modulations of the input-output characteristics, layer-dependent VASO responses were investigated in the ipsilateral hemisphere during unilateral finger tapping. Positive activation in ipsilateral primary motor cortex and negative activation in ipsilateral primary sensory cortex were observed. This feature is only visible in high-resolution fMRI where opposing sides of a sulcus can be investigated independently because of a lack of partial volume effects. Based on the results presented here, we conclude that VASO offers good reproducibility, high sensitivity and lower sensitivity than GE-BOLD to changes in larger vessels, making it a valuable tool for layer-dependent fMRI studies in humans

Transcranial direct current stimulation with functional magnetic resonance imaging: a detailed validation and operational guide [version 1; peer review: 1 approved with reservations]

Introduction: Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation technique used to modulate human brain and behavioural function in both research and clinical interventions. The combination of functional magnetic resonance imaging (fMRI) with tDCS enables researchers to directly test causal contributions of stimulated brain regions, answering questions about the physiology and neural mechanisms underlying behaviour. Despite the promise of the technique, advances have been hampered by technical challenges and methodological variability between studies, confounding comparability/replicability. / Methods: Here tDCS-fMRI at 3T was developed for a series of experiments investigating language recovery after stroke. To validate the method, one healthy volunteer completed an fMRI paradigm with three conditions: (i) No-tDCS, (ii) Sham-tDCS, (iii) 2mA Anodal-tDCS. MR data were analysed in SPM12 with region-of-interest (ROI) analyses of the two electrodes and reference sites. / Results: Quality assessment indicated no visible signal dropouts or distortions introduced by the tDCS equipment. After modelling scanner drift, motion-related variance, and temporal autocorrelation, we found no field inhomogeneity in functional sensitivity metrics across conditions in grey matter and in the three ROIs. / Discussion: Key safety factors and risk mitigation strategies that must be taken into consideration when integrating tDCS into an fMRI environment are outlined. To obtain reliable results, we provide practical solutions to technical challenges and complications of the method. It is hoped that sharing these data and SOP will promote methodological replication in future studies, enhancing the quality of tDCS-fMRI application, and improve the reliability of scientific results in this field. / Conclusions: The method and data provided here provide a technically safe, reliable tDCS-fMRI procedure to obtain high quality MR data. The detailed framework of the Standard Operation Procedure SOP (https://doi.org/10.5281/zenodo.4606564) systematically reports the technical and procedural elements of our tDCS-fMRI approach, which we hope can be adopted and prove useful in future studies