Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring

Practical considerations of diffusion-weighted MRS with ultra-strong diffusion gradients

Introduction

Diffusion-weighted magnetic resonance spectroscopy (DW-MRS) offers improved cellular specificity to microstructure-compared to water-based methods alone-but spatial resolution and SNR is severely reduced and slow-diffusing metabolites necessitate higher b-values to accurately characterize their diffusion properties. Ultra-strong gradients allow access to higher b-values per-unit time, higher SNR for a given b-value, and shorter diffusion times, but introduce additional challenges such as eddy-current artefacts, gradient non-uniformity, and mechanical vibrations.

Methods

In this work, we present initial DW-MRS data acquired on a 3T Siemens Connectom scanner equipped with ultra-strong (300 mT/m) gradients. We explore the practical issues associated with this manner of acquisition, the steps that may be taken to mitigate their impact on the data, and the potential benefits of ultra-strong gradients for DW-MRS. An in-house DW-PRESS sequence and data processing pipeline were developed to mitigate the impact of these confounds. The interaction of TE, b-value, and maximum gradient amplitude was investigated using simulations and pilot data, whereby maximum gradient amplitude was restricted. Furthermore, two DW-MRS voxels in grey and white matter were acquired using ultra-strong gradients and high b-values.

Results

Simulations suggest T2-based SNR gains that are experimentally confirmed. Ultra-strong gradient acquisitions exhibit similar artefact profiles to those of lower gradient amplitude, suggesting adequate performance of artefact mitigation strategies. Gradient field non-uniformity influenced ADC estimates by up to 4% when left uncorrected. ADC and Kurtosis estimates for tNAA, tCho, and tCr align with previously published literature.

Discussion

In conclusion, we successfully implemented acquisition and data processing strategies for ultra-strong gradient DW-MRS and results indicate that confounding effects of the strong gradient system can be ameliorated, while achieving shorter diffusion times and improved metabolite SNR.

Characterization and correction of diffusion gradient-induced eddy currents in second-order motion-compensated echo-planar and spiral cardiac DTI

PURPOSE

Very high gradient amplitudes played out over extended time intervals as required for second-order motion-compensated cardiac DTI may violate the assumption of a linear time-invariant gradient system model. The aim of this work was to characterize diffusion gradient-related system nonlinearity and propose a correction approach for echo-planar and spiral spin-echo motion-compensated cardiac DTI.

METHODS

Diffusion gradient-induced eddy currents of 9 diffusion directions were characterized at b values of 150 s/mm² and 450 s/mm² for a 1.5 Tesla system and used to correct phantom, ex vivo, and in vivo motion-compensated cardiac DTI data acquired with echo-planar and spiral trajectories. Predicted trajectories were calculated using gradient impulse response function and diffusion gradient strength- and direction-dependent zeroth- and first-order eddy current responses. A reconstruction method was implemented using the predicted <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mi>k</mml:mi></mml:mrow> <mml:annotation>$$ k $$</mml:annotation></mml:semantics> </mml:math> -space trajectories to additionally include off-resonances and concomitant fields. Resulting images were compared to a reference reconstruction omitting diffusion gradient-induced eddy current correction.

RESULTS

Diffusion gradient-induced eddy currents exhibited nonlinear effects when scaling up the gradient amplitude and could not be described by a 3D basis alone. This indicates that a gradient impulse response function does not suffice to describe diffusion gradient-induced eddy currents. Zeroth- and first-order diffusion gradient-induced eddy current effects of up to -1.7 rad and -16 to +12 rad/m, respectively, were identified. Zeroth- and first-order diffusion gradient-induced eddy current correction yielded improved image quality upon image reconstruction.

CONCLUSION

The proposed approach offers correction of diffusion gradient-induced zeroth- and first-order eddy currents, reducing image distortions to promote improvements of second-order motion-compensated spin-echo cardiac DTI.

Comparison of interpolation methods of predominant cardiomyocyte orientation from in vivo and ex vivo cardiac diffusion tensor imaging data

