Split Bregman multicoil accelerated reconstruction technique: A new framework for rapid reconstruction of cardiac perfusion MRI

Signal intensity informed multi-coil encoding operator for physics-guided deep learning reconstruction of highly accelerated myocardial perfusion CMR

PURPOSE

To develop a physics-guided deep learning (PG-DL) reconstruction strategy based on a signal intensity informed multi-coil (SIIM) encoding operator for highly-accelerated simultaneous multislice (SMS) myocardial perfusion cardiac MRI (CMR).

METHODS

First-pass perfusion CMR acquires highly-accelerated images with dynamically varying signal intensity/SNR following the administration of a gadolinium-based contrast agent. Thus, using PG-DL reconstruction with a conventional multi-coil encoding operator leads to analogous signal intensity variations across different time-frames at the network output, creating difficulties in generalization for varying SNR levels. We propose to use a SIIM encoding operator to capture the signal intensity/SNR variations across time-frames in a reformulated encoding operator. This leads to a more uniform/flat contrast at the output of the PG-DL network, facilitating generalizability across time-frames. PG-DL reconstruction with the proposed SIIM encoding operator is compared to PG-DL with conventional encoding operator, split slice-GRAPPA, locally low-rank (LLR) regularized reconstruction, low-rank plus sparse (L + S) reconstruction, and regularized ROCK-SPIRiT.

RESULTS

Results on highly accelerated free-breathing first pass myocardial perfusion CMR at three-fold SMS and four-fold in-plane acceleration show that the proposed method improves upon the reconstruction methods use for comparison. Substantial noise reduction is achieved compared to split slice-GRAPPA, and aliasing artifacts reduction compared to LLR regularized reconstruction, L + S reconstruction and PG-DL with conventional encoding. Furthermore, a qualitative reader study indicated that proposed method outperformed all methods.

CONCLUSION

PG-DL reconstruction with the proposed SIIM encoding operator improves generalization across different time-frames /SNRs in highly accelerated perfusion CMR.

Quantitative susceptibility mapping using plug-and-play alternating direction method of multipliers

Quantitative susceptibility mapping employs regularization to reduce artifacts, yet many recent denoisers are unavailable for reconstruction. We developed a plug-and-play approach to QSM reconstruction (PnP QSM) and show its flexibility using several patch-based denoisers. We developed PnP QSM using alternating direction method of multiplier framework and applied collaborative filtering denoisers. We apply the technique to the 2016 QSM Challenge and in 10 glioblastoma multiforme datasets. We compared its performance with four published QSM techniques and a multi-orientation QSM method. We analyzed magnetic susceptibility accuracy using brain region-of-interest measurements, and image quality using global error metrics. Reconstructions on glioblastoma data were analyzed using ranked and semiquantitative image grading by three neuroradiologist observers to assess image quality (IQ) and sharpness (IS). PnP-BM4D QSM showed good correlation (β = 0.84, R² = 0.98, p < 0.05) with COSMOS and no significant bias (bias = 0.007 ± 0.012). PnP-BM4D QSM achieved excellent quality when assessed using structural similarity index metric (SSIM = 0.860), high frequency error norm (HFEN = 58.5), cross correlation (CC = 0.804), and mutual information (MI = 0.475) and also maintained good conspicuity of fine features. In glioblastoma datasets, PnP-BM4D QSM showed higher performance (IQ_Grade = 2.4 ± 0.4, IS_Grade = 2.7 ± 0.3, IQ_Rank = 3.7 ± 0.3, IS_Rank = 3.9 ± 0.3) compared to MEDI (IQ_Grade = 2.1 ± 0.5, IS_Grade = 2.1 ± 0.6, IQ_Rank = 2.4 ± 0.6, IS_Rank = 2.9 ± 0.2) and FANSI-TGV (IQ_Grade = 2.2 ± 0.6, IS_Grade = 2.1 ± 0.6, IQ_Rank = 2.7 ± 0.3, IS_Rank = 2.2 ± 0.2). We illustrated the modularity of PnP QSM by interchanging two additional patch-based denoisers. PnP QSM reconstruction was feasible, and its flexibility was shown using several patch-based denoisers. This technique may allow rapid prototyping and validation of new denoisers for QSM reconstruction for an array of useful clinical applications.

