Functional MRI of the motor cortex using a conventional gradient system: comparison of FLASH and EPI techniques

Gradient echo (GE) and echo planar imaging (EPI) techniques are two different approaches to functional MRI (fMRI). In contrast to GE sequences, the ultra short EPI technique facilitates fMRI experiments with high spatial and temporal resolution or mapping of the whole brain. Although it has become the method of choice for fMRI, EPI is generally restricted to modern scanners with a strong gradient system. The aim of our study was to evaluate the applicability of EPI for fMRI of the motor cortex using a 1.5 T scanner with a conventional gradient system of 10 mT/m (rise time: 1 ms). Therefore, EPI was compared with a well-established high resolution fast low angle shot (FLASH) technique (matrix size 128²). The FLASH technique was applied additionally with a 64² matrix size to exclude influences caused by different spatial resolution, because the EPI sequence was restricted to a 64² matrix size. A total of 35 healthy volunteers were included in this study. The task consisted of clenching and spreading of the right hand. FLASH and EPI techniques were compared regarding geometric distortions as well as qualitative and quantitative fMRI criteria: Mean signal increase between activation and rest and the area of activation were measured within the contralateral, ipsilateral, and supplementary motor cortex. The quality of subtraction images between activation and rest, as well as the quality of z-maps and time course within activated regions of interest, was evaluated visually. EPI revealed significant distortions of the anterior and postior brain margins; lateral distortions (relevant for the motor cortex) could be neglected in most cases. The mean signal increase was significantly higher using FLASH 128² compared to FLASH 64² and EPI 64², whereas the activated areas proved to be smaller in FLASH 128² functional images. Both results can be explained by well-documented partial volume effects, caused by different voxel size. Similar quality of the subtraction images and of the time courses in different regions of interest were found for all techniques under investigation, but slightly reduced quality of z-map in FLASH 128². Within the limits of reproducibility and measurement accuracy, the location of contralateral activation was similar using FLASH and EPI sequences. In conclusion, EPI proved to be a reliable technique for fMRI of the motor cortex, even on an MR scanner with a conventional gradient system.     

Elimination of k-space spikes in fMRI data

The subtle signal changes in functional magnetic resonance imaging (fMRI) can be easily overwhelmed by noise of various origins. Spikes in the collected fMRI raw data often arise from high-duty usage of the scanner hardware and can introduce significant noise in the image and thereby in the image time series. Consequently, the spikes will corrupt the functional data and degrade the result of functional mapping. In this work, a simple method based on processing the time course of the k-space data are introduced and implemented to remove the spikes in the acquired data. Application of the method to experimental data shows that the methods are robust and effective for eliminating of spike-related noise in fMRI time series.     

Concurrent correction of geometric distortion and motion using the map-slice-to-volume method in echo-planar imaging

The accuracy of measuring voxel intensity changes between stimulus and rest images in fMRI echo-planar imaging (EPI) data is severely degraded in the presence of head motion. In addition, EPI is sensitive to susceptibility-induced geometric distortions. Head motion causes image shifts and associated field map changes that induce different geometric distortion at different time points. Conventionally, geometric distortion is "corrected" with a static field map independently of image registration. That approach ignores all field map changes induced by head motion. This work evaluates the improved motion correction capability of mapping slice to volume with concurrent iterative field corrected reconstruction using updated field maps derived from an initial static field map that has been spatially transformed and resampled. It accounts for motion-induced field map changes for translational and in-plane rotation motion. The results from simulated EPI time series data, in which motion, image intensity and activation ground truths are available, show improved accuracy in image registration, field corrected image reconstruction and activation detection.     

Surface-based functional magnetic resonance imaging analysis of partial brain echo planar imaging data at 1.5 T