Cardiac electrophysiology and cardiac mechanics both depend on the average cardiomyocyte long-axis orientation. In the realm of personalized medicine, knowledge of the patient-specific changes in cardiac microstructure plays a crucial role. Patient-specific computational modelling has emerged as a tool to better understand disease progression. In vivo cardiac diffusion tensor imaging (cDTI) is a vital tool to non-destructively measure the average cardiomyocyte long-axis orientation in the heart. However, cDTI suffers from long scan times, rendering volumetric, high-resolution acquisitions challenging. Consequently, interpolation techniques are needed to populate bio-mechanical models with patient-specific average cardiomyocyte long-axis orientations. In this work, we compare five interpolation techniques applied to in vivo and ex vivo porcine input data. We compare two tensor interpolation approaches, one rule-based approximation, and two data-driven, low-rank models. We demonstrate the advantage of tensor interpolation techniques, resulting in lower interpolation errors than do low-rank models and rule-based methods adapted to cDTI data. In an ex vivo comparison, we study the influence of three imaging parameters that can be traded off against acquisition time: in-plane resolution, signal to noise ratio, and number of acquired short-axis imaging slices.

Characterization and correction of time-varying eddy currents for diffusion MRI

PURPOSE

To develop and test a method for reducing artifacts due to time-varying eddy currents in oscillating gradient spin-echo (OGSE) diffusion images.

METHODS

An in-house algorithm (TVEDDY), that for the first time retrospectively models eddy current decay, was tested on pulsed gradient spin echo and OGSE brain images acquired at 7 T. Image pairs were acquired using opposite polarity diffusion gradients. A three-parameter exponential decay model (two amplitudes and a time constant) was used to characterize and correct eddy current distortions by minimizing the intensity difference between image pairs. Correction performance was compared with conventional correction methods by evaluating the mean squared error (MSE) between diffusion-weighted images acquired with opposite polarity diffusion gradients. As a ground-truth comparison, images were corrected using field dynamics up to third order in space, measured using a field monitoring system.

RESULTS

Time-varying eddy currents were observed for OGSE, which introduced blurring that was not reduced using the traditional approach but was diminished considerably with TVEDDY and field monitoring-informed model-based reconstruction. No MSE difference was observed between the conventional approach and TVEDDY for pulsed gradient spin echo, but for OGSE TVEDDY resulted in significantly lower MSE than the conventional approach. The field-monitoring reconstruction had the lowest MSE for both pulsed gradient spin echo and OGSE.

CONCLUSION

This work establishes that it is possible to estimate time-varying eddy currents from the actual diffusion data, which provides substantial image-quality improvements for gradient-intensive diffusion MRI acquisitions like OGSE.

Widespread effects of dMRI data quality on diffusion measures in children

Diffusion magnetic resonance imaging (dMRI) datasets are susceptible to several confounding factors related to data quality, which is especially true in studies involving young children. With the recent trend of large‐scale multicenter studies, it is more critical to be aware of the varied impacts of data quality on measures of interest. Here, we investigated data quality and its effect on different diffusion measures using a multicenter dataset. dMRI data were obtained from 691 participants (5–17 years of age) from six different centers. Six data quality metrics—contrast to noise ratio, outlier slices, and motion (absolute, relative, translation, and rotational)—and four diffusion measures—fractional anisotropy, mean diffusivity, tract density, and length—were computed for each of 36 major fiber tracts for all participants. The results indicated that four out of six data quality metrics (all except absolute and translation motion) differed significantly between centers. Associations between these data quality metrics and the diffusion measures differed significantly across the tracts and centers. Moreover, these effects remained significant after applying recently proposed harmonization algorithms that purport to remove unwanted between‐site variation in diffusion data. These results demonstrate the widespread impact of dMRI data quality on diffusion measures. These tracts and measures have been routinely associated with individual differences as well as group‐wide differences between neurotypical populations and individuals with neurological or developmental disorders. Accordingly, for analyses of individual differences or group effects (particularly in multisite dataset), we encourage the inclusion of data quality metrics in dMRI analysis.

