Accelerated MR diffusion tensor imaging using distributed compressed sensing

Continuous diffusion spectrum computation for diffusion-weighted magnetic resonance imaging of the kidney tubule system

Background

The use of rigid multi-exponential models (with a priori predefined numbers of components) is common practice for diffusion-weighted MRI (DWI) analysis of the kidney. This approach may not accurately reflect renal microstructure, as the data are forced to conform to the a priori assumptions of simplified models. This work examines the feasibility of less constrained, data-driven non-negative least squares (NNLS) continuum modelling for DWI of the kidney tubule system in simulations that include emulations of pathophysiological conditions.

Methods

Non-linear least squares (LS) fitting was used as reference for the simulations. For performance assessment, a threshold of 5% or 10% for the mean absolute percentage error (MAPE) of NNLS and LS results was used. As ground truth, a tri-exponential model using defined volume fractions and diffusion coefficients for each renal compartment (tubule system: D_tubules , f_tubules ; renal tissue: D_tissue , f_tissue ; renal blood: D_blood , f_blood ;) was applied. The impact of: (I) signal-to-noise ratio (SNR) =40-1,000, (II) number of b-values (n=10-50), (III) diffusion weighting (b-range_small =0-800 up to b-range_large =0-2,180 s/mm²), and (IV) fixation of the diffusion coefficients D_tissue and D_blood was examined. NNLS was evaluated for baseline and pathophysiological conditions, namely increased tubular volume fraction (ITV) and renal fibrosis (10%: grade I, mild) and 30% (grade II, moderate).

Results

NNLS showed the same high degree of reliability as the non-linear LS. MAPE of the tubular volume fraction (f_tubules ) decreased with increasing SNR. Increasing the number of b-values was beneficial for f_tubules precision. Using the b-range_large led to a decrease in MAPE _ftubules compared to b-range_small. The use of a medium b-value range of b=0-1,380 s/mm² improved f_tubules precision, and further b_max increases beyond this range yielded diminishing improvements. Fixing D_blood and D_tissue significantly reduced MAPE _ftubules and provided near perfect distinction between baseline and ITV conditions. Without constraining the number of renal compartments in advance, NNLS was able to detect the (fourth) fibrotic compartment, to differentiate it from the other three diffusion components, and to distinguish between 10% vs. 30% fibrosis.

Conclusions

This work demonstrates the feasibility of NNLS modelling for DWI of the kidney tubule system and shows its potential for examining diffusion compartments associated with renal pathophysiology including ITV fraction and different degrees of fibrosis.

Optimal acceleration factor for image acquisition in turbo spin echo: diffusion-weighted imaging with compressed sensing

In this study, the change in the image quality and apparent diffusion coefficient (ADC) with increase in the acceleration factor (AF) was analyzed and the most optimal AF was determined to reduce the scan time while preserving the image quality. The AF was changed from 2 to 20 in the MR acquisitions. The similarities between the accelerated and reference images were determined based on the structural similarity (SSIM) index for DWI image and coefficient of variation (%CV) for ADC. The SSIM index decreased significantly when the AF ≥ 8 compared with when the AF = 2 (p < 0.05). In the reference image, the %CV of the ADC increased significantly when the AF ≥ 10 (p < 0.01). In conclusion, a remarkable decrease in the image quality and ADC was observed when the AF was > 8. Thus, an AF < 8 would be optimal for reducing the scan time while preserving the image quality.

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.

Acceleration of three-dimensional diffusion magnetic resonance imaging using a kernel low-rank compressed sensing method

