Accelerated Cardiac Diffusion Tensor Imaging Using Joint Low-Rank and Sparsity Constraints

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.

Second-Order Motion-Compensated Echo-Planar Cardiac Diffusion-Weighted MRI: Usefulness of Compressed Sensitivity Encoding

BACKGROUND

Cardiac diffusion-weighted imaging (DWI) using second-order motion-compensated spin echo (M2C) can provide noninvasive in-vivo microstructural assessment, but limited by relatively low signal-to-noise ratio (SNR). Echo-planar imaging (EPI) with compressed sensitivity encoding (EPICS) could address these issues.

PURPOSE

To combine M2C DWI and EPCIS (M2C EPICS DWI), and compare image quality for M2C DWI.

STUDY TYPE

Prospective.

POPULATION

Ten ex-vivo hearts, 10 healthy volunteers (females, 5 [50%]; mean ± SD of age, 25 ± 4 years), and 12 patients with diseased hearts (female, 1 [8.3%]; mean ± SD of age, 44 ± 16 years; including coronary artery heart disease, congenital heart disease, dilated cardiomyopathy, amyloidosis, and myocarditis).

FIELD STRENGTH/SEQUENCE

3-T, M2C EPICS DWI, and M2C DWI.

ASSESSMENT

The apparent SNR (aSNR) and the rating scores were used to evaluate and compared image quality of all three groups. The aSNR was calculated usingaSNR = Mean intensity myocardium / Standard deviation myocardium $$ \mathrm{aSNR}={\mathrm{Mean}\ \mathrm{intensity}}_{\mathrm{myocardium}}/{\mathrm{Standard}\ \mathrm{deviation}}_{\mathrm{myocardium}} $$ , and the myocardium was segmented manually. Three observers independently rated subjective image quality using a 5-point Likert scale.

STATISTICAL TESTS

Bland-Altman analysis and paired t-tests. The threshold for statistical significance was set at P < 0.05.

RESULTS

In healthy volunteers, the aSNR with a b-value of 450 s/mm2 acquired by M2C EPICS DWI was significantly higher than M2C DWI at in-plane resolutions of 3.0 × 3.0, 2.5 × 2.5, and 2.0 × 2.0 mm2. In patients with diseased hearts, the aSNR ofM2C EPICS DWI was also significantly higher than that for M2C DWI (bias of M2C EPICS-M2C = 1.999, 95% limits of agreement, 0.362 to 3.636; mean ± SD, 7.80 ± 1.37 vs. 5.80 ± 0.81). The ADC values of M2C EPICS was significantly higher than M2C DWI in in-vivo hearts. Over 80% of the images with rating scores for M2C EPICS DWI were higher than M2C DWI in in-vivo hearts.

DATA CONCLUSION

Cardiac imaging by M2C EPICS DWI may demonstrate better overall image quality and higher aSNR than M2C DWI.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Deep learning-based diffusion tensor cardiac magnetic resonance reconstruction: a comparison study

In vivo cardiac diffusion tensor imaging (cDTI) is a promising Magnetic Resonance Imaging (MRI) technique for evaluating the microstructure of myocardial tissue in living hearts, providing insights into cardiac function and enabling the development of innovative therapeutic strategies. However, the integration of cDTI into routine clinical practice poses challenging due to the technical obstacles involved in the acquisition, such as low signal-to-noise ratio and prolonged scanning times. In this study, we investigated and implemented three different types of deep learning-based MRI reconstruction models for cDTI reconstruction. We evaluated the performance of these models based on the reconstruction quality assessment, the diffusion tensor parameter assessment as well as the computational cost assessment. Our results indicate that the models discussed in this study can be applied for clinical use at an acceleration factor (AF) of × 2 and × 4 , with the D5C5 model showing superior fidelity for reconstruction and the SwinMR model providing higher perceptual scores. There is no statistical difference from the reference for all diffusion tensor parameters at AF × 2 or most DT parameters at AF × 4 , and the quality of most diffusion tensor parameter maps is visually acceptable. SwinMR is recommended as the optimal approach for reconstruction at AF × 2 and AF × 4 . However, we believe that the models discussed in this study are not yet ready for clinical use at a higher AF. At AF × 8 , the performance of all models discussed remains limited, with only half of the diffusion tensor parameters being recovered to a level with no statistical difference from the reference. Some diffusion tensor parameter maps even provide wrong and misleading information.

