Correction of eddy current distortions in high angular resolution diffusion imaging

Purpose:

To correct distortions caused by eddy currents induced by large diffusion gradients during high angular resolution diffusion imaging without any auxiliary reference scans.

Materials and Methods:

Image distortion parameters were obtained by image coregistration, performed only between diffusion‐weighted images with close diffusion gradient orientations. A linear model that describes distortion parameters (translation, scale, and shear) as a function of diffusion gradient directions was numerically computed to allow individualized distortion correction for every diffusion‐weighted image.

Results:

The assumptions of the algorithm were successfully verified in a series of experiments on phantom and human scans. Application of the proposed algorithm in high angular resolution diffusion images markedly reduced eddy current distortions when compared to results obtained with previously published methods.

Conclusion:

The method can correct eddy current artifacts in the high angular resolution diffusion images, and it avoids the problematic procedure of cross‐correlating images with significantly different contrasts resulting from very different gradient orientations or strengths. J. Magn. Reson. Imaging 2013;37:1460–1467. © 2012 Wiley Periodicals, Inc.     

Eddy current correction in diffusion-weighted imaging using pairs of images acquired with opposite diffusion gradient polarity.

In echo-planar-based diffusion-weighted imaging (DWI) and diffusion tensor imaging (DTI), the evaluation of diffusion parameters such as apparent diffusion coefficients and anisotropy indices is affected by image distortions that arise from residual eddy currents produced by the diffusion-sensitizing gradients. Correction methods that coregister diffusion-weighted and non-diffusion-weighted images suffer from the different contrast properties inherent in these image types. Here, a postprocessing correction scheme is introduced that makes use of the inverse characteristics of distortions generated by gradients with reversed polarity. In this approach, only diffusion-weighted images with identical contrast are included for correction. That is, non-diffusion-weighted images are not needed as a reference for registration. Furthermore, the acquisition of an additional dataset with moderate diffusion-weighting as suggested by Haselgrove and Moore (Magn Reson Med 1996;36:960-964) is not required. With phantom data it is shown that the theoretically expected symmetry of distortions is preserved in the images to a very high degree, demonstrating the practicality of the new method. Results from human brain images are also presented.     

k-space correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging.

This paper describes a method for correcting eddy-current (EC)-induced distortions in diffusion-weighted echo-planar imaging (DW-EPI). First, reference measurements of EC fields within the EPI acquisition window are performed for DW gradient pulses applied separately along each physical axis of the gradient set and for a range of gradient amplitudes. EC fields caused by the DW gradients of the DW-MRI protocol are then calculated using the reference EC measurements. Finally, these calculated fields are used to correct the respective DW-EPI raw (k-space) data during image reconstruction. The technique was implemented in a small-bore MRI scanner with no digital preemphasis. It corrected EC-induced image distortions in both phantom and in vivo brain diffusion tensor imaging (DTI) data more effectively than commonly used image-based techniques. The method did not increase imaging time, since the same reference EC measurements were used to correct data acquired from different phantoms, subjects, and DTI protocols. Because of the simplicity of the reference EC measurements, the method can easily be implemented in clinical scanners.     

Correction of high-order eddy current induced geometric distortion in diffusion-weighted echo-planar images.

Diffusion-weighted images acquired with the echo-planar imaging technique are highly sensitive to eddy current induced geometric distortions that vary with the magnitude and direction of the diffusion sensitizing gradients. Such distortions cause misalignment of images acquired with different diffusion strengths and orientations. This in turn can result in errors when calculating maps of the apparent diffusion coefficient and diffusion tensor. Previous correction methods either require separate calibration data or only deal with low-order errors. In this study, we demonstrate a method that can correct for higher-order errors. The method relies on collecting pairs of images with diffusion sensitizing gradients reversed. This paired data are first corrected for shifts and linear distortion and then combined to cancel higher-order errors. All acquired data contribute to the final results. The method has been tested by simulation, on phantoms, on adult volunteers, and on neonatal brain examinations.     

Echoplanar diffusion tensor imaging of the lower leg musculature using eddy current nulled stimulated echo preparation.

