Free-running 3D whole heart myocardial T1 mapping with isotropic spatial resolution

Free-breathing 3D whole-heart joint T1/T2 mapping and water/fat imaging at 0.55 T

PURPOSE

To develop and validate a highly efficient motion compensated free-breathing isotropic resolution 3D whole-heart joint T1/T2 mapping sequence with anatomical water/fat imaging at 0.55 T.

METHODS

The proposed sequence takes advantage of shorter T1 at 0.55 T to acquire three interleaved water/fat volumes with inversion-recovery preparation, no preparation, and T2 preparation, respectively. Image navigators were used to facilitate nonrigid motion-compensated image reconstruction. T1 and T2 maps were jointly calculated by a dictionary matching method. Validations were performed with simulation, phantom, and in vivo experiments on 10 healthy volunteers and 1 patient. The performance of the proposed sequence was compared with conventional 2D mapping sequences including modified Look-Locker inversion recovery and T2-prepared balanced steady-SSFP sequence.

RESULTS

The proposed sequence has a good T1 and T2 encoding sensitivity in simulation, and excellent agreement with spin-echo reference T1 and T2 values was observed in a standardized T1/T2 phantom (R2 = 0.99). In vivo experiments provided good-quality co-registered 3D whole-heart T1 and T2 maps with 2-mm isotropic resolution in a short scan time of about 7 min. For healthy volunteers, left-ventricle T1 mean and SD measured by the proposed sequence were both comparable with those of modified Look-Locker inversion recovery (640 ± 35 vs. 630 ± 25 ms [p = 0.44] and 49.9 ± 9.3 vs. 54.4 ± 20.5 ms [p = 0.42]), whereas left-ventricle T2 mean and SD measured by the proposed sequence were both slightly lower than those of T2-prepared balanced SSFP (53.8 ± 5.5 vs. 58.6 ± 3.3 ms [p < 0.01] and 5.2 ± 0.9 vs. 6.1 ± 0.8 ms [p = 0.03]). Myocardial T1 and T2 in the patient measured by the proposed sequence were in good agreement with conventional 2D sequences and late gadolinium enhancement.

CONCLUSION

The proposed sequence simultaneously acquires 3D whole-heart T1 and T2 mapping with anatomical water/fat imaging at 0.55 T in a fast and efficient 7-min scan. Further investigation in patients with cardiovascular disease is now warranted.

Free-breathing, fat-corrected T1 mapping of the liver with stack-of-stars MRI, and joint estimation of T1, PDFF, R 2 * , and B 1 +

PURPOSE

Quantitative T1 mapping has the potential to replace biopsy for noninvasive diagnosis and quantitative staging of chronic liver disease. Conventional T1 mapping methods are confounded by fat andB 1 + $$ {B}_1^{+} $$ inhomogeneities, resulting in unreliable T1 estimations. Furthermore, these methods trade off spatial resolution and volumetric coverage for shorter acquisitions with only a few images obtained within a breath-hold. This work proposes a novel, volumetric (3D), free-breathing T1 mapping method to account for multiple confounding factors in a single acquisition.

THEORY AND METHODS

Free-breathing, confounder-corrected T1 mapping was achieved through the combination of non-Cartesian imaging, magnetization preparation, chemical shift encoding, and a variable flip angle acquisition. A subspace-constrained, locally low-rank image reconstruction algorithm was employed for image reconstruction. The accuracy of the proposed method was evaluated through numerical simulations and phantom experiments with a T1/proton density fat fraction phantom at 3.0 T. Further, the feasibility of the proposed method was investigated through contrast-enhanced imaging in healthy volunteers, also at 3.0 T.

RESULTS

The method showed excellent agreement with reference measurements in phantoms across a wide range of T1 values (200 to 1000 ms, slope = 0.998 (95% confidence interval (CI) [0.963 to 1.035]), intercept = 27.1 ms (95% CI [0.4 54.6]), r2 = 0.996), and a high level of repeatability. In vivo imaging studies demonstrated moderate agreement (slope = 1.099 (95% CI [1.067 to 1.132]), intercept = -96.3 ms (95% CI [-82.1 to -110.5]), r2 = 0.981) compared to saturation recovery-based T1 maps.

CONCLUSION

The proposed method produces whole-liver, confounder-corrected T1 maps through simultaneous estimation of T1, proton density fat fraction, andB 1 + $$ {B}_1^{+} $$ in a single, free-breathing acquisition and has excellent agreement with reference measurements in phantoms.

3D whole heart k-space-based super-resolution cardiac T1 mapping using rotated stacks

Objective. To provide three-dimensional (3D) whole-heart high-resolution isotropic cardiac T1 maps using a k-space-based through-plane super-resolution reconstruction (SRR) with rotated multi-slice stacks.Approach. Due to limited SNR and cardiac motion, often only 2D T1 maps with low through-plane resolution (4-8 mm) can be obtained. Previous approaches used SRR to calculate 3D high-resolution isotropic cardiac T1 maps. However, they were limited to the ventricles. The proposed approach acquires rotated stacks in long-axis orientation with high in-plane resolution but low through-plane resolution. This results in radially overlapping stacks from which high-resolution T1 maps of the whole heart are reconstructed using a k-space-based SRR framework considering the complete acquisition model. Cardiac and residual respiratory motion between different breath holds is estimated and incorporated into the reconstruction. The proposed approach was evaluated in simulations and phantom experiments and successfully applied to ten healthy subjects.Main results. 3D T1 maps of the whole heart were obtained in the same acquisition time as previous methods covering only the ventricles. T1 measurements were possible even for small structures, such as the atrial wall. The proposed approach provided accurate (P> 0.4;R2> 0.99) and precise T1 values (SD of 64.32 ± 22.77 ms in the proposed approach, 44.73 ± 31.9 ms in the reference). The edge sharpness of the T1 maps was increased by 6.20% and 4.73% in simulation and phantom experiments, respectively. Contrast-to-noise ratios between the septum and blood pool increased by 14.50% inin vivomeasurements with a k-space compared to an image-space-based SRR.Significance. The proposed approach provided whole-heart high-resolution 1.3 mm isotropic T1 maps in an overall acquisition time of approximately three minutes. Small structures, such as the atrial and right ventricular walls, could be visualized in the T1 maps.

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

Free-breathing myocardial T1 mapping using inversion-recovery radial FLASH and motion-resolved model-based reconstruction