Accelerated cardiac diffusion tensor imaging using deep neural network

Cardiac diffusion tensor imaging (DTI) is a noninvasive method for measuring the microstructure of the myocardium. However, its long scan time significantly hinders its wide application. In this study, we developed a deep learning framework to obtain high-quality DTI parameter maps from six diffusion-weighted images (DWIs) by combining deep-learning-based image generation and tensor fitting, and named the new framework FG-Net. In contrast to frameworks explored in previous deep-learning-based fast DTI studies, FG-Net generates inter-directional DWIs from six input DWIs to supplement the loss information and improve estimation accuracy for DTI parameters. FG-Net was evaluated using two datasets ofex vivohuman hearts. The results showed that FG-Net can generate fractional anisotropy, mean diffusivity maps, and helix angle maps from only six raw DWIs, with a quantification error of less than 5%. FG-Net outperformed conventional tensor fitting and black-box network fitting in both qualitative and quantitative metrics. We also demonstrated that the proposed FG-Net can achieve highly accurate fractional anisotropy and helix angle maps in DWIs with differentb-values.

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.

WA, Boyle J, Coll‐Font J, Dall'Armellina E, Ennis DB, Ferreira PF, Kalra P, Kolipaka A, Kozerke S, Lohr D, Mongeon F, Moulin K, Nguyen C, Nielles‐Vallespin S, Raterman B, Schreiber LM, Scott AD, Sosnovik DE, Stoeck CT, Tous C, Tunnicliffe EM, Weng AM, Croisille P, Viallon M, Schneider JE. Validation of cardiac diffusion tensor imaging sequences: A multicentre test-retest phantom study

Cardiac diffusion tensor imaging (DTI) is an emerging technique for the in vivo characterisation of myocardial microstructure, and there is a growing need for its validation and standardisation. We sought to establish the accuracy, precision, repeatability and reproducibility of state-of-the-art pulse sequences for cardiac DTI among 10 centres internationally. Phantoms comprising 0%-20% polyvinylpyrrolidone (PVP) were scanned with DTI using a product pulsed gradient spin echo (PGSE; N = 10 sites) sequence, and a custom motion-compensated spin echo (SE; N = 5) or stimulated echo acquisition mode (STEAM; N = 5) sequence suitable for cardiac DTI in vivo. A second identical scan was performed 1-9 days later, and the data were analysed centrally. The average mean diffusivities (MDs) in 0% PVP were (1.124, 1.130, 1.113) x 10^-3  mm² /s for PGSE, SE and STEAM, respectively, and accurate to within 1.5% of reference data from the literature. The coefficients of variation in MDs across sites were 2.6%, 3.1% and 2.1% for PGSE, SE and STEAM, respectively, and were similar to previous studies using only PGSE. Reproducibility in MD was excellent, with mean differences in PGSE, SE and STEAM of (0.3 ± 2.3, 0.24 ± 0.95, 0.52 ± 0.58) x 10^-5  mm² /s (mean ± 1.96 SD). We show that custom sequences for cardiac DTI provide accurate, precise, repeatable and reproducible measurements. Further work in anisotropic and/or deforming phantoms is warranted.

Accelerated in Vivo Cardiac Diffusion-Tensor MRI Using Residual Deep Learning-based Denoising in Participants with Obesity

Purpose

To develop and assess a residual deep learning algorithm to accelerate in vivo cardiac diffusion-tensor MRI (DT-MRI) by reducing the number of averages while preserving image quality and DT-MRI parameters.

Materials and Methods

