A method for the dynamic correction of B0-related distortions in single-echo EPI at 7 T

We propose a method to calculate field maps from the phase of each EPI in an fMRI time series. These field maps can be used to correct the corresponding magnitude images for distortion caused by inhomogeneity in the static magnetic field. In contrast to conventional static distortion correction, in which one 'snapshot’ field map is applied to all subsequent fMRI time points, our method also captures dynamic changes to B0which arise due to motion and respiration. The approach is based on the assumption that the non-B0-related contribution to the phase measured by each radio-frequency coil, which is dominated by the coil sensitivity, is stable over time and can therefore be removed to yield a field map from EPI. Our solution addresses imaging with multi-channel coils at ultra-high field (7 T), where phase offsets vary rapidly in space, phase processing is non-trivial and distortions are comparatively large. We propose using dual-echo gradient echo reference scan for the phase offset calculation, which yields estimates with high signal-to-noise ratio. An extrapolation method is proposed which yields reliable estimates for phase offsets even where motion is large and a tailored phase unwrapping procedure for EPI is suggested which gives robust results in regions with disconnected tissue or strong signal decay. Phase offsets are shown to be stable during long measurements (40 min) and for large head motions. The dynamic distortion correction proposed here is found to work accurately in the presence of large motion (up to 8.1°), whereas a conventional method based on single field map fails to correct or even introduces distortions (up to 11.2 mm). Finally, we show that dynamic unwarping increases the temporal stability of EPI in the presence of motion. Our approach can be applied to any EPI measurements without the need for sequence modification

Deconvolution‐based distortion correction of EPI using analytic single‐voxel point‐spread functions

Purpose To develop a postprocessing algorithm that corrects geometric distortions due to spatial variations of the static magnetic field amplitude, B0, and effects from relaxation during signal acquisition in EPI. Theory and Methods An analytic, complex point‐spread function is deduced for k‐space trajectories of EPI variants and applied to corresponding acquisitions in a resolution phantom and in human volunteers at 3 T. With the analytic point‐spread function and experimental maps of B0 (and, optionally, the effective transverse relaxation time, urn:x-wiley:07403194:media:mrm28591:mrm28591-math-0004) as input, a point‐spread function matrix operator is devised for distortion correction by a Thikonov‐regularized deconvolution in image space. The point‐spread function operator provides additional information for an appropriate correction of the signal intensity distribution. A previous image combination algorithm for acquisitions with opposite phase blip polarities is adapted to the proposed method to recover destructively interfering signal contributions. Results Applications of the proposed deconvolution‐based distortion correction (“DecoDisCo”) algorithm demonstrate excellent distortion corrections and superior performance regarding the recovery of an undistorted intensity distribution in comparison to a multifrequency reconstruction. Examples include full and partial Fourier standard EPI scans as well as double‐shot center‐out trajectories. Compared with other distortion‐correction approaches, DecoDisCo permits additional deblurring to obtain sharper images in cases of significant urn:x-wiley:07403194:media:mrm28591:mrm28591-math-0005 effects. Conclusion Robust distortion corrections in EPI acquisitions are feasible with high quality by regularized deconvolution with an analytic point‐spread function. The general algorithm, which is publicly released on GitHub, can be straightforwardly adapted for specific EPI variants or other acquisition schemes

Joint B0 and image estimation integrated with model based reconstruction for field map update and distortion correction in prostate diffusion MRI

In prostate Diffusion Weighted MRI, differences in susceptibility values exist at the interface between the prostate and rectal-air. This can result in off-resonance magnetic field leading to geometric distortions including signal stretching and signal pile-up in the reconstructed images. Using a set of EPI data acquired with blip-up and blip-down phase encoding gradient directions, model based reconstruction has recently been proposed that can correct these distortions by using a B0 field estimated from a separate B0 scan. However, change in the size of the rectal air region across time can occur that can result in a mismatch of the B0 field to the EPI scan. Also, the measured B0 field itself can be erroneous in regions of low Signal to Noise ratio around the prostate rectal air interface. In this work, using a set of single shot EPI data acquired with blip-up and blip-down phase encoding gradient directions, a novel joint model based reconstruction is proposed that can account for changes in the off resonance effects between the B0 and EPI scans. For ten prostate patients, using a measured B0 field as an initial B0 estimate, on a 5-point scale (1-5) image quality scores evaluated by an experienced radiologist, the proposed framework achieved scores of 3.50+/-0.85 and 3.40+/-0.51 for bvalues of 0 and 500 s/mm2, respectively compared to 3.40+/-0.70 and 3.30+/-0.67 for model based reconstruction. The proposed framework is also capable of estimating a distortion corrected EPI image even without an initial B0 field estimate in situations where a separate B0 scan cannot be obtained due to time constraint

Simultaneous multislice acquisition with multi-contrast segmented EPI for separation of signal contributions in dynamic contrast-enhanced imaging