PURPOSE

To develop a free-breathing myocardialT 1 $$ {\mathrm{T}}_1 $$ mapping technique using inversion-recovery (IR) radial fast low-angle shot (FLASH) and calibrationless motion-resolved model-based reconstruction.

METHODS

Free-running (free-breathing, retrospective cardiac gating) IR radial FLASH is used for data acquisition at 3T. First, to reduce the waiting time between inversions, an analytical formula is derived that takes the incompleteT 1 $$ {\mathrm{T}}_1 $$ recovery into account for an accurateT 1 $$ {\mathrm{T}}_1 $$ calculation. Second, the respiratory motion signal is estimated from the k-space center of the contrast varying acquisition using an adapted singular spectrum analysis (SSA-FARY) technique. Third, a motion-resolved model-based reconstruction is used to estimate both parameter and coil sensitivity maps directly from the sorted k-space data. Thus, spatiotemporal total variation, in addition to the spatial sparsity constraints, can be directly applied to the parameter maps. Validations are performed on an experimental phantom, 11 human subjects, and a young landrace pig with myocardial infarction.

RESULTS

In comparison to an IR spin-echo reference, phantom results confirm goodT 1 $$ {\mathrm{T}}_1 $$ accuracy, when reducing the waiting time from 5 s to 1 s using the new correction. The motion-resolved model-based reconstruction further improvesT 1 $$ {\mathrm{T}}_1 $$ precision compared to the spatial regularization-only reconstruction. Aside from showing that a reliable respiratory motion signal can be estimated using modified SSA-FARY, in vivo studies demonstrate that dynamic myocardialT 1 $$ {\mathrm{T}}_1 $$ maps can be obtained within 2 min with good precision and repeatability.

CONCLUSION

Motion-resolved myocardialT 1 $$ {\mathrm{T}}_1 $$ mapping during free-breathing with good accuracy, precision and repeatability can be achieved by combining inversion-recovery radial FLASH, self-gating and a calibrationless motion-resolved model-based reconstruction.

Simultaneous T1, T2, and T2* Mapping of Carotid Plaque: The SIMPLE* Technique

Background Quantitative T1, T2, and T2* measurements of carotid atherosclerotic plaque are important in evaluating plaque vulnerability and monitoring its progression. Purpose To develop a sequence to simultaneously quantify T1, T2, and T2* of carotid plaque. Materials and Methods The simultaneous T1, T2, and T2* mapping of carotid plaque (SIMPLE*) sequence is composed of three modules with different T2 preparation pulses, inversion-recovery pulses, and acquisition schemas. Single-echo data were used for T1 and T2 quantification, while the multiecho (ME) data were used for T2* quantification. The quantitative accuracy of SIMPLE* was tested in a phantom study by comparing its measurements with those of reference standard sequences. In vivo feasibility of the technique was prospectively evaluated between November 2020 and February 2022 in healthy volunteers and participants with carotid atherosclerotic plaque. The Pearson or Spearman correlation test, Student t test, and Wilcoxon rank-sum test were used. Results T1, T2, and T2* estimated with SIMPLE* strongly correlated with inversion-recovery spin-echo (SE) (correlation coefficient [r] = 0.99), ME-SE (r = 0.99), and ME gradient-echo (r = 0.99) sequences in the phantom study. In five healthy volunteers (mean age, 25 years ± 3 [SD]; three women), measurements were similar between SIMPLE* and modified Look-Locker inversion recovery, or MOLLI (1151 msec ± 71 vs 1098 msec ± 64; P = .14), ME turbo SE (31 msec ± 1 vs 31 msec ± 1; P = .32), and ME turbo field echo (24 msec ± 2 vs 25 msec ± 2; P = .19). In 18 participants with carotid plaque (mean age, 65 years ± 9; 16 men), quantitative T1, T2, and T2* of plaque components were consistent with their signal characteristics on multicontrast images. Conclusion A quantitative technique for simultaneous T1, T2, and T2* mapping of carotid plaque with 100-mm³ coverage and 0.8-mm³ resolution was developed using the proposed SIMPLE* sequence and demonstrated high accuracy and in vivo feasibility. © RSNA, 2023 Supplemental material is available for this article.

Free-breathing simultaneous water-fat separation and T1 mapping of the whole liver (SWALI) with isotropic resolution using 3D golden-angle radial trajectory

Background

Conventional liver T1 mapping techniques are typically performed under breath-holding conditions; they have limited slice coverage and often rely on multiple acquisitions. Furthermore, liver fat affects the accuracy of T1 quantification. Therefore, we aim to propose a free-breathing technique for simultaneous water-fat separation and T1 mapping of the whole liver (SWALI) in a single scan.

Methods

The proposed SWALI sequence included an inversion recovery (IR) preparation pulse followed by a series of multiecho three-dimensional (3D) golden-angle radial acquisitions. For each echo time (TE), a series of images containing a mix of water and fat were reconstructed using a sliding window method. For each inversion time (TI), water and fat were separated, and then water and fat T1 estimation was conducted. The fat fraction (FF) was calculated based on the last TI image. The FF and water T1 quantification accuracy were compared with the gold standard sequences in the phantom. The in vivo feasibility was tested in 9 healthy volunteers, 2 patients with fatty liver, and 3 patients with hepatocellular carcinoma (HCC). The reproducibility was evaluated in the patients with fatty liver and in the healthy volunteers.

Results

The mean FF and the mean water T1 values obtained by the SWALI sequence showed good agreements with chemical shift-encoded magnetic resonance imaging (CSE-MRI; r=0.998; P<0.001) and fat-suppressed (FS) IR-spin echo (SE; r=0.997; P<0.001) in the phantom. For the patients with fatty liver and the healthy volunteers, the SWALI sequence showed no significant difference with CSE-MRI in FF quantification (P=0.53). In T1 quantification, comparable T1 values were obtained with the SWALI sequence and modified Look-Locker inversion recovery (MOLLI; P=0.10) in healthy volunteers, while the water T1 estimated by the SWALI sequence was significantly lower than the water-fat compound T1 estimated by MOLLI (P<0.001) in patients with fatty liver. In the reproducibility study, the intraclass correlation coefficients (ICCs) for the estimated FF and water T1 were 0.997 and 0.943, respectively. Water T1 of the patients with HCC calculated using the SWALI sequence showed a significant reduction after the contrast administration (P<0.001).