Diffusion Magnetic Resonance Imaging (dMRI) has shown great potential in probing tissue microstructure and structural connectivity in the brain but is often limited by the lengthy scan time needed to sample the diffusion profile by acquiring multiple diffusion weighted images (DWIs). Although parallel imaging technique has improved the speed of dMRI acquisition, attaining high resolution three dimensional (3D) dMRI on preclinical MRI systems remained still time consuming. In this paper, kernel principal component analysis, a machine learning approach, was employed to estimate the correlation among DWIs. We demonstrated the feasibility of such correlation estimation from low-resolution training DWIs and used the correlation as a constraint to reconstruct high-resolution DWIs from highly under-sampled k-space data, which significantly reduced the scan time. Using full k-space 3D dMRI data of post-mortem mouse brains, we retrospectively compared the performance of the so-called kernel low rank (KLR) method with a conventional compressed sensing (CS) method in terms of image quality and ability to resolve complex fiber orientations and connectivity. The results demonstrated that the KLR-CS method outperformed the conventional CS method for acceleration factors up to 8 and was likely to enhance our ability to investigate brain microstructure and connectivity using high-resolution 3D dMRI.

Molecular Subtype Classification in Lower-Grade Glioma with Accelerated DTI

BACKGROUND AND PURPOSE

Image-based classification of lower-grade glioma molecular subtypes has substantial prognostic value. Diffusion tensor imaging has shown promise in lower-grade glioma subtyping but currently requires lengthy, nonstandard acquisitions. Our goal was to investigate lower-grade glioma classification using a machine learning technique that estimates fractional anisotropy from accelerated diffusion MR imaging scans containing only 3 diffusion-encoding directions.

MATERIALS AND METHODS

Patients with lower-grade gliomas (n = 41) (World Health Organization grades II and III) with known isocitrate dehydrogenase (IDH) mutation and 1p/19q codeletion status were imaged preoperatively with DTI. Whole-tumor volumes were autodelineated using conventional anatomic MR imaging sequences. In addition to conventional ADC and fractional anisotropy reconstructions, fractional anisotropy estimates were computed from 3-direction DTI subsets using DiffNet, a neural network that directly computes fractional anisotropy from raw DTI data. Differences in whole-tumor ADC, fractional anisotropy, and estimated fractional anisotropy were assessed between IDH-wild-type and IDH-mutant lower-grade gliomas with and without 1p/19q codeletion. Multivariate classification models were developed using whole-tumor histogram and texture features from ADC, ADC + fractional anisotropy, and ADC + estimated fractional anisotropy to identify the added value provided by fractional anisotropy and estimated fractional anisotropy.

RESULTS

ADC (P = .008), fractional anisotropy (P < .001), and estimated fractional anisotropy (P < .001) significantly differed between IDH-wild-type and IDH-mutant lower-grade gliomas. ADC (P < .001) significantly differed between IDH-mutant gliomas with and without codeletion. ADC-only multivariate classification predicted IDH mutation status with an area under the curve of 0.81 and codeletion status with an area under the curve of 0.83. Performance improved to area under the curve = 0.90/0.94 for the ADC + fractional anisotropy classification and to area under the curve = 0.89/0.89 for the ADC + estimated fractional anisotropy classification.

CONCLUSIONS

Fractional anisotropy estimates made from accelerated 3-direction DTI scans add value in classifying lower-grade glioma molecular status.

Common Information Enhanced Reconstruction for Accelerated High-resolution Multi-shot Diffusion Imaging

PURPOSE

Multi-shot technique can effectively achieve high-resolution diffusion weighted images, but the acquisition time of multi-shot technique is prolonged, especially for multiple direction diffusion encoding. Thus, increasing acquisition efficiency is highly desirable for high-resolution diffusion tensor imaging (DTI). In this study, based on the assumption that different diffusion directions share the common information, image ratio constrained reconstruction (IRCR) combined with iterative self-consistent parallel imaging reconstruction (SPIRiT) is proposed to improve data sampling efficiency and image reconstruction fidelity for high-resolution DTI.

THEORY AND METHODS

The proposed reconstruction framework is named Common Information Enhanced Reconstruction (CIER). Inter-image correlation among different direction diffusion-weighted images is used through common information, which is an isotropic component and structure, for improving the performance of reconstruction. The framework consists of three steps. (i) Pre-processing: three intermediate multi-shot images, low-resolution composite image, high-resolution composite image and low-resolution diffusion weighted image, are generated based on the SPIRiT method. (ii) IRCR: the initial high-resolution diffusion weighted image is calculated from the images in step (i) based on that the ratio map between high-resolution images is approximated by the ratio map between the corresponding low-resolution images. (iii) Final SPIRiT reconstruction: the final image is generated with the image from IRCR as initialization by considering data consistency only in the SPIRiT calculation. A specific implementation based on multishot variable density spiral (VDS) DTI is used to demonstrate the method.

