Analysis and correction of gradient nonlinearity bias in apparent diffusion coefficient measurements

Velocity-compensated intravoxel incoherent motion of the human calf muscle

PURPOSE

To determine whether intravoxel incoherent motion (IVIM) describes the blood perfusion in muscles better, assuming pseudo diffusion (Bihan Model 1) or ballistic motion (Bihan Model 2).

METHODS

IVIM parameters were measured in 18 healthy subjects with three different diffusion gradient time profiles (bipolar with two diffusion times and one with velocity compensation) and 17 b-values (0-600 s/mm2) at rest and after muscle activation. The diffusion coefficient, perfusion fraction, and pseudo-diffusion coefficient were estimated with a segmented fit in the gastrocnemius medialis (GM) and tibialis anterior (TA) muscles.

RESULTS

Velocity-compensated gradients resulted in a decreased perfusion fraction (6.9% ± 1.4% vs. 4.4% ± 1.3% in the GM after activation) and pseudo-diffusion coefficient (0.069 ± 0.046 mm2/s vs. 0.014 ± 0.006 in the GM after activation) compared to the bipolar gradients with the longer diffusion encoding time. Increased diffusion coefficients, perfusion fractions, and pseudo-diffusion coefficients were observed in the GM after activation for all gradient profiles. However, the increase was significantly smaller for the velocity-compensated gradients. A diffusion time dependence was found for the pseudo-diffusion coefficient in the activated muscle.

CONCLUSION

Velocity-compensated diffusion gradients significantly suppress the IVIM effect in the calf muscle, indicating that the ballistic limit is mostly reached, which is supported by the time dependence of the pseudo-diffusion coefficient.

The effect of elimination of gibbs ringing, noise and systematic errors on the DTI metrics and tractography in a rat brain

Diffusion tensor imaging (DTI) metrics and tractography can be biased due to low signal-to-noise ratio (SNR) and systematic errors resulting from image artifacts and imperfections in magnetic field gradients. The imperfections include non-uniformity and nonlinearity, effects caused by eddy currents, and the influence of background and imaging gradients. We investigated the impact of systematic errors on DTI metrics of an isotropic phantom and DTI metrics and tractography of a rat brain measured at high resolution. We tested denoising and Gibbs ringing removal methods combined with the B matrix spatial distribution (BSD) method for magnetic field gradient calibration. The results showed that the performance of the BSD method depends on whether Gibbs ringing is removed and the effectiveness of stochastic error removal. Region of interest (ROI)-based analysis of the DTI metrics showed that, depending on the size of the ROI and its location in space, correction methods can remove systematic bias to varying degrees. The preprocessing pipeline proposed and dedicated to this type of data together with the BSD method resulted in an even - 90% decrease in fractional anisotropy (FA) (globally and locally) in the isotropic phantom and - 45% in the rat brain. The largest global changes in the rat brain tractogram compared to the standard method without preprocessing (sDTI) were noticed after denoising. The direction of the first eigenvector obtained from DTI after denoising, Gibbs ringing removal and BSD differed by an average of 56 and 10 degrees in the ROI from sDTI and from sDTI after denoising and Gibbs ringing removal, respectively. The latter can be identified with the amount of improvement in tractography due to the elimination of systematic errors related to imperfect magnetic field gradients. Based on the results, the systematic bias for high resolution data mainly depended on SNR, but the influence of non-uniform gradients could also be seen. After denoising, the BSD method was able to further correct both the metrics and tractography of the diffusion tensor in the rat brain by taking into account the actual distribution of magnetic field gradients independent of the examined object and uniquely dependent on the scanner and sequence. This means that in vivo studies are also subject to this type of errors, which should be taken into account when processing such data.

Frequency-dependent diffusion kurtosis imaging in the human brain using an oscillating gradient spin echo sequence and a high-performance head-only gradient

Measuring the time/frequency dependence of diffusion MRI is a promising approach to distinguish between the effects of different tissue microenvironments, such as membrane restriction, tissue heterogeneity, and compartmental water exchange. In this study, we measure the frequency dependence of diffusivity (D) and kurtosis (K) with oscillating gradient diffusion encoding waveforms and a diffusion kurtosis imaging (DKI) model in human brains using a high-performance, head-only MAGNUS gradient system, with a combination of b-values, oscillating frequencies (f), and echo time that has not been achieved in human studies before. Frequency dependence of diffusivity and kurtosis are observed in both global and local white matter (WM) and gray matter (GM) regions and characterized with a power-law model ∼Λ*fθ. The frequency dependences of diffusivity and kurtosis (including changes between fmin and fmax, Λ, and θ) vary over different WM and GM regions, indicating potential microstructural differences between regions. A trend of decreasing kurtosis over frequency in the short-time limit is successfully captured for in vivo human brains. The effects of gradient nonlinearity (GNL) on frequency-dependent diffusivity and kurtosis measurements are investigated and corrected. Our results show that the GNL has prominent scaling effects on the measured diffusivity values (3.5∼5.5% difference in the global WM and 6∼8% difference in the global cortex) and subsequently affects the corresponding power-law parameters (Λ, θ) while having a marginal influence on the measured kurtosis values (<0.05% difference) and power-law parameters (Λ, θ). This study expands previous OGSE studies and further demonstrates the translatability of frequency-dependent diffusivity and kurtosis measurements to human brains, which may provide new opportunities to probe human brain microstructure in health and disease.

Efficient approximate signal reconstruction for correction of gradient nonlinearities in diffusion-weighted imaging

In diffusion weighted MRI (DW-MRI), hardware nonlinearities lead to spatial variations in the orientation and magnitude of diffusion weighting. While the correction of these spatial distortions has been well established for analyses of DW-MRI, the existing voxel-wise empirical correction for gradient nonlinearities requires reimplementation of existing models, as the resultant gradients vary by voxel. Herein, we propose a two-step signal approximation after voxel-wise correction of gradient nonlinearity effects in DW-MRI. The proposed technique (1) scales the diffusion signal and (2) resamples the gradient orientations. This results in uniform gradients across the corrected image and provides the key advantage of seamless integration into current diffusion workflows. We investigated the validity of our technique by fitting a multi-compartment neurite orientation dispersion and density imaging (NODDI) model to the empirical correction and proposed approximation in five subjects from the MASiVar pediatric dataset. We evaluated intra-cellular volume fraction (iVF), CSF volume fraction (cVF), and orientation dispersion index (ODI) from NODDI. The Cohen's d of iVF, cVF and ODI between the techniques was <0.2 indicating the proposed technique does not exhibit significant differences from the voxel-wise correction technique. Our two-step signal approximation is an efficient representation of the voxel-wise gradient table correction. Using this approximation, correction of gradient nonlinearities can be easily incorporated into existing diffusion preprocessing pipelines and is implemented in "PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images".

Mapping the impact of nonlinear gradient fields with noise on diffusion MRI

In diffusion MRI, gradient nonlinearities cause spatial variations in the magnitude and direction of diffusion gradients. Studies have shown artifacts from these distortions can results in biased diffusion tensor information and tractography. Here, we investigate the impact of gradient nonlinearity correction in the presence of noise. We introduced empirically derived gradient nonlinear fields at different signal-to-noise ratio (SNR) levels in two experiments: tensor field simulation and simulation of the brain. For each experiment, this work compares two techniques empirically: voxel-wise gradient table correction and approximate correction by scaling the signal directly. The impact was assessed through diffusion metrics including mean diffusivity (MD), fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD), and principal eigen vector (V1). The study shows (1) the correction of gradient nonlinearities will not lead to substantively incorrect estimation of diffusion metrics in a linear system, (2) gradient nonlinearity correction does not interact adversely with noise, (3) nonlinearity correction suppresses the impact of nonlinearities in typical SNR data, (4) for SNR below 30, the performance of both the gradient nonlinearity correction techniques were similar, and (5) larger impacts are seen in regions where the gradient nonlinearities are distinct. Thus, this study suggests that there were greater beneficial effects than adverse effects due to the correction of nonlinearities. Additionally, correction of nonlinearities is recommended when region of interests are in areas with pronounced nonlinearities.

Evaluation of Apparent Diffusion Coefficient Repeatability and Reproducibility for Preclinical MRIs Using Standardized Procedures and a Diffusion-Weighted Imaging Phantom