Surface-based functional magnetic resonance imaging (fMRI) analysis is more sensitive and accurate than volume-based analysis for detecting neural activation. However, these advantages are less important in practical fMRI experiments with commonly used 1.5-T magnetic resonance devices because of the resolution gap between the echo planar imaging data and the cortical surface models. We expected high-resolution segmented partial brain echo planar imaging (EPI) data to overcome this problem, and the activation patterns of the high-resolution data could be different from the low-resolution data. For the practical applications of surface-based fMRI analysis using segmented EPI techniques, the effects of some important factors (e.g., activation patterns, registration and local distortions) should be intensively evaluated because the results of surface-based fMRI analyses could be influenced by them. In this study, we demonstrated the difference between activations detected from low-resolution EPI data, which were covering whole brain, and high-resolution segmented EPI data covering partial brain by volume- and surface-based analysis methods. First, we compared the activation maps of low- and high-resolution EPI datasets detected by volume- and surface-based analyses, with the spatial patterns of activation clusters, and analyzed the distributions of activations in occipital lobes. We also analyzed the high-resolution EPI data covering motor areas and fusiform gyri of human brain, and presented the differences of activations detected by volume- and surface-based methods.     

Semiautomated segmentation of the human spine based on echoplanar images

The number of functional magnetic resonance imaging (fMRI) studies performed on the human spinal cord (SC) has considerably increased in recent years. The lack of a validated processing pipeline is, however, a significant obstacle to the spread of SC fMRI. One component likely to be involved in any such pipeline is the process of SC masking, analogous to brain extraction in cerebral fMRI. In general, SC masking has been performed manually, with the incumbent costs of being very time consuming and operator dependent. To overcome these drawbacks, we have developed a tailored semiautomatic method for segmenting echoplanar images (EPI) of human spine that is able to identify the spinal canal and the SC. The method exploits both temporal and spatial features of the EPI series and was tested and optimized on EPI images of cervical spine acquired at 3 T. The dependence of algorithm performance on the degree of EPI image distortion was assessed by computing the displacement warping field that best matched the EPI to the corresponding high-resolution T(2) images. Segmentation accuracy was above 80%, a significant improvement over values obtained with similar approaches, but not exploiting temporal information. Geometric distortion was found to explain about 50% of the variance of algorithm classification efficiency.     

Functional magnetic resonance imaging of awake monkeys: some approaches for improving imaging quality

Functional magnetic resonance imaging (fMRI) at high magnetic field strength can suffer from serious degradation of image quality because of motion and physiological noise, as well as spatial distortions and signal losses due to susceptibility effects. Overcoming such limitations is essential for sensitive detection and reliable interpretation of fMRI data. These issues are particularly problematic in studies of awake animals. As part of our initial efforts to study functional brain activations in awake, behaving monkeys using fMRI at 4.7 T, we have developed acquisition and analysis procedures to improve image quality with encouraging results.We evaluated the influence of two main variables on image quality. First, we show how important the level of behavioral training is for obtaining good data stability and high temporal signal-to-noise ratios. In initial sessions, our typical scan session lasted 1.5 h, partitioned into short (<10 min) runs. During reward periods and breaks between runs, the monkey exhibited movements resulting in considerable image misregistrations. After a few months of extensive behavioral training, we were able to increase the length of individual runs and the total length of each session. The monkey learned to wait until the end of a block for fluid reward, resulting in longer periods of continuous acquisition. Each additional 60 training sessions extended the duration of each session by 60 min, culminating, after about 140 training sessions, in sessions that last about 4 h. As a result, the average translational movement decreased from over 500 μm to less than 80 μm, a displacement close to that observed in anesthetized monkeys scanned in a 7-T horizontal scanner.Another major source of distortion at high fields arises from susceptibility variations. To reduce such artifacts, we used segmented gradient-echo echo-planar imaging (EPI) sequences. Increasing the number of segments significantly decreased susceptibility artifacts and image distortion. Comparisons of images from functional runs using four segments with those using a single-shot EPI sequence revealed a roughly twofold improvement in functional signal-to-noise-ratio and 50% decrease in distortion. These methods enabled reliable detection of neural activation and permitted blood-oxygenation-level-dependent-based mapping of early visual areas in monkeys using a volume coil.In summary, both extensive behavioral training of monkeys and application of segmented gradient-echo EPI sequence improved signal-to-noise ratio and image quality. Understanding the effects these factors have is important for the application of high field imaging methods to the detection of submillimeter functional structures in the awake monkey brain.     

High-resolution segmented EPI in a motor task fMRI study