Closed-loop control of k-space sampling via physiologic feedback for cine MRI

BACKGROUND

Segmented cine cardiac MRI combines data from multiple heartbeats to achieve high spatiotemporal resolution cardiac images, yet predefined k-space segmentation trajectories can lead to suboptimal k-space sampling. In this work, we developed and evaluated an autonomous and closed-loop control system for radial k-space sampling (ARKS) to increase sampling uniformity.

METHODS

The closed-loop system autonomously selects radial k-space sampling trajectory during live segmented cine MRI and attempts to optimize angular sampling uniformity by selecting views in regions of k-space that were not previously well-sampled. Sampling uniformity and the ability to detect cardiac phase in vivo was assessed using ECG data acquired from 10 normal subjects in an MRI scanner. The approach was then implemented with a fast gradient echo sequence on a whole-body clinical MRI scanner and imaging was performed in 4 healthy volunteers. The closed-loop k-space trajectory was compared to random, uniformly distributed and golden angle view trajectories via measurement of k-space uniformity and the point spread function. Lastly, an arrhythmic dataset was used to evaluate a potential application of the approach.

RESULTS

The autonomous trajectory increased k-space sampling uniformity by 15±7%, main lobe point spread function (PSF) signal intensity by 6±4%, and reduced ringing relative to golden angle sampling. When implemented, the autonomous pulse sequence prescribed radial view angles faster than the scan TR (0.98 ± 0.01 ms, maximum = 1.38 ms) and increased k-space sampling mean uniformity by 10±11%, decreased uniformity variability by 44±12%, and increased PSF signal ratio by 6±6% relative to golden angle sampling.

CONCLUSION

The closed-loop approach enables near-uniform radial sampling in a segmented acquisition approach which was higher than predetermined golden-angle radial sampling. This can be utilized to increase the sampling or decrease the temporal footprint of an acquisition and the closed-loop framework has the potential to be applied to patients with complex heart rhythms.

Quantitative 3D myocardial perfusion with an efficient arterial input function

PURPOSE

The purpose of this study was to further develop and combine several innovative sequence designs to achieve quantitative 3D myocardial perfusion. These developments include an optimized 3D stack-of-stars readout (150 ms per beat), efficient acquisition of a 2D arterial input function, tailored saturation pulse design, and potential whole heart coverage during quantitative stress perfusion.

THEORY AND METHODS

All studies were performed free-breathing on a Prisma 3T MRI scanner. Phantom validation was used to verify sequence accuracy. A total of 21 subjects (3 patients with known disease) were scanned, 12 with a rest only protocol and 9 with both stress (regadenoson) and rest protocols. First pass quantitative perfusion was performed with gadoteridol (0.075 mmol/kg).

RESULTS

Implementation and quantitative perfusion results are shown for healthy subjects and subjects with known coronary disease. Average rest perfusion for the 15 included healthy subjects was 0.79 ± 0.19 mL/g/min, the average stress perfusion for 6 healthy subject studies was 2.44 ± 0.61 mL/g/min, and the average global myocardial perfusion reserve ratio for 6 healthy subjects was 3.10 ± 0.24. Perfusion deficits for 3 patients with ischemia are shown. Average resting heart rate was 59 ± 7 bpm and the average stress heart rate was 81 ± 10 bpm.

CONCLUSION

This work demonstrates that a quantitative 3D myocardial perfusion sequence with the acquisition of a 2D arterial input function is feasible at high stress heart rates such as during stress. T1 values and gadolinium concentrations of the sequence match the reference standard well in a phantom, and myocardial rest and stress perfusion and myocardial perfusion reserve values are consistent with those published in literature.

Data-Driven Quantitative Susceptibility Mapping Using Loss Adaptive Dipole Inversion (LADI)

BACKGROUND

Quantitative susceptibility mapping (QSM) uses prior information to reconstruct maps, but prior information may not show pathology and introduce inconsistencies with susceptibility maps, degrade image quality and inadvertently smoothing image features.