A sequence for echoplanar diffusion tensor imaging of musculature was developed using a stimulated echo preparation. The strategy was optimized in order to obtain reliable diffusion tensor data in a short measuring time. Image distortion problems due to eddy currents arising from long-lasting diffusion sensitizing gradients could be overcome by insertion of additional gradient pulses in the TM interval of the stimulated echo preparation. In contrast to former approaches with similar intention, the proposed strategy does not influence the stimulated echo signal itself and does not lead to prolonged echo time as in the case of spin echo methods. Phantom measurements were performed to compare eddy current induced distortion effects in diffusion weighted images. The diffusion tensor in the musculature of the lower leg was investigated in four healthy subjects and maps of the trace and the three eigenvalues of the diffusion tensor, fractional anisotropy maps, and angle maps were calculated.     

Assessing atrophy of the major white matter fiber bundles of the brain from diffusion tensor MRI data.

Brain atrophy is a typical feature of many neurological conditions. Therefore, quantitative evaluation and spatial characterization of atrophy are potentially useful for monitoring the evolution of central nervous system (CNS) disorders. In this study, a method for measuring atrophy of the major white matter (WM) fiber bundles in the brain using diffusion tensor (DT) MRI data is developed. To this end, an atlas was created from sets of diffusion anisotropy images from normal subjects, and the deformations necessary to match single subject anisotropy images to this atlas were then computed. Because diffusion anisotropy images were used, this approach should be sensitive to fiber bundle volume changes in the same way that using T1-weighted images allows gray matter volume changes to be measured. The Jacobian determinant of the deformation field for each subject was then used as a measure of contraction or expansion of the tissue at each image voxel. An overview of the nonlinear registration problem is given; then an optimization of the parameters for the chosen algorithm is performed and the method for producing the atlas is described. The effectiveness of the method was then tested on data from five patients with multiple sclerosis (MS) and two patients with amyotrophic lateral sclerosis (ALS).     

Stimulated echo induced misestimates on diffusion tensor indices and its remedy

Purpose:

To report possible erroneous estimates of diffusion parameters in the twice‐refocused spin‐echo (TRSE) technique, proposed to eliminate eddy‐current‐induced geometric distortions in diffusion‐weighted echo‐planar imaging, when stimulated echo signals are inappropriately included.

Materials and Methods:

Eleven subjects were included for imaging experiments on two 1.5 Tesla systems using the TRSE sequence. Three versions, two with unbalanced crusher gradients inserted to dephase the stimulated echo from the b = 0 images and one with balanced crusher gradients, were implemented. The apparent diffusion coefficients (ADC) and fractional anisotropy (FA) were derived and compared.

Results:

The ADCs obtained with unbalanced crusher gradients were closer to values reported in the literature. Stimulated echo led to ADC over‐estimations by 34.2%, 50.4%, 54.0%, 51.5%, 24.0%, and 41.9% in the genu of corpus callosum, splenium of corpus callosum, bilateral corona radiata, internal capsule, mediofrontal gyrus, and the cuneus, respectively (P < 0.01), with concomitant reduction in FA in highly anisotropic regions. Over‐estimations of diffusion coefficients were found to be roughly equal along all directions.

Conclusion:

Formation of stimulated echo in the TRSE technique can lead to erroneous estimations of the diffusion parameters, even if no prominent morphological artifacts are seen. J. Magn. Reson. Imaging 2010;31:1522–1529. © 2010 Wiley‐Liss, Inc.     

Slab scan diffusion imaging.

For maximum robustness of a diffusion-weighted MR imaging sequence, it is desirable to use a single-shot imaging method. This article introduces a new single-shot imaging approach that combines the advantages of multiple spin-echoes with the technique of line scan diffusion imaging. A slab volume, which can be spatially encoded with fewer phase encodes than a regular field of view, is selected with 2D selective pulses. With the shorter echo train, the sensitivity to field inhomogeneities and chemical shift is thus greatly diminished. Further reduction is achieved by interleaving short gradient echo trains with refocusing spin-echo pulses. Optimized slice-selective RF pulses that produce flip angles close to 180 degrees are used to minimize the stimulated echo component. Motion-related phase shifts, which change polarity with each spin-echo excitation, will give rise to artifacts that are avoidable by processing even and odd spin-echoes separately. As with line scan diffusion imaging, the complete field of view is acquired by sequential scanning. Since with each shot several lines of data are collected, a considerable improvement over line scan diffusion imaging in terms of scanning speed is achieved. Diffusion data obtained in phantoms and normal subjects demonstrate the feasibility of this novel approach.     