Correcting Susceptibility Artifacts of MRI Sensors in Brain Scanning: A 3D Anatomy-Guided Deep Learning Approach

Echo planar imaging (EPI), a fast magnetic resonance imaging technique, is a powerful tool in functional neuroimaging studies. However, susceptibility artifacts, which cause misinterpretations of brain functions, are unavoidable distortions in EPI. This paper proposes an end-to-end deep learning framework, named TS-Net, for susceptibility artifact correction (SAC) in a pair of 3D EPI images with reversed phase-encoding directions. The proposed TS-Net comprises a deep convolutional network to predict a displacement field in three dimensions to overcome the limitation of existing methods, which only estimate the displacement field along the dominant-distortion direction. In the training phase, anatomical T1-weighted images are leveraged to regularize the correction, but they are not required during the inference phase to make TS-Net more flexible for general use. The experimental results show that TS-Net achieves favorable accuracy and speed trade-off when compared with the state-of-the-art SAC methods, i.e., TOPUP, TISAC, and S-Net. The fast inference speed (less than a second) of TS-Net makes real-time SAC during EPI image acquisition feasible and accelerates the medical image-processing pipelines.

Diffusion-PEPTIDE: Distortion- and blurring-free diffusion imaging with self-navigated motion-correction and relaxometry capabilities

PURPOSE

To implement the time-resolved relaxometry PEPTIDE technique into a diffusion acquisition to provide self-navigated, distortion- and blurring-free diffusion imaging that is robust to motion, while simultaneously providing T₂ and T 2 ∗ mapping.

THEORY AND METHODS

The PEPTIDE readout was implemented into a spin-echo diffusion acquisition, enabling reconstruction of a time-series of T₂ - and T 2 ∗ -weighted images, free from conventional echo planar imaging (EPI) distortion and blurring, for each diffusion-encoding. Robustness of PEPTIDE to motion and shot-to-shot phase variation was examined through a deliberate motion-corrupted diffusion experiment. Two diffusion-relaxometry in vivo brain protocols were also examined: (1)1 × 1 × 3 mm³ across 32 diffusion directions in 20 min, (2)1.5 × 1.5 × 3.0 mm³ across 6 diffusion-weighted images in 3.4 min. T₂ , T 2 ∗ , and diffusion parameter maps were calculated from these data. As initial exploration of the rich diffusion-relaxometry data content for use in multi-compartment modeling, PEPTIDE data were acquired of a gadolinium-doped asparagus phantom. These datasets contained two compartments with different relaxation parameters and different diffusion orientation properties, and T₂ relaxation variations across these diffusion directions were explored.

RESULTS

Diffusion-PEPTIDE showed the capability to provide high quality diffusion images and T₂ and T 2 ∗ maps from both protocols. The reconstructions were distortion-free, avoided potential resolution losses exceeding 100% in equivalent EPI acquisitions, and showed tolerance to nearly 30° of rotational motion. Expected variation in T₂ values as a function of diffusion direction was observed in the two-compartment asparagus phantom (P < .01), demonstrating potential to explore diffusion-PEPTIDE data for multi-compartment modeling.

CONCLUSIONS

Diffusion-PEPTIDE provides highly robust diffusion and relaxometry data and offers potential for future applications in diffusion-relaxometry multi-compartment modeling.

ADC measurements on the Unity MR-linac - A recommendation on behalf of the Elekta Unity MR-linac consortium

BACKGROUND AND PURPOSE

Diffusion-weighted imaging (DWI) for treatment response monitoring is feasible on hybrid magnetic resonance linear accelerator (MR-linac) systems. The MRI scanner of the Elekta Unity system has an adjusted design compared to diagnostic scanners. We investigated its impact on measuring the DWI-derived apparent diffusion coefficient (ADC) regarding three aspects: the choice of b-values, the spatial variation of the ADC, and scanning during radiation treatment. The aim of this study is to give recommendations for accurate ADC measurements on Unity systems.