RESULTS

The proposed CIER method was compared with the traditional reconstruction methods, conjugate gradient SENSE (CG-SENSE), L1-regularized SPIRiT (L1-SPIRiT), and anisotropic-sparsity SPIRiT (AS-SPIRiT) in brain DTI at acceleration factors of 3 to 7. CIER provided better diffusion image quality than other methods shown by both qualitative and quantitative results, especially at higher undersampling acceleration factors.

CONCLUSION

CIER offers better diffusion image quality at higher undersampling acceleration factors for high-resolution DTI. Both qualitative and quantitative results prove that common information can be used to improve sampling efficiency and maintain the image quality of diffusion-weighted images.

Highly accelerated, model-free diffusion tensor MRI reconstruction using neural networks

PURPOSE

The purpose of this study was to develop a neural network that accurately performs diffusion tensor imaging (DTI) reconstruction from highly accelerated scans.

MATERIALS AND METHODS

This retrospective study was conducted using data acquired between 2013 and 2018 and was approved by the local institutional review board. DTI acquired in healthy volunteers (N = 10) was used to train a neural network, DiffNet, to reconstruct fractional anisotropy (FA) and mean diffusivity (MD) maps from small subsets of acquired DTI data with between 3 and 20 diffusion-encoding directions. FA and MD maps were then reconstructed in volunteers and in patients with glioblastoma multiforme (GBM, N = 12) using both DiffNet and conventional reconstructions. Accuracy and precision were quantified in volunteer scans and compared between reconstructions. The accuracy of tumor delineation was compared between reconstructed patient data by evaluating agreement between DTI-derived tumor volumes and volumes defined by contrast-enhanced T1-weighted MRI. Comparisons were performed using areas under the receiver operating characteristic curves (AUC).

RESULTS

DiffNet FA reconstructions were more accurate and precise compared with conventional reconstructions for all acceleration factors. DiffNet permitted reconstruction with only three diffusion-encoding directions with significantly lower bias than the conventional method using six directions (0.01 ± 0.01 vs 0.06 ± 0.01, P < 0.001). While MD-based tumor delineation was not substantially different with DiffNet (AUC range: 0.888-0.902), DiffNet FA had higher AUC than conventional reconstructions for fixed scan time and achieved similar performance with shorter scans (conventional, six directions: AUC = 0.926, DiffNet, three directions: AUC = 0.920).

CONCLUSION

DiffNet improved DTI reconstruction accuracy, precision, and tumor delineation performance in GBM while permitting reconstruction from only three diffusion-encoding directions.&!#6.

Cardiac diffusion tensor imaging based on compressed sensing using joint sparsity and low-rank approximation

Diffusion tensor magnetic resonance (DTMR) imaging and diffusion tensor imaging (DTI) have been widely used to probe noninvasively biological tissue structures. However, DTI suffers from long acquisition times, which limit its practical and clinical applications. This paper proposes a new Compressed Sensing (CS) reconstruction method that employs joint sparsity and rank deficiency to reconstruct cardiac DTMR images from undersampled k-space data. Diffusion-weighted images acquired in different diffusion directions were firstly stacked as columns to form the matrix. The matrix was row sparse in the transform domain and had a low rank. These two properties were then incorporated into the CS reconstruction framework. The underlying constrained optimization problem was finally solved by the first-order fast method. Experiments were carried out on both simulation and real human cardiac DTMR images. The results demonstrated that the proposed approach had lower reconstruction errors for DTI indices, including fractional anisotropy (FA) and mean diffusivities (MD), compared to the existing CS-DTMR image reconstruction techniques.

Prospective acceleration of diffusion tensor imaging with compressed sensing using adaptive dictionaries

PURPOSE