PROPELLER EPI: an MRI technique suitable for diffusion tensor imaging at high field strength with reduced geometric distortions.

A technique suitable for diffusion tensor imaging (DTI) at high field strengths is presented in this work. The method is based on a periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) k-space trajectory using EPI as the signal readout module, and hence is dubbed PROPELLER EPI. The implementation of PROPELLER EPI included a series of correction schemes to reduce possible errors associated with the intrinsically higher sensitivity of EPI to off-resonance effects. Experimental results on a 3.0 Tesla MR system showed that the PROPELLER EPI images exhibit substantially reduced geometric distortions compared with single-shot EPI, at a much lower RF specific absorption rate (SAR) than the original version of the PROPELLER fast spin-echo (FSE) technique. For DTI, the self-navigated phase-correction capability of the PROPELLER EPI sequence was shown to be effective for in vivo imaging. A higher signal-to-noise ratio (SNR) compared to single-shot EPI at an identical total scan time was achieved, which is advantageous for routine DTI applications in clinical practice.     

Unifying the analyses of anatomical and diffusion tensor images using volume-preserved warping

PURPOSE: To introduce a framework that automatically identifies regions of anatomical abnormality within anatomical MR images and uses those regions in hypothesis-driven selection of seed points for fiber tracking with diffusion tensor (DT) imaging (DTI). MATERIALS AND METHODS: Regions of interest (ROIs) are first extracted from MR images using an automated algorithm for volume-preserved warping (VPW) that identifies localized volumetric differences across groups. ROIs then serve as seed points for fiber tracking in coregistered DT images. Another algorithm automatically clusters and compares morphologies of detected fiber bundles. We tested our framework using datasets from a group of patients with Tourette's syndrome (TS) and normal controls. RESULTS: Our framework automatically identified regions of localized volumetric differences across groups and then used those regions as seed points for fiber tracking. In our applied example, a comparison of fiber tracts in the two diagnostic groups showed that most fiber tracts failed to correspond across groups, suggesting that anatomical connectivity was severely disrupted in fiber bundles leading from regions of known anatomical abnormality. CONCLUSION: Our framework automatically detects volumetric abnormalities in anatomical MRIs to aid in generating a priori hypotheses concerning anatomical connectivity that then can be tested using DTI. Additionally, automation enhances the reliability of ROIs, fiber tracking, and fiber clustering.     

Diffusion tensor imaging (DTI) of rodent brains in vivo using a 1.5T clinical MR scanner

PURPOSE: To evaluate the feasibility of using a clinical 1.5T MR scanner to perform magnetic resonance (MR) diffusion tensor imaging (DTI) on in vivo rodent brains and to trace major rodent neuronal bundles with anatomical correlation. MATERIALS AND METHODS: Two normal adult Sprague Dawley (SD) rats were anesthetized and imaged in a 1.5T MR scanner with a microscopic coil. DTI was performed at a resolution of 0.94 mm x 0.94 mm x 0.5 mm (reconstructed to 0.47 mm x 0.47 mm x 0.5 mm, with b-factors of 600 seconds/mm2 and 1000 seconds/mm2) and a higher resolution of 0.63 mm x 0.63 mm x 0.5 mm (reconstructed to 0.235 mm x 0.235 mm x 0.5 mm, with a b-factor of 1500 seconds/mm2). The fiber-tracking results were correlated with corresponding anatomical sections stained to visualize neuronal fibers. The apparent diffusion coefficient (ADC) and fractional anisotropy (FA) of the neuronal fibers were measured and compared with results in published reports. RESULTS: Several major neuronal fiber tracts, including the corticospinal cord, corpus callosum, and anterior commissure, were identified in all DTI data sets. Stained anatomical sections obtained from the rats confirmed the location of these fibers. The ADC values (0.6-0.8 +/- 10(-3) mm2/second) of the fibers were similar to published figures. However, the FA values (0.3-0.35) were lower than those obtained in previous studies of white matter in rodent spinal cord. CONCLUSION: We have demonstrated the feasibility of using a 1.5T clinical MR scanner for neuronal fiber tracking in rodent brains. The technique will be useful in rodent neuroanatomy studies. Further investigation is encouraged to verify the FA values generated by DTI with such techniques.     