MATERIALS AND METHODS

Signal-to-noise ratio (SNR) measurements with increasing b-values were done to determine the highest bvalue that can be measured reliably. The spatial variation of the ADC was assessed on six Unity systems with a cylindrical phantom of 40 cm diameter. The influence of gantry rotation and irradiation was investigated by acquiring DWI images before and during treatment of 11 prostate cancer patients.

RESULTS

On the Unity system, a maximum b-value of 500 s/mm2 should be used for ADC quantification, as a trade-off between SNR and diffusion weighting. Accurate ADC values were obtained within 7 cm from the iso-center, while outside this region ADC values deviated more than 5%. The ADC was not influenced by the rotating linac or irradiation during treatment.

CONCLUSION

We provide Unity system specific recommendations for measuring the ADC. This will increase the consistency of ADC values acquired in different centers on the Unity system, enabling large cohort studies for biomarker discovery and treatment response monitoring.

Strong diffusion gradients allow the separation of intra- and extra-axonal gradient-echo signals in the human brain

The quantification of brain white matter properties is a key area of application of Magnetic Resonance Imaging (MRI), with much effort focused on using MR techniques to quantify tissue microstructure. While diffusion MRI probes white matter (WM) microstructure by characterising the sensitivity of Brownian motion of water molecules to anisotropic structures, susceptibility-based techniques probe the tissue microstructure by observing the effect of interaction between the tissue and the magnetic field. Here, we unify these two complementary approaches by combining ultra-strong (300 mT/m) gradients with a novel Diffusion-Filtered Asymmetric Spin Echo (D-FASE) technique. Using D-FASE we can separately assess the evolution of the intra- and extra-axonal signals under the action of susceptibility effects, revealing differences in the behaviour in different fibre tracts. We observed that the effective relaxation rate of the ASE signal in the corpus callosum decreases with increasing b-value in all subjects (from 17.1 ± 0.7 s⁻¹ at b = 0 s/mm² to 14.6 ± 0.7 s⁻¹ at b = 4800 s/mm²), while this dependence on b in the corticospinal tract is less pronounced (from 12.0± 1.1 s⁻¹ at b = 0s/mm² to 10.7 ± 0.5 s⁻¹ at b = 4800 s/mm²). Voxelwise analysis of the signal evolution with respect to b-factor and acquisition delay using a microscopic model demonstrated differences in gradient echo signal evolution between the intra- and extra-axonal pools.

Motion and eddy current-induced signal dephasing in in vivo cardiac DTI

PURPOSE

To address motion in cardiac DWI, stimulated-echo acquisition mode (STEAM) and second-order motion-compensated spin-echo (SE) sequences have been proposed. Despite applying motion-compensation strategies, residual motion can cause misleading signal attenuation. The purpose of this study is to estimate the motion-induced error in both sequences by analysis of image phase.

METHODS

Diffusion-weighted motion-compensated SE sequences and STEAM imaging was applied in vivo with diffusion encoding along 3 orthogonal directions. A b-value range of 100 to 600 s/mm² and trigger delays of 25%, 50%, and 75% of end systole and middiastole were used. Eddy-current contributions were obtained from phantom measurements. After computation of motion-induced phase maps, the amount of signal dephasing was computed from phase gradients, and the resulting errors in diffusion tensor parameters were calculated.

RESULTS

Motion-induced dephasing from the STEAM sequence showed less dependency on the b-value and no dependency on the heart phase, whereas SE imaging performed best at 75% end systole followed by 50% end systole and middiastole. For a typical experimental setting, errors of 3.3%/3.0% mean diffusivity, 4.9%/4.8% fractional anisotropy, 2.9º/3.2º helix angulation, 0.8º/0.7º transverse angulation, and 9.9º/10.0º sheet angulation (SE/STEAM) were calculated.

CONCLUSION

Image phase contains valuable information regarding uncompensated motion and eddy currents in cardiac DTI. Although the trigger delay window for SE is narrower compared with the STEAM-based approach, imaging in both systole and diastole is feasible and both sequences perform similarly if the trigger delays are selected carefully with SE.

