Acceleration Strategies for MR-STAT: Achieving High-Resolution Reconstructions on a Desktop PC Within 3 Minutes

GPU-accelerated Bloch simulations and MR-STAT reconstructions using the Julia programming language

PURPOSE

MR-STAT is a relatively new multiparametric quantitative MRI technique in which quantitative paramater maps are obtained by solving a large-scale nonlinear optimization problem. Managing reconstruction times is one of the main challenges of MR-STAT. In this work we leverage GPU hardware to reduce MR-STAT reconstruction times. A highly optimized, GPU-compatible Bloch simulation toolbox is developed as part of this work that can be utilized for other quantitative MRI techniques as well.

METHODS

The Julia programming language was used to develop a flexible yet highly performant and GPU-compatible Bloch simulation toolbox called BlochSimulators.jl. The runtime performance of the toolbox is benchmarked against other Bloch simulation toolboxes. Furthermore, a (partially matrix-free) modification of a previously presented (matrix-free) MR-STAT reconstruction algorithm is proposed and implemented using the Julia language on GPU hardware. The proposed algorithm is combined with BlochSimulators.jl and the resulting MR-STAT reconstruction times on GPU hardware are compared to previously presented MR-STAT reconstruction times.

RESULTS

The BlochSimulators.jl package demonstrates superior runtime performance on both CPU and GPU hardware when compared to other existing Bloch simulation toolboxes. The GPU-accelerated partially matrix-free MR-STAT reconstruction algorithm, which relies on BlochSimulators.jl, allows for reconstructions of 68 seconds per two-dimensional (2D slice).

CONCLUSION

By combining the proposed Bloch simulation toolbox and the partially matrix-free reconstruction algorithm, 2D MR-STAT reconstructions can be performed in the order of one minute on a modern GPU card. The Bloch simulation toolbox can be utilized for other quantitative MRI techniques as well, for example for online dictionary generation for MR Fingerprinting.

A three-dimensional Magnetic Resonance Spin Tomography in Time-domain protocol for high-resolution multiparametric quantitative magnetic resonance imaging

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) is a multiparametric quantitative MR framework, which allows for simultaneously acquiring quantitative tissue parameters such as T1, T2, and proton density from one single short scan. A typical two-dimensional (2D) MR-STAT acquisition uses a gradient-spoiled, gradient-echo sequence with a slowly varying RF flip-angle train and Cartesian readouts, and the quantitative tissue maps are reconstructed by an iterative, model-based optimization algorithm. In this work, we design a three-dimensional (3D) MR-STAT framework based on previous 2D work, in order to achieve better image signal-to-noise ratio, higher though-plane resolution, and better tissue characterization. Specifically, we design a 7-min, high-resolution 3D MR-STAT sequence, and the corresponding two-step reconstruction algorithm for the large-scale dataset. To reduce the long acquisition time, Cartesian undersampling strategies such as SENSE are adopted in our transient-state quantitative framework. To reduce the computational burden, a data-splitting scheme is designed for decoupling the 3D reconstruction problem into independent 2D reconstructions. The proposed 3D framework is validated by numerical simulations, phantom experiments, and in vivo experiments. High-quality knee quantitative maps with 0.8 × 0.8 × 1.5 mm3 resolution and bilateral lower leg maps with 1.6 mm isotropic resolution can be acquired using the proposed 7-min acquisition sequence and the 3-min-per-slice decoupled reconstruction algorithm. The proposed 3D MR-STAT framework could have wide clinical applications in the future.

Synthetic MRI with Magnetic Resonance Spin TomogrAphy in Time-Domain (MR-STAT): Results from a Prospective Cross-Sectional Clinical Trial

BACKGROUND

Magnetic Resonance Spin TomogrAphy in Time-domain (MR-STAT) can reconstruct whole-brain multi-parametric quantitative maps (eg, T₁ , T₂ ) from a 5-minute MR acquisition. These quantitative maps can be leveraged for synthetization of clinical image contrasts.

PURPOSE

The objective was to assess image quality and overall diagnostic accuracy of synthetic MR-STAT contrasts compared to conventional contrast-weighted images.

STUDY TYPE

Prospective cross-sectional clinical trial.

POPULATION

Fifty participants with a median age of 45 years (range: 21-79 years) consisting of 10 healthy participants and 40 patients with neurological diseases (brain tumor, epilepsy, multiple sclerosis or stroke).

FIELD STRENGTH/SEQUENCE

3T/Conventional contrast-weighted imaging (T₁ /T₂ weighted, proton density [PD] weighted, and fluid-attenuated inversion recovery [FLAIR]) and a MR-STAT acquisition (2D Cartesian spoiled gradient echo with varying flip angle preceded by a non-selective inversion pulse).

ASSESSMENT

Quantitative T₁ , T₂ , and PD maps were computed from the MR-STAT acquisition, from which synthetic contrasts were generated. Three neuroradiologists blinded for image type and disease randomly and independently evaluated synthetic and conventional datasets for image quality and diagnostic accuracy, which was assessed by comparison with the clinically confirmed diagnosis.

STATISTICAL TESTS

Image quality and consequent acceptability for diagnostic use was assessed with a McNemar's test (one-sided α = 0.025). Wilcoxon signed rank test with a one-sided α = 0.025 and a margin of Δ = 0.5 on the 5-level Likert scale was used to assess non-inferiority.

RESULTS

All data sets were similar in acceptability for diagnostic use (≥3 Likert-scale) between techniques (T₁ w:P = 0.105, PDw:P = 1.000, FLAIR:P = 0.564). However, only the synthetic MR-STAT T₂ weighted images were significantly non-inferior to their conventional counterpart; all other synthetic datasets were inferior (T₁ w:P = 0.260, PDw:P = 1.000, FLAIR:P = 1.000). Moreover, true positive/negative rates were similar between techniques (conventional: 88%, MR-STAT: 84%).

DATA CONCLUSION

MR-STAT is a quantitative technique that may provide radiologists with clinically useful synthetic contrast images within substantially reduced scan time.

EVIDENCE LEVEL

1 Technical Efficacy: Stage 2.