A high-resolution gradient echo, multi-slice segmented echo planar imaging method was used for functional MRI (fMRI) using a motor task at 1.5 Tesla. Functional images with an in-plane resolution of 1 mm and slice thickness of 4 mm were obtained with good white-gray matter contrast. The multi-shot approach, combined with a short total readout period of 82 ms, limits blurring effects for short T(2)(*) tissues (such as gray matter), assuring truly high-resolution images. In all subjects, motor functions were clearly depicted in the contralateral central sulcus over several slices and sometimes activation was detected in the supplementary motor area and/or ipsilateral central sulcus. The average signal change of 11+/-3% was much higher than in standard low-resolution fMRI EPI experiments, as a result of larger relative blood fractions.     

K-space trajectory mapping and its application for ultrashort Echo time imaging

MR images are affected by system delays and gradient field imperfections which induce discrepancies between prescribed and actual k-space trajectories. This could be even more critical for non-Cartesian data acquisitions where even a small deviation from the assumed k-space trajectory results in severe image degradation and artifacts. Knowledge of the actual k-space trajectories is therefore crucial and can be incorporated in the reconstruction of high quality non-Cartesian images. A novel MR method for the calibration of actual gradient waveforms was developed using a combination of phase encoding increments and subsequent detection of the exact time point at which the corresponding trajectory is crossing the k-space origin. The measured sets of points were fitted to a parametrical model to calculate the complete actual acquisition trajectory. Measurements performed on phantoms and volunteers, positioned both in- and off-isocenter of the magnet, clearly demonstrate the improvement in reconstructed ultrashort echo time (UTE) images, when information from calibration of k-space sampling trajectories is employed in the MR image reconstruction procedure. The unique feature of the proposed method is its robustness and simple experimental setup, making it suitable for quick acquisition trajectory calibration procedures e.g. for non-Cartesian radial fast imaging.     

k-space inherited parallel acquisition (KIPA): application on dynamic magnetic resonance imaging thermometry

In this study, a novel method for dynamic parallel image acquisition and reconstruction is presented. In this method, called k-space inherited parallel acquisition (KIPA), localized reconstruction coefficients are used to achieve higher reduction factors, and lower noise and artifact levels compared to that of generalized autocalibrating partially parallel acquisition (GRAPPA) reconstruction. In KIPA, the full k-space for the first frame and the partial k-space for later frames are required to reconstruct a whole series of images. Reconstruction coefficients calculated for different segments of k-space from the first frame data set are used to estimate missing k-space lines in corresponding k-space segments of other frames. The local determination of KIPA reconstruction coefficients is essential to adjusting them according to the local signal-to-noise ratio characteristics of k-space data. The proposed algorithm is applicable to dynamic imaging with arbitrary k-space sampling trajectories. Simulations of magnetic resonance thermometry using the KIPA method with a reduction factor of 6 and using dynamic imaging studies of human subjects with reduction factors of 4 and 6 have been performed to prove the feasibility of our method and to show apparent improvement in image quality in comparison with GRAPPA for dynamic imaging.     

Improved B0 field map estimation for high field EPI

Echo planar imaging (EPI) is an ultrafast magnetic resonance imaging (MRI) technique that allows one to acquire a 2D image in about 100 ms. Unfortunately, the standard EPI images suffer from substantial geometric distortions, mainly originating from susceptibility differences in adjacent tissues. To reduce EPI distortions, correction methods based on a field map, which is a map of the off-resonance frequencies, have been developed. In this work, a nonlinear least squares estimator is used to optimize the estimation of the field map of the B₀ field. The model of the EPI and reference data includes parameters for the phase evolution, the complex magnitude, the relaxation of the MRI signal and the EPI-specific phase difference between odd and even echoes, and from these parameters, additional corrections might be computed. The reference data required to estimate the field map can be acquired with a modified EPI-sequence. The proposed method is tested on simulated as well as experimental data and proves to be significantly more robust against noise, compared to the previously suggested method.     

Single-shot fast spin-echo diffusion tensor imaging of the brain and spine with head and phased array coils at 1.5 T and 3.0 T

