Motion-compensated gradient waveforms for tensor-valued diffusion encoding by constrained numerical optimization

In vivo diffusion MRI of the human heart using a 300 mT/m gradient system

PURPOSE

This work reports for the first time on the implementation and application of cardiac diffusion-weighted MRI on a Connectom MR scanner with a maximum gradient strength of 300 mT/m. It evaluates the benefits of the increased gradient performance for the investigation of the myocardial microstructure.

METHODS

Cardiac diffusion-weighted imaging (DWI) experiments were performed on 10 healthy volunteers using a spin-echo sequence with up to second- and third-order motion compensation (M 2 $$ {M}_2 $$ andM 3 $$ {M}_3 $$ ) andb = 100 ,   450 $$ b=100,450 $$ , and 1000s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ (twice theb max $$ {b}_{\mathrm{max}} $$ commonly used on clinical scanners). Mean diffusivity (MD), fractional anisotropy (FA), helix angle (HA), and secondary eigenvector angle (E2A) were calculated for b = [100, 450]s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ and b = [100, 1000]s / m m 2 $$ \mathrm{s}/\mathrm{m}{\mathrm{m}}^2 $$ for bothM 2 $$ {M}_2 $$ andM 3 $$ {M}_3 $$ .

RESULTS

The MD values withM 3 $$ {M}_3 $$ are slightly higher than withM 2 $$ {M}_2 $$ withΔ MD = 0 . 05 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 5 ) $$ \Delta \mathrm{MD}=0.05\pm 0.05\kern0.3em \left[\times 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-5\right) $$ forb max = 450 s / mm 2 $$ {b}_{\mathrm{max}}=450\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ andΔ MD = 0 . 03 ± 0 . 03 [ × 1 0 - 3 mm 2 / s ] ( p = 4 e - 4 ) $$ \Delta \mathrm{MD}=0.03\pm 0.03\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=4e-4\right) $$ forb max = 1000 s / mm 2 $$ {b}_{\mathrm{max}}=1000\kern0.3em \mathrm{s}/{\mathrm{mm}}^2 $$ . A reduction in MD is observed by increasing theb max $$ {b}_{\mathrm{max}} $$ from 450 to 1000s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ (Δ MD = 0 . 06 ± 0 . 04 [ × 1 0 - 3 mm 2 / s ] ( p = 1 . 6 e - 9 ) $$ \Delta \mathrm{MD}=0.06\pm 0.04\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1.6e-9\right) $$ forM 2 $$ {M}_2 $$ andΔ MD = 0 . 08 ± 0 . 05 [ × 1 0 - 3 mm 2 / s ] ( p = 1 e - 9 ) $$ \Delta \mathrm{MD}=0.08\pm 0.05\kern0.3em \left[\times \kern0.3em 1{0}^{-3}\kern0.3em {\mathrm{mm}}^2/\mathrm{s}\right]\kern0.3em \left(p=1e-9\right) $$ forM 3 $$ {M}_3 $$ ). The difference between FA, E2A, and HA was not significant in different schemes (p > 0 . 05 $$ p>0.05 $$ ).

CONCLUSION

This work demonstrates cardiac DWI in vivo with higher b-value and higher order of motion compensated diffusion gradient waveforms than is commonly used. Increasing the motion compensation order fromM 2 $$ {M}_2 $$ toM 3 $$ {M}_3 $$ and the maximum b-value from 450 to 1000  s / mm 2 $$ \mathrm{s}/{\mathrm{mm}}^2 $$ affected the MD values but FA and the angular metrics (HA and E2A) remained unchanged. Our work paves the way for cardiac DWI on the next-generation MR scanners with high-performance gradient systems.

Motion-induced phase-corrected homodyne reconstruction for partial Fourier single-shot diffusion-weighted echo planar imaging of the liver