Diffusion MRI requires acquisition of multiple diffusion-weighted images, resulting in long scan times. Here, we investigate combining compressed sensing and a fast imaging sequence to dramatically reduce acquisition times in cardiac diffusion MRI.

METHODS

Fully sampled and prospectively undersampled diffusion tensor imaging data were acquired in five rat hearts at acceleration factors of between two and six using a fast spin echo (FSE) sequence. Images were reconstructed using a compressed sensing framework, enforcing sparsity by means of decomposition by adaptive dictionaries. A tensor was fit to the reconstructed images and fiber tractography was performed.

RESULTS

Acceleration factors of up to six were achieved, with a modest increase in root mean square error of mean apparent diffusion coefficient (ADC), fractional anisotropy (FA), and helix angle. At an acceleration factor of six, mean values of ADC and FA were within 2.5% and 5% of the ground truth, respectively. Marginal differences were observed in the fiber tracts.

CONCLUSION

We developed a new k-space sampling strategy for acquiring prospectively undersampled diffusion-weighted data, and validated a novel compressed sensing reconstruction algorithm based on adaptive dictionaries. The k-space undersampling and FSE acquisition each reduced acquisition times by up to 6× and 8×, respectively, as compared to fully sampled spin echo imaging. Magn Reson Med 76:248-258, 2016. © 2015 Wiley Periodicals, Inc.

Fast iterative algorithm for the reconstruction of multishot non-cartesian diffusion data

PURPOSE

To accelerate the motion-compensated iterative reconstruction of multishot non-Cartesian diffusion data.

METHOD

The motion-compensated recovery of multishot non-Cartesian diffusion data is often performed using a modified iterative sensitivity-encoded algorithm. Specifically, the encoding matrix is replaced with a combination of nonuniform Fourier transforms and composite sensitivity functions, which account for the motion-induced phase errors. The main challenge with this scheme is the significantly increased computational complexity, which is directly related to the total number of composite sensitivity functions (number of shots × number of coils). The dimensionality of the composite sensitivity functions and hence the number of Fourier transforms within each iteration is reduced using a principal component analysis-based scheme. Using a Toeplitz-based conjugate gradient approach in combination with an augmented Lagrangian optimization scheme, a fast algorithm is proposed for the sparse recovery of diffusion data.

RESULTS

The proposed simplifications considerably reduce the computation time, especially in the recovery of diffusion data from under-sampled reconstructions using sparse optimization. By choosing appropriate number of basis functions to approximate the composite sensitivities, faster reconstruction (close to 9 times) with effective motion compensation is achieved.

CONCLUSION

The proposed enhancements can offer fast motion-compensated reconstruction of multishot diffusion data for arbitrary k-space trajectories.

Parallel imaging and compressed sensing combined framework for accelerating high-resolution diffusion tensor imaging using inter-image correlation

PURPOSE

Increasing acquisition efficiency is always a challenge in high-resolution diffusion tensor imaging (DTI), which has low signal-to-noise ratio and is sensitive to reconstruction artifacts. In this study, a parallel imaging (PI) and compressed sensing (CS) combined framework is proposed, which features motion error correction, PI calibration, and sparsity model using inter-image correlation tailored for high-resolution DTI.

THEORY AND METHODS

The proposed method, named anisotropic sparsity SPIRiT, consists of three steps: (i) motion-induced phase error estimation, (ii) initial CS reconstruction and PI kernel calibration, and (iii) final reconstruction combining PI and CS. Inter-image correlation of diffusion-weighted images are used through anisotropic signals for improved sparsity. A specific implementation based on multishot variable density spiral DTI is used to demonstrate the method.

RESULTS

The proposed reconstruction method was compared with CG-SENSE, CS-based joint reconstruction, and PI and CS combined methods with L1 and joint sparsity regularization, in brain DTI experiments at acceleration factors of 3 to 5. Both qualitative and quantitative results demonstrated that the proposed method resulted in better preserved image quality and more accurate DTI parameters than other methods.

CONCLUSION

The proposed method can accelerate high-resolution DTI acquisition effectively by using the sharable information among different diffusion encoding directions.