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, R2 = 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 (IQGrade = 2.4 ± 0.4, ISGrade = 2.7 ± 0.3, IQRank = 3.7 ± 0.3, ISRank = 3.9 ± 0.3) compared to MEDI (IQGrade = 2.1 ± 0.5, ISGrade = 2.1 ± 0.6, IQRank = 2.4 ± 0.6, ISRank = 2.9 ± 0.2) and FANSI-TGV (IQGrade = 2.2 ± 0.6, ISGrade = 2.1 ± 0.6, IQRank = 2.7 ± 0.3, ISRank = 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.

Endogenous T1ρ cardiovascular magnetic resonance in hypertrophic cardiomyopathy.

BACKGROUND: Hypertrophic cardiomyopathy (HCM) is characterized by increased left ventricular wall thickness, cardiomyocyte hypertrophy, and fibrosis. Adverse cardiac risk characterization has been performed using late gadolinium enhancement (LGE), native T1, and extracellular volume (ECV). Relaxation time constants are affected by background field inhomogeneity. T1ρ utilizes a spin-lock pulse to decrease the effect of unwanted relaxation. The objective of this study was to study T1ρ as compared to T1, ECV, and LGE in HCM patients. METHODS: HCM patients were recruited as part of the Novel Markers of Prognosis in Hypertrophic Cardiomyopathy study, and healthy controls were matched for comparison. In addition to cardiac functional imaging, subjects underwent T1 and T1ρ cardiovascular magnetic resonance imaging at short-axis positions at 1.5T. Subjects received gadolinium and underwent LGE imaging 15-20 min after injection covering the entire heart. Corresponding basal and mid short axis LGE slices were selected for comparison with T1 and T1ρ. Full-width half-maximum thresholding was used to determine the percent enhancement area in each LGE-positive slice by LGE, T1, and T1ρ. Two clinicians independently reviewed LGE images for presence or absence of enhancement. If in agreement, the image was labeled positive (LGE + +) or negative (LGE --); otherwise, the image was labeled equivocal (LGE + -). RESULTS: In 40 HCM patients and 10 controls, T1 percent enhancement area (Spearman's rho = 0.61, p < 1e-5) and T1ρ percent enhancement area (Spearman's rho = 0.48, p < 0.001e-3) correlated with LGE percent enhancement area. T1 and T1ρ percent enhancement areas were also correlated (Spearman's rho = 0.28, p = 0.047). For both T1 and T1ρ, HCM patients demonstrated significantly longer relaxation times compared to controls in each LGE category (p < 0.001 for all). HCM patients also showed significantly higher ECV compared to controls in each LGE category (p < 0.01 for all), and LGE -- slices had lower ECV than LGE + + (p = 0.01). CONCLUSIONS: Hyperenhancement areas as measured by T1ρ and LGE are moderately correlated. T1, T1ρ, and ECV were elevated in HCM patients compared to controls, irrespective of the presence of LGE. These findings warrant additional studies to investigate the prognostic utility of T1ρ imaging in the evaluation of HCM patients.

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.

MR Intracranial Vessel Wall Imaging: A Systematic Review.

The purpose of this systematic review is to identify trends and extent of variability in intracranial vessel wall MR imaging (VWI) techniques and protocols. Although variability in selection of protocol design and pulse sequence type is known, data on what and how protocols vary are unknown. Three databases were searched to identify publications using intracranial VWI. Publications were screened by predetermined inclusion/exclusion criteria. Technical development publications were scored for completeness of reporting using a modified Nature Reporting Summary Guideline to assess reproducibility. From 2,431 articles, 122 met the inclusion criteria. Trends over the last 23 years (1995-2018) show increased use of 3-Tesla MR (P < .001) and 3D volumetric T1-weighted acquisitions (P < .001). Most (65%) clinical VWI publications report achieving a noninterpolated in-plane spatial resolution of ≤.55 mm. In the last decade, an increasing number of technical development (n = 20) and 7 Tesla (n = 12) publications have been published, focused on pulse sequence development, improving cerebrospinal fluid suppression, scan efficiency, and imaging ex vivo specimen for histologic validation. Mean Reporting Summary Score for the technical development publications was high (.87, range: .63-1.0) indicating strong scientific technical reproducibility. Innovative work continues to emerge to address implementation challenges. Gradual adoption into the research and scientific community was suggested by a shift in the name in the literature from "high-resolution MR" to "vessel wall imaging," specifying diagnostic intent. Insight into current practices and identifying the extent of technical variability in the literature will help to direct future clinical and technical efforts to address needs for implementation.

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, R2 = 0.96), morphology-enabled dipole inversion (MEDI) (ß = 0.91, R2 = 0.97), and fast nonlinear susceptibility inversion (FANSI) (ß = 0.81, R2 = 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 (IQGrade = 2.6 ± 0.4, ISGrade = 2.6 ± 0.3, IQRank = 3.5 ± 0.4, ISRank = 3.9 ± 0.2) compared with MEDI (IQGrade = 2.1 ± 0.3, ISGrade = 2 ± 0.5, IQRank = 2.4 ± 0.5, ISRank = 2.8 ± 0.2) and FANSI (IQGrade = 2.2 ± 0.5, ISGrade = 2 ± 0.4, IQRank = 2.8 ± 0.3, ISRank = 2.1 ± 0.2). Dark artifact visible near the infarcted region in MEDI (InfMEDI = -0.27 ± 0.06 ppm) was better mitigated by FANSI (InfFANSI-TGV = -0.17 ± 0.05 ppm) and LADI (InfLADI = -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.

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.

Technical Note: Evaluation of pre-reconstruction interpolation methods for iterative reconstruction of radial k-space data.

PURPOSE: To evaluate the use of three different pre-reconstruction interpolation methods to convert non-Cartesian k-space data to Cartesian samples such that iterative reconstructions can be performed more simply and more rapidly. METHODS: Phantom as well as cardiac perfusion radial datasets were reconstructed by four different methods. Three of the methods used pre-reconstruction interpolation once followed by a fast Fourier transform (FFT) at each iteration. The methods were: bilinear interpolation of nearest-neighbor points (BINN), 3-point interpolation, and a multi-coil interpolator called GRAPPA Operator Gridding (GROG). The fourth method performed a full non-Uniform FFT (NUFFT) at each iteration. An iterative reconstruction with spatiotemporal total variation constraints was used with each method. Differences in the images were quantified and compared. RESULTS: The GROG multicoil interpolation, the 3-point interpolation, and the NUFFT-at-each-iteration approaches produced high quality images compared to BINN, with the GROG-derived images having the fewest streaks among the three preinterpolation approaches. However, all reconstruction methods produced approximately equal results when applied to perfusion quantitation tasks. Pre-reconstruction interpolation gave approximately an 83% reduction in reconstruction time. CONCLUSION: Image quality suffers little from using a pre-reconstruction interpolation approach compared to the more accurate NUFFT-based approach. GROG-based pre-reconstruction interpolation appears to offer the best compromise by using multicoil information to perform the interpolation to Cartesian sample points prior to image reconstruction. Speed gains depend on the implementation and relatively standard optimizations on a MATLAB platform result in preinterpolation speedups of ~ 6 compared to using NUFFT at every iteration, reducing the reconstruction time from around 42 min to 7 min.

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

PURPOSE: Rapid reconstruction of undersampled multicoil MRI data with iterative constrained reconstruction method is a challenge. The authors sought to develop a new substitution based variable splitting algorithm for faster reconstruction of multicoil cardiac perfusion MRI data. METHODS: The new method, split Bregman multicoil accelerated reconstruction technique (SMART), uses a combination of split Bregman based variable splitting and iterative reweighting techniques to achieve fast convergence. Total variation constraints are used along the spatial and temporal dimensions. The method is tested on nine ECG-gated dog perfusion datasets, acquired with a 30-ray golden ratio radial sampling pattern and ten ungated human perfusion datasets, acquired with a 24-ray golden ratio radial sampling pattern. Image quality and reconstruction speed are evaluated and compared to a gradient descent (GD) implementation and to multicoil k-t SLR, a reconstruction technique that uses a combination of sparsity and low rank constraints. RESULTS: Comparisons based on blur metric and visual inspection showed that SMART images had lower blur and better texture as compared to the GD implementation. On average, the GD based images had an ∼18% higher blur metric as compared to SMART images. Reconstruction of dynamic contrast enhanced (DCE) cardiac perfusion images using the SMART method was ∼6 times faster than standard gradient descent methods. k-t SLR and SMART produced images with comparable image quality, though SMART was ∼6.8 times faster than k-t SLR. CONCLUSIONS: The SMART method is a promising approach to reconstruct good quality multicoil images from undersampled DCE cardiac perfusion data rapidly.

Compressed sensing for rapid late gadolinium enhanced imaging of the left atrium: A preliminary study.

Current late gadolinium enhancement (LGE) imaging of left atrial (LA) scar or fibrosis is relatively slow and requires 5-15min to acquire an undersampled (R=1.7) 3D navigated dataset. The GeneRalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) based parallel imaging method is the current clinical standard for accelerating 3D LGE imaging of the LA and permits an acceleration factor ~R=1.7. Two compressed sensing (CS) methods have been developed to achieve higher acceleration factors: a patch based collaborative filtering technique tested with acceleration factor R~3, and a technique that uses a 3D radial stack-of-stars acquisition pattern (R~1.8) with a 3D total variation constraint. The long reconstruction time of these CS methods makes them unwieldy to use, especially the patch based collaborative filtering technique. In addition, the effect of CS techniques on the quantification of percentage of scar/fibrosis is not known. We sought to develop a practical compressed sensing method for imaging the LA at high acceleration factors. In order to develop a clinically viable method with short reconstruction time, a Split Bregman (SB) reconstruction method with 3D total variation (TV) constraints was developed and implemented. The method was tested on 8 atrial fibrillation patients (4 pre-ablation and 4 post-ablation datasets). Blur metric, normalized mean squared error and peak signal to noise ratio were used as metrics to analyze the quality of the reconstructed images, Quantification of the extent of LGE was performed on the undersampled images and compared with the fully sampled images. Quantification of scar from post-ablation datasets and quantification of fibrosis from pre-ablation datasets showed that acceleration factors up to R~3.5 gave good 3D LGE images of the LA wall, using a 3D TV constraint and constrained SB methods. This corresponds to reducing the scan time by half, compared to currently used GRAPPA methods. Reconstruction of 3D LGE images using the SB method was over 20 times faster than standard gradient descent methods.

Edge-enhanced spatiotemporal constrained reconstruction of undersampled dynamic contrast-enhanced radial MRI.

Dynamic contrast-enhanced magnetic resonance imaging (MRI) is a technique used to study and track contrast kinetics in an area of interest in the body over time. Reconstruction of images with high contrast and sharp edges from undersampled data is a challenge. While good results have been reported using a radial acquisition and a spatiotemporal constrained reconstruction (STCR) method, we propose improvements from using spatially adaptive weighting and an additional edge-based constraint. The new method uses intensity gradients from a sliding window reference image to improve the sharpness of edges in the reconstructed image. The method was tested on eight radial cardiac perfusion data sets with 24 rays and compared to the STCR method. The reconstructions showed that the new method, termed edge-enhanced spatiotemporal constrained reconstruction, was able to reconstruct images with sharper edges, and there were a 36%±13.7% increase in contrast-to-noise ratio and a 24%±11% increase in contrast near the edges when compared to STCR. The novelty of this paper is the combination of spatially adaptive weighting for spatial total variation (TV) constraint along with a gradient matching term to improve the sharpness of edges. The edge map from a reference image allows the reconstruction to trade-off between TV and edge enhancement, depending on the spatially varying weighting provided by the edge map.