DTI常用扫描序列原理及比较

磁共振弥散张量成像技术是利用水分子的弥散各向异性进行成像,可用于脑白质纤维研究,常用扫描技术包括单次激发平面回波成像(EPI),线阵扫描弥散成像,导航自旋回波弥散加权成像(LSDI),半傅立叶探测单发射快速自旋回波成像等。每种成像技术各有其优缺点,EPI扫描时间短,图像信噪比高,但存在化学位移伪影、磁敏感性伪影、几何变形;LSDI精确度高,几乎无伪影及变形,但扫描时间过长;导航自旋回波弥散加权成像运动伪影少,但扫描时间长;半傅立叶探测单发射快速自旋回波成像扫描时间短,但图像模糊。综合比较,单次激发平面回波成像是用于临床研究较适宜的方法。     

Contiguous‐slice zonally oblique multislice (CO‐ZOOM) diffusion tensor imaging: Examples of in vivo spinal cord and optic nerve applications

Purpose

To describe and demonstrate a new technique that allows diffusion tensor imaging of small structures such as the spinal cord (SC) and optic nerve (ON) with contiguous slices and reduced image distortions using a narrow field of view (FOV).

Materials and Methods

Images were acquired with a modified single‐shot echo‐planar imaging (EPI) sequence that contains a refocusing radio frequency (RF) pulse in the presence of the phase‐encoding (rather than slice‐select) gradient. As a result, only a narrow volume may be both excited and refocused, removing the problem of signal aliasing for narrow FOVs. Two variants of this technique were developed: cardiac gating is included in the study of the SC to reduce pulsation artifacts, whereas inversion‐recovery (IR) cerebrospinal fluid (CSF) suppression is utilized in the study of the ON to eliminate partial volume effects. The technique was evaluated with phantoms, and mean diffusivity (MD) and fractional anisotropy (FA) measurements were made in the SC and ON of two healthy volunteers.

Results

The technique provides contiguous‐slice, reduced‐FOV images that do not suffer from aliasing and have reduced magnetic susceptibility artifacts. MD and FA values determined here lie within the ranges quoted in the literature.

Conclusion