Partial Fourier encoding is popular in single-shot (ss) diffusion-weighted (DW) echo planar imaging (EPI) because it enables a shorter echo time (TE) and, hence, improves the signal-to-noise-ratio. Motion during diffusion encoding causes k-space shifting and dispersion, which compromises the quality of the homodyne reconstruction. This work provides a comprehensive understanding of the artifacts in homodyne reconstruction of partial Fourier ss-DW-EPI data in the presence of motion-induced phase and proposes the motion-induced phase-corrected homodyne (mpc-hdyne) reconstruction method to ameliorate these artifacts. Simulations with different types of motion-induced phase were performed to provide an understanding of the potential artifacts that occur in the homodyne reconstruction of partial Fourier ss-DW-EPI data. To correct for the artifacts, the mpc-hdyne reconstruction is proposed. The algorithm recenters k-space, updates the partial Fourier factor according to detected global k-space shifts, and removes low-resolution nonlinear phase before the conventional homodyne reconstruction. The mpc-hdyne reconstruction is tested on both simulation and in vivo data. Motion-induced phase can cause signal overestimation, worm artifacts, and signal loss in partial Fourier ss-DW-EPI data with the conventional homodyne reconstruction. Simulation and in vivo data showed that the proposed mpc-hdyne reconstruction ameliorated artifacts, yielding higher quality DW images compared with conventional homodyne reconstruction. Based on the understanding of the artifacts in homodyne reconstruction of partial Fourier ss-DW-EPI data, the mpc-hdyne reconstruction was proposed and showed superior performance compared with the conventional homodyne reconstruction on both simulation and in vivo data.

Partial Fourier in the presence of respiratory motion in prostate diffusion-weighted echo planar imaging

PURPOSE

To investigate the effect of respiratory motion in terms of signal loss in prostate diffusion-weighted imaging (DWI), and to evaluate the usage of partial Fourier in a free-breathing protocol in a clinically relevant b-value range using both single-shot and multi-shot acquisitions.

METHODS

A controlled breathing DWI acquisition was first employed at 3 T to measure signal loss from deep breathing patterns. Single-shot and multi-shot (2-shot) acquisitions without partial Fourier (no pF) and with partial Fourier (pF) factors of 0.75 and 0.65 were employed in a free-breathing protocol. The apparent SNR and ADC values were evaluated in 10 healthy subjects to measure if low pF factors caused low apparent SNR or overestimated ADC.

RESULTS

Controlled breathing experiments showed a difference in signal coefficient of variation between shallow and deep breathing. In free-breathing single-shot acquisitions, the pF 0.65 scan showed a significantly (p < 0.05) higher apparent SNR than pF 0.75 and no pF in the peripheral zone (PZ) of the prostate. In the multi-shot acquisitions in the PZ, pF 0.75 had a significantly higher apparent SNR than 0.65 pF and no pF. The single-shot pF 0.65 scan had a significantly lower ADC than single-shot no pF.

CONCLUSION

Deep breathing patterns can cause intravoxel dephasing in prostate DWI. For single-shot acquisitions at a b-value of 800 s/mm2, any potential risks of motion-related artefacts at low pF factors (pF 0.65) were outweighed by the increase in signal from a lower TE, as shown by the increase in apparent SNR. In multi-shot acquisitions however, the minimum pF factor should be larger, as shown by the lower apparent SNR at low pF factors.

Comparison of methods for intravoxel incoherent motion parameter estimation in the brain from flow-compensated and non-flow-compensated diffusion-encoded data

PURPOSE

Joint analysis of flow-compensated (FC) and non-flow-compensated (NC) diffusion MRI (dMRI) data has been suggested for increased robustness of intravoxel incoherent motion (IVIM) parameter estimation. For this purpose, a set of methods commonly used or previously found useful for IVIM analysis of dMRI data obtained with conventional diffusion encoding were evaluated in healthy human brain.

METHODS

Five methods for joint IVIM analysis of FC and NC dMRI data were compared: (1) direct non-linear least squares fitting, (2) a segmented fitting algorithm with estimation of the diffusion coefficient from higher b-values of NC data, (3) a Bayesian algorithm with uniform prior distributions, (4) a Bayesian algorithm with spatial prior distributions, and (5) a deep learning-based algorithm. Methods were evaluated on brain dMRI data from healthy subjects and simulated data at multiple noise levels. Bipolar diffusion encoding gradients were used with b-values 0-200 s/mm2 and corresponding flow weighting factors 0-2.35 s/mm for NC data and by design 0 for FC data. Data were acquired twice for repeatability analysis.