Relevant to co-clinical trials, the goal of this work was to assess repeatability, reproducibility, and bias of the apparent diffusion coefficient (ADC) for preclinical MRIs using standardized procedures for comparison to performance of clinical MRIs. A temperature-controlled phantom provided an absolute reference standard to measure spatial uniformity of these performance metrics. Seven institutions participated in the study, wherein diffusion-weighted imaging (DWI) data were acquired over multiple days on 10 preclinical scanners, from 3 vendors, at 6 field strengths. Centralized versus site-based analysis was compared to illustrate incremental variance due to processing workflow. At magnet isocenter, short-term (intra-exam) and long-term (multiday) repeatability were excellent at within-system coefficient of variance, wCV [±CI] = 0.73% [0.54%, 1.12%] and 1.26% [0.94%, 1.89%], respectively. The cross-system reproducibility coefficient, RDC [±CI] = 0.188 [0.129, 0.343] µm2/ms, corresponded to 17% [12%, 31%] relative to the reference standard. Absolute bias at isocenter was low (within 4%) for 8 of 10 systems, whereas two high-bias (>10%) scanners were primary contributors to the relatively high RDC. Significant additional variance (>2%) due to site-specific analysis was observed for 2 of 10 systems. Base-level technical bias, repeatability, reproducibility, and spatial uniformity patterns were consistent with human MRIs (scaled for bore size). Well-calibrated preclinical MRI systems are capable of highly repeatable and reproducible ADC measurements.

Mapping the Impact of Approximate Gradient Nonlinearity Fields Correction on Tractography

Nonlinear gradients impact diffusion weighted MRI by introducing spatial variation in estimated diffusion tensors. Recent studies have shown that increasing signal-to-noise ratios and the use of ultra-strong gradients may lead to clinically significant impacts on analyses due to these nonlinear gradients in microstructural measures. These effects can potentially bias tractography results and cause misinterpretation of data. Herein, we characterize the impact of an "approximate" gradient nonlinearity correction technique in tractography using empirically derived gradient nonlinear fields. This technique scales the diffusion signal by the change in magnitude due to the gradient nonlinearities, without concomitant correction of gradient direction errors. The impact of this correction on tractography is assessed through white matter bundle segmentation and connectomics via bundle-wise volume, fractional anisotropy, mean diffusivity, radial diffusivity, axial diffusivity, primary eigenvector, and length; as well as the modularity, global efficiency, and characteristic path length connectomics graph measures. We investigate the differences between (1) these measures directly and (2) the within session variability of these measures before and after approximate correction in 61 subjects from the MASiVar pediatric reproducibility dataset. We find approximate correction results is little to no differences on the population level, but large differences on the subject-specific level for both the measures directly and their within session variability. Thus, this study suggests though approximate correction of gradient nonlinearities may not change tractography findings on the population level, subject-specific interpretations may exhibit large fluctuations. A limitation is the lack of comparison with the empirical voxel-wise gradient table correction.

Diffusion-Weighted Imaging of the Abdomen: Correction for Gradient Nonlinearity Bias in Apparent Diffusion Coefficient

BACKGROUND

Gradient nonlinearity (GNL) introduces spatial nonuniformity bias in apparent diffusion coefficient (ADC) measurements, especially at large offsets from the magnet isocenter.

PURPOSE

To investigate the effects of GNL in abdominal ADC measurements and to develop an ADC bias correction procedure.

STUDY TYPE

Retrospective.

PHANTOM/POPULATION

Two homemade ultrapure water phantoms/25 patients with histologically confirmed pancreatic ductal adenocarcinoma (PDAC).

FIELD STRENGTH/SEQUENCE

A 3.0 T/diffusion-weighted imaging (DWI) with single-shot echo-planar imaging sequence.

ASSESSMENT

ADC bias was computed in the three orthogonal directions at different offset locations. The spatial-dependent correctors of ADC bias were generated from the ADCs of phantom 1. The ADCs were estimated before and after corrections for the phantom 1 with both the proposed approach and the theoretical GNL correction method. For the patients, ADCs were measured in abdominal tissues including left and right liver lobes, PDAC, spleen, bilateral kidneys, and bilateral paraspinal muscles.

STATISTICAL TEST

Friedman tests and Wilcoxon tests.

RESULTS

The ADC bias measured by phantom 1 was 9.7% and 12.6% higher in the right-left and anterior-posterior directions and 9.2% lower in the superior-inferior direction at the 150 mm offsets from the magnetic isocenter. The corrected vs. the uncorrected ADCs measurements (median: 2.20 × 10^-3  mm² /sec for both the proposed method and the theoretical GNL method vs. 2.31 × 10^-3  mm² /sec, respectively) and their relative ADC errors (0.014, 0.016, and 0.054, respectively) were lower in the phantom 1. The relative ADC errors substantially decreased after correction in the phantom 2 (median: 0.048 and -0.008, respectively). The ADCs of all the abdominal tissues were lower after correction except for the left liver lobes (P = 0.13).

DATA CONCLUSION

GNL bias in abdominal ADC can be measured by a DWI phantom. The proposed correction procedure was successfully applied for the bias correction in abdominal ADC.

EVIDENCE LEVEL

3.

TECHNICAL EFFICACY

Stage 1.

Contrastive semi-supervised harmonization of single-shell to multi-shell diffusion MRI

Diffusion weighted MRI (DW-MRI) harmonization is necessary for multi-site or multi-acquisition studies. Current statistical methods address the need to harmonize from one site to another, but do not simultaneously consider the use of multiple datasets which are comprised of multiple sites, acquisitions protocols, and age demographics. This work explores deep learning methods which can generalize across these variations through semi-supervised and unsupervised learning while also learning to estimate multi-shell data from single-shell data using the Multi-shell Diffusion MRI Harmonization Challenge (MUSHAC) and Baltimore Longitudinal Study on Aging (BLSA) datasets. We compare disentanglement harmonization models, which seek to encode anatomy and acquisition in separate latent spaces, and a CycleGAN harmonization model, which uses generative adversarial networks (GAN) to perform style transfer between sites, to the baseline preprocessing and to SHORE interpolation. We find that the disentanglement models achieve superior performance in harmonizing all data while at the same transforming the input data to a single target space across several diffusion metrics (fractional anisotropy, mean diffusivity, mean kurtosis, primary eigenvector).

Clinical translation of quantitative magnetic resonance imaging biomarkers - An overview and gap analysis of current practice

PURPOSE

This overview of the current landscape of quantitative magnetic resonance imaging biomarkers (qMR IBs) aims to support the standardisation of academic IBs to assist their translation to clinical practice.

METHODS

We used three complementary approaches to investigate qMR IB use and quality management practices within the UK: 1) a literature search of qMR and quality management terms during 2011-2015 and 2016-2020; 2) a database search for clinical research studies using qMR IBs during 2016-2020; and 3) a survey to ascertain the current availability and quality management practices for clinical MRI scanners and associated equipment at research institutions across the UK.

RESULTS

The analysis showed increased use of all qMR methods between the periods 2011-2015 and 2016-2020 and diffusion-tensor MRI and volumetry to be popular methods. However, the "translation ratio" of journal articles to clinical research studies was higher for qMR methods that have evidence of clinical translation via a commercial route, such as fat fraction and T2 mapping. The number of journal articles citing quality management terms doubled between the periods 2011-2015 and 2016-2020; although, its proportion relative to all journal articles only increased by 3.0%. The survey suggested that quality assurance (QA) and quality control (QC) of data acquisition procedures are under-reported in the literature and that QA/QC of acquired data/data analysis are under-developed and lack consistency between institutions.

CONCLUSIONS

We summarise current attempts to standardise and translate qMR IBs, and conclude by outlining the ideal quality management practices and providing a gap analysis between current practice and a metrological standard.

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Technical Note: Temperature and concentration dependence of water diffusion in polyvinylpyrrolidone solutions

Objective

The goal of this work is to provide temperature and concentration calibration of water diffusivity in polyvinylpyrrolidone (PVP) solutions used in phantoms to assess system bias and linearity in apparent diffusion coefficient (ADC) measurements.

Method

ADC measurements were performed for 40 kDa (K40) PVP of six concentrations (0%, 10%, 20%, 30%, 40%, and 50% by weight) at three temperatures (19.5°C, 22.5°C, and 26.4°C), with internal phantom temperature monitored by optical thermometer (±0.2°C). To achieve ADC measurement and fit accuracy of better than 0.5%, three orthogonal diffusion gradients were calibrated using known water diffusivity at 0°C and system gradient nonlinearity maps. Noise‐floor fit bias was also controlled by limiting the maximum b‐value used for ADC calculation of each sample. The ADC temperature dependence was modeled by Arrhenius functions of each PVP concentration. The concentration dependence was modeled by quadratic function for ADC normalized by the theoretical water diffusion values. Calibration coefficients were obtained from linear regression model fits.

Results

Measured phantom ADC values increased with temperature and decreasing PVP concentration, [PVP]. The derived Arrhenius model parameters for [PVP] between 0% and 50%, are reported and can be used for K40 ADC temperature calibration with absolute ADC error within ±0.016 μm²/ms. Arrhenius model fit parameters normalized to water value scaled with [PVP] between 10% and 40%, and proportional change in activation energy increased faster than collision frequency. ADC normalization by water diffusivity, D_W, from the Speedy–Angell relation accounted for the bulk of temperature dependence (±0.035 μm²/ms) and yielded quadratic calibration for ADC_PVP/D_W = (12.5 ± 0.7) ·10⁻⁵·[PVP]² − (23.2 ± 0.3)·10⁻³·[PVP]+1, nearly independent of PVP molecular weight and temperature.