In this prospective study, a denoising convolutional neural network (DnCNN) for DT-MRI was developed; a total of 26 participants, including 20 without obesity (body mass index [BMI] < 30 kg/m²; mean age, 28 years ± 3 [standard deviation]; 11 women) and six with obesity (BMI ≥ 30 kg/m²; mean age, 48 years ± 11; five women), were recruited from June 19, 2019, to July 29, 2020. DT-MRI data were constructed at four averages (4Av), two averages (2Av), and one average (1Av) without and with the application of the DnCNN (4Av_DnCNN, 2Av_DnCNN, 1Av_DnCNN). All data were compared against the reference DT-MRI data constructed at eight averages (8Av). Image quality, characterized by using the signal-to-noise ratio (SNR) and structural similarity index (SSIM), and the DT-MRI parameters of mean diffusivity (MD), fractional anisotropy (FA), and helix angle transmurality (HAT) were quantified.

Results

No differences were found in image quality or DT-MRI parameters between the accelerated 4Av_DnCNN DT-MRI and the reference 8Av DT-MRI data for the SNR (29.1 ± 2.7 vs 30.5 ± 2.9), SSIM (0.97 ± 0.01), MD (1.3 µm²/msec ± 0.1 vs 1.31 µm²/msec ± 0.11), FA (0.32 ± 0.05 vs 0.30 ± 0.04), or HAT (1.10°/% ± 0.13 vs 1.11°/% ± 0.09). The relationship of a higher MD and lower FA and HAT in individuals with obesity compared with individuals without obesity in reference 8Av DT-MRI measurements was retained in 4Av_DnCNN and 2Av_DnCNN DT-MRI measurements but was not retained in 4Av or 2Av DT-MRI measurements.

Conclusion

Cardiac DT-MRI can be performed at an at least twofold-accelerated rate by using DnCNN to preserve image quality and DT-MRI parameter quantification.Keywords: Adults, Cardiac, Obesity, Technology Assessment, MR-Diffusion Tensor Imaging, Heart, Tissue CharacterizationSupplemental material is available for this article.© RSNA, 2021.

Fast submillimeter diffusion MRI using gSlider-SMS and SNR-enhancing joint reconstruction

PURPOSE

We evaluate a new approach for achieving diffusion MRI data with high spatial resolution, large volume coverage, and fast acquisition speed.

THEORY AND METHODS

A recent method called gSlider-SMS enables whole-brain submillimeter diffusion MRI with high signal-to-noise ratio (SNR) efficiency. However, despite the efficient acquisition, the resulting images can still suffer from low SNR due to the small size of the imaging voxels. This work proposes to mitigate the SNR problem by combining gSlider-SMS with a regularized SNR-enhancing reconstruction approach.

RESULTS

Illustrative results show that, from gSlider-SMS data acquired over a span of only 25 minutes on a 3T scanner, the proposed method is able to produce 71 MRI images (64 diffusion encoding orientations with b = 1500 s/ mm 2 , and 7 images without diffusion weighting) of the entire in vivo human brain with nominal 0.66 mm spatial resolution. Using data acquired from 75 minutes of acquisition as a gold standard reference, we demonstrate that the proposed SNR-enhancement procedure leads to substantial improvements in estimated diffusion parameters compared to conventional gSlider reconstruction. Results also demonstrate that the proposed method has advantages relative to denoising methods based on low-rank matrix modeling. A theoretical analysis of the trade-off between spatial resolution and SNR suggests that the proposed approach has high efficiency.

CONCLUSIONS

The combination of gSlider-SMS with advanced regularized reconstruction enables high-resolution quantitative diffusion MRI from a relatively fast acquisition.

Three-dimensional simultaneous brain T1 , T2 , and ADC mapping with MR Multitasking

PURPOSE

To develop a simultaneous T₁ , T₂ , and ADC mapping method that provides co-registered, distortion-free images and enables multiparametric quantification of 3D brain coverage in a clinically feasible scan time with the MR Multitasking framework.

METHODS

The T₁ /T₂ /diffusion weighting was generated by a series of T₂ preparations and diffusion preparations. The underlying multidimensional image containing 3 spatial dimensions, 1 T₁ weighting dimension, 1 T₂ -preparation duration dimension, 1 b-value dimension, and 1 diffusion direction dimension was modeled as a 5-way low-rank tensor. A separate real-time low-rank model incorporating time-resolved phase correction was also used to compensate for both inter-shot and intra-shot phase inconsistency induced by physiological motion. The proposed method was validated on both phantom and 16 healthy subjects. The quantification of T₁ /T₂ /ADC was evaluated for each case. Three post-surgery brain tumor patients were scanned for demonstration of clinical feasibility.