Diffusion MRI in the brain - Theory and concepts

Over the past two decades, diffusion MRI has become an essential tool in neuroimaging investigations. This is due to its sensitivity to the motion of water molecules as they diffuse through the microstructural environment, allowing diffusion MRI to be used as a 'probe' of tissue microstructure. Furthermore, this sensitivity is strongly direction-dependent, notably in brain white matter, due to the alignment of structures that restrict or hinder the motion of water molecules, notably axonal membranes. This provides a means of inferring the orientation of fibres in vivo, and by use of appropriate fibre-tracking algorithms, of delineating the path of white matter tracts in the brain. The ability to perform so-called tractography in humans in vivo non-invasively is unique to diffusion MRI, and is now used in applications such as neurosurgery planning and more broadly within investigations of brain connectomics. This review describes the theory and concepts of diffusion MRI and describes its most important areas of application in the brain, with a strong focus on tractography.

Huntington's disease: Brain imaging in Huntington's disease

Huntington's disease (HD) gene-carriers show prominent neuronal loss by end-stage disease, and the use of magnetic resonance imaging (MRI) has been increasingly used to quantify brain changes during earlier stages of the disease. MRI offers an in vivo method of measuring structural and functional brain change. The images collected via MRI are processed to measure different anatomical features, such as brain volume, macro- and microstructural changes within white matter and functional brain activity. Structural imaging has demonstrated significant volume loss across multiple white and gray matter regions in HD, particularly within subcortical structures. There also appears to be increasing disorganization of white matter tracts and between-region connectivity with increasing disease progression. Finally, functional changes are thought to represent changes in brain activity underlying compensatory mechanisms in HD. This chapter will provide an overview of the principles of MRI and practicalities associated with using MRI in HD studies, and summarize findings from MRI studies investigating brain structure and function in HD.

The quest for high spatial resolution diffusion-weighted imaging of the human brain in vivo

Diffusion-weighted imaging, a contrast unique to MRI, is used for assessment of tissue microstructure in vivo. However, this exquisite sensitivity to finer scales far above imaging resolution comes at the cost of vulnerability to errors caused by sources of motion other than diffusion motion. Addressing the issue of motion has traditionally limited diffusion-weighted imaging to a few acquisition techniques and, as a consequence, to poorer spatial resolution than other MRI applications. Advances in MRI imaging methodology have allowed diffusion-weighted MRI to push to ever higher spatial resolution. In this review we focus on the pulse sequences and associated techniques under development that have pushed the limits of image quality and spatial resolution in diffusion-weighted MRI.

Diffusion-prepared stimulated-echo turbo spin echo (DPsti-TSE): An eddy current-insensitive sequence for three-dimensional high-resolution and undistorted diffusion-weighted imaging

In this study, we present a new three-dimensional (3D), diffusion-prepared turbo spin echo sequence based on a stimulated-echo read-out (DPsti-TSE) enabling high-resolution and undistorted diffusion-weighted imaging (DWI). A dephasing gradient in the diffusion preparation module and rephasing gradients in the turbo spin echo module create stimulated echoes, which prevent signal loss caused by eddy currents. Near to perfect agreement of apparent diffusion coefficient (ADC) values between DPsti-TSE and diffusion-weighted echo planar imaging (DW-EPI) was demonstrated in both phantom transient signal experiments and phantom imaging experiments. High-resolution and undistorted DPsti-TSE was demonstrated in vivo in prostate and carotid vessel wall. 3D whole-prostate DWI was achieved with four b values in only 6 min. Undistorted ADC maps of the prostate peripheral zone were obtained at low and high imaging resolutions with no change in mean ADC values [(1.60 ± 0.10) × 10^-3 versus (1.60 ± 0.02) × 10^-3  mm² /s]. High-resolution 3D DWI of the carotid vessel wall was achieved in 12 min, with consistent ADC values [(1.40 ± 0.23) × 10^-3  mm² /s] across different subjects, as well as slice locations through the imaging volume. This study shows that DPsti-TSE can serve as a robust 3D diffusion-weighted sequence and is an attractive alternative to the traditional two-dimensional DW-EPI approaches.