Conclusion

The study provides ground‐truth ADC values for K40 PVP solutions commonly used in diffusion phantoms for scanning at ambient room temperature. The described procedures and the reported calibration can be used for quality control and standardization of measured ADC values of PVP at different concentrations and temperatures.

Long-Term Stability of Gradient Characteristics Warrants Model-Based Correction of Diffusion Weighting Bias

The study aims to test the long-term stability of gradient characteristics for model-based correction of diffusion weighting (DW) bias in an apparent diffusion coefficient (ADC) for multisite imaging trials. Single spin echo (SSE) DWI of a long-tube ice-water phantom was acquired quarterly on six MR scanners over two years for individual diffusion gradient channels, along with B0 mapping, as a function of right-left (RL) and superior-inferior (SI) offsets from the isocenter. Additional double spin-echo (DSE) DWI was performed on two systems. The offset dependences of derived ADC were fit to 4th-order polynomials. Chronic shim gradients were measured from spatial derivatives of B0 maps along the tube direction. Gradient nonlinearity (GNL) was modeled using vendor-provided gradient field descriptions. Deviations were quantified by root-mean-square differences (RMSD), normalized to reference ice-water ADC, between the model and reference (RMSDREF), measurement and model (RMSDEXP), and temporal measurement variations (RMSDTMP). Average RMSDREF was 4.9 ± 3.2 (%RL) and –14.8 ± 3.8 (%SI), and threefold larger than RMSDEXP. RMSDTMP was close to measurement errors (~3%). GNL-induced bias across gradient systems varied up to 20%, while deviation from the model accounted at most for 6.5%, and temporal variation for less than 3% of ADC reproducibility error. Higher SSE RMSDEXP = 7.5–11% was reduced to 2.5–4.8% by DSE, consistent with the eddy current origin. Measured chronic shim gradients below 0.1 mT/m had a minor contribution to ADC bias. The demonstrated long-term stability of spatial ADC profiles and consistency with system GNL models justifies retrospective and prospective DW bias correction based on system gradient design models. Residual errors due to eddy currents and shim gradients should be corrected independent of GNL.

Mapping the Impact of Non-Linear Gradient Fields on Diffusion MRI Tensor Estimation

Non-linear gradients impact diffusion weighted (DW) MRI by corrupting the experimental setup and lead to problems during image encoding including the effects in-plane distortion, in-plane shifts, intensity modulations and phase errors. Recent studies have been shown this may present significant complication in the interpretation of results and conclusion while studying tractography and tissue microstructure in data. To interpret the degree in consequences of gradient non-linearities between the desired and achieved gradients, we introduced empirically derived gradient nonlinear fields at different orientations and different tensor properties. The impact is assessed through diffusion tensor properties including mean diffusivity (MD), fractional anisotropy (FA) and principal eigen vector (PEV). The study shows lower FA are more susceptible to LR fields and LR fields with determinant <1 or >1 corrupt tensor more. The corruption can result in significantly different FA based on true-FA and LR field. Apparent MD decreases for negative determinant, on the other hand positive determinant shows the opposite effect. LR field have a larger impact on PEV when FA value is small. The results are dependent on the underlying orientation, non-linear field corruption can cause both increase and decrease of estimated FA, MD and PEV value. This work provides insight into characterizing the non-linear gradient error and aid in selecting correction techniques to address the inaccuracies in b-values.

Volumetric and multispectral DWI near metallic implants using a non-linear phase Carr-Purcell-Meiboom-Gill diffusion preparation

PURPOSE

DWI near metal implants has not been widely explored due to substantial challenges associated with through-slice and in-plane distortions, the increased encoding requirement of different spectral bins, and limited SNR. There is no widely adopted clinical protocol for DWI near metal since the commonly used EPI trajectory fails completely due to distortion from extreme off-resonance ranging from 2 to 20 kHz. We present a sequence that achieves DWI near metal with moderate b-values (400-500 s/mm² ) and volumetric coverage in clinically feasible scan times.

THEORY AND METHODS

Multispectral excitation with Cartesian sampling, view angle tilting, and kz phase encoding reduce in-plane and through-plane off-resonance artifacts, and Carr-Purcell-Meiboom-Gill (CPMG) spin-echo refocusing trains counteract T2* effects. The effect of random phase on the refocusing train is eliminated using a stimulated echo diffusion preparation. Root-flipped Shinnar-Le Roux refocusing pulses permits preparation of a high spectral bandwidth, which improves imaging times by reducing the number of excitations required to cover the desired spectral range. B₁ sensitivity is reduced by using an excitation that satisfies the CPMG condition in the preparation. A method for ADC quantification insensitive to background gradients is presented.

RESULTS

Non-linear phase refocusing pulses reduces the peak B₁ by 46% which allows RF bandwidth to be doubled. Simulations and phantom experiments show that a non-linear phase CPMG pulse pair reduces B₁ sensitivity. Application in vivo demonstrates complementary contrast to conventional multispectral acquisitions and improved visualization compared to DW-EPI.

CONCLUSION

Volumetric and multispectral DW imaging near metal can be achieved with a 3D encoded sequence.

Gradient nonlinearity correction in liver DWI using motion-compensated diffusion encoding waveforms

OBJECTIVE : To experimentally characterize the effectiveness of a gradient nonlinearity correction method in removing ADC bias for different motion-compensated diffusion encoding waveforms.

METHODS

The diffusion encoding waveforms used were the standard monopolar Stejskal-Tanner pulsed gradient spin echo (pgse) waveform, the symmetric bipolar velocity-compensated waveform (sym-vc), the asymmetric bipolar velocity-compensated waveform (asym-vc) and the asymmetric bipolar partial velocity-compensated waveform (asym-pvc). The effectiveness of the gradient nonlinearity correction method using the spherical harmonic expansion of the gradient coil field was tested with the aforementioned waveforms in a phantom and in four healthy subjects.

RESULTS

The gradient nonlinearity correction method reduced the ADC bias in the phantom experiments for all used waveforms. The range of the ADC values over a distance of ± 67.2 mm from isocenter reduced from 1.29 × 10^-4 to 0.32 × 10^-4 mm²/s for pgse, 1.04 × 10^-4 to 0.22 × 10^-4 mm²/s for sym-vc, 1.22 × 10^-4 to 0.24 × 10^-4 mm²/s for asym-vc and 1.07 × 10^-4 to 0.11 × 10^-4 mm²/s for asym-pvc. The in vivo results showed that ADC overestimation due to motion or bright vessels can be increased even further by the gradient nonlinearity correction.

CONCLUSION

The investigated gradient nonlinearity correction method can be used effectively with various motion-compensated diffusion encoding waveforms. In coronal liver DWI, ADC errors caused by motion and residual vessel signal can be increased even further by the gradient nonlinearity correction.

Tractography methods and findings in brain tumors and traumatic brain injury

White matter fiber tracking using diffusion magnetic resonance imaging (dMRI) provides a noninvasive approach to map brain connections, but improving anatomical accuracy has been a significant challenge since the birth of tractography methods. Utilizing tractography in brain studies therefore requires understanding of its technical limitations to avoid shortcomings and pitfalls. This review explores tractography limitations and how different white matter pathways pose different challenges to fiber tracking methodologies. We summarize the pros and cons of commonly-used methods, aiming to inform how tractography and its related analysis may lead to questionable results. Extending these experiences, we review the clinical utilization of tractography in patients with brain tumors and traumatic brain injury, starting from tensor-based tractography to more advanced methods. We discuss current limitations and highlight novel approaches in the context of these two conditions to inform future tractography developments.

Mapping gradient nonlinearity and miscalibration using diffusion-weighted MR images of a uniform isotropic phantom

PURPOSE

To use diffusion measurements to map the spatial dependence of the magnetic field produced by the gradient coils of an MRI scanner with sufficient accuracy to correct errors in quantitative diffusion MRI (DMRI) caused by gradient nonlinearity and gradient amplifier miscalibration.

THEORY AND METHODS

The field produced by the gradient coils is expanded in regular solid harmonics. The expansion coefficients are found by fitting a model to a minimum set of diffusion-weighted images of an isotropic diffusion phantom. The accuracy of the resulting gradient coil field maps is evaluated by using them to compute corrected b-matrices that are then used to process a multi-shell diffusion tensor imaging (DTI) dataset with 32 diffusion directions per shell.

RESULTS

The method substantially reduces both the spatial inhomogeneity of the computed mean diffusivities (MD) and the computed values of the fractional anisotropy (FA), as well as virtually eliminating any artifactual directional bias in the tensor field secondary to gradient nonlinearity. When a small scaling miscalibration was purposely introduced in the x, y, and z, the method accurately detected the amount of miscalibration on each gradient axis.

CONCLUSION