RESULTS

Multitasking T₁ /T₂ /ADC maps were perfectly co-registered and free from image distortion. Phantom studies showed substantial quantitative agreement ( R 2 = 0.999 ) with reference protocols for T₁ /T₂ /ADC. In vivo studies showed nonsignificant T₁ (P = .248), T₂ (P = .97), ADC (P = .328) differences among the frontal, parietal, and occipital regions. Although Multitasking showed significant differences of T₁ (P = .03), T₂ (P < .001), and ADC (P = .001) biases against the references, the mean bias estimates were small (ΔT₁ % < 5%, ΔT₂ % < 7%, ΔADC% < 5%), with all intraclass correlation coefficients greater than 0.82 indicating "excellent" agreement. Patient studies showed that Multitasking T₁ /T₂ /ADC maps were consistent with the clinical qualitative images.

CONCLUSION

The Multitasking approach simultaneously quantifies T₁ /T₂ /ADC with substantial agreement with the references and is promising for clinical applications.

Multi-shot diffusion-weighted MRI reconstruction with magnitude-based spatial-angular locally low-rank regularization (SPA-LLR)

PURPOSE

To resolve the motion-induced phase variations in multi-shot multi-direction diffusion-weighted imaging (DWI) by applying regularization to magnitude images.

THEORY AND METHODS

A nonlinear model was developed to estimate phase and magnitude images separately. A locally low-rank regularization (LLR) term was applied to the magnitude images from all diffusion-encoding directions to exploit the spatial and angular correlation. In vivo experiments with different resolutions and b-values were performed to validate the proposed method.

RESULTS

The proposed method significantly reduces the noise level compared to the conventional reconstruction method and achieves submillimeter (0.8mm and 0.9mm isotropic resolutions) DWI with a b-value of 1,000  s / mm 2 and 1-mm isotropic DWI with a b-value of 2,000  s / mm 2 without modification of the sequence.

CONCLUSIONS

A joint reconstruction method with spatial-angular LLR regularization on magnitude images substantially improves multi-direction DWI reconstruction, simultaneously removes motion-induced phase artifacts, and denoises images.

Accelerating cardiac diffusion tensor imaging combining local low-rank and 3D TV constraint

OBJECTIVE

Diffusion tensor magnetic resonance imaging (DT-MRI, or DTI) is a promising technique for invasively probing biological tissue structures. However, DTI is known to suffer from much longer acquisition time with respect to conventional MRI and the problem is worsened when dealing with in vivo acquisitions. Therefore, faster DTI for both ex vivo and in vivo scans is highly desired.

MATERIALS AND METHODS

This paper proposes a new compressed sensing (CS) reconstruction method that employs local low-rank (LLR) model and three-dimensional (3D) total variation (TV) constraint to reconstruct cardiac diffusion-weighted (DW) images from highly undersampled k-space data. The LLR model takes the set of DW images corresponding to different diffusion gradient directions as a 3D image volume and decomposes the latter into overlapping 3D blocks. Then, the 3D blocks are stacked as two-dimensional (2D) matrix. Finally, low-rank property is applied to each block matrix and the 3D TV constraint to the 3D image volume. The underlying constrained optimization problem is finally solved using the first-order fast method. The proposed method is evaluated on real ex vivo cardiac DTI data as a prerequisite to in vivo cardiac DTI applications.

RESULTS

The results on real human ex vivo cardiac DTI images demonstrate that the proposed method exhibits lower reconstruction errors for DTI indices, including fractional anisotropy (FA), mean diffusivities (MD), transverse angle (TA), and helix angle (HA), compared to existing CS-based DTI image reconstruction techniques.

CONCLUSION

The proposed method provides better reconstruction quality and more accurate DTI indices in comparison with the state-of-the-art CS-based DW image reconstruction methods.