RESULTS

Measurement repeatability as well as estimation bias and variability were at similar levels or better with the Bayesian algorithm with spatial prior distributions and the deep learning-based algorithm for IVIM parametersD $$ D $$ andf $$ f $$ , and for the Bayesian algorithm only forv d $$ {v}_d $$ , relative to the other methods.

CONCLUSION

A Bayesian algorithm with spatial prior distributions is preferable for joint IVIM analysis of FC and NC dMRI data in the healthy human brain, but deep learning-based algorithms appear promising.

Feature-guided deep learning reduces signal loss and increases lesion CNR in diffusion-weighted imaging of the liver

PURPOSE

This research aims to develop a feature-guided deep learning approach and compare it with an optimized conventional post-processing algorithm in order to enhance the image quality of diffusion-weighted liver images and, in particular, to reduce the pulsation-induced signal loss occurring predominantly in the left liver lobe.

METHODS

Data from 40 patients with liver lesions were used. For the conventional approach, the best-suited out of five examined algorithms was chosen. For the deep learning approach, a U-Net was trained. Instead of learning "gold-standard" target images, the network was trained to optimize four image features (lesion CNR, vessel darkness, data consistency, and pulsation artifact reduction), which could be assessed quantitatively using manually drawn ROIs. A quality score was calculated from these four features. As an additional quality assessment, three radiologists rated different features of the resulting images.

RESULTS

The conventional approach could substantially increase the lesion CNR and reduce the pulsation-induced signal loss. However, the vessel darkness was reduced. The deep learning approach increased the lesion CNR and reduced the signal loss to a slightly lower extent, but it could additionally increase the vessel darkness. According to the image quality score, the quality of the deep-learning images was higher than that of the images obtained using the conventional approach. The radiologist ratings were mostly consistent with the quantitative scores, but the overall quality ratings differed among the readers.

CONCLUSION

Unlike the conventional algorithm, the deep-learning algorithm increased the vessel darkness. Therefore, it may be a viable alternative to conventional algorithms.

Self-calibrated subspace reconstruction for multidimensional MR fingerprinting for simultaneous relaxation and diffusion quantification

PURPOSE

To propose a new reconstruction method for multidimensional MR fingerprinting (mdMRF) to address shading artifacts caused by physiological motion-induced measurement errors without navigating or gating.

METHODS

The proposed method comprises two procedures: self-calibration and subspace reconstruction. The first procedure (self-calibration) applies temporally local matrix completion to reconstruct low-resolution images from a subset of under-sampled data extracted from the k-space center. The second procedure (subspace reconstruction) utilizes temporally global subspace reconstruction with pre-estimated temporal subspace from low-resolution images to reconstruct aliasing-free, high-resolution, and time-resolved images. After reconstruction, a customized outlier detection algorithm was employed to automatically detect and remove images corrupted by measurement errors. Feasibility, robustness, and scan efficiency were evaluated through in vivo human brain imaging experiments.

RESULTS

The proposed method successfully reconstructed aliasing-free, high-resolution, and time-resolved images, where the measurement errors were accurately represented. The corrupted images were automatically and robustly detected and removed. Artifact-free T1, T2, and ADC maps were generated simultaneously. The proposed reconstruction method demonstrated robustness across different scanners, parameter settings, and subjects. A high scan efficiency of less than 20 s per slice has been achieved.

CONCLUSION

The proposed reconstruction method can effectively alleviate shading artifacts caused by physiological motion-induced measurement errors. It enables simultaneous and artifact-free quantification of T1, T2, and ADC using mdMRF scans without prospective gating, with robustness and high scan efficiency.

Impact of velocity- and acceleration-compensated encodings on signal dropout and black-blood state in diffusion-weighted magnetic resonance liver imaging at clinical TEs

PURPOSE

The study aims to develop easy-to-implement concomitant field-compensated gradient waveforms with varying velocity-weighting (M1) and acceleration-weighting (M2) levels and to evaluate their efficacy in correcting signal dropouts and preserving the black-blood state in liver diffusion-weighted imaging. Additionally, we seek to determine an optimal degree of compensation that minimizes signal dropouts while maintaining blood signal suppression.