The method presented detects and corrects the effects of gradient nonlinearity and gradient gain miscalibration using a simple isotropic diffusion phantom. The correction would improve the accuracy of DMRI measurements in the brain and other organs for both DTI and higher order diffusion analysis. In particular, it would allow calibration of MRI systems, improving data harmony in multicenter studies.

Whole-body diffusion-weighted MRI of normal lymph nodes: prospective apparent diffusion coefficient histogram and nodal distribution analysis in a healthy cohort

BACKGROUND

Whole body DWI (WB-DWI) enables the identification of lymph nodes for disease evaluation. However, quantitative data of benign lymph nodes across the body are lacking to allow meaningful comparison of diseased states. We evaluated apparent diffusion coefficient (ADC) histogram parameters of all visible lymph nodes in healthy volunteers on WB-DWI and compared differences in nodal ADC values between anatomical regions.

METHODS

WB-DWI was performed on a 1.5 T MR system in 20 healthy volunteers (7 female, 13 male, mean age 35 years). The b900 images were evaluated by two radiologists and all visible nodes from the neck to groin areas were segmented and individual nodal median ADC recorded. All segmented nodes in a patient were summated to generate the total nodal volume. Descriptors of the global ADC histogram, derived from individual node median ADCs, including mean, median, skewness and kurtosis were obtained for the global volume and each nodal region per patient. ADC values between nodal regions were compared using one-way ANOVA with Bonferroni post hoc tests and a p-value ≤0.05 was deemed statistically significant.

RESULTS

One thousand sixty-seven lymph nodes were analyzed. The global mean and median ADC of all lymph nodes were 1.12 ± 0.27 (10^- 3 mm²/s) and 1.09 (10^- 3 mm²/s). The average median ADC skewness was 0.25 ± 0.02 and average median ADC kurtosis was 0.34 ± 0.04. The ADC values of intrathoracic, portal and retroperitoneal nodes were significantly higher (1.53 × 10^- 3, 1.75 × 10^- 3 and 1.58 × 10^- 3 mm²/s respectively) than in other regions. Intrathoracic, portal and mesenteric nodes were relatively uncommon, accounting for only 3% of the total nodes segmented.

CONCLUSIONS

The global mean and median ADC of all lymph nodes were 1.12 ± 0.27 (10^- 3 mm²/s) and 1.09 (10^- 3 mm²/s). Intrathoracic, portal and retroperitoneal nodes display significantly higher ADCs. Normal intrathoracic, portal and mesenteric nodes are infrequently visualized on WB-DWI of healthy individuals.

TRIAL REGISTRATION

Royal Marsden Hospital committee for clinical research registration number 09/H0801/86, 19.10.2009.

Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group

On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T₁ . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.

Cross-scanner reproducibility and harmonization of a diffusion MRI structural brain network: A traveling subject study of multi-b acquisition