Conclusions

Free-breathing water-fat separation and T1 mapping of the whole liver with 2.5 mm isotropic spatial resolution were achieved simultaneously using the SWALI sequence in a 5-min scan.

Free-running 3D whole-heart T1 and T2 mapping and cine MRI using low-rank reconstruction with non-rigid cardiac motion correction

PURPOSE

To introduce non-rigid cardiac motion correction into a novel free-running framework for the simultaneous acquisition of 3D whole-heart myocardial <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images, enabling a <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3-min scan.

METHODS

Data were acquired using a free-running 3D golden-angle radial readout interleaved with inversion recovery and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> -preparation pulses. After correction for translational respiratory motion, non-rigid cardiac-motion-corrected reconstruction with dictionary-based low-rank compression and patch-based regularization enabled 3D whole-heart <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> mapping at any given cardiac phase as well as whole-heart cardiac cine imaging. The framework was validated and compared with established methods in 11 healthy subjects.

RESULTS

Good quality 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images were reconstructed for all subjects. Septal <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> values using the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1200</mml:mn> <mml:mo>±</mml:mo> <mml:mn>50</mml:mn></mml:mrow> <mml:annotation>$$ 1200\pm 50 $$</mml:annotation></mml:semantics> </mml:math> ms) were higher than those from a 2D MOLLI sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1063</mml:mn> <mml:mo>±</mml:mo> <mml:mn>33</mml:mn></mml:mrow> <mml:annotation>$$ 1063\pm 33 $$</mml:annotation></mml:semantics> </mml:math> ms), which is known to underestimate <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> , while <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> values from the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>4</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 4 $$</mml:annotation></mml:semantics> </mml:math> ms) were in good agreement with those from a 2D GraSE sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 2 $$</mml:annotation></mml:semantics> </mml:math> ms).

CONCLUSION

The proposed technique provides 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images with isotropic spatial resolution in a single <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3.3-min scan.

3D model-based super-resolution motion-corrected cardiac T1 mapping

Objective. To provide 3D high-resolution cardiac T1 maps using model-based super-resolution reconstruction (SRR).Approach. Due to signal-to-noise ratio limitations and the motion of the heart during imaging, often 2D T1 maps with only low through-plane resolution (i.e. slice thickness of 6-8 mm) can be obtained. Here, a model-based SRR approach is presented, which combines multiple stacks of 2D acquisitions with 6-8 mm slice thickness and generates 3D high-resolution T1 maps with a slice thickness of 1.5-2 mm. Every stack was acquired in a different breath hold (BH) and any misalignment between BH was corrected retrospectively. The novelty of the proposed approach is the BH correction and the application of model-based SRR on cardiac T1 Mapping. The proposed approach was evaluated in numerical simulations and phantom experiments and demonstrated in four healthy subjects.Main results. Alignment of BH states was essential for SRR even in healthy volunteers. In simulations, respiratory motion could be estimated with an RMS error of 0.18 ± 0.28 mm. SRR improved the visualization of small structures. High accuracy and precision (average standard deviation of 69.62 ms) of the T1 values was ensured by SRR while the detectability of small structures increased by 40%.Significance. The proposed SRR approach provided T1 maps with high in-plane and high through-plane resolution (1.3 × 1.3 × 1.5-2 mm³). The approach led to improvements in the visualization of small structures and precise T1 values.

Estimation of Cardiac Short Axis Slice Levels with a Cascaded Deep Convolutional and Recurrent Neural Network Model

Automatic identification of short axis slice levels in cardiac magnetic resonance imaging (MRI) is important in efficient and precise diagnosis of cardiac disease based on the geometry of the left ventricle. We developed a combined model of convolutional neural network (CNN) and recurrent neural network (RNN) that takes a series of short axis slices as input and predicts a series of slice levels as output. Each slice image was labeled as one of the following five classes: out-of-apical, apical, mid, basal, and out-of-basal levels. A variety of multi-class classification models were evaluated. When compared with the CNN-alone models, the cascaded CNN-RNN models resulted in higher mean F1-score and accuracy. In our implementation and testing of four different baseline networks with different combinations of RNN modules, MobileNet as the feature extractor cascaded with a two-layer long short-term memory (LSTM) network produced the highest scores in four of the seven evaluation metrics, i.e., five F1-scores, area under the curve (AUC), and accuracy. Our study indicates that the cascaded CNN-RNN models are superior to the CNN-alone models for the classification of short axis slice levels in cardiac cine MR images.

Free-breathing, non-ECG, simultaneous myocardial T1 , T2 , T2 *, and fat-fraction mapping with motion-resolved cardiovascular MR multitasking

PURPOSE

To develop a free-breathing, non-electrocardiogram technique for simultaneous myocardial T1 , T2 , T2 *, and fat-fraction (FF) mapping in a single scan.

METHODS

The MR Multitasking framework is adapted to quantify T1 , T2 , T2 *, and FF simultaneously. A variable TR scheme is developed to preserve temporal resolution and imaging efficiency. The underlying high-dimensional image is modeled as a low-rank tensor, which allows accelerated acquisition and efficient reconstruction. The accuracy and/or repeatability of the technique were evaluated on static and motion phantoms, 12 healthy volunteers, and 3 patients by comparing to the reference techniques.

RESULTS

In static and motion phantoms, T1 /T2 /T2 */FF measurements showed substantial consistency (R > 0.98) and excellent agreement (intraclass correlation coefficient > 0.93) with reference measurements. In human subjects, the proposed technique yielded repeatable T1 , T2 , T2 *, and FF measurements that agreed with those from references.

CONCLUSIONS

The proposed free-breathing, non-electrocardiogram, motion-resolved Multitasking technique allows simultaneous quantification of myocardial T1 , T2 , T2 *, and FF in a single 2.5-min scan.

SAturation-recovery and Variable-flip-Angle-based three-dimensional free-breathing cardiovascular magnetic resonance T1 mapping at 3 T