We present a method to efficiently separate signal in magnetic resonance imaging (MRI) into a base signal S0, representing the mainly T1-weighted component without T2*-relaxation, and its T2*-weighted counterpart by the rapid acquisition of multiple contrasts for advanced pharmacokinetic modelling. This is achieved by incorporating simultaneous multislice (SMS) imaging into a multi-contrast, segmented echo planar imaging (EPI) sequence to allow extended spatial coverage, which covers larger body regions without time penalty. Simultaneous acquisition of four slices was combined with segmented EPI for fast imaging with three gradient echo times in a preclinical perfusion study. Six female domestic pigs, German-landrace or hybrid-form, were scanned for 11 minutes respectively during administration of gadolinium-based contrast agent. Influences of reconstruction methods and training data were investigated. The separation into T1- and T2*-dependent signal contributions was achieved by fitting a standard analytical model to the acquired multi-echo data. The application of SMS yielded sufficient temporal resolution for the detection of the arterial input function in major vessels, while anatomical coverage allowed perfusion analysis of muscle tissue. The separation of the MR signal into T1- and T2*-dependent components allowed the correction of susceptibility related changes. We demonstrate a novel sequence for dynamic contrast-enhanced MRI that meets the requirements of temporal resolution (Δt < 1.5 s) and image quality. The incorporation of SMS into multi-contrast, segmented EPI can overcome existing limitations of dynamic contrast enhancement and dynamic susceptibility contrast methods, when applied separately. The new approach allows both techniques to be combined in a single acquisition with a large spatial coverage

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

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

Impact of susceptibility-induced distortion correction on perfusion imaging by pCASL with a segmented 3D GRASE readout

Purpose: The consensus for the clinical implementation of arterial spin labeling (ASL) perfusion imaging recommends a segmented 3D Gradient and Spin-Echo (GRASE) readout for optimal signal-to-noise-ratio (SNR). The correction of the associated susceptibility-induced geometric distortions has been shown to improve diagnostic precision, but its impact on ASL data has not been systematically assessed and it is not consistently part of pre-processing pipelines. Here, we investigate the effects of susceptibility-induced distortion correction on perfusion imaging by pseudo-continuous ASL (pCASL) with a segmented 3D GRASE readout. Methods: Data acquired from 28 women using pCASL with 3D GRASE at 3T was analyzed using three pre-processing options: without distortion correction, with distortion correction, and with spatial smoothing (without distortion correction) matched to control for blurring effects induced by distortion correction. Maps of temporal SNR (tSNR) and relative perfusion were analyzed in eight regions-of-interest (ROIs) across the brain. Results: Distortion correction significantly affected tSNR and relative perfusion across the brain. Increases in tSNR were like those produced by matched spatial smoothing in most ROIs, indicating that they were likely due to blurring effects. However, that was not the case in the frontal and temporal lobes, where we also found increased relative perfusion with distortion correction even compared with matched spatial smoothing. These effects were found in both controls and patients, with no interactions with the participant group. Conclusion: Correction of susceptibility-induced distortions significantly impacts ASL perfusion imaging using a segmented 3D GRASE readout, and this step should therefore be considered in ASL pre-processing pipelines. This is of special importance in clinical studies, reporting perfusion across ROIs defined on relatively undistorted images and when conducting group analyses requiring the alignment of images across different subjects.info:eu-repo/semantics/publishedVersio

Interleaved EPI based fMRI improved by multiplexed sensitivity encoding (MUSE) and simultaneous multi-band imaging

© 2014 Chang et al. Functional magnetic resonance imaging (fMRI) is a non-invasive and powerful imaging tool for detecting brain activities. The majority of fMRI studies are performed with single-shot echo-planar imaging (EPI) due to its high temporal resolution. Recent studies have demonstrated that, by increasing the spatial-resolution of fMRI, previously unidentified neuronal networks can be measured. However, it is challenging to improve the spatial resolution of conventional single-shot EPI based fMRI. Although multi-shot interleaved EPI is superior to single-shot EPI in terms of the improved spatial-resolution, reduced geometric distortions, and sharper point spread function (PSF), interleaved EPI based fMRI has two main limitations: 1) the imaging throughput is lower in interleaved EPI; 2) the magnitude and phase signal variations among EPI segments (due to physiological noise, subject motion, and B0 drift) are translated to significant in-plane aliasing artifact across the field of view (FOV). Here we report a method that integrates multiple approaches to address the technical limitations of interleaved EPI-based fMRI. Firstly, the multiplexed sensitivity-encoding (MUSE) post-processing algorithm is used to suppress in-plane aliasing artifacts resulting from time-domain signal instabilities during dynamic scans. Secondly, a simultaneous multi-band interleaved EPI pulse sequence, with a controlled aliasing scheme incorporated, is implemented to increase the imaging throughput. Thirdly, the MUSE algorithm is then generalized to accommodate fMRI data obtained with our multi-band interleaved EPI pulse sequence, suppressing both in-plane and through-plane aliasing artifacts. The blood-oxygenation-level-dependent (BOLD) signal detectability and the scan throughput can be significantly improved for interleaved EPI-based fMRI. Our human fMRI data obtained from 3 Tesla systems demonstrate the effectiveness of the developed methods. It is expected that future fMRI studies requiring high spatial-resolvability and fidelity will largely benefit from the reported techniques.published_or_final_versio