Contiguous‐slice zonally orthogonal multislice (CO‐ZOOM‐EPI is a new technique for diffusion‐weighted imaging of small structures such as the ON and SC with high resolution and reduced distortions due to susceptibility variations. This technique is able to acquire contiguous slices that may allow further nerve‐tracking analyses. J. Magn. Reson. Imaging 2009;29:454–460. © 2009 Wiley‐Liss, Inc.     

Correction for eddy current induced Bo shifts in diffusion‐weighted echo‐planar imaging

The importance of diffusion‐weighted MRI in the assessment of acute stroke is well‐recognized, and quantitative maps of the apparent diffusion coefficient (ADC) are now widely used. Echo‐planar imaging provides a robust method of acquiring diffusion‐weighted images free of motion artifact. However, initial experience with clinical MRI systems indicates that calculation of artifact‐free ADC maps from a series of echo‐planar diffusion‐weighted images is not necessarily straight‐forward. One of the problems is that frequency shifts resulting from eddy currents can cause misregistration of base diffusion‐weighted images. In this study, an on‐line correction method that overcomes this problem is described, and phantom and human images that demonstrate the validity of the technique are presented. The method uses a non‐phase‐encoded reference scan to correct the phase of each echo in the echo train, and can provide ADC maps that are free of misregistration artifacts, without the need for off‐line postprocessing. Magn Reson Med 41:95‐102, 1999. © 1999 Wiley‐Liss, Inc.     

Phantom-based evaluation of geometric distortions in functional magnetic resonance and diffusion tensor imaging.

Phantom-based evaluation of geometric distortions in functional MRI and diffusion tensor imaging (DTI) was investigated. An acrylic water-filled phantom with a grid structure was designed and manufactured to provide accurate geometric information over the volume measured in human brain imaging. The grid structures were well detected in data acquired using a 3-T MRI scanner with echo-planar imaging (EPI) sequences commonly applied in functional MRI and DTI. A method for quantifying distortions in the phantom data was presented and applied for the images. The validity of the phantom for EPI was evaluated by quantitatively comparing the distortions present in and induced by the phantom and a human brain when imaged under identical conditions. The results suggest that the new phantom can reveal geometric distortions easily undermined by standard MRI phantoms. For example, prominent variability in the distortions was found as a function of the orientation of the diffusion-sensitizing gradient. Possible future applications for this type of phantom include quality assurance and calibration of the hardware and software used in EPI-based functional MRI and DTI.     

Differentiation of fibroblastic meningiomas from other benign subtypes using diffusion tensor imaging

PURPOSE: To differentiate fibroblastic meningiomas, usually considered to be of a hard consistency, from other benign subtypes using diffusion tensor imaging (DTI). MATERIALS AND METHODS: From DTI data sets of 30 patients with benign meningiomas, we calculated diffusion tensors and mean diffusivity (MD) and fractional anisotropy (FA) maps as well as barycentric maps representing the geometrical shape of the tensors. The findings were compared to postoperative histology. The study was approved by the local ethics committee, and informed consent was given by the patients. RESULTS: According to one-way analysis of variance (ANOVA), FA was the best parameter to differentiate between the subtypes (F=32.2; p<0.0001). Regarding tensor shape, endothelial meningiomas were represented by spherical tensors (80%) corresponding to isotropic diffusion, whereas the fibroblastic meningiomas showed a high percentage (43%) of nonspherical tensors, indicating planar or longitudinal diffusion. The difference was highly significant (F=28.4; p<0.0001) and may be due to the fascicular arrangement of long spindle-shaped tumor cells and the high content of intra- and interfascicular fibers as shown in the histology. In addition, a capsule-like rim of the in-plane diffusion surrounded most meningiomas irrespective of their histological type. CONCLUSION: If these results correlate to the intraoperative findings of meningioma consistency, DTI-based measurement of FA and analysis of the shape of the diffusion tensor is a promising method to differentiate between fibroblastic and other subtypes of benign meningiomas in order to get information about their "hard" or "soft" consistency prior to removal.     

Generalized diffusion tensor imaging and analytical relationships between diffusion tensor imaging and high angular resolution diffusion imaging.

A new method for mapping diffusivity profiles in tissue is presented. The Bloch-Torrey equation is modified to include a diffusion term with an arbitrary rank Cartesian tensor. This equation is solved to give the expression for the generalized Stejskal-Tanner formula quantifying diffusive attenuation in complicated geometries. This makes it possible to calculate the components of higher-rank tensors without using the computationally-difficult spherical harmonic transform. General theoretical relations between the diffusion tensor (DT) components measured by traditional (rank-2) DT imaging (DTI) and 3D distribution of diffusivities, as measured by high angular resolution diffusion imaging (HARDI) methods, are derived. Also, the spherical tensor components from HARDI are related to the rank-2 DT. The relationships between higher- and lower-rank Cartesian DTs are also presented. The inadequacy of the traditional rank-2 tensor model is demonstrated with simulations, and the method is applied to excised rat brain data collected in a spin-echo HARDI experiment.     

Efficient measurement and calculation of MR diffusion anisotropy images using the Platonic variance method.

The Platonic variance method produces MR diffusion anisotropy (DA) images with a minimum amount of computational effort. It can be programmed in a self-contained MR sequence, thus eliminating the need for postprocessing on a separate workstation. The method uses gradient acquisition schemes, based on Platonic solids: the "icosahedric" scheme (N = 6), the "dodecahedric" scheme (N = 10), and combinations thereof. For these schemes the average of the diffusion tensor eigenvalues equals the average of the measured apparent diffusion coefficients (ADCs), and the variance of the eigenvalues equals 5/2 times the variance of the diffusion coefficients. This results in compact expressions for anisotropy measures, directly in terms of the acquired images, i.e., without calculating the eigenvalues or even the tensor elements. The resulting anisotropy images were shown to be identical to the ones traditionally derived. It is expected that this method will considerably promote the routine use of DA imaging.