In this study, we investigated the use of a single-shot fast spin-echo-based sequence to perform diffusion tensor imaging (DTI) with improved anatomic fidelity through the entire brain and the cervical spine. Traditionally, diffusion tensor images have been acquired by single-shot echo-planar imaging (EPI) methods in which large distortions result from magnetic susceptibility effects, especially near air-tissue interfaces. These distortions can be problematic, especially in anterior and inferior portions of the brain, and they also can severely limit applications in the spine. At higher magnetic fields these magnetic susceptibility artifacts are increased. The single-shot fast spin-echo (SSFSE) method used in this study utilizes radiofrequency rephasing in the transverse plane and thus provides diffusion images with negligible distortion even at 3 Tesla. In addition, the SSFSE sequence does not require multiple fast-receivers, which are not available on many magnetic resonance (MR) systems. Phased array coils were used to increase the signal-to-noise ratio of the images, offering a major inherent advantage in diffusion tensor imaging of the spine and brain. The mean diffusion measurements obtained with the SSFSE acquisition were not statistically different (p > 0.05) from EPI-based acquisitions. Compared to routine T(2)-weighted MR images, the DTI-EPI sequence showed up to 20% in elongation of the brain in the anterior-posterior direction on a sagittal image due to magnetic susceptibility distortions, whereas in the DTI-SSFSE, the image distortions were negligible. The diffusion tensor SSFSE method was also able to assess diffusion abnormalities in a brain stem hemorrhage, unaffected by the spatial distortions that limited conventional EPI acquisition.     

Gradient induced artifacts in simultaneous EEG-fMRI: Effect of synchronization on spiral and EPI k-space trajectories

The nature of the gradient induced electroencephalography (EEG) artifact is analyzed and compared for two functional magnetic resonance imaging (fMRI) pulse sequences with different k-space trajectories: echo planar imaging (EPI) and spiral. Furthermore, the performance of the average artifact subtraction algorithm (AAS) to remove the gradient artifact for both sequences is evaluated. The results show that the EEG gradient artifact for spiral sequences is one order of magnitude higher than for EPI sequences due to the chirping spectrum of the spiral sequence and the dB/dt of its crusher gradients. However, in the presence of accurate synchronization, the use of AAS yields the same artifact suppression efficiency for both pulse sequences below 80 Hz. The quality of EEG signal after AAS is demonstrated for phantom and human data. EEG spectrogram and visual evoked potential (VEP) are compared outside the scanner and use both EPI and spiral pulse sequences. MR related artifact residues affect the spectra over 40 Hz (less than 0.2 μV up to 120 Hz) and modify the amplitude of P1, N2 and P300 in the VEP. These modifications in the EEG signal have to be taken into account when interpreting EEG data acquired in simultaneous EEG-fMRI experiments.     

Quantitative analysis of spatial distortions of diffusion techniques at 3T

Diffusion has been widely adopted in the clinical setting to study the microstructural tissue changes in conjunction with anatomic imaging and metabolic imaging to offer insights on the status of the tissue injury or lesion. However, geometric distortions caused by magnetic susceptibility effects, eddy currents and gradient imperfections greatly affect the clinical utility of the diffusion images. Several diffusion methods have been proposed in the recent years to obtain diffusion parameters with increased accuracy. In most cases, the comparisons to the clinical standard echo-planar imaging (EPI) diffusion are done visually without quantitative measurements. In this study, we present three simple, complementary quantitative methods of nonrigid image registration and shape analyses for evaluating spatial distortions on magnetic resonance images with application in comparing single-shot fast spin-echo (SSFSE) and EPI based diffusion measurements. These methods have confirmed the SSFSE based diffusion method is less distorted than the EPI based one, which is generally accepted through visual inspection.     

Optimized distortion correction of epi-based statistical parametrical maps for stereotactic neurosurgery

A new approach is presented for unwarping geometrical distortions in echo planar imaging (EPI)-based statistical parametrical maps (SPM). With the developed method, a substantial reduction of geometrical distortions in the order of 5-15 mm down to a residue of approximately 1-2 mm was possible, thus allowing reliable activation localization and, hence, interpretation after fusion with undistorted conventional anatomic images. The unwarping procedure is an optimized combination of two distortion correction methods, the "inhomogeneity mapping" and the coregistration approach. Local and global distortions, the first one object-dependent, the second one caused by imperfect gradient waveforms, are reduced. With several examples of patients and volunteers, the applicability for stereotactic neurosurgery is demonstrated.     