Concurrent recording of RF pulses and gradient fields - comprehensive field monitoring for MRI

Reconstruction of MRI data is based on exact knowledge of all magnetic field dynamics, since the interplay of RF and gradient pulses generates the signal, defines the contrast and forms the basis of resolution in spatial and spectral dimensions. Deviations caused by various sources, such as system imperfections, delays, eddy currents, drifts or externally induced fields, can therefore critically limit the accuracy of MRI examinations. This is true especially at ultra-high fields, because many error terms scale with the main field strength, and higher available SNR renders even smaller errors relevant. Higher baseline field also often requires higher acquisition bandwidths and faster signal encoding, increasing hardware demands and the severity of many types of hardware imperfection. To address field imperfections comprehensively, in this work we propose to expand the concept of magnetic field monitoring to also encompass the recording of RF fields. In this way, all dynamic magnetic fields relevant for spin evolution are covered, including low- to audio-frequency magnetic fields as produced by main magnets, gradients and shim systems, as well as RF pulses generated with single- and multiple-channel transmission systems. The proposed approach permits field measurements concurrently with actual MRI procedures on a strict common time base. The combined measurement is achieved with an array of miniaturized field probes that measure low- to audio-frequency fields via (19) F NMR and simultaneously pick up RF pulses in the MRI system's (1) H transmit band. Field recordings can form the basis of system calibration, retrospective correction of imaging data or closed-loop feedback correction, all of which hold potential to render MRI more robust and relax hardware requirements. The proposed approach is demonstrated for a range of imaging methods performed on a 7 T human MRI system, including accelerated multiple-channel RF pulses. Copyright © 2015 John Wiley & Sons, Ltd.

Demonstration of Non-Gaussian Restricted Diffusion in Tumor Cells Using Diffusion Time-Dependent Diffusion-Weighted Magnetic Resonance Imaging Contrast

The diffusion-weighted magnetic resonance imaging (DWI) technique enables quantification of water mobility for probing microstructural properties of biological tissue and has become an effective tool for collecting information about the underlying pathology of cancerous tissue. Measurements using multiple b-values have indicated biexponential signal attenuation, ascribed to “fast” (high ADC) and “slow” (low ADC) diffusion components. In this empirical study, we investigate the properties of the diffusion time (Δ)-dependent components of the diffusion-weighted (DW) signal in a constant b-value experiment. A xenograft gliobastoma mouse was imaged using Δ = 11 ms, 20 ms, 40 ms, 60 ms, and b = 500–4000 s/mm² in intervals of 500 s/mm². Data were corrected for EPI distortions, and the Δ-dependence on the DW-signal was measured within three regions of interest [intermediate- and high-density tumor regions and normal-appearing brain (NAB) tissue regions]. In this study, we verify the assumption that the slow decaying component of the DW-signal is non-Gaussian and dependent on Δ, consistent with restricted diffusion of the intracellular space. As the DW-signal is a function of Δ and is specific to restricted diffusion, manipulating Δ at constant b-value (cb) provides a complementary and direct approach for separating the restricted from the hindered diffusion component. We found that Δ-dependence is specific to the tumor tissue signal. Based on an extended biexponential model, we verified the interpretation of the diffusion time-dependent contrast and successfully estimated the intracellular restricted ADC, signal volume fraction, and cell size within each ROI.

An integrated approach to correction for off-resonance effects and subject movement in diffusion MR imaging