Characterization of brain networks by diffusion MRI (dMRI) has rapidly evolved, and there are ongoing movements toward data sharing and multi-center studies. To extract meaningful information from multi-center data, methods to correct for the bias caused by scanner differences, that is, harmonization, are urgently needed. In this work, we report the cross-scanner differences in structural network analyses using data from nine traveling subjects (four males and five females, 21-49 years-old) who underwent scanning using four 3T scanners (public database available from the Brain/MINDS Beyond Human Brain MRI project (http://mriportal.umin.jp/)). The reliability and reproducibility were compared to those of data from another set of four subjects (all males, 29-42 years-old) who underwent scan-rescan (interval, 105-147 days) with the same scanner as well as scan-rescan data from the Human Connectome Project database. The results demonstrated that the reliability of the edge weights and graph theory metrics was lower for data including different scanners, compared to the scan-rescan with the same scanner. Besides, systematic differences between scanners were observed, indicating the risk of bias in comparing networks obtained from different scanners directly. We further demonstrate that it is feasible to reduce inter-scanner variabilities while preserving the inter-subject differences among healthy individuals by modeling the scanner effects at the level of network matrices, when traveling-subject data are available for calibration between scanners. The present data and results are expected to serve as a basis for developing and evaluating novel harmonization methods.

Quantitative diffusion MRI of the abdomen and pelvis

Diffusion MRI has enormous potential and utility in the evaluation of various abdominal and pelvic disease processes including cancer and noncancer imaging of the liver, prostate, and other organs. Quantitative diffusion MRI is based on acquisitions with multiple diffusion encodings followed by quantitative mapping of diffusion parameters that are sensitive to tissue microstructure. Compared to qualitative diffusion-weighted MRI, quantitative diffusion MRI can improve standardization of tissue characterization as needed for disease detection, staging, and treatment monitoring. However, similar to many other quantitative MRI methods, diffusion MRI faces multiple challenges including acquisition artifacts, signal modeling limitations, and biological variability. In abdominal and pelvic diffusion MRI, technical acquisition challenges include physiologic motion (respiratory, peristaltic, and pulsatile), image distortions, and low signal-to-noise ratio. If unaddressed, these challenges lead to poor technical performance (bias and precision) and clinical outcomes of quantitative diffusion MRI. Emerging and novel technical developments seek to address these challenges and may enable reliable quantitative diffusion MRI of the abdomen and pelvis. Through systematic validation in phantoms, volunteers, and patients, including multicenter studies to assess reproducibility, these emerging techniques may finally demonstrate the potential of quantitative diffusion MRI for abdominal and pelvic imaging applications.

Retrospective Correction of ADC for Gradient Nonlinearity Errors in Multicenter Breast DWI Trials: ACRIN6698 Multiplatform Feasibility Study

The presented analysis of multisite, multiplatform clinical oncology trial data sought to enhance quantitative utility of the apparent diffusion coefficient (ADC) metric, derived from diffusion-weighted magnetic resonance imaging, by reducing technical interplatform variability owing to systematic gradient nonlinearity (GNL). This study tested the feasibility and effectiveness of a retrospective GNL correction (GNC) implementation for quantitative quality control phantom data, as well as in a representative subset of 60 subjects from the ACRIN 6698 breast cancer therapy response trial who were scanned on 6 different gradient systems. The GNL ADC correction based on a previously developed formalism was applied to trace-DWI using system-specific gradient-channel fields derived from vendor-provided spherical harmonic tables. For quantitative DWI phantom images acquired in typical breast imaging positions, the GNC improved interplatform accuracy from a median of 6% down to 0.5% and reproducibility of 11% down to 2.5%. Across studied trial subjects, GNC increased low ADC (<1 µm²/ms) tumor volume by 16% and histogram percentiles by 5%–8%, uniformly shifting percentile-dependent ADC thresholds by ∼0.06 µm²/ms. This feasibility study lays the grounds for retrospective GNC implementation in multiplatform clinical imaging trials to improve accuracy and reproducibility of ADC metrics used for breast cancer treatment response prediction.

Empirical validation of gradient field models for an accurate ADC measured on clinical 3T MR systems in body oncologic applications

PURPOSE

To empirically corroborate vendor-provided gradient nonlinearity (GNL) characteristics and demonstrate efficient GNL bias correction for human brain apparent diffusion coefficient (ADC) across 3T MR systems and spatial locations.

METHODS

Spatial distortion vector fields (DVF) were mapped in 3D using a surface fiducial array phantom for individual gradient channels on three 3T MR platforms from different vendors. Measured DVF were converted into empirical 3D GNL tensors and compared with their theoretical counterparts derived from vendor-provided spherical harmonic (SPH) coefficients. To illustrate spatial impact of GNL on ADC, diffusion weighted imaging using three orthogonal gradient directions was performed on a volunteer brain positioned at isocenter (as a reference) and offset superiorly by 10-17 cm (>10% predicted GNL bias). The SPH tensor-based GNL correction was applied to individual DWI gradient directions, and derived ADC was compared with low-bias reference for human brain white matter (WM) ROIs.

RESULTS

Empiric and predicted GNL errors were comparable for all three studied 3T MR systems, with <1.0% differences in the median and width of spatial histograms for individual GNL tensor elements. Median (±width) of ADC (10-3mm2/s) histograms measured at isocenter in WM reference ROIs from three MR systems were: 0.73 ± 0.11, 0.71 ± 0.14, 0.74 ± 0.17, and at off-isocenters (before versus after GNL correction) were respectively 0.63 ± 0.14 versus 0.72 ± 0.11, 0.53 ± 0.16 versus 0.74 ± 0.18, and 0.65 ± 0.16 versus 0.76 ± 0.18.

CONCLUSION

The phantom-based spatial distortion measurements validated vendor-provided gradient fields, and accurate WM ADC was recovered regardless of spatial locations and clinical MR platforms using system-specific tensor-based GNL correction for routine DWI.

In vivo human whole-brain Connectom diffusion MRI dataset at 760 µm isotropic resolution

We present a whole-brain in vivo diffusion MRI (dMRI) dataset acquired at 760 μm isotropic resolution and sampled at 1260 q-space points across 9 two-hour sessions on a single healthy participant. The creation of this benchmark dataset is possible through the synergistic use of advanced acquisition hardware and software including the high-gradient-strength Connectom scanner, a custom-built 64-channel phased-array coil, a personalized motion-robust head stabilizer, a recently developed SNR-efficient dMRI acquisition method, and parallel imaging reconstruction with advanced ghost reduction algorithm. With its unprecedented resolution, SNR and image quality, we envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance the understanding of human brain structures and connectivity. This comprehensive dataset can also be used as a test bed for new modeling, sub-sampling strategies, denoising and processing algorithms, potentially providing a common testing platform for further development of in vivo high resolution dMRI techniques. Whole brain anatomical T1-weighted and T2-weighted images at submillimeter scale along with field maps are also made available.

The effect of gradient nonlinearities on fiber orientation estimates from spherical deconvolution of diffusion magnetic resonance imaging data

Gradient nonlinearities in magnetic resonance imaging (MRI) cause spatially varying mismatches between the imposed and the effective gradients and can cause significant biases in rotationally invariant diffusion MRI measures derived from, for example, diffusion tensor imaging. The estimation of the orientational organization of fibrous tissue, which is nowadays frequently performed with spherical deconvolution techniques ideally using higher diffusion weightings, can likewise be biased by gradient nonlinearities. We explore the sensitivity of two established spherical deconvolution approaches to gradient nonlinearities, namely constrained spherical deconvolution (CSD) and damped Richardson-Lucy (dRL). Additionally, we propose an extension of dRL to take into account gradient imperfections, without the need of data interpolation. Simulations show that using the effective b-matrix can improve dRL fiber orientation estimation and reduces angular deviations, while CSD can be more robust to gradient nonlinearity depending on the implementation. Angular errors depend on a complex interplay of many factors, including the direction and magnitude of gradient deviations, underlying microstructure, SNR, anisotropy of the effective response function, and diffusion weighting. Notably, angular deviations can also be observed at lower b-values in contrast to the perhaps common assumption that only high b-value data are affected. In in vivo Human Connectome Project data and acquisitions from an ultrastrong gradient (300 mT/m) scanner, angular differences are observed between applying and not applying the effective gradients in dRL estimation. As even small angular differences can lead to error propagation during tractography and as such impact connectivity analyses, incorporating gradient deviations into the estimation of fiber orientations should make such analyses more reliable.

Empirical field mapping for gradient nonlinearity correction of multi-site diffusion weighted MRI

BACKGROUND

Achieving inter-site / inter-scanner reproducibility of diffusion weighted magnetic resonance imaging (DW-MRI) metrics has been challenging given differences in acquisition protocols, analysis models, and hardware factors.

PURPOSE

Magnetic field gradients impart scanner-dependent spatial variations in the applied diffusion weighting that can be corrected if the gradient nonlinearities are known. However, retrieving manufacturer nonlinearity specifications is not well supported and may introduce errors in interpretation of units or coordinate systems. We propose an empirical approach to mapping the gradient nonlinearities with sequences that are supported across the major scanner vendors.

STUDY TYPE

Prospective observational study.

SUBJECTS

A spherical isotropic diffusion phantom, and a single human control volunteer.

FIELD STRENGTH/SEQUENCE

3 T (two scanners). Stejskal-Tanner spin echo sequence with b-values of 1000, 2000 s/mm² with 12, 32, and 384 diffusion gradient directions per shell.

ASSESSMENT

We compare the proposed correction with the prior approach using manufacturer specifications against typical diffusion pre-processing pipelines (i.e., ignoring spatial gradient nonlinearities). In phantom data, we evaluate metrics against the ground truth. In human and phantom data, we evaluate reproducibility across scans, sessions, and hardware.

STATISTICAL TESTS

Wilcoxon rank-sum test between uncorrected and corrected data.

RESULTS

In phantom data, our correction method reduces variation in mean diffusivity across sessions over uncorrected data (p < 0.05). In human data, we show that this method can also reduce variation in mean diffusivity across scanners (p < 0.05).

CONCLUSION

Our method is relatively simple, fast, and can be applied retroactively. We advocate incorporating voxel-specific b-value and b-vector maps should be incorporated in DW-MRI harmonization preprocessing pipelines to improve quantitative accuracy of measured diffusion parameters.

Probing neural tissues at small scales: Recent progress of oscillating gradient spin echo (OGSE) neuroimaging in humans

The detection sensitivity of diffusion MRI (dMRI) is dependent on diffusion times. A shorter diffusion time can increase the sensitivity to smaller length scales. However, the conventional dMRI uses the pulse gradient spin echo (PGSE) sequence that probes relatively long diffusion times only. To overcome this, the oscillating gradient spin echo (OGSE) sequence has been developed to probe much shorter diffusion times with hardware limitations on preclinical and clinical MRI systems. The OGSE sequence has been previously used on preclinical animal MRI systems. Recently, several studies have translated the OGSE sequence to humans on clinical MRI systems and achieved new information that is invisible using conventional PGSE sequence. This paper provides an overview of the recent progress of the OGSE neuroimaging in humans, including the technical improvements in the translation of the OGSE sequence to human imaging and various applications in different neurological disorders and stroke. Some possible future directions of the OGSE sequence are also discussed.

A comparative study of gradient nonlinearity correction strategies for processing diffusion data obtained with ultra-strong gradient MRI scanners

Purpose

The analysis of diffusion data obtained under large gradient nonlinearities necessitates corrections during data reconstruction and analysis. While two such preprocessing pipelines have been proposed, no comparative studies assessing their performance exist. Furthermore, both pipelines neglect the impact of subject motion during acquisition, which, in the presence of gradient nonlinearities, induces spatio‐temporal B‐matrix variations. Here, spatio‐temporal B‐matrix tracking (STB) is proposed and its performance compared to established pipelines.

Methods

Diffusion tensor MRI (DT‐MRI) was performed using a 300 mT/m gradient system. Data were acquired with volunteers positioned in regions with pronounced gradient nonlinearities, and used to compare the performance of six different processing pipelines, including STB.

Results

Up to 30% errors were observed in DT‐MRI parameter estimates when neglecting gradient nonlinearities. Moreover, the order in which B0 inhomogeneity, eddy current and gradient nonlinearity corrections were performed was found to impact the consistency of parameter estimates significantly. Although, no pipeline emerged as a clear winner, the STB approach seemed to yield the most consistent parameter estimates under large gradient nonlinearities.

Conclusions

Under large gradient nonlinearities, the choice of preprocessing pipeline significantly impacts the estimated diffusion parameters. Motion‐induced spatio‐temporal B‐matrix variations can lead to systematic bias in the parameter estimates, that can be ameliorated using the proposed STB framework.

Click here for author‐reader discussions

Assessment of the systematic errors caused by diffusion gradient inhomogeneity in DTI-computer simulations

Diffusion tensor imaging (DTI) is a powerful MRI modality that allows the investigation of the microstructure of tissues both in vivo and noninvasively. Its reliability is strictly dependent on the performance of diffusion-sensitizing gradients, of which spatial nonuniformity is a known issue in the case of virtually all clinical MRI scanners. The influence of diffusion gradient inhomogeneity on the accuracy of the diffusion tensor imaging was investigated by means of computer simulations supported by an MRI experiment performed at the isocenter and 15 cm away. The DTI measurements of two diffusion phantoms were simulated assuming a nonuniform diffusion-sensitizing gradient and various levels of noise. Thereafter, the tensors were calculated by two methods: (i) assuming a spatially constant b-matrix (standard DTI) and (ii) applying the b-matrix spatial distribution in the DTI (BSD-DTI) technique, a method of indicating the b-matrix for each voxel separately using an anisotropic phantom as a standard of diffusion. The average eigenvalues and fractional anisotropy across the homogeneous region of interest were calculated and compared with the expected values. Diffusion gradient inhomogeneity leads to overestimation of the largest eigenvalue, underestimation of the smallest one and thus overestimation of fractional anisotropy. The effect is similar to that caused by noise; however, it could not be corrected by increasing SNR. The MRI measurements, performed using a 3 T clinical scanner, revealed that the split of the eigenvalues measured 15 cm away from the isocenter is significant (up to 25%). The BSD-DTI calibration allowed the reduction of the measured fractional anisotropy of the isotropic medium from 0.174 to 0.031, suggesting that gradient inhomogeneity was the main cause of this error. For the phantom measured at the isocenter, however, the split was almost not observed; the average eigenvalues were shifted from the expected value by ~ 5%.

Diffusion gradient nonlinearity bias correction reduces bias of breast cancer bone metastasis ADC values

CONTRACT GRANT SPONSOR

Health Research Fund of Central Denmark Region.

BACKGROUND

Diffusion gradient nonlinearity (DGNL) bias causes apparent diffusion coefficient (ADC) values to drop with increasing superior-inferior (SI) isocenter offset. This is a concern when performing quantitative diffusion-weighted imaging (DWI).

PURPOSE/HYPOTHESIS

To investigate if DGNL ADC bias can be corrected in breast cancer bone metastases using a clinical DWI protocol and an online correction algorithm.

STUDY TYPE

Prospective.

SUBJECTS/PHANTOM

A diffusion phantom (Model 128, High Precision Devices, Boulder, CO) was used for in vitro validation. Twenty-three women with bone-metastasizing breast cancer were enrolled to assess DGNL correction in vivo.

FIELD STRENGTH/SEQUENCE

DWI was performed on a 1.5T MRI system as single-shot, spin-echo, echo-planar imaging with short-tau inversion recovery (STIR) fat-saturation. ADC maps with and without DGNL correction were created from the b50 and b800 images.

ASSESSMENT

Uncorrected and DGNL-corrected ADC values were measured in phantom and bone metastases by placing regions of interest on b800 images and copying them to the ADC map. The SI offset was recorded.

STATISTICAL TESTS

In all, 79 bone metastases were assessed. ADC values with and without DGNL correction were compared at 14 cm SI offset using a two-tailed t-test.

RESULTS

In the diffusion phantom, DGNL correction increased SI offset, where ADC bias was lower than 5%, from 7.3-13.8 cm. Of the 23 patients examined, six had no metastases in the covered regions. In the remaining patients, bias of uncorrected bone metastasis ADC values was 19.1% (95% confidence interval [CI]: 15.4-22.9%) at 14 cm SI offset. After DGNL correction, ADC bias was significantly reduced to 3.5% (95% CI: 0.7-6.3%, P < 0.001), thus reducing bias due to DGNL by 82%.

DATA CONCLUSION

Online DGNL correction corrects DGNL ADC value bias and allows increased station lengths in the SI direction.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:904-911.

Distortion-free imaging: A double encoding method (DIADEM) combined with multiband imaging for rapid distortion-free high-resolution diffusion imaging on a compact 3T with high-performance gradients

BACKGROUND

Distortion-free, high-resolution diffusion imaging using DIADEM (Distortion-free Imaging: A Double Encoding Method), proposed recently, has great potential for clinical applications. However, it can suffer from prolonged scan times and its reliability for quantitative diffusion imaging has not been evaluated.

PURPOSE

To investigate the clinical feasibility of DIADEM-based high-resolution diffusion imaging on a novel compact 3T (C3T) by evaluating the reliability of quantitative diffusion measurements and utilizing both the high-performance gradients (80 mT/m, 700 T/m/s) and the sequence optimization with the navigator acquisition window reduction and simultaneous multislice (multiband) imaging.

STUDY TYPE

Prospective feasibility study.

PHANTOM/SUBJECTS

Diffusion quality control phantom scans to evaluate the reliability of quantitative diffusion measurements; 36 normal control scans for B₀ -field mapping; six healthy and two patient subject scans with a brain tumor for comparisons of diffusion and anatomical imaging.

FIELD STRENGTH/SEQUENCE

3T; the standard single-shot echo-planar-imaging (EPI), multishot DIADEM diffusion, and anatomical (2D-FSE [fast-spin-echo], 2D-FLAIR [fluid-attenuated-inversion-recovery], and 3D-MPRAGE [magnetization prepared rapid acquisition gradient echo]) imaging.

ASSESSMENT

The scan time reduction, the reliability of quantitative diffusion measurements, and the clinical efficacy for high-resolution diffusion imaging in healthy control and brain tumor volunteers.

STATISTICAL TEST

Bland-Altman analysis.

RESULTS

The scan time for high in-plane (0.86 mm² ) resolution, distortion-free, and whole brain diffusion imaging were reduced from 10 to 5 minutes with the sequence optimizations. All of the mean apparent diffusion coefficient (ADC) values in phantom were within the 95% confidence interval in the Bland-Altman plot. The proposed acquisition with a total off-resonance coverage of 597.2 Hz wider than the expected bandwidth of 500 Hz in human brain could yield a distortion-free image without foldover artifacts. Compared with EPI, therefore, this approach allowed direct image matching with the anatomical images and enabled improved delineation of the tumor boundaries.

DATA CONCLUSION

The proposed high-resolution diffusion imaging approach is clinically feasible on C3T due to a combination of hardware and sequence improvements.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1 J. Magn. Reson. Imaging 2020;51:296-310.

Improving apparent diffusion coefficient accuracy on a compact 3T MRI scanner using gradient nonlinearity correction

BACKGROUND

Gradient nonlinearity (GNL) leads to biased apparent diffusion coefficients (ADCs) in diffusion-weighted imaging. A gradient nonlinearity correction (GNLC) method has been developed for whole body systems, but is yet to be tested for the new compact 3T (C3T) scanner, which exhibits more complex GNL due to its asymmetrical design.

PURPOSE

To assess the improvement of ADC quantification with GNLC for the C3T scanner.

STUDY TYPE

Phantom measurements and retrospective analysis of patient data.

PHANTOM/SUBJECTS

A diffusion quality control phantom with vials containing 0-30% polyvinylpyrrolidone in water was used. For in vivo data, 12 patient exams were analyzed (median age, 33).

FIELD STRENGTH/SEQUENCE

Imaging was performed on the C3T and two commercial 3T scanners. A clinical DWI (repetition time [TR] = 10,000 msec, echo time [TE] = minimum, b = 1000 s/mm2 ) sequence was used for phantom imaging and 10 patient cases and a clinical DTI (TR = 6000-10,000 msec, TE = minimum, b = 1000 s/mm2 ) sequence was used for two patient cases.

ASSESSMENT

The 0% vial was measured along three orthogonal axes, and at two different temperatures. The ADC for each concentration was compared between the C3T and two whole-body scanners. Cerebrospinal fluid and white matter ADCs were quantified for each patient and compared to values in literature.

STATISTICAL TESTS

Paired t-test and two-way analysis of variance (ANOVA).

RESULTS

For all PVP concentrations, the corrected ADC was within 2.5% of the reference ADC. On average, the ADC of cerebrospinal fluid and white matter post-GNLC were within 1% and 6%, respectively, of values reported in the literature and were significantly different from the uncorrected data (P < 0.05).

DATA CONCLUSION

This study demonstrated that GNL effects were more severe for the C3T due to the asymmetric gradient design, but our implementation of a GNLC compensated for these effects, resulting in ADC values that are in good agreement with values from the literature.

LEVEL OF EVIDENCE

4 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;48:1498-1507.

Microstructural imaging of the human brain with a 'super-scanner': 10 key advantages of ultra-strong gradients for diffusion MRI

The key component of a microstructural diffusion MRI 'super-scanner' is a dedicated high-strength gradient system that enables stronger diffusion weightings per unit time compared to conventional gradient designs. This can, in turn, drastically shorten the time needed for diffusion encoding, increase the signal-to-noise ratio, and facilitate measurements at shorter diffusion times. This review, written from the perspective of the UK National Facility for In Vivo MR Imaging of Human Tissue Microstructure, an initiative to establish a shared 300 mT/m-gradient facility amongst the microstructural imaging community, describes ten advantages of ultra-strong gradients for microstructural imaging. Specifically, we will discuss how the increase of the accessible measurement space compared to a lower-gradient systems (in terms of Δ, b-value, and TE) can accelerate developments in the areas of 1) axon diameter distribution mapping; 2) microstructural parameter estimation; 3) mapping micro-vs macroscopic anisotropy features with gradient waveforms beyond a single pair of pulsed-gradients; 4) multi-contrast experiments, e.g. diffusion-relaxometry; 5) tractography and high-resolution imaging in vivo and 6) post mortem; 7) diffusion-weighted spectroscopy of metabolites other than water; 8) tumour characterisation; 9) functional diffusion MRI; and 10) quality enhancement of images acquired on lower-gradient systems. We finally discuss practical barriers in the use of ultra-strong gradients, and provide an outlook on the next generation of 'super-scanners'.