METHODS

Numerically optimized gradient waveforms were adapted using a novel method that allows for the simultaneous tuning of M1- and M2-weighting by changing only one timing variable. Seven healthy volunteers underwent diffusion-weighted magnetic resonance imaging (DWI) with five diffusion encoding schemes (monopolar, velocity-compensated (M1 = 0), acceleration-compensated (M1 = M2 = 0), 84%-M1-M2-compensated, 67%-M1-M2-compensated) at b-values of 50 and 800 s/mm2 at a constant echo time of 70 ms. Signal dropout correction and apparent diffusion coefficients (ADCs) were quantified using regions of interest in the left and right liver lobe. The blood appearance was evaluated using two five-point Likert scales.

RESULTS

Signal dropout was more pronounced in the left lobe (19%-42% less signal than in the right lobe with monopolar scheme) and best corrected by acceleration-compensation (8%-10% less signal than in the right lobe). The black-blood state was best with monopolar encodings and decreased significantly (p < 0.001) with velocity- and/or acceleration-compensation. The partially M1-M2-compensated encoding schemes could restore the black-blood state again. Strongest ADC bias occurred for monopolar encodings (difference between left/right lobe of 0.41 μm2/ms for monopolar vs. < 0.12 μm2/ms for the other encodings).

CONCLUSION

All of the diffusion encodings used in this study demonstrated suitability for routine DWI application. The results indicate that a perfect value for the level of M1-M2-compensation does not exist. However, among the examined encodings, the 84%-M1-M2-compensated encodings provided a suitable tradeoff.

Cardiac q-space trajectory imaging by motion-compensated tensor-valued diffusion encoding in human heart in vivo

PURPOSE

Tensor-valued diffusion encoding can probe more specific features of tissue microstructure than what is available by conventional diffusion weighting. In this work, we investigate the technical feasibility of tensor-valued diffusion encoding at high b-values with q-space trajectory imaging (QTI) analysis, in the human heart in vivo.

METHODS

Ten healthy volunteers were scanned on a 3T scanner. We designed time-optimal gradient waveforms for tensor-valued diffusion encoding (linear and planar) with second-order motion compensation. Data were analyzed with QTI. Normal values and repeatability were investigated for the mean diffusivity (MD), fractional anisotropy (FA), microscopic FA (μFA), isotropic, anisotropic and total mean kurtosis (MKi, MKa, and MKt), and orientation coherence (Cc ). A phantom, consisting of two fiber blocks at adjustable angles, was used to evaluate sensitivity of parameters to orientation dispersion and diffusion time.

RESULTS

QTI data in the left ventricular myocardium were MD = 1.62 ± 0.07 μm2 /ms, FA = 0.31 ± 0.03, μFA = 0.43 ± 0.07, MKa = 0.20 ± 0.07, MKi = 0.13 ± 0.03, MKt = 0.33 ± 0.09, and Cc  = 0.56 ± 0.22 (mean ± SD across subjects). Phantom experiments showed that FA depends on orientation dispersion, whereas μFA was insensitive to this effect.

CONCLUSION

We demonstrated the first tensor-valued diffusion encoding and QTI analysis in the heart in vivo, along with first measurements of myocardial μFA, MKi, MKa, and Cc . The methodology is technically feasible and provides promising novel biomarkers for myocardial tissue characterization.

MR Fingerprinting with b-Tensor Encoding for Simultaneous Quantification of Relaxation and Diffusion in a Single Scan

PURPOSE

Although both relaxation and diffusion imaging are sensitive to tissue microstructure, studies have reported limited sensitivity and robustness of using relaxation or conventional diffusion alone to characterize tissue microstructure. Recently, it has been shown that tensor-valued diffusion encoding and joint relaxation-diffusion quantification enable more reliable quantification of compartment-specific microstructural properties. However, scan times to acquire such data can be prohibitive. Here, we aim to simultaneously quantify relaxation and diffusion using MR fingerprinting (MRF) and b-tensor encoding in a clinically feasible time.

METHODS