In this paper we describe a method for retrospective estimation and correction of eddy current (EC)-induced distortions and subject movement in diffusion imaging. In addition a susceptibility-induced field can be supplied and will be incorporated into the calculations in a way that accurately reflects that the two fields (susceptibility- and EC-induced) behave differently in the presence of subject movement. The method is based on registering the individual volumes to a model free prediction of what each volume should look like, thereby enabling its use on high b-value data where the contrast is vastly different in different volumes. In addition we show that the linear EC-model commonly used is insufficient for the data used in the present paper (high spatial and angular resolution data acquired with Stejskal-Tanner gradients on a 3T Siemens Verio, a 3T Siemens Connectome Skyra or a 7T Siemens Magnetome scanner) and that a higher order model performs significantly better. The method is already in extensive practical use and is used by four major projects (the WU-UMinn HCP, the MGH HCP, the UK Biobank and the Whitehall studies) to correct for distortions and subject movement.

Eddy current compensated double diffusion encoded (DDE) MRI

PURPOSE

Eddy currents might lead to image distortions in diffusion-weighted echo planar imaging. A method is proposed to reduce their effects on double diffusion encoding (DDE) MRI experiments and the thereby derived microscopic fractional anisotropy (μFA).

METHODS

The twice-refocused spin echo scheme was adapted for DDE measurements. To assess the effect of individual diffusion encodings on the image distortions, measurements of a grid of plastic rods in water were performed. The effect of eddy current compensation on μFA measurements was evaluated in the brains of six healthy volunteers.

RESULTS

The use of an eddy current compensation reduced the signal variation. As expected, the distortions caused by the second encoding were larger than those of the first encoding, entailing a stronger need to compensate for them. For an optimal result, however, both encodings had to be compensated. The artifact reduction strongly improved the measurement of the μFA in ventricles and gray matter by reducing the overestimation. An effect of the compensation on absolute μFA values in white matter was not observed.

CONCLUSION

It is advisable to compensate both encodings in DDE measurements for eddy currents. Magn Reson Med 77:328-335, 2017. © 2015 Wiley Periodicals, Inc.

Non-invasive pressure difference estimation from PC-MRI using the work-energy equation

Pressure difference is an accepted clinical biomarker for cardiovascular disease conditions such as aortic coarctation. Currently, measurements of pressure differences in the clinic rely on invasive techniques (catheterization), prompting development of non-invasive estimates based on blood flow. In this work, we propose a non-invasive estimation procedure deriving pressure difference from the work-energy equation for a Newtonian fluid. Spatial and temporal convergence is demonstrated on in silico Phase Contrast Magnetic Resonance Image (PC-MRI) phantoms with steady and transient flow fields. The method is also tested on an image dataset generated in silico from a 3D patient-specific Computational Fluid Dynamics (CFD) simulation and finally evaluated on a cohort of 9 subjects. The performance is compared to existing approaches based on steady and unsteady Bernoulli formulations as well as the pressure Poisson equation. The new technique shows good accuracy, robustness to noise, and robustness to the image segmentation process, illustrating the potential of this approach for non-invasive pressure difference estimation.

Spin echo versus stimulated echo diffusion tensor imaging of the in vivo human heart

Purpose

To compare signal‐to‐noise ratio (SNR) efficiency and diffusion tensor metrics of cardiac diffusion tensor mapping using acceleration‐compensated spin‐echo (SE) and stimulated echo acquisition mode (STEAM) imaging.

Methods

Diffusion weighted SE and STEAM sequences were implemented on a clinical 1.5 Tesla MR system. The SNR efficiency of SE and STEAM was measured (b = 50–450 s/mm²) in isotropic agar, anisotropic diffusion phantoms and the in vivo human heart. Diffusion tensor analysis was performed on mean diffusivity, fractional anisotropy, helix and transverse angles.

Results

In the isotropic phantom, the ratio of SNR efficiency for SE versus STEAM, SNR_t(SE/STEAM), was 2.84 ± 0.08 for all tested b‐values. In the anisotropic diffusion phantom the ratio decreased from 2.75 ± 0.05 to 2.20 ± 0.13 with increasing b‐value, similar to the in vivo decrease from 2.91 ± 0.43 to 2.30 ± 0.30. Diffusion tensor analysis revealed reduced deviation of helix angles from a linear transmural model and reduced transverse angle standard deviation for SE compared with STEAM. Mean diffusivity and fractional anisotropy were measured to be statistically different (P < 0.001) between SE and STEAM.