Echo-shifted multislice EPI for high-speed fMRI

The advantages of event-related functional Magnetic Resonance Imaging (fMRI) and the increasing use of fMRI in cognitive experiments are both driving the development of techniques that allow images sensitive to the blood oxygen level-dependent effect to be acquired at ever-higher temporal resolution. Here, we present a technique based on the use of echo shifting (ES) in conjunction with a multislice (MS) echo planar imaging (EPI) readout, which allows T2*-weighted images to be generated with a repetition time per slice that is less than the echo time (TE). Using this ES-MS-EPI approach, it is shown that images with a TE of 40 ms can be acquired with an acquisition time per slice of only 27 ms. The utility of the MS-ES-EPI sequence is demonstrated in a visual-motor, event-related fMRI study in which nine-slice image volumes are acquired continuously at a rate of 4.1 Hz. The sequence is shown to produce reliable activation associated with both visual stimuli and motor actions.     

Reduced distortion artifact whole brain CBF mapping using blip-reversed non-segmented 3D echo planar imaging with pseudo-continuous arterial spin labeling

PurposeTo implement and evaluate interleaved blip-up, blip-down, non-segmented 3D echo planar imaging (EPI) with pseudo-continuous arterial spin labeling (pCASL) and post-processing for reduced susceptibility artifact cerebral blood flow (CBF) maps.Materials and methods3D EPI non-segmented acquisition with a pCASL labeling sequence was modified to include alternating k-space coverage along phase encoding direction (referred to as “blip-reversed”) for alternating dynamic acquisitions of control and label pairs. Eight volunteers were imaged on a 3T scanner. Images were corrected for distortion using spatial shifting transformation of the underlying field map. CBF maps were calculated and compared with maps obtained without blip reversal using matching gray matter (GM) images from a high resolution 3D scan. Additional benefit of using the correction for alternating blip-up and blip-down acquisitions was assessed by comparing to corrected blip-up only and corrected blip-down only CBF maps. Matched Student t-test of overlapping voxels for the eight volunteers was done to ascertain statistical improvement in distortion.ResultsMean CBF value in GM for the eight volunteers from distortion corrected CBF maps was 50.8 ± 9.9 ml/min/100 gm tissue. Corrected CBF maps had 6.3% and 4.1% more voxels in GM when compared with uncorrected blip up (BU) and blip down (BD) images, respectively. Student t-test showed significant reduction in distortion when compared with blip-up images and blip-down images (p < 0.001). When compared with corrected BU and corrected BD only CBF maps, BU and BD corrected maps had 2.3% and 1% more voxels (p = 0.006 and 0.04, respectively).ConclusionPseudo-continuous arterial spin labeling with non-segmented 3D EPI acquisition using alternating blip-reversed k-space traversal and distortion correction provided significantly better matching GM CBF maps. In addition, employing alternating blip-reversed acquisitions during pCASL acquisition resulted in statistically significant improvement over corrected blip-up and blip-down CBF maps.     

An unsupervised deep learning technique for susceptibility artifact correction in reversed phase-encoding EPI images

Echo planar imaging (EPI) is a fast and non-invasive magnetic resonance imaging technique that supports data acquisition at high spatial and temporal resolutions. However, susceptibility artifacts, which cause the misalignment to the underlying structural image, are unavoidable distortions in EPI. Traditional susceptibility artifact correction (SAC) methods estimate the displacement field by optimizing an objective function that involves one or more pairs of reversed phase-encoding (PE) images. The estimated displacement field is then used to unwarp the distorted images and produce the corrected images. Since this conventional approach is time-consuming, we propose an end-to-end deep learning technique, named S-Net, to correct the susceptibility artifacts the reversed-PE image pair. The proposed S-Net consists of two components: (i) a convolutional neural network to map a reversed-PE image pair to the displacement field; and (ii) a spatial transform unit to unwarp the input images and produce the corrected images. The S-Net is trained using a set of reversed-PE image pairs and an unsupervised loss function, without ground-truth data. For a new image pair of reversed-PE images, the displacement field and corrected images are obtained simultaneously by evaluating the trained S-Net directly. Evaluations on three different datasets demonstrate that S-Net can correct the susceptibility artifacts in the reversed-PE images. Compared with two state-of-the-art SAC methods (TOPUP and TISAC), the proposed S-Net runs significantly faster: 20 times faster than TISAC and 369 times faster than TOPUP, while achieving a similar correction accuracy. Consequently, S-Net accelerates the medical image processing pipelines and makes the real-time correction for MRI scanners feasible. Our proposed technique also opens up a new direction in learning-based SAC.     