We developed multidimensional MRF scans (mdMRF) with linear and spherical b-tensor encoding (LTE and STE) to simultaneously quantify T1, T2, and ADC maps from a single scan. The image quality, accuracy, and scan efficiency were compared between the mdMRF using LTE and STE. Moreover, we investigated the robustness of different sequence designs to signal errors and their impact on the maps.

RESULTS

T1 and T2 maps derived from the mdMRF scans have consistently high image quality, while ADC maps are sensitive to different sequence designs. Notably, the fast imaging steady state precession (FISP)-based mdMRF scan with peripheral pulse gating provides the best ADC maps that are free of image distortion and shading artifacts.

CONCLUSION

We demonstrated the feasibility of quantifying T1, T2, and ADC maps simultaneously from a single mdMRF scan in around 24 s/slice. The map quality and quantitative values are consistent with the reference scans.

Motion compensated renal diffusion weighted imaging

PURPOSE

To assess the effect of respiratory motion and cardiac driven pulsation in renal DWI and to examine asymmetrical velocity-compensated diffusion encoding waveforms for robust ADC mapping in the kidneys.

METHODS

The standard monopolar Stejskal-Tanner pulsed gradient spin echo (pgse) and the asymmetric bipolar velocity-compensated (asym-vc) diffusion encoding waveforms were used for coronal renal DWI at 3T. The robustness of the ADC quantification in the kidneys was tested with the aforementioned waveforms in respiratory-triggered and breath-held cardiac-triggered scans at different trigger delays in 10 healthy subjects.

RESULTS

The pgse waveform showed higher ADC values in the right kidney at short trigger delays in comparison to longer trigger delays in the respiratory triggered scans when the diffusion gradient was applied in the feet-head (FH) direction. The coefficient of variation over all respiratory trigger delays, averaged over all subjects was 0.15 for the pgse waveform in the right kidney when diffusion was measured in the FH direction; the corresponding coefficient of variation for the asym-vc waveform was 0.06. The effect of cardiac driven pulsation was found to be small in comparison to the effect of respiratory motion.

CONCLUSION

Short trigger delays in respiratory-triggered scans can cause higher ADC values in comparison to longer trigger delays in renal DWI, especially in the right kidney when diffusion is measured in the FH direction. The asym-vc waveform can reduce ADC variation due to respiratory motion in respiratory-triggered scans at the cost of reduced SNR compared to the pgse waveform.

Reduction of the cardiac pulsation artifact and improvement of lesion conspicuity in flow‐compensated diffusion images in the liver—A quantitative evaluation of postprocessing algorithms

PURPOSE

To enhance image quality of flow-compensated diffusion-weighted liver MRI data by increasing the lesion conspicuity and reducing the cardiac pulsation artifact using postprocessing algorithms.

METHODS

Diffusion-weighted image data of 40 patients with liver lesions had been acquired at 1.5 T. These data were postprocessed with 5 different algorithms (weighted averaging, p-mean, percentile, outlier exclusion, and exception set). Four image properties of the postprocessed data were evaluated for optimizing the algorithm parameters. These properties were the lesion to tissue contrast-to-noise ratio (CNR), the reduction of the cardiac pulsation artifact, the data consistency, and the vessel darkness. They were combined into a total quality score ( Q total , $$ {Q}_{\mathrm{total}}, $$ set to 1 for the trace-weighted reference image), which was used to rate the image quality objectively.

RESULTS

The weighted averaging algorithm performed best according to the total quality score ( Q total = 1.111 ± 0.067 $$ {Q}_{\mathrm{total}}=1.111\pm 0.067 $$ ). The further ranking was outlier exclusion algorithm ( Q total = 1.086 ± 0.061 $$ {Q}_{\mathrm{total}}=1.086\pm 0.061 $$ ), p-mean algorithm ( Q total = 1.045 ± 0.049 $$ {Q}_{\mathrm{total}}=1.045\pm 0.049 $$ ), percentile algorithm ( Q total = 1.012 ± 0.049 $$ {Q}_{\mathrm{total}}=1.012\pm 0.049 $$ ), and exception set algorithm ( Q total = 0.957 ± 0.027 $$ {Q}_{\mathrm{total}}=0.957\pm 0.027 $$ ). All optimized algorithms except for the exception set algorithm corrected the pulsation artifact and increased the lesion CNR. Changes in Q total $$ {Q}_{\mathrm{total}} $$ were significant for all optimized algorithms except for the percentile algorithm. Liver ADC was significantly reduced (except for the exception set algorithm), particularly in the left lobe.