Conclusion

Cardiac DTI using motion‐compensated SE yields a 2.3–2.9× increase in SNR efficiency relative to STEAM and improved accuracy of tensor metrics. The SE method hence presents an attractive alternative to STEAM based approaches. Magn Reson Med 76:862–872, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Correction of MRI-induced geometric distortions in whole-body small animal PET-MRI

PURPOSE

The fusion of positron emission tomography (PET) and magnetic resonance imaging (MRI) data can be a challenging task in whole-body PET-MRI. The quality of the registration between these two modalities in large field-of-views (FOV) is often degraded by geometric distortions of the MRI data. The distortions at the edges of large FOVs mainly originate from MRI gradient nonlinearities. This work describes a method to measure and correct for these kind of geometric distortions in small animal MRI scanners to improve the registration accuracy of PET and MRI data.

METHODS

The authors have developed a geometric phantom which allows the measurement of geometric distortions in all spatial axes via control points. These control points are detected semiautomatically in both PET and MRI data with a subpixel accuracy. The spatial transformation between PET and MRI data is determined with these control points via 3D thin-plate splines (3D TPS). The transformation derived from the 3D TPS is finally applied to real MRI mouse data, which were acquired with the same scan parameters used in the phantom data acquisitions. Additionally, the influence of the phantom material on the homogeneity of the magnetic field is determined via field mapping.

RESULTS

The spatial shift according to the magnetic field homogeneity caused by the phantom material was determined to a mean of 0.1 mm. The results of the correction show that distortion with a maximum error of 4 mm could be reduced to less than 1 mm with the proposed correction method. Furthermore, the control point-based registration of PET and MRI data showed improved congruence after correction.

CONCLUSIONS

The developed phantom has been shown to have no considerable negative effect on the homogeneity of the magnetic field. The proposed method yields an appropriate correction of the measured MRI distortion and is able to improve the PET and MRI registration. Furthermore, the method is applicable to whole-body small animal imaging routines including different standard MRI sequences.

The Age-ility Project (Phase 1): Structural and functional imaging and electrophysiological data repository

Our understanding of the complex interplay between structural and functional organisation of brain networks is being advanced by the development of novel multi-modal analyses approaches. The Age-ility Project (Phase 1) data repository offers open access to structural MRI, diffusion MRI, and resting-state fMRI scans, as well as resting-state EEG recorded from the same community participants (n=131, 15-35 y, 66 male). Raw imaging and electrophysiological data as well as essential demographics are made available via the NITRC website. All data have been reviewed for artifacts using a rigorous quality control protocol and detailed case notes are provided.

Comparisons of Reproducibility and Mean Values of Diffusion Tensor Imaging-Derived Indices between Unipolar and Bipolar Diffusion Pulse Sequences

Eddy current distortion is an important issue that may influence the quantitative measurements of diffusion tensor imaging (DTI). The corrections of eddy current artifacts could be performed using bipolar diffusion gradients or unipolar gradients with affine registration. Whether the diffusion pulse sequence affects the quantification of DTI indices and the technique that produces more reliable DTI indices in terms of reproducibility both remain unclear. Therefore, the purpose of this study was to compare the reproducibility and mean values of DTI-derived indices between unipolar and bipolar diffusion pulse sequences based on actual human brain data. Five repeated datasets of unipolar and bipolar DTI were acquired from 10 healthy subjects at different echo times (TEs). The reproducibility and mean values of DTI indices were assessed by calculating the coefficient of variation and mean values of the 5 repeated measurements. The results revealed that the reproducibility and mean values of DTI indices were significantly affected by the pulse sequence. Unipolar DTI exhibited significantly higher reproducibility than bipolar DTI even at the same TE, and the mean values of DTI indices were significantly different between them. Therefore, we concluded that the reproducibility and mean values of DTI indices were significantly influenced by diffusion pulse sequences.