Functional neuroimaging of inner fields-of-view with 2D-selective RF excitations

Echo-planar imaging is widely used in functional neuroimaging but suffers from its pronounced sensitivity to field inhomogeneities that cause geometric distortions and image blurring which both limit the effective in-plane resolution achievable. In this work, it is shown how inner-field-of-view techniques based on 2D-selective RF excitations (2DRF) can be applied to reduce the field-of-view in the phase-encoding direction without aliasing and increase the in-plane resolution accordingly. Free-induction-decay (FID) EPI and echo-train-shifted (T2*-weighted) and standard (T2-weighted) spin-echo (SE) EPI with in-plane resolutions of up to 0.5×1.0 mm² (slice thickness 5 mm) were acquired at 3 T. Unwanted signal contributions of 2DRF side excitations were shifted out of the object (FID-EPI) or of the refocusing plane by tilting the excitation plane (SE-EPI). Brain activation in healthy volunteers was investigated with checkerboard and finger-tapping block-design paradigms. Brain activation could be detected with all sequences and contrasts, most reliably with FID-EPI due to its higher signal amplitude and the longer 2DRF excitation that are more sensitive to magnetic field inhomogeneities. In conclusion, inner-FOV EPI based on 2DRF excitations could help to improve the spatial resolution of fMRI of focal target regions, e.g. for applications in the spinal cord.     

Motion artifact reduction technique for dual-contrast FSE imaging

There is considerable similarity between proton density-weighted (PDw) and T2-weighted (T2w) images acquired by dual-contrast fast spin-echo (FSE) sequences. The similarity manifests itself in image space as consistency between the phases of PDw and T2w images and in k-space as correspondence between PDw and T2w k-space data. A method for motion artifact reduction for dual-contrast FSE imaging has been developed. The method uses projection onto convex sets (POCS) formalism and is based on image space phase consistency and the k-space similarity between PDw and T2w images. When coupled with a modified dual-contrast FSE phase encoding scheme the method can yield considerable artifact reduction, as long as less than half of the acquired data is corrupted by motion. The feasibility and efficiency of the developed method were demonstrated using phantom and human MRI data.     

Reduced aliasing artifacts using shaking projection k-space sampling trajectory

Radial imaging techniques, such as projection-reconstruction （PR）, are used in magnetic resonance imaging （MRI） for dynamic imaging, angiography, and short-T2 imaging. They are less sensitive to flow and motion artifacts, and support fast imaging with short echo times. However, aliasing and streaking artifacts are two main sources which degrade radial imaging quality. For a given fixed number of k-space projections, data distributions along radial and angular directions will influence the level of aliasing and streaking artifacts. Conventional radial k-space sampling trajectory introduces an aliasing artifact at the first principal ring of point spread function （PSF）. In this paper, a shaking projection （SP） k-space sampling trajectory was proposed to reduce aliasing artifacts in MR images. SP sampling trajectory shifts the projection alternately along the k-space center, which separates k-space data in the azimuthal direction. Simulations based on conventional and SP sampling trajectories were compared with the same number projections. A significant reduction of aliasing artifacts was observed using the SP sampling trajectory. These two trajectories were also compared with different sampling frequencies. ASP trajectory has the same aliasing character when using half sampling frequency （or half data） for reconstruction. SNR comparisons with different white noise levels show that these two trajectories have the same SNR character. In conclusion, the SP trajectory can reduce the aliasing artifact without decreasing SNR and also provide a way for undersampling recon- struction. Furthermore, this method can be applied to three-dimensional （3D） hybrid or spherical radial k-space sampling for a more efficient reduction of aliasing artifacts.