1955

Self-gated 5-minute whole-heart 4D flow imaging
Aaron Pruitt1, Adam Rich1, Yingmin Liu2, Ning Jin3, Lee Potter4, Orlando Simonetti1,2,5, and Rizwan Ahmad1,2,4

1Biomedical Engineering, The Ohio State University, Columbus, OH, United States, 2Davis Heart and Lung Research Institute, The Ohio State University, Columbus, OH, United States, 3Siemens Medical Solutions, USA, Columbus, OH, United States, 4Electrical and Computer Engineering, The Ohio State University, Columbus, OH, United States, 5Internal Medicine, The Ohio State University, Columbus, OH, United States

Synopsis

4D flow imaging can provide comprehensive hemodynamical analysis of blood flow through the heart and great vessels; however, acquiring 4D flow images with whole-heart coverage is prohibitively time-consuming. In this work we describe a highly accelerated and fully self-gated whole-heart 4D flow acquisition and reconstruction methodology. Additionally, we show proof-of-concept of a fully self-gated 5-minute whole-heart 4D flow exam and demonstrate good agreement in aortic flow quantification compared to conventional 2D phase contrast.

Introduction

Volumetric and temporally-resolved 3-directional velocity (4D flow) MRI provides comprehensive hemodynamical analysis of blood flow through the chambers and great vessels of the heart.1 Although whole-heart coverage can be particularly valuable for the evaluation of congenital defects characterized by pathological flow,2 the extensive acquisition time required precludes its use in common clinical practice, leading to recent efforts to accelerate 4D flow acquisition.3 Furthermore, self-gating surrogate signals have been used for both respiratory and cardiac motion compensation4 in accelerated whole-heart 4D flow, obviating the need for ECG or time-consuming respiratory navigator-based approaches to further reduce the overhead of setup time in 4D flow imaging. In this work, we extend these concepts to propose a highly accelerated and fully self-gated 4D flow acquisition and reconstruction methodology to explore the feasibility of a 5-minute whole-heart 4D flow exam in a healthy subject.

Methods

One healthy subject was recruited for this study in accordance with the local ethics board. Highly accelerated 4D flow images were acquired with a 1.5T clinical scanner (MAGNETOM Avanto, Siemens Healthcare, Germany) for 5-minutes. Incoherent k-space sampling was achieved using a variable density, cartesian sampling pattern in the ky-kz plane (Figure 1D). Starting from the first sample at the center of the k-space, the location of each subsequent sample is advanced angularly by the golden angle and radially by another irrational number, e.g., $$$\sqrt[3](35)$$$ . The mutual irrationality of the golden angle and $$$\sqrt[3](35)$$$ leads to a k-space coverage that is approximately uniform over any arbitrary, non-contiguous time duration. A sagittal slab covering the heart and aorta was prescribed with parameters conforming to the guidelines outlined in the consensus statement (Table 1A).1 Flow-compensated readout lines through the center of k-space were interleaved every 9 TR’s to obtain a self-gating signal (SG in Figure 1A); bandpass filtering followed by principal-component-analysis yielded self-gating surrogate signals for respiratory (Figure 1B) and cardiac (Figure 1C) motion. k-Space data were binned into 3 respiratory and 20 cardiac phases. Following binning, 4D flow images were reconstructed using the previously proposed ReVEAL4D algorithm5 from the respiratory expiration bin. Note that the 5-minute acquisition time included both the time to acquire the self-gating lines, which were used solely for binning, and data from respiratory phases, which were not included in the final images. Unoptimized reconstruction time for the respiratory expiration bin was approximately 6 hours (12-core Intel i7 CPU). A summary of the reconstruction pipeline is given in Figure 2. Volumetric flow rate as a function of cardiac phase, stroke volume, and peak velocity were computed in the ascending aorta. Standard breath-held 2D-phase contrast cines transecting the ascending aorta were acquired for reference Parameters for 2D-phase contrast were matched as closely as possible to the whole-heart 4D flow acquisition.

Results