CONCLUSION

Postprocessing algorithms should be used for flow-compensated liver DWI. The proposed weighted averaging algorithm seems to be suited best to increase the image quality of artifact-corrupted flow-compensated diffusion-weighted liver data.

A review on investigation of the basic contrast mechanism underlying multidimensional diffusion MRI in assessment of neurological disorders

INTRODUCTION

Multidimensional diffusion MRI (MDD MRI) is a novel diffusion technique that uses advanced gradient waveforms for microstructural tissue characterization to provide information about average rate, anisotropy and orientation of the diffusion and to disentangle the signal fraction from specific cell types i.e., elongated cells, isotropic cells and free water.

AIM

To review the diagnostic potential of MDD MRI in the clinical setting for microstructural tissue characterization in patients with neurological disorders to aid in patient care and treatment.

METHOD

A scoping review on the clinical applications of MDD MRI was conducted from original articles published in PubMed and Scopus from 2015 to 2021 using the keywords "Multidimensional diffusion MRI" OR "diffusion tensor distribution" OR "Tensor-Valued Diffusion" OR "b-tensor encoding" OR "microscopic diffusion anisotropy" OR "microscopic anisotropy" OR "microscopic fractional anisotropy" OR "double diffusion encoding" OR "triple diffusion encoding" OR "double pulsed field gradients" OR "double wave vector" OR "correlation tensor imaging" AND "brain" OR "axons".

RESULTS

Initially 145 articles were screened and after applying inclusion and exclusion criteria, nine articles were included in the final analysis. In most of these studies, microscopic diffusion anisotropy within the lesion showed deviation from the normal-appearing tissue.

CONCLUSION

Multidimensional diffusion MRI can provide better quantification and visualization of tissue microstructure than conventional diffusion MRI and can be used in the clinical setting for diagnosis of neurological disorders.

Flow-compensated diffusion encoding in MRI for improved liver metastasis detection

Magnetic resonance (MR) diffusion-weighted imaging (DWI) is often used to detect focal liver lesions (FLLs), though DWI image quality can be limited in the left liver lobe owing to the pulsatile motion of the nearby heart. Flow-compensated (FloCo) diffusion encoding has been shown to reduce this pulsation artifact. The purpose of this prospective study was to intra-individually compare DWI of the liver acquired with conventional monopolar and FloCo diffusion encoding for assessing metastatic FLLs in non-cirrhotic patients. Forty patients with known or suspected multiple metastatic FLLs were included and measured at 1.5 T field strength with a conventional (monopolar) and a FloCo diffusion encoding EPI sequence (single refocused; b-values, 50 and 800 s/mm²). Two board-certified radiologists analyzed the DWI images independently. They issued Likert-scale ratings (1 = worst, 5 = best) for pulsation artifact severity and counted the difference of lesions visible at b = 800 s/mm² separately for small and large FLLs (i.e., < 1 cm or > 1 cm) and separately for left and right liver lobe. Differences between the two diffusion encodings were assessed with the Wilcoxon signed-rank test. Both readers found a reduction in pulsation artifact in the liver with FloCo encoding (p < 0.001 for both liver lobes). More small lesions were detected with FloCo diffusion encoding in both liver lobes (left lobe: six and seven additional lesions by readers 1 and 2, respectively; right lobe: five and seven additional lesions for readers 1 and 2, respectively). Both readers found one additional large lesion in the left liver lobe. Thus, flow-compensated diffusion encoding appears more effective than monopolar diffusion encoding for the detection of liver metastases.

Separating Glioma Hyperintensities From White Matter by Diffusion-Weighted Imaging With Spherical Tensor Encoding

Background