Analysis and correction of errors in DTI-based tractography due to diffusion gradient inhomogeneity

The DTI-based tractography, despite its restrictions, is the most widely utilized fiber tracking method in clinical practice. Its fidelity is strictly dependent on the precision and accuracy of the DTI measurement, which in turn is limited by the linearity of the diffusion sensitizing gradient. The influence of the gradient distortions on the differences between the real and measured orientation of fibers was investigated by computer simulations. In addition, the potential of the b-matrix Spatial Distribution in DTI (BSD-DTI) technique in correcting such kind of errors was demonstrated experimentally. The simulations revealed that the diffusion gradient inhomogeneity, if not corrected, leads to the erroneous indication of the fiber direction. The average and maximum deviations were about 1° and 15°, respectively. Remarkably, the deviation between the real and measured orientation of fibers is directionally dependent, what was confirmed in MRI measurement. The deviation errors can be effectively corrected by preceding the DTI measurement with the BSD-DTI calibration.

The generalized Stejskal-Tanner equation for non-uniform magnetic field gradients

The intensity of the diffusion weighted NMR signal is described by the Stejskal-Tanner equation, which was derived under the assumption that the gradients are uniform throughout the sample. Nevertheless, it has been demonstrated numerous times that this condition is not fulfilled in the cases of virtually any clinical or research MRI scanners. Therefore, technically, the Stejskal-Tanner equation is valid only for a very specific case of homogeneous gradients. In this paper the Stejskal-Tanner equation was derived for the general case on non-uniform diffusion gradients. To this end, the magnetic field was expressed as linear in a curvilinear coordinate system defined by a vector function p(r). Thereafter, the expression for the non-linear magnetic field was put into the Bloch-Torrey equation and solved. Moreover, the meaning of so-called coil tensor, which is used for the gradients inhomogeneity correction, was explained. It was proven that in the case of the spin echo-based sequences, the Stejskal-Tenner equation is still valid, even if the diffusion gradients are non-uniform. However, in such a case, the b-matrix should be derived for each voxel separately. For other sequence, the derived relation possesses an imaginary term, which corresponds do the phase shift of the diffusion weighted signal.