PURPOSE

To develop a local field data-driven QSM reconstruction that does not depend on spatial edge prior information.

STUDY TYPE

Retrospective.

SUBJECTS, ANIMAL MODELS

A dataset from 2016 ISMRM QSM Challenge, 11 patients with glioblastoma, a patient with microbleeds and porcine heart.

SEQUENCE/FIELD STRENGTH

3D gradient echo sequence on 3T and 7T scanners.

ASSESSMENT

Accuracy was compared to Calculation of Susceptibility through Multiple Orientation Sampling (COSMOS), and several published techniques using region of interest (ROI) measurements, root-mean-squared error (RMSE), structural similarity index metric (SSIM), and high-frequency error norm (HFEN). Numerical ranking and semiquantitative image grading was performed by three expert observers to assess overall image quality (IQ) and image sharpness (IS).

STATISTICAL TESTS

Bland-Altman, Friedman test, and Conover multiple comparisons.

RESULTS

Loss adaptive dipole inversion (LADI) (β = 0.82, R² = 0.96), morphology-enabled dipole inversion (MEDI) (β = 0.91, R² = 0.97), and fast nonlinear susceptibility inversion (FANSI) (β = 0.81, R² = 0.98) had excellent correlation with COSMOS and no bias was detected (bias = 0.006 ± 0.014, P < 0.05). In glioblastoma patients, LADI showed consistently better performance (IQ_Grade = 2.6 ± 0.4, IS_Grade = 2.6 ± 0.3, IQ_Rank = 3.5 ± 0.4, IS_Rank = 3.9 ± 0.2) compared with MEDI (IQ_Grade = 2.1 ± 0.3, IS_Grade = 2 ± 0.5, IQ_Rank = 2.4 ± 0.5, IS_Rank = 2.8 ± 0.2) and FANSI (IQ_Grade = 2.2 ± 0.5, IS_Grade = 2 ± 0.4, IQ_Rank = 2.8 ± 0.3, IS_Rank = 2.1 ± 0.2). Dark artifact visible near the infarcted region in MEDI (Inf_MEDI = -0.27 ± 0.06 ppm) was better mitigated by FANSI (Inf_FANSI-TGV = -0.17 ± 0.05 ppm) and LADI (Inf_LADI = -0.18 ± 0.05 ppm).

CONCLUSION

For neuroimaging applications, LADI preserved image sharpness and fine features in glioblastoma and microbleed patients. LADI performed better at mitigating artifacts in cardiac QSM.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY STAGE: 1 J. Magn. Reson. Imaging 2020;52:823-835.

Efficient Dynamic Parallel MRI Reconstruction for the Low-Rank Plus Sparse Model

The low-rank plus sparse (L+S) decomposition model enables the reconstruction of under-sampled dynamic parallel magnetic resonance imaging (MRI) data. Solving for the low-rank and the sparse components involves non-smooth composite convex optimization, and algorithms for this problem can be categorized into proximal gradient methods and variable splitting methods. This paper investigates new efficient algorithms for both schemes. While current proximal gradient techniques for the L+S model involve the classical iterative soft thresholding algorithm (ISTA), this paper considers two accelerated alternatives, one based on the fast iterative shrinkage-thresholding algorithm (FISTA), and the other with the recent proximal optimized gradient method (POGM). In the augmented Lagrangian (AL) framework, we propose an efficient variable splitting scheme based on the form of the data acquisition operator, leading to simpler computation than the conjugate gradient (CG) approach required by existing AL methods. Numerical results suggest faster convergence of the efficient implementations for both frameworks, with POGM providing the fastest convergence overall and the practical benefit of being free of algorithm tuning parameters.

Accelerated free-breathing 3D T1ρ cardiovascular magnetic resonance using multicoil compressed sensing

BACKGROUND

Endogenous contrast T1ρ cardiovascular magnetic resonance (CMR) can detect scar or infiltrative fibrosis in patients with ischemic or non-ischemic cardiomyopathy. Existing 2D T1ρ techniques have limited spatial coverage or require multiple breath-holds. The purpose of this project was to develop an accelerated, free-breathing 3D T1ρ mapping sequence with whole left ventricle coverage using a multicoil, compressed sensing (CS) reconstruction technique for rapid reconstruction of undersampled k-space data.