Tumor-related hyperintensities in high b-value diffusion-weighted imaging (DWI) are radiologically important in the workup of gliomas. However, the white matter may also appear as hyperintense, which may conflate interpretation.

Purpose

To investigate whether DWI with spherical b-tensor encoding (STE) can be used to suppress white matter and enhance the conspicuity of glioma hyperintensities unrelated to white matter.

Materials and Methods

Twenty-five patients with a glioma tumor and at least one pathology-related hyperintensity on DWI underwent conventional MRI at 3 T. The DWI was performed both with linear and spherical tensor encoding (LTE-DWI and STE-DWI). The LTE-DWI here refers to the DWI obtained with conventional diffusion encoding and averaged across diffusion-encoding directions. Retrospectively, the differences in contrast between LTE-DWI and STE-DWI, obtained at a b-value of 2,000 s/mm², were evaluated by comparing hyperintensities and contralateral normal-appearing white matter (NAWM) both visually and quantitatively in terms of the signal intensity ratio (SIR) and contrast-to-noise ratio efficiency (CNR_eff).

Results

The spherical tensor encoding DWI was more effective than LTE-DWI at suppressing signals from white matter and improved conspicuity of pathology-related hyperintensities. The median SIR improved in all cases and on average by 28%. The median (interquartile range) SIR was 1.9 (1.6 – 2.1) for STE and 1.4 (1.3 – 1.7) for LTE, with a significant difference of 0.4 (0.3 –0.5) (p < 10^–4, paired U-test). In 40% of the patients, the SIR was above 2 for STE-DWI, but with LTE-DWI, the SIR was below 2 for all patients. The CNR_eff of STE-DWI was significantly higher than of LTE-DWI: 2.5 (2 – 3.5) vs. 2.3 (1.7 – 3.1), with a significant difference of 0.4 (−0.1 –0.6) (p < 10^–3, paired U-test). The STE improved CNR_eff in 70% of the cases. We illustrate the benefits of STE-DWI in three patients, where STE-DWI may facilitate an improved radiological description of tumor-related hyperintensity, including one case that could have been missed out if only LTE-DWI was inspected.

Conclusion

The contrast mechanism of high b-value STE-DWI results in a stronger suppression of white matter than conventional LTE-DWI, and may, therefore, be more sensitive and specific for assessment of glioma tumors and DWI-hyperintensities.

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.

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.

Multi-tissue spherical deconvolution of tensor-valued diffusion MRI

Multi-tissue constrained spherical deconvolution (MT-CSD) leverages the characteristic b-value dependency of each tissue type to estimate both the apparent tissue densities and the white matter fiber orientation distribution function from diffusion MRI data. In this work, we generalize MT-CSD to tensor-valued diffusion encoding with arbitrary b-tensor shapes. This enables the use of data encoded with mixed b-tensors, rather than being limited to the subset of linear (conventional) b-tensors. Using the complete set of data, including all b-tensor shapes, provides a categorical improvement in the estimation of apparent tissue densities, fiber ODF, and resulting tractography. Furthermore, we demonstrate that including multiple b-tensor shapes in the analysis provides improved contrast between tissue types, in particular between gray matter and white matter. We also show that our approach provides high-quality apparent tissue density maps and high-quality fiber tracking from data, even with sparse sampling across b-tensors that yield whole-brain coverage at 2 mm isotropic resolution in approximately 5:15 min.

Cross-term-compensated gradient waveform design for tensor-valued diffusion MRI