Dependence of apparent diffusion coefficient measurement on diffusion gradient direction and spatial position - A quality assurance intercomparison study of forty-four scanners for quantitative diffusion-weighted imaging

PURPOSE

To propose an MRI quality assurance procedure that can be used for routine controls and multi-centre comparison of different MR-scanners for quantitative diffusion-weighted imaging (DWI).

MATERIALS AND METHODS

44 MR-scanners with different field strengths (1 T, 1.5 T and 3 T) were included in the study. DWI acquisitions (b-value range 0-1000 s/mm²), with three different orthogonal diffusion gradient directions, were performed for each MR-scanner. All DWI acquisitions were performed by using a standard spherical plastic doped water phantom. Phantom solution ADC value and its dependence with temperature was measured using a DOSY sequence on a 600 MHz NMR spectrometer. Apparent diffusion coefficient (ADC) along each diffusion gradient direction and mean ADC were estimated, both at magnet isocentre and in six different position 50 mm away from isocentre, along positive and negative AP, RL and HF directions.

RESULTS

A good agreement was found between the nominal and measured mean ADC at isocentre: more than 90% of mean ADC measurements were within 5% from the nominal value, and the highest deviation was 11.3%. Away from isocentre, the effect of the diffusion gradient direction on ADC estimation was larger than 5% in 47% of included scanners and a spatial non uniformity larger than 5% was reported in 13% of centres.

CONCLUSION

ADC accuracy and spatial uniformity can vary appreciably depending on MR scanner model, sequence implementation (i.e. gradient diffusion direction) and hardware characteristics. The DWI quality assurance protocol proposed in this study can be employed in order to assess the accuracy and spatial uniformity of estimated ADC values, in single- as well as multi-centre studies.

Phantom-based field maps for gradient nonlinearity correction in diffusion imaging

Gradient coils in magnetic resonance imaging do not produce perfectly linear gradient fields. For diffusion imaging, the field nonlinearities cause the amplitude and direction of the applied diffusion gradients to vary over the field of view. This leads to site- and scan-specific systematic errors in estimated diffusion parameters such as diffusivity and anisotropy, reducing reliability especially in studies that take place over multiple sites. These errors can be substantially reduced if the actual scanner-specific gradient coil magnetic fields are known. The nonlinearity of the coil fields is measured by scanner manufacturers and used internally for geometric corrections, but obtaining and using the information for a specific scanner may be impractical for many sites that operate without special-purpose local engineering and research support. We have implemented an empirical field-mapping procedure using a large phantom combined with a solid harmonic approximation to the coil fields that is simple to perform and apply. Here we describe the accuracy and precision of the approach in reproducing manufacturer gold standard field maps and in reducing spatially varying errors in quantitative diffusion imaging for a specific scanner. Before correction, median B value error ranged from 33 - 41 relative to manufacturer specification at 100 mm from isocenter; correction reduced this to 0 - 4. On-axis spatial variation in the estimated mean diffusivity of an isotropic phantom was 2.2% - 4.1% within 60 mm of isocenter before correction, 0.5% - 1.6% after. Expected fractional anisotropy in the phantom was 0; highest estimated fractional anisotropy within 60 mm of isocenter was reduced from 0.024 to 0.012 in the phase encoding direction (48% reduction) and from 0.020 to 0.006 in the frequency encoding direction (72% reduction).

How reliable are ADC measurements? A phantom and clinical study of cervical lymph nodes

Objective

To assess the reliability of ADC measurements in vitro and in cervical lymph nodes of healthy volunteers.

Methods

We used a GE 1.5 T MRI scanner and a first ice-water phantom according to recommendations released by the Quantitative Imaging Biomarker Alliance (QIBA) for assessing ADC against reference values. We analysed the target size effect by using a second phantom made of six inserted spheres with diameters ranging from 10 to 37 mm. Thirteen healthy volunteers were also scanned to assess the inter- and intra-observer reproducibility of volumetric ADC measurements of cervical lymph nodes.

Results

On the ice-water phantom, the error in ADC measurements was less than 4.3 %. The spatial bias due to the non-linearity of gradient fields was found to be 24 % at 8 cm from the isocentre. ADC measure reliability decreased when addressing small targets due to partial volume effects (up to 12.8 %). The mean ADC value of cervical lymph nodes was 0.87.10^-3 ± 0.12.10^-3 mm²/s with a good intra-observer reliability. Inter-observer reproducibility featured a bias of -5.5 % due to segmentation issues.

Conclusion

ADC is a potentially important imaging biomarker in oncology; however, variability issues preclude its broader adoption. Reliable use of ADC requires technical advances and systematic quality control.

Key Points

• ADC is a promising quantitative imaging biomarker.

• ADC has a fair inter-reader variability and good intra-reader variability.

• Partial volume effect, post-processing software and non-linearity of scanners are limiting factors.

• No threshold values for detecting cervical lymph node malignancy can be drawn.

Electronic supplementary material

The online version of this article (10.1007/s00330-017-5265-2) contains supplementary material, which is available to authorized users.

Stability of Gradient Field Corrections for Quantitative Diffusion MRI

In magnetic resonance diffusion imaging, gradient nonlinearity causes significant bias in the estimation of quantitative diffusion parameters such as diffusivity, anisotropy, and diffusion direction in areas away from the magnet isocenter. This bias can be substantially reduced if the scanner- and coil-specific gradient field nonlinearities are known. Using a set of field map calibration scans on a large (29 cm diameter) phantom combined with a solid harmonic approximation of the gradient fields, we predicted the obtained b-values and applied gradient directions throughout a typical field of view for brain imaging for a typical 32-direction diffusion imaging sequence. We measured the stability of these predictions over time. At 80 mm from scanner isocenter, predicted b-value was 1-6% different than intended due to gradient nonlinearity, and predicted gradient directions were in error by up to 1 degree. Over the course of one month the change in these quantities due to calibration-related factors such as scanner drift and variation in phantom placement was <0.5% for b-values, and <0.5 degrees for angular deviation. The proposed calibration procedure allows the estimation of gradient nonlinearity to correct b-values and gradient directions ahead of advanced diffusion image processing for high angular resolution data, and requires only a five-minute phantom scan that can be included in a weekly or monthly quality assurance protocol.

QIN DAWG Validation of Gradient Nonlinearity Bias Correction Workflow for Quantitative Diffusion-Weighted Imaging in Multicenter Trials

Previous research has shown that system-dependent gradient nonlinearity (GNL) introduces a significant spatial bias (nonuniformity) in apparent diffusion coefficient (ADC) maps. Here, the feasibility of centralized retrospective system-specific correction of GNL bias for quantitative diffusion-weighted imaging (DWI) in multisite clinical trials is demonstrated across diverse scanners independent of the scanned object. Using corrector maps generated from system characterization by ice-water phantom measurement completed in the previous project phase, GNL bias correction was performed for test ADC measurements from an independent DWI phantom (room temperature agar) at two offset locations in the bore. The precomputed three-dimensional GNL correctors were retrospectively applied to test DWI scans by the central analysis site. The correction was blinded to reference DWI of the agar phantom at magnet isocenter where the GNL bias is negligible. The performance was evaluated from changes in ADC region of interest histogram statistics before and after correction with respect to the unbiased reference ADC values provided by sites. Both absolute error and nonuniformity of the ADC map induced by GNL (median, 12%; range, -35% to +10%) were substantially reduced by correction (7-fold in median and 3-fold in range). The residual ADC nonuniformity errors were attributed to measurement noise and other non-GNL sources. Correction of systematic GNL bias resulted in a 2-fold decrease in technical variability across scanners (down to site temperature range). The described validation of GNL bias correction marks progress toward implementation of this technology in multicenter trials that utilize quantitative DWI.