Representative magnitude and phase images in the sagittal, coronal, and transverse orientations at systole are shown in Figure 3 for the reconstructed self-gated 5-minute 4D flow images. The net acceleration rate for these images is approximately 20.6. Volumetric flow rate curves as a function of cardiac phase are given in Figure 4. The self-gated 5-minute whole-heart 4D flow acquisition yields good agreement with 2D phase contrast in the ascending aorta with mean flow rates of 63.6 mL/s and 64.3 mL/s, respectively. Likewise, stroke volume and peak velocity calculated from the self-gated 5-minute whole-heart 4D flow images agree well with values derived from 2D phase contrast, with 76.4 mL and 79.3 mL, respectively for stroke volume, and 94.3 cm/s and 88.8 cm/s, respectively, for peak velocity (Table 1B).

Conclusions

In this work, we have demonstrated the feasibility of acquiring diagnostic quality whole-heart 4D flow images in 5-minutes relying only on self-gating surrogate signals for motion compensation without the need of ECG or respiratory navigators. While our proposed whole-heart 4D flow acquisition and reconstruction strategy produced good agreement with conventional 2D-phase contrast on several clinical metrics, further validation and optimization is required in a larger cohort. In the future, we will also consider implementing a “soft-gating” strategy6 to leverage additional information from multiple respiratory states.

Acknowledgements

This work was partially funded by NIH grants R21EB021655 and R21EB022277.

References

[1] Dyverfeldt et al. 4D flow cardiovascular magnetic resonance consensus statement. 2015; JCMR 17:72.

[2] Vasanawala et al. Congenital heart disease assessment with 4D flow MRI. JMRI 2015; 42(4):870-86.

[3] Bollache et al. k-t accelerated aortic 4D flow MRI in under two minutes: feasibility and impact of resolution, k-space sampling patterns, and respiratory navigator gating on hemodynamic measurements. MRM 2018;79(1):195-207.

[4] Bastkowski et al. Self-gated golden-angle spiral 4D flow MRI. MRM 2018;80(3):904-13.

[5] Rich et al. A Bayesian model for highly accelerated phase-contrast MRI. MRM 2016;76(2):689-701.

[6] Han et al. Respiratory motion-resolved, self-gated 4D-MRI using rotating cartesian k-space (ROCK). Med Phy 2017;44(4):1359-68.

Figures

Figure 1. (A) Self-gating encoding strategy for whole-heart 4D flow acquisition. Self-gating (SG) lines through the center of k-space are interleaved every 9 TR’s. (B) Respiratory and (C) cardiac surrogate signals after filtering and principal component analysis of the self-gating signal. Cardiac trigger points are marked in red. (D) Example k-space lines sampled for 5, 10, and 20-minute acquisition for one encoding and cardiac phase at end-expiration. Although only a 5-minute acquisition is included in this study, the resulting coverage of k-space is approximately uniform regardless of acquisition time (and acceleration rate).

Figure 2. Flow chart summarizing the processing pipeline for highly accelerated, self-gated 4D flow after data acquisition.

Figure 3. Representative magnitude and phase 4D flow images from a sagittal, coronal, and transverse slice at systole from the self-gated 5-minute whole-heart acquisition.

Figure 4. Calculated volumetric flow rate in the ascending aorta as a function of cardiac phase in the healthy subject for the 5-minute whole-heart acquisition compared with the mean of three 2D phase contrast replicants; error bars represent the minimum and maximum. Here, the self-gated 5-minute whole-heart 4D flow acquisition produces comparable volumetric flow rate to 2D phase contrast.

Table 1. (A) Acquisition parameters for self-gated 5-minute whole-heart 4D flow. (B) Flow quantification metrics - stroke volume and peak velocity - calculated for 2D phase contrast and the self-gated 5-minute whole-heart 4D flow acquisition in the ascending aorta. Note that for 2D phase contrast, values are stated as mean [min, max] over three replicates performed over the course of the exam.

Proc. Intl. Soc. Mag. Reson. Med. 27 (2019)
1955