Diffusion MRI uses magnetic field gradients to sensitize the signal to the random motion of spins. In addition to the prescribed gradient waveforms, background field gradients contribute to the diffusion weighting and thereby cause an error in the measured signal and consequent parameterization. The most prominent contribution to the error comes from so-called 'cross-terms.' In this work we present a novel gradient waveform design that enables diffusion encoding that cancels such cross-terms and yields a more accurate measurement. This is achieved by numerical optimization that maximizes encoding efficiency with a simultaneous constraint on the 'cross-term sensitivity' (c = 0). We found that the optimized cross-term-compensated waveforms were superior to previous cross-term-compensated designs for a wide range of waveform types that yield linear, planar, and spherical b-tensor encoding. The efficacy of the proposed design was also demonstrated in practical experiments using a clinical MRI system. The sensitivity to cross-terms was evaluated in a water phantom with a folded surface which provoked strong internal field gradients. In every comparison, the cross-term-compensated waveforms were robust to the effects of background gradients, whereas conventional designs were not. We also propose a method to measure background gradients from diffusion-weighted data, and show that cross-term-compensated waveforms produce parameters that are markedly less dependent on the background compared to non-compensated designs. Finally, we also used simulations to show that the proposed cross-term compensation was robust to background gradients in the interval 0 to 3 mT/m, whereas non-compensated designs were impacted in terms of a severe signal and parameter bias. In conclusion, we have proposed and demonstrated a waveform design that yields efficient cross-term compensation and facilitates accurate diffusion MRI in the presence of static background gradients regardless of their amplitude and direction. The optimization framework is compatible with arbitrary spin-echo sequence timing and RF events, b-tensor shapes, suppression of concomitant gradient effects and motion encoding, and is shared in open source.

Mapping prostatic microscopic anisotropy using linear and spherical b-tensor encoding: A preliminary study

Purpose:

Tensor-valued diffusion encoding provides more specific information than conventional diffusion-weighted imaging (DWI), but has mainly been applied in neuroimaging studies. This study aimed to assess its potential for the imaging of prostate cancer (PCa).

Methods:

Seventeen patients with histologically proven PCa were enrolled. DWI of the prostate was performed with linear and spherical tensor encoding using a maximal b-value of 1.5 ms/μm² and a voxel size of 3 × 3 × 4 mm³. The gamma-distribution model was used to estimate the mean diffusivity (MD), the isotropic kurtosis (MK_I), and the anisotropic kurtosis (MK_A). Regions of interest were placed in MR-defined cancerous tissues, as well as in apparently healthy tissues in the peripheral and transitional zones (PZs and TZs).

Results:

DWI with linear and spherical encoding yielded different image contrasts at high b-values, which enabled the estimation of MK_A and MK_I. Compared with healthy tissue (PZs and TZs combined) the cancers displayed a significantly lower MD (P < .05), higher MK_I (P < 10⁻⁵), and lower MK_A (P < .05). Compared with the TZ, tissue in the PZ showed lower MD (P < 10⁻³) and higher MK_A (P < 10⁻³). No significant differences were found between cancers of different Gleason scores, possibly because of the limited sample size.

Conclusion:

Tensor-valued diffusion encoding enabled mapping of MK_A and MK_I in the prostate. The elevated MK_I in PCa compared with normal tissues suggests an elevated heterogeneity in the cancers. Increased in-plane resolution could improve tumor delineation in future studies.

Motion-compensated gradient waveforms for tensor-valued diffusion encoding by constrained numerical optimization

PURPOSE

Diffusion-weighted MRI is sensitive to incoherent tissue motion, which may confound the measured signal and subsequent analysis. We propose a "motion-compensated" gradient waveform design for tensor-valued diffusion encoding that negates the effects bulk motion and incoherent motion in the ballistic regime.

METHODS

Motion compensation was achieved by constraining the magnitude of gradient waveform moment vectors. The constraint was incorporated into a numerical optimization framework, along with existing constraints that account for b-tensor shape, hardware restrictions, and concomitant field gradients. We evaluated the efficacy of encoding and motion compensation in simulations, and we demonstrated the approach by linear and planar b-tensor encoding in a healthy heart in vivo.

RESULTS

The optimization framework produced asymmetric motion-compensated waveforms that yielded b-tensors of arbitrary shape with improved efficiency compared with previous designs for tensor-valued encoding, and equivalent efficiency to previous designs for linear (conventional) encoding. Technical feasibility was demonstrated in the heart in vivo, showing vastly improved data quality when using motion compensation. The optimization framework is available online in open source.

CONCLUSION

Our gradient waveform design is both more flexible and efficient than previous methods, facilitating tensor-valued diffusion encoding in tissues in which motion would otherwise confound the signal. The proposed design exploits asymmetric encoding times, a single refocusing pulse or multiple refocusing pulses, and integrates compensation for concomitant gradient effects throughout the imaging volume.