Correction of Gradient Nonlinearity Bias in Quantitative Diffusion Parameters of Renal Tissue with Intra Voxel Incoherent Motion

Spatially nonuniform diffusion weighting bias as a result of gradient nonlinearity (GNL) causes substantial errors in apparent diffusion coefficient (ADC) maps for anatomical regions imaged distant from the magnet isocenter. Our previously described approach effectively removed spatial ADC bias from 3 orthogonal diffusion-weighted imaging (DWI) measurements for monoexponential media of arbitrary anisotropy. This work evaluates correction feasibility and performance for quantitative diffusion parameters of the 2-component intravoxel incoherent motion (IVIM) model for well-perfused and nearly isotropic renal tissue. Sagittal kidney DWI scans of a volunteer were performed on a clinical 3T magnetic resonance imaging scanner near isocenter and offset superiorly. Spatially nonuniform diffusion weighting caused by GNL resulted both in shifting and broadening of perfusion-suppressed ADC histograms for off-center DWI relative to unbiased measurements close to the isocenter. Direction-average diffusion weighting bias correctors were computed based on the known gradient design provided by the vendor. The computed bias maps were empirically confirmed by coronal DWI measurements for an isotropic gel-flood phantom. Both phantom and renal tissue ADC bias for off-center measurements was effectively removed by applying precomputed 3D correction maps. Comparable ADC accuracy was achieved for corrections of both b maps and DWI intensities in the presence of IVIM perfusion. No significant bias impact was observed for the IVIM perfusion fraction.

Practical estimate of gradient nonlinearity for implementation of apparent diffusion coefficient bias correction

PURPOSE

To describe an efficient procedure to empirically characterize gradient nonlinearity and correct for the corresponding apparent diffusion coefficient (ADC) bias on a clinical magnetic resonance imaging (MRI) scanner.

MATERIALS AND METHODS

Spatial nonlinearity scalars for individual gradient coils along superior and right directions were estimated via diffusion measurements of an isotropicic e-water phantom. Digital nonlinearity model from an independent scanner, described in the literature, was rescaled by system-specific scalars to approximate 3D bias correction maps. Correction efficacy was assessed by comparison to unbiased ADC values measured at isocenter.

RESULTS

Empirically estimated nonlinearity scalars were confirmed by geometric distortion measurements of a regular grid phantom. The applied nonlinearity correction for arbitrarily oriented diffusion gradients reduced ADC bias from 20% down to 2% at clinically relevant offsets both for isotropic and anisotropic media. Identical performance was achieved using either corrected diffusion-weighted imaging (DWI) intensities or corrected b-values for each direction in brain and ice-water. Direction-average trace image correction was adequate only for isotropic medium.

CONCLUSION

Empiric scalar adjustment of an independent gradient nonlinearity model adequately described DWI bias for a clinical scanner. Observed efficiency of implemented ADC bias correction quantitatively agreed with previous theoretical predictions and numerical simulations. The described procedure provides an independent benchmark for nonlinearity bias correction of clinical MRI scanners.

Demonstration of nonlinearity bias in the measurement of the apparent diffusion coefficient in multicenter trials

PURPOSE

Characterize system-specific bias across common magnetic resonance imaging (MRI) platforms for quantitative diffusion measurements in multicenter trials.

METHODS

Diffusion weighted imaging (DWI) was performed on an ice-water phantom along the superior-inferior (SI) and right-left (RL) orientations spanning ± 150 mm. The same scanning protocol was implemented on 14 MRI systems at seven imaging centers. The bias was estimated as a deviation of measured from known apparent diffusion coefficient (ADC) along individual DWI directions. The relative contributions of gradient nonlinearity, shim errors, imaging gradients, and eddy currents were assessed independently. The observed bias errors were compared with numerical models.

RESULTS

The measured systematic ADC errors scaled quadratically with offset from isocenter, and ranged between -55% (SI) and 25% (RL). Nonlinearity bias was dependent on system design and diffusion gradient direction. Consistent with numerical models, minor ADC errors (± 5%) due to shim, imaging and eddy currents were mitigated by double echo DWI and image coregistration of individual gradient directions.

CONCLUSION

The analysis confirms gradient nonlinearity as a major source of spatial DW bias and variability in off-center ADC measurements across MRI platforms, with minor contributions from shim, imaging gradients and eddy currents. The developed protocol enables empiric description of systematic bias in multicenter quantitative DWI studies.

Gradient nonlinearity correction to improve apparent diffusion coefficient accuracy and standardization in the american college of radiology imaging network 6698 breast cancer trial

PURPOSE

To evaluate a gradient nonlinearity correction (GNC) program for quantitative apparent diffusion coefficient (ADC) measurements on phantom and human subject diffusion-weighted (DW) magnetic resonance imaging (MRI) scans in a multicenter breast cancer treatment response study

MATERIALS AND METHODS

A GNC program using fifth-order spherical harmonics for gradient modeling was applied retrospectively to qualification phantom and human subject scans. Ice-water phantoms of known diffusion coefficient were scanned at five different study centers with different scanners and receiver coils. Human in vivo data consisted of baseline and early-treatment exams on 54 patients from four sites. ADC maps were generated with and without GNC. Regions of interest were defined to quantify absolute errors and changes with GNC over breast imaging positions.

RESULTS

Phantom ADC errors varied with region of interest (ROI) position and scanner configuration; the mean error by configuration ranged from 1.4% to 19.9%. GNC significantly reduced the overall mean error for all sites from 9.9% to 0.6% (P = 0.016). Spatial dependence of GNC was highest in the right-left (RL) and anterior-posterior (AP) directions. Human subject mean tumor ADC was reduced 0.2 to 12% by GNC at different sites. By regression, every 1-cm change in tumor ROI position between baseline and follow-up visits resulted in an estimated change of 2.4% in the ADC early-treatment response measurement.

CONCLUSION

GNC is effective for removing large, system-dependent errors in quantitative breast DWI. GNC may be important in ensuring reproducibility in multicenter studies and in reducing errors in longitudinal treatment response measures arising from spatial variations in tumor position between visits.

Pediatric brain tumor consortium multisite assessment of apparent diffusion coefficient z-axis variation assessed with an ice-water phantom

RATIONALE AND OBJECTIVES

Magnetic resonance diffusion imaging can characterize physiologic characteristics of pediatric brain tumors used to assess therapy response. The purpose of this study was to assess the variability of the apparent diffusion coefficient (ADC) along z-axis of scanners in the multicenter Pediatric Brain Tumor Consortium (PBTC).

MATERIALS AND METHODS

Ice-water diffusion phantoms for each PBTC site were distributed with a specific diffusion imaging protocol. The phantom was scanned four successive times to 1) confirm water in the tube reached thermal equilibrium and 2) allow for assessment of intra-examination ADC repeatability. ADC profiles across slice positions for each vendor and institution combination were characterized using linear regression modeling with a quadratic fit.

RESULTS

Eleven sites collected data with a high degree of compliance to the diffusion protocol for each scanner. The mean ADC value at slice position zero for vendor A was 1.123 × 10(-3) mm(2)/s, vendor B was 1.0964 × 10(-3) mm(2)/s, and vendor C was 1.110 × 10(-3) mm(2)/s. The percentage coefficient of variation across all sites was 0.309% (standard deviation = 0.322). The ADC values conformed well to a second-order polynomial along the z-axis, (ie, following a linear model pattern with quadratic fit) for vendor-institution combinations and across vendor-institution combinations as shown in the longitudinal model.

CONCLUSIONS

Assessment of the variability of diffusion metrics is essential for establishing the validity of using these quantitative metrics in multicenter trials. The low variability in ADC values across vendors and institutions and validates the use of ADC as a quantitative tumor marker in pediatric multicenter trials.

Image registration for quantitative parametric response mapping of cancer treatment response

Imaging biomarkers capable of early quantification of tumor response to therapy would provide an opportunity to individualize patient care. Image registration of longitudinal scans provides a method of detecting treatment associated changes within heterogeneous tumors by monitoring alterations in the quantitative value of individual voxels over time, which is unattainable by traditional volumetric-based histogram methods. The concepts involved in the use of image registration for tracking and quantifying breast cancer treatment response using parametric response mapping (PRM), a voxel-based analysis of diffusion-weighted magnetic resonance imaging (DW-MRI) scans, are presented. Application of PRM to breast tumor response detection is described, wherein robust registration solutions for tracking small changes in water diffusivity in breast tumors during therapy are required. Methodologies that employ simulations are presented for measuring expected statistical accuracy of PRM for response assessment. Test-retest clinical scans are used to yield estimates of system noise to indicate significant changes in voxel-based changes in water diffusivity. Overall, registration-based PRM image analysis provides significant opportunities for voxel-based image analysis to provide the required accuracy for early assessment of response to treatment in breast cancer patients receiving neoadjuvant chemotherapy.