METHODS

We developed a cardiac- and respiratory-gated, free-breathing 3D T1ρ sequence and acquired data using a variable-density k-space sampling pattern (A = 3). The effect of the transient magnetization trajectory, incomplete recovery of magnetization between T1ρ-preparations (heart rate dependence), and k-space sampling pattern on T1ρ relaxation time error and edge blurring was analyzed using Bloch simulations for normal and chronically infarcted myocardium. Sequence accuracy and repeatability was evaluated using MnCl2 phantoms with different T1ρ relaxation times and compared to 2D measurements. We further assessed accuracy and repeatability in healthy subjects and compared these results to 2D breath-held measurements.

RESULTS

The error in T1ρ due to incomplete recovery of magnetization between T1ρ-preparations was T1ρhealthy = 6.1% and T1ρinfarct = 10.8% at 60 bpm and T1ρhealthy = 13.2% and T1ρinfarct = 19.6% at 90 bpm. At a heart rate of 60 bpm, error from the combined effects of readout-dependent magnetization transients, k-space undersampling and reordering was T1ρhealthy = 12.6% and T1ρinfarct = 5.8%. CS reconstructions had improved edge sharpness (blur metric = 0.15) compared to inverse Fourier transform reconstructions (blur metric = 0.48). There was strong agreement between the mean T1ρ estimated from the 2D and accelerated 3D data (R2 = 0.99; P < 0.05) acquired on the MnCl2 phantoms. The mean R1ρ estimated from the accelerated 3D sequence was highly correlated with MnCl2 concentration (R2 = 0.99; P < 0.05). 3D T1ρ acquisitions were successful in all human subjects. There was no significant bias between undersampled 3D T1ρ and breath-held 2D T1ρ (mean bias = 0.87) and the measurements had good repeatability (COV2D = 6.4% and COV3D = 7.1%).

CONCLUSIONS

This is the first report of an accelerated, free-breathing 3D T1ρ mapping of the left ventricle. This technique may improve non-contrast myocardial tissue characterization in patients with heart disease in a scan time appropriate for patients.

Self-gated MRI of multiple beat morphologies in the presence of arrhythmias

PURPOSE

Develop self-gated MRI for distinct heartbeat morphologies in subjects with arrhythmias.

METHODS

Golden angle radial data was obtained in seven sinus and eight arrhythmias subjects. An image-based cardiac navigator was derived from single-shot images, distinct beat types were identified, and images were reconstructed for repeated morphologies. Image sharpness, contrast, and volume variation were quantified and compared with self-gated MRI. Images were scored for image quality and artifacts. Hemodynamic parameters were computed for each distinct beat morphology in bigeminy and trigeminy subjects and for sinus beats in patients with infrequent premature ventricular contractions.

RESULTS

Images of distinct beat types were reconstructed except for two patients with infrequent premature ventricular contractions. Image contrast and sharpness were similar to sinus self-gated images (contrast = 0.45 ± 0.13 and 0.43 ± 0.15; sharpness = 0.21 ± 0.11 and 0.20 ± 0.05). Visual scoring was highest in self-gated images (4.1 ± 0.3) compared with real-time (3.9 ± 0.4) and ECG-gated cine (3.4 ± 1.5). ECG-gated cine had less artifacts than self-gating (2.3 ± 0.7 and 2.1 ± 0.2), but was affected by misgating in two subjects. Among arrhythmia subjects, post-extrasystole/sinus (58.1 ± 8.6 mL) and interrupted sinus (61.4 ± 5.9 mL) stroke volume was higher than extrasystole (32.0 ± 16.5 mL; P < 0.02).

CONCLUSION

Self-gated imaging can reconstruct images during ectopy and allowed for quantification of hemodynamic function of different beat morphologies. Magn Reson Med 78:678-688, 2017. © 2016 International Society for Magnetic Resonance in Medicine.