The purpose of the current study was to develop and validate a three-dimensional (3D) free-breathing cardiac T₁ -mapping sequence using SAturation-recovery and Variable-flip-Angle (SAVA). SAVA sequentially acquires multiple electrocardiogram-triggered volumes using a multishot spoiled gradient-echo sequence. The first volume samples the equilibrium signal of the longitudinal magnetization, where a flip angle of 2° is used to reduce the time for the magnetization to return to equilibrium. The succeeding three volumes are saturation prepared with variable delays, and are acquired using a 15° flip angle to maintain the signal-to-noise ratio. A diaphragmatic navigator is used to compensate the respiratory motion. T₁ is calculated using a saturation-recovery model that accounts for the flip angle. We validated SAVA by simulations, phantom, and human subject experiments at 3 T. SAVA was compared with modified Look-Locker inversion recovery (MOLLI) and saturation-recovery single-shot acquisition (SASHA) in vivo. In phantoms, T₁ by SAVA had good agreement with the reference (R²  = 0.99). In vivo 3D T₁ mapping by SAVA could achieve an imaging resolution of 1.25 × 1.25 × 8 mm³ . Both global and septal T₁ values by SAVA (1347 ± 37 and 1332 ± 42 ms) were in between those by SASHA (1612 ± 63 and 1618 ± 51 ms) and MOLLI (1143 ± 59 and 1188 ± 65 ms). According to the standard deviation (SD) and coefficient of variation (CV), T₁ precision measured by SAVA (SD: 99 ± 14 and 60 ± 8 ms; CV: 7.4% ± 0.9% and 4.5% ± 0.6%) was comparable with MOLLI (SD: 99 ± 25 and 46 ± 12 ms; CV: 8.8% ± 2.5% and 3.9% ± 1.1%) and superior to SASHA (SD: 222 ± 89 and 132 ± 33 ms; CV: 13.8% ± 5.5% and 8.1% ± 2.0%). It was concluded that the proposed free-breathing SAVA sequence enables more efficient 3D whole-heart T₁ estimation with good accuracy and precision.

Free-breathing self-gated continuous-IR spiral T1 mapping: Comparison of dual flip-angle and Bloch-Siegert B1-corrected techniques

Purpose

To develop a B1‐corrrected single flip‐angle continuous acquisition strategy with free‐breathing and cardiac self‐gating for spiral T1 mapping, and compare it to a previous dual flip‐angle technique.

Methods

Data were continuously acquired using a spiral‐out trajectory, rotated by the golden angle in time. During the first 2 s, off‐resonance Fermi RF pulses were applied to generate a Bloch‐Siegert shift B1 map, and the subsequent data were acquired with an inversion RF pulse applied every 4 s to create a T1* map. The final T1 map was generated from the B1 and the T1* maps by using a look‐up table that accounted for slice profile effects, yielding more accurate T1 values. T1 values were compared to those from inversion recovery (IR) spin echo (phantom only), MOLLI, SAturation‐recovery single‐SHot Acquisition (SASHA), and previously proposed dual flip‐angle results. This strategy was evaluated in a phantom and 25 human subjects.

Results

The proposed technique showed good agreement with IR spin‐echo results in the phantom experiment. For in‐vivo studies, the proposed technique and the previously proposed dual flip‐angle method were more similar to SASHA results than to MOLLI results.

Conclusions

B1‐corrected single flip‐angle T1 mapping successfully acquired B1 and T1 maps in a free‐breathing, continuous‐IR spiral acquisition, providing a method with improved accuracy to measure T1 using a continuous Look‐Locker acquisition, as compared to the previously proposed dual excitation flip‐angle technique.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Myocardial T1, T2, T2*, and fat fraction quantification via low-rank motion-corrected cardiac MR fingerprinting

Purpose

Develop a novel 2D cardiac MR fingerprinting (MRF) approach to enable simultaneous T1, T2, T2*, and fat fraction (FF) myocardial tissue characterization in a single breath‐hold scan.

Methods

Simultaneous, co‐registered, multi‐parametric mapping of T1, T2, and FF has been recently achieved with cardiac MRF. Here, we further incorporate T2* quantification within this approach, enabling simultaneous T1, T2, T2*, and FF myocardial tissue characterization in a single breath‐hold scan. T2* quantification is achieved with an eight‐echo readout that requires a long cardiac acquisition window. A novel low‐rank motion‐corrected (LRMC) reconstruction is exploited to correct for cardiac motion within the long acquisition window. The proposed T1/T2/T2*/FF cardiac MRF was evaluated in phantom and in 10 healthy subjects in comparison to conventional mapping techniques.

Results

The proposed approach achieved high quality parametric mapping of T1, T2, T2*, and FF with corresponding normalized RMS error (RMSE) T1 = 5.9%, T2 = 9.6% (T2 values <100 ms), T2* = 3.3% (T2* values <100 ms), and FF = 0.8% observed in phantom scans. In vivo, the proposed approach produced higher left‐ventricular myocardial T1 values than MOLLI (1148 vs 1056 ms), lower T2 values than T2‐GraSE (42.8 vs 50.6 ms), lower T2* values than eight‐echo gradient echo (GRE) (35.0 vs 39.4 ms), and higher FF values than six‐echo GRE (0.8 vs 0.3 %) reference techniques. The proposed approach achieved considerable reduction in motion artifacts compared to cardiac MRF without motion correction, improved spatial uniformity, and statistically higher apparent precision relative to conventional mapping for all parameters.

Conclusion

The proposed cardiac MRF approach enables simultaneous, co‐registered mapping of T1, T2, T2*, and FF in a single breath‐hold for comprehensive myocardial tissue characterization, achieving higher apparent precision than conventional methods.

Dual flip-angle IR-FLASH with spin history mapping for B1+ corrected T1 mapping: Application to T1 cardiovascular magnetic resonance multitasking

PURPOSE

To develop a single-scan method for B 1 + -corrected T1 mapping and apply it for free-breathing (FB) cardiac MR multitasking without electrocardiogram (ECG) triggering.

METHODS

One dual flip-angle (2FA) inversion recovery (IR)-FLASH scan provides two observations of T 1 ∗ (apparent T1 ) corresponding to two distinct combinations of the nominal FA α and B 1 + . Spatiotemporally coregistered T1 and B 1 + spin history maps are obtained by fitting the 2FA signal model. T1 estimate accuracy and repeatability for single flip-angle (1FA) and 2FA IR-FLASH sequence MR multitasking were evaluated at 3T. A T1 phantom was first imaged on the scanner table, then on two human subjects' thoraxes in both breath-hold (BH) and FB conditions. IR-turbo spin echo (IR-TSE) static phantom T1 measurements served as reference. In 10 healthy subjects, myocardial T1 was evaluated with ECG-free, FB multitasking sequences alongside ECG-triggered BH MOLLI.

RESULTS

For phantom-on-table T1 estimates, 2FA agreed better with IR-TSE (intraclass correlation coefficient [ICC] = 0.996, mean error ± SD = -1.6% ± 1.9%) than did 1FA (ICC = 0.922; mean error ± SD = -4.3% ± 12%). For phantom-on-thorax, 2FA was more repeatable and robust to respiration than 1FA (coefficient of variation [CoV] = 1.2% 2FA, = 11.3% 1FA). In vivo, in intrasession T1 repeatability, 2FA (septal CoV = 2.4%, six-segment CoV = 4.4%) outperformed 1FA (septal CoV = 3.1%, six-segment CoV = 5.5%). In six-segment T1 homogeneity, 2FA (CoV = 7.9%) also outperformed 1FA (CoV = 11.1%).

CONCLUSION

The 2FA IR-FLASH improves T1 estimate accuracy and repeatability over 1FA IR-FLASH, and enables single-scan B 1 + -corrected T1 mapping without BHs or ECG when used with MR multitasking.

Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group

On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T₁ . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.

Free-breathing 3D cardiac T1 mapping with transmit B1 correction at 3T

PURPOSE

To develop a cardiac T₁ mapping method for free-breathing 3D T₁ mapping of the whole heart at 3 T with transmit B₁ ( B 1 + ) correction.

METHODS

A free-breathing, electrocardiogram-gated inversion-recovery sequence with spoiled gradient-echo readout was developed and optimized for cardiac T₁ mapping at 3 T. High-frame-rate dynamic images were reconstructed from sparse (k,t)-space data acquired along a stack-of-stars trajectory using a subspace-based method for accelerated imaging. Joint T₁ and flip-angle estimation was performed in T₁ mapping to improve its robustness to B 1 + inhomogeneity. Subject-specific timing of data acquisition was used in the estimation to account for natural heart-rate variations during the imaging experiment.

RESULTS

Simulations showed that accuracy and precision of T₁ mapping can be improved with joint T₁ and flip-angle estimation and optimized electrocardiogram-gated spoiled gradient echo-based inversion-recovery acquisition scheme. The phantom study showed good agreement between the T₁ maps from the proposed method and the reference method. Three-dimensional cardiac T₁ maps (40 slices) were obtained at a 1.9-mm in-plane and 4.5-mm through-plane spatial resolution from healthy subjects (n = 6) with an average imaging time of 14.2 ± 1.6 minutes (heartbeat rate: 64.2 ± 7.1 bpm), showing myocardial T₁ values comparable to those obtained from modified Look-Locker inversion recovery. The proposed method generated B 1 + maps with spatially smooth variation showing 21%-32% and 11%-15% variations across the septal-lateral and inferior-anterior regions of the myocardium in the left ventricle.

CONCLUSION

The proposed method allows free-breathing 3D T₁ mapping of the whole heart with transmit B₁ correction in a practical imaging time.

Generalized low-rank nonrigid motion-corrected reconstruction for MR fingerprinting

PURPOSE

Develop a novel low-rank motion-corrected (LRMC) reconstruction for nonrigid motion-corrected MR fingerprinting (MRF).

METHODS

Generalized motion-corrected (MC) reconstructions have been developed for steady-state imaging. Here we extend this framework to enable nonrigid MC for transient imaging applications with varying contrast, such as MRF. This is achieved by integrating low-rank dictionary-based compression into the generalized MC model to reconstruct MC singular images, reducing motion artifacts in the resulting parametric maps. The proposed LRMC reconstruction was applied for cardiac motion correction in 2D myocardial MRF (T₁ and T₂ ) with extended cardiac acquisition window (~450 ms) and for respiratory MC in free-breathing 3D myocardial and 3D liver MRF. Experiments were performed in phantom and 22 healthy subjects. The proposed approach was compared with reference spin echo (phantom) and with 2D electrocardiogram-triggered/breath-hold MOLLI and T₂ gradient-and-spin echo conventional maps (in vivo 2D and 3D myocardial MRF).

RESULTS

Phantom results were in general agreement with reference spin-echo measurements, presenting relative errors of approximately 5.4% and 5.5% for T₁ and short T₂ (<100 ms), respectively. The proposed LRMC MRF reduced residual blurring artifacts with respect to no MC for cardiac or respiratory motion in all cases (2D and 3D myocardial, 3D abdominal). In 2D myocardial MRF, left-ventricle T₁ values were 1150 ± 41 ms for LRMC MRF and 1010 ± 56 ms for MOLLI; T₂ values were 43.8 ± 2.3 ms for LRMC MRF and 49.5 ± 4.5 ms for T₂ gradient and spin echo. Corresponding measurements for 3D myocardial MRF were 1085 ± 30 ms and 1062 ± 29 ms for T₁ , and 43.5 ± 1.9 ms and 51.7 ± 1.7 ms for T₂ . For 3D liver, LRMC MRF measured liver T₁ at 565 ± 44 ms and liver T₂ at 35.4 ± 2.4 ms.

CONCLUSION

The proposed LRMC reconstruction enabled generalized (nonrigid) MC for 2D and 3D MRF, both for cardiac and respiratory motion. The proposed approach reduced motion artifacts in the MRF maps with respect to no motion compensation and achieved good agreement with reference measurements.

Deep Learning-Based ECG-Free Cardiac Navigation for Multi-Dimensional and Motion-Resolved Continuous Magnetic Resonance Imaging

For the clinical assessment of cardiac vitality, time-continuous tomographic imaging of the heart is used. To further detect e.g., pathological tissue, multiple imaging contrasts enable a thorough diagnosis using magnetic resonance imaging (MRI). For this purpose, time-continous and multi-contrast imaging protocols were proposed. The acquired signals are binned using navigation approaches for a motion-resolved reconstruction. Mostly, external sensors such as electrocardiograms (ECG) are used for navigation, leading to additional workflow efforts. Recent sensor-free approaches are based on pipelines requiring prior knowledge, e.g., typical heart rates. We present a sensor-free, deep learning-based navigation that diminishes the need for manual feature engineering or the necessity of prior knowledge compared to previous works. A classifier is trained to estimate the R-wave timepoints in the scan directly from the imaging data. Our approach is evaluated on 3-D protocols for continuous cardiac MRI, acquired in-vivo and free-breathing with single or multiple imaging contrasts. We achieve an accuracy of > 98% on previously unseen subjects, and a well comparable image quality with the state-of-the-art ECG-based reconstruction. Our method enables an ECG-free workflow for continuous cardiac scans with simultaneous anatomic and functional imaging with multiple contrasts. It can be potentially integrated without adapting the sampling scheme to other continuous sequences by using the imaging data for navigation and reconstruction.

Magnetization-prepared GRASP MRI for rapid 3D T1 mapping and fat/water-separated T1 mapping

PURPOSE

This study aimed to (i) develop Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI using a stack-of-stars trajectory for rapid free-breathing T1 mapping and (ii) extend MP-GRASP to multi-echo acquisition (MP-Dixon-GRASP) for fat/water-separated (water-specific) T1 mapping.

METHODS

An adiabatic non-selective 180° inversion-recovery pulse was added to a gradient-echo-based golden-angle stack-of-stars sequence for magnetization-prepared 3D single-echo or 3D multi-echo acquisition. In combination with subspace-based GRASP-Pro reconstruction, the sequence allows for standard T1 mapping (MP-GRASP) or fat/water-separated T1 mapping (MP-Dixon-GRASP), respectively. The accuracy of T1 mapping using MP-GRASP was evaluated in a phantom and volunteers (brain and liver) against clinically accepted reference methods. The repeatability of T1 estimation was also assessed in the phantom and volunteers. The performance of MP-Dixon-GRASP for water-specific T1 mapping was evaluated in a fat/water phantom and volunteers (brain and liver).

RESULTS

ROI-based mean T1 values are correlated between the references and MP-GRASP in the phantom (R2 = 1.0), brain (R2 = 0.96), and liver (R2 = 0.73). MP-GRASP achieved good repeatability of T1 estimation in the phantom (R2 = 1.0), brain (R2 = 0.99), and liver (R2 = 0.82). Water-specific T1 is different from in-phase and out-of-phase composite T1 (composite T1 when fat and water signal are mixed in phase or out of phase) both in the phantom and volunteers.

CONCLUSION

This work demonstrated the initial performance of MP-GRASP and MP-Dixon-GRASP MRI for rapid 3D T1 mapping and 3D fat/water-separated T1 mapping in the brain (without motion) and in the liver (during free breathing). With fat/water-separated T1 estimation, MP-Dixon-GRASP could be potentially useful for imaging patients with fatty-liver diseases.

Dual-excitation flip-angle simultaneous cine and T1 mapping using spiral acquisition with respiratory and cardiac self-gating

PURPOSE

To develop a free-breathing cardiac self-gated technique that provides cine images and B1+ slice profile-corrected T₁ maps from a single acquisition.

METHODS

Without breath-holding or electrocardiogram gating, data were acquired continuously on a 3T scanner using a golden-angle gradient-echo spiral pulse sequence, with an inversion RF pulse applied every 4 seconds. Flip angles of 3° and 15° were used for readouts after the first four and second four inversions. Self-gating cardiac triggers were extracted from heart image navigators, and respiratory motion was corrected by rigid registration on each heartbeat. Cine images were reconstructed from the steady-state portion of 15° readouts using a low-rank plus sparse reconstruction. The T₁ maps were fit using a projection onto convex sets approach from images reconstructed using slice profile-corrected dictionary learning. This strategy was evaluated in a phantom and 14 human subjects.

RESULTS

The self-gated signal demonstrated close agreement with the acquired electrocardiogram signal. The image quality for the proposed cine images and standard clinical balanced SSFP images were 4.31 (±0.50) and 4.65 (±0.30), respectively. The slice profile-corrected T₁ values were similar to those of the inversion-recovery spin-echo technique in phantom, and had a higher global T₁ value than that of MOLLI in human subjects.

CONCLUSIONS

Cine and T₁ mapping using spiral acquisition with respiratory and cardiac self-gating successfully acquired T₁ maps and cine images in a single acquisition without the need for electrocardiogram gating or breath-holding. This dual-excitation flip-angle approach provides a novel approach for measuring T₁ while accounting for B1+ and slice profile effect on the apparent T1∗ .

Synergistic multi-contrast cardiac magnetic resonance image reconstruction

Cardiac magnetic resonance imaging (CMR) is an important tool for the non-invasive diagnosis of a variety of cardiovascular diseases. Parametric mapping with multi-contrast CMR is able to quantify tissue alterations in myocardial disease and promises to improve patient care. However, magnetic resonance imaging is an inherently slow imaging modality, resulting in long acquisition times for parametric mapping which acquires a series of cardiac images with different contrasts for signal fitting or dictionary matching. Furthermore, extra efforts to deal with respiratory and cardiac motion by triggering and gating further increase the scan time. Several techniques have been developed to speed up CMR acquisitions, which usually acquire less data than that required by the Nyquist-Shannon sampling theorem, followed by regularized reconstruction to mitigate undersampling artefacts. Recent advances in CMR parametric mapping speed up CMR by synergistically exploiting spatial-temporal and contrast redundancies. In this article, we will review the recent developments in multi-contrast CMR image reconstruction for parametric mapping with special focus on low-rank and model-based reconstructions. Deep learning-based multi-contrast reconstruction has recently been proposed in other magnetic resonance applications. These developments will be covered to introduce the general methodology. Current technical limitations and potential future directions are discussed. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Deep learning-enhanced T1 mapping with spatial-temporal and physical constraint

PURPOSE

To propose a reconstruction framework to generate accurate T₁ maps for a fast MR T₁ mapping sequence.

METHODS

A deep learning-enhanced T₁ mapping method with spatial-temporal and physical constraint (DAINTY) was proposed. This method explicitly imposed low-rank and sparsity constraints on the multiframe T₁ -weighted images to exploit the spatial-temporal correlation. A deep neural network was used to efficiently perform T₁ mapping as well as denoise and reduce undersampling artifacts. Additionally, the physical constraint was used to build a bridge between low-rank and sparsity constraint and deep learning prior, so the benefits of constrained reconstruction and deep learning can be both available. The DAINTY method was trained on simulated brain data sets, but tested on real acquired phantom, 6 healthy volunteers, and 7 atherosclerosis patients, compared with the narrow-band k-space-weighted image contrast filter conjugate-gradient SENSE (NK-CS) method, kt-sparse-SENSE (kt-SS) method, and low-rank plus sparsity (L+S) method with least-squares T₁ fitting and direct deep learning mapping.

RESULTS

The DAINTY method can generate more accurate T₁ maps and higher-quality T₁ -weighted images compared with other methods. For atherosclerosis patients, the intraplaque hemorrhage can be successfully detected. The computation speed of DAINTY was 10 times faster than traditional methods. Meanwhile, DAINTY can reconstruct images with comparable quality using only 50% of k-space data.

CONCLUSION

The proposed method can provide accurate T₁ maps and good-quality T₁ -weighted images with high efficiency.

Cardiac magnetic resonance fingerprinting: Trends in technical development and potential clinical applications

Quantitative cardiac magnetic resonance has emerged in recent years as an approach for evaluating a range of cardiovascular conditions, with T1 and T2 mapping at the forefront of these developments. Cardiac Magnetic Resonance Fingerprinting (cMRF) provides a rapid and robust framework for simultaneous quantification of myocardial T1 and T2 in addition to other tissue properties. Since the advent of cMRF, a number of technical developments and clinical validation studies have been reported. This review provides an overview of cMRF, recent technical developments, healthy subject and patient studies, anticipated technical improvements, and potential clinical applications. Recent technical developments include slice profile and pulse efficiency corrections, improvements in image reconstruction, simultaneous multislice imaging, 3D whole-ventricle imaging, motion-resolved imaging, fat-water separation, and machine learning for rapid dictionary generation. Future technical developments in cMRF, such as B0 and B1 field mapping, acceleration of acquisition and reconstruction, imaging of patients with implanted devices, and quantification of additional tissue properties are also described. Potential clinical applications include characterization of infiltrative, inflammatory, and ischemic cardiomyopathies, tissue characterization in the left atrium and right ventricle, post-cardiac transplantation assessment, reduction of contrast material, pre-procedural planning for electrophysiology interventions, and imaging of patients with implanted devices.

Cardiac MRI-Update 2020

PURPOSE

To review emerging techniques in cardiac magnetic resonance imaging (CMR) and their clinical applications with a special emphasis on new technologies, recent trials, and updated guidelines.

TECHNOLOGICAL INNOVATIONS

The utility of CMR has expanded with the development of new MR sequences, postprocessing techniques, and artificial intelligence-based technologies, which have substantially increased the spectrum, quality, and reliability of information that can be obtained by CMR.

ESTABLISHED AND EMERGING INDICATIONS

The CMR modality has become an irreplaceable tool for diagnosis, treatment guidance and follow-up of patients with ischemic heart disease, myocarditis, and cardiomyopathies. Its role has been further strengthened by recent trials and guidelines. Quantitative mapping techniques are increasingly used for tissue characterization and detection of diffuse myocardial changes including myocardial storage diseases.

PRACTICAL RECOMMENDATIONS

With state-of-the-art CMR sequences, postprocessing techniques and understanding of their interpretation, CMR makes invaluable contributions to provide state-of-the-art diagnostics and care for cardiac patients in a multidisciplinary team.

Fully self-gated free-running 3D Cartesian cardiac CINE with isotropic whole-heart coverage in less than 2 min

PURPOSE

To develop a novel fast water-selective free-breathing 3D Cartesian cardiac CINE scan with full self-navigation and isotropic whole-heart (WH) coverage.

METHODS

A free-breathing 3D Cartesian cardiac CINE scan with a water-selective balanced steady-state free precession and a continuous (non-ECG-gated) variable-density Cartesian sampling with spiral profile ordering, out-inward sampling and acquisition-adaptive alternating tiny golden and golden angle increment between spiral arms is proposed. Data is retrospectively binned based on respiratory and cardiac self-navigation signals. A translational respiratory-motion-corrected and cardiac-motion-resolved image is reconstructed with a multi-bin patch-based low-rank reconstruction (MB-PROST) within about 15 min. A respiratory-motion-resolved approach is also investigated. The proposed 3D Cartesian cardiac CINE is acquired in sagittal orientation in 1 min 50 s for 1.9 mm³ isotropic WH coverage. Left ventricular (LV) function parameters and image quality derived from a blinded reading of the proposed 3D CINE framework are compared against conventional multi-slice 2D CINE imaging in 10 healthy subjects and 10 patients with suspected cardiovascular disease.

RESULTS

The proposed framework provides free-breathing 3D cardiac CINE images with 1.9 mm³ spatial and about 45 ms temporal resolution in a short acquisition time (<2 min). LV function parameters derived from 3D CINE were in good agreement with 2D CINE (10 healthy subjects and 10 patients). Bias and confidence intervals were obtained for end-systolic volume, end-diastolic volume and ejection fraction of 0.1 ± 3.5 mL, -0.6 ± 8.2 mL and -0.1 ± 2.2%, respectively.

CONCLUSION

The proposed framework enables isotropic 3D Cartesian cardiac CINE under free breathing for fast assessment of cardiac anatomy and function.

3D free-breathing cardiac magnetic resonance fingerprinting

PURPOSE

To develop a novel respiratory motion compensated three-dimensional (3D) cardiac magnetic resonance fingerprinting (cMRF) approach for whole-heart myocardial T₁ and T₂ mapping from a free-breathing scan.

METHODS

Two-dimensional (2D) cMRF has been recently proposed for simultaneous, co-registered T₁ and T₂ mapping from a breath-hold scan; however, coverage is limited. Here we propose a novel respiratory motion compensated 3D cMRF approach for whole-heart myocardial T₁ and T₂ tissue characterization from a free-breathing scan. Variable inversion recovery and T₂ preparation modules are used for parametric encoding, respiratory bellows driven localized autofocus is proposed for beat-to-beat translation motion correction and a subspace regularized reconstruction is employed to accelerate the scan. The proposed 3D cMRF approach was evaluated in a standardized T₁ /T₂ phantom in comparison with reference spin echo values and in 10 healthy subjects in comparison with standard 2D MOLLI, SASHA and T₂ -GraSE mapping techniques at 1.5 T.

RESULTS

3D cMRF T₁ and T₂ measurements were generally in good agreement with reference spin echo values in the phantom experiments, with relative errors of 2.9% and 3.8% for T₁ and T₂ (T₂ < 100 ms), respectively. in vivo left ventricle (LV) myocardial T₁ values were 1054 ± 19 ms for MOLLI, 1146 ± 20 ms for SASHA and 1093 ± 24 ms for the proposed 3D cMRF; corresponding T₂ values were 51.8 ± 1.6 ms for T2-GraSE and 44.6 ± 2.0 ms for 3D cMRF. LV coefficients of variation were 7.6 ± 1.6% for MOLLI, 12.1 ± 2.7% for SASHA and 5.8 ± 0.8% for 3D cMRF T₁ , and 10.5 ± 1.4% for T2-GraSE and 11.7 ± 1.6% for 3D cMRF T₂ .

CONCLUSION

The proposed 3D cMRF can provide whole-heart, simultaneous and co-registered T₁ and T₂ maps with accuracy and precision comparable to those of clinical standards in a single free-breathing scan of about 7 min.

Prospective correction of patient-specific respiratory motion in myocardial T1 and T2 mapping

PURPOSE

Respiratory motion in cardiovascular MRI presents a challenging problem with many potential solutions. Current approaches require breath-holds, apply retrospective image registration, or significantly increase scan time by respiratory gating. Myocardial T₁ and T₂ mapping techniques are particularly sensitive to motion as they require multiple source images to be accurately aligned prior to the estimation of tissue relaxation. We propose a patient-specific prospective motion correction (PROCO) strategy that corrects respiratory motion on the fly with the goal of reducing the spatial variation of myocardial parametric mapping techniques.

METHODS

A rapid, patient-specific training scan was performed to characterize respiration-induced motion of the heart relative to a diaphragmatic navigator, and a parametric mapping pulse sequence utilized the resulting motion model to prospectively update the scan plane in real-time. Midventricular short-axis T₁ and T₂ maps were acquired under breath-hold or free-breathing conditions with and without PROCO in 7 healthy volunteers and 3 patients. T₁ and T₂ were measured in 6 segments and compared to reference standard breath-hold measurements using Bland-Altman analysis.

RESULTS

PROCO significantly reduced the spatial variation of parametric maps acquired during free-breathing, producing limits of agreement of -47.16 to 30.98 ms (T₁ ) and -1.35 to 4.02 ms (T₂ ), compared to -67.77 to 74.34 ms (T₁ ) and -2.21 to 5.62 ms (T₂ ) for free-breathing acquisition without PROCO.

CONCLUSION

Patient-specific respiratory PROCO method significantly reduced the spatial variation of myocardial T₁ and T₂ mapping, while allowing for 100% efficient free-breathing acquisitions.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

Isotropic 3D Cartesian single breath-hold CINE MRI with multi-bin patch-based low-rank reconstruction

PURPOSE

To develop a novel acquisition and reconstruction framework for isotropic 3D Cartesian cardiac CINE within a single breath-hold for left ventricle (LV) and whole-heart coverage.

METHODS

A variable-density Cartesian acquisition with spiral profile ordering, out-inward sampling, and acquisition-adaptive alternating tiny golden/golden angle increment between spiral arms is proposed to provide incoherent and nonredundant sampling within and among cardiac phases. A novel multi-bin patch-based low-rank reconstruction, named MB-PROST, is proposed to exploit redundant information on a local (within a patch), nonlocal (similar patches within a spatial neighborhood), and temporal (among all cardiac phases) scale with an implicit motion alignment among patches. The proposed multi-bin patch-based low-rank reconstruction reconstruction is compared against compressed sensing reconstruction, whereas LV function parameters derived from the proposed 3D CINE framework are compared against those estimated from conventional multislice 2D CINE imaging in 10 healthy subjects and 15 patients.

RESULTS

The proposed framework provides 3D cardiac CINE images with high spatial (1.9 mm³ ) and temporal resolution (˜50 ms) in a single breath-hold of ˜20 s for LV and ˜26 s for whole-heart coverage in healthy subjects. Shorter breath-hold durations of ˜13 to 15 s are feasible for LV coverage with slightly anisotropic resolution (1.9 × 1.9 × 2.5 mm) in patients. LV function parameters derived from 3D CINE were in good agreement with 2D CINE, with a bias of -0.1 mL/0.1 mL, -0.9 mL/-1.0 mL, -0.1%/-0.8%; and confidence intervals of ±1.7 mL/±3.7 mL, ±1.2 mL/±2.6 mL, and ±1.2%/±3.6% (10 healthy subjects/15 patients) for end-systolic volume, end-diastolic volume, and ejection fraction, respectively.

CONCLUSION

The proposed framework enables 3D isotropic cardiac CINE in a single breath-hold scan of ˜20 s/˜26 s for LV/whole-heart coverage, showing good agreement with clinical 2D CINE scans in terms of LV functional assessment.

Free-running 3D whole heart myocardial T1 mapping with isotropic spatial resolution

PURPOSE

To develop a free-running (free-breathing, retrospective cardiac gating) 3D myocardial T1 mapping with isotropic spatial resolution.

METHODS

The free-running sequence is inversion recovery (IR)-prepared followed by continuous 3D golden angle radial data acquisition. 1D respiratory motion signal is extracted from the k-space center of all spokes and used to bin the k-space data into different respiratory states, enabling estimation and correction of 3D translational respiratory motion, whereas cardiac motion is recorded using electrocardiography and synchronized with data acquisition. 3D translational respiratory motion compensated T1 maps at diastole and systole were generated with 1.5 mm isotropic spatial resolution with low-rank inversion and high-dimensionality patch-based undersampled reconstruction. The technique was validated against conventional methods in phantom and 9 healthy subjects.

RESULTS

Phantom results demonstrated good agreement (R2 = 0.99) of T1 estimation with reference method. Homogeneous systolic and diastolic 3D T1 maps were reconstructed from the proposed technique. Diastolic septal T1 estimated with the proposed method (1140 ± 36 ms) was comparable to the saturation recovery single-shot acquisition (SASHA) sequence (1153 ± 49 ms), but was higher than the modified Look-Locker inversion recovery (MOLLI) sequence (1037 ± 33 ms). Precision of the proposed method (42 ± 8 ms) was comparable to MOLLI (41 ± 7 ms) and improved with respect to SASHA (87 ± 19 ms).

CONCLUSIONS

The proposed free-running whole heart T1 mapping method allows for reconstruction of isotropic resolution 3D T1 maps at different cardiac phases, serving as a promising tool for whole heart myocardial tissue characterization.