Skip to main content

Advertisement

Accelerated aortic 4D flow cardiovascular magnetic resonance using compressed sensing: applicability, validation and clinical integration

Article metrics

Abstract

Background

Three-dimensional time-resolved phase-contrast cardiovascular magnetic resonance (4D flow CMR) enables the quantification and visualisation of blood flow, but its clinical applicability remains hampered by its long scan time. The aim of this study was to evaluate the use of compressed sensing (CS) with on-line reconstruction to accelerate the acquisition and reconstruction of 4D flow CMR of the thoracic aorta.

Methods

4D flow CMR of the thoracic aorta was acquired in 20 healthy subjects using CS with acceleration factors ranging from 4 to 10. As a reference, conventional parallel imaging (SENSE) with acceleration factor 2 was used. Flow curves, net flows, peak flows and peak velocities were extracted from six contours along the aorta. To measure internal data consistency, a quantitative particle trace analysis was performed. Additionally, scan-rescan, inter- and intraobserver reproducibility were assessed. Subsequently, 4D flow CMR with CS factor 6 was acquired in 3 patients with differing aortopathies. The flow patterns resulting from particle trace visualisation were qualitatively analysed.

Results

All collected data were successfully acquired and reconstructed on-line. The average acquisition time including respiratory navigator efficiency with CS factor 6 was 5:02 ± 2:23 min while reconstruction took approximately 9 min. For CS factors of 8 or less, mean differences in net flow, peak flow and peak velocity as compared to SENSE were below 2.2 ± 7.8 ml/cycle, 4.6 ± 25.2 ml/s and − 7.9 ± 13.0 cm/s, respectively. For a CS factor of 10 differences reached 5.4 ± 8.0 ml/cycle, 14.4 ± 28.3 ml/s and − 4.0 ± 12.2 cm/s. Scan-rescan analysis yielded mean differences in net flow of − 0.7 ± 4.9 ml/cycle for SENSE and − 0.2 ± 8.5 ml/cycle for CS factor of 6.

Conclusions

A six- to eightfold acceleration of 4D flow CMR using CS is feasible. Up to a CS acceleration rate of 6, no statistically significant differences in measured flow parameters could be observed with respect to the reference technique. Acquisitions in patients with aortopathies confirm the potential to integrate the proposed method in a clinical routine setting, whereby its main benefits are scan-time savings and direct on-line reconstruction.

Background

The assessment of cardiovascular haemodynamics is important for the diagnosis and monitoring of many cardiovascular diseases including valvular heart diseases, congenital heart defects and aortopathies. In addition to echocardiography, phase-contrast (PC) cardiovascular magnetic resonance imaging (CMR) provides a tool to non-invasively quantify blood flow [1]. 4D flow CMR allows the assessment of the 3D, three-directional, time-resolved blood flow [2]. As compared to 2D PC-CMR, planning is simplified and a retrospective quantification of flow parameters in any desired imaging plane within the acquired volume makes 4D flow CMR less operator dependent [2]. However, as compared to 2D PC-CMR which is widely available and used in standardised CMR protocols [3, 4], 4D flow CMR currently remains technically challenging and is rarely used in clinical routine settings. This is mainly due to the inherently long scan time which is further increased by the need for respiratory motion compensation [5]. An aortic acquisition with parameters as described in the 4D flow CMR consensus statement [6] takes approximately 15 min by using standard parallel imaging techniques on a typical clinical CMR system. Therefore, a range of acceleration techniques that exploit spatio-temporal correlations have been applied to 4D flow CMR of the heart and/or surrounding vessels including k-t BLAST [7], k-t GRAPPA [8, 9] or k-t PCA [10, 11] allowing for acceleration factors in the order of 5 to 8. Acceleration techniques based on compressed sensing (CS) − often combined with non-Cartesian sampling, with and without exploitation of spatio-temporal correlations − have also been applied to 4D flow CMR achieving similar acceleration rates [12,13,14,15]. One major disadvantage of most proposed CS techniques is their relatively long off-line reconstruction time, reportedly being in the range of 45 to 60 min [12, 13, 16].

The aim of this study was to systematically evaluate the use of a CS technique, combining CS and parallel imaging, including on-line reconstruction to accelerate the acquisition and reconstruction of 4D flow CMR of the thoracic aorta. The impact of different CS acceleration factors on the measurement of direct haemodynamic parameters including net flow (NF), peak flow (PF) and peak velocity (PV) was assessed in healthy subjects. As a reference, conventional parallel imaging 4D flow CMR was used. A scan-rescan and inter- and intraobserver analysis were additionally performed on these parameters.

Finally, the proposed CS 4D flow CMR technique was applied to three patients with bicuspid aortic valves or Marfan syndrome to show feasibility in a clinical routine CMR examination of diseases in which aortic haemodynamic changes are relevant [17,18,19,20].

Methods

Study population

Between June 2018 and January 2019, 20 healthy subjects (27.1 years, range: 21–41; 11 women) without known cardiovascular diseases and 3 patients with aortopathies (further characterised in Table 1) were recruited. The study was approved by the local ethics committee and written informed consent was obtained from all participants.

Table 1 Patient characteristics

CMR imaging

Healthy subjects

4D flow CMR in healthy subjects was acquired on a clinical 3 T CMR system (Ingenia; Philips Healthcare, Best, The Netherlands) with a 28-channel coil array. A sagittal volume covering the thoracic aorta was acquired. No contrast agent was used. Acquisition parameters were chosen in accordance with the consensus statement [6] and are listed in Table 2. A non-symmetric four-point flow encoding scheme was used, generating one flow-compensated and three flow-encoded acquisitions [21]. All acceleration factors were defined as scan time reductions with respect to a fully sampled dataset. The number of healthy subjects was divided into two sub-groups, group 1 consisting of 15 subjects and group 2 consisting of 5 subjects.

Table 2 Acquisition parameters

In group 1, five or six 4D flow CMR acquisitions were performed per subject: four CS (acceleration factors R = 4, 6, 8 and 10 - termed CS4, CS6, CS8 and CS10) and one parallel imaging (SENSE acceleration factor R = 2 - termed S2) sequences were acquired. In addition, 9 subjects further underwent an acquisition with a SENSE acceleration factor of R = 6 (S6). To decrease the probability of systematic effects due to long examination lengths such as heart rate variations and motion, acquisitions were performed in a randomised order for each subject (except for the S6 scan, which was always acquired last).

In group 2, the scan-rescan reproducibility of the reference technique (S2) and of CS6 and CS8 was assessed using identical acquisition parameters as in group 1. The scans were also performed in a randomised order for each subject and repeated in the same order after a short break including repositioning and replanning.

Patients

Three patients were acquired on a clinical 1.5 T CMR system (Ingenia, Philips Healthcare) with a 28-channel coil array. The patients underwent a clinically indicated standard-of-care CMR examination, complemented by a CS6 4D flow measurement. As all patients underwent a contrast-enhanced angiography, the flip angle of the subsequent 4D flow CMR acquisition was increased to account for reduced longitudinal relaxation times in the blood pool (Table 2).

Coil compression [22] to 4 virtual coils was used to accelerate reconstruction times. For the S6 scan, coil compression to 12 virtual coils was used.

All acquisitions were triggered retrospectively to the electrocardiogram (ECG) and a pencil-beam respiratory navigator was placed on the right diaphragm-liver interface. The acceptance window for respiratory gating was set to 6 mm.

The CS technique used in this work combines the compressed sensing and the parallel imaging (SENSE) approaches. The CS as well as the used SENSE technique were provided by the manufacturer as product sequences (Compressed SENSE/ SENSE, Philips Healthcare). In line with most CS techniques, a variable density incoherent undersampling pattern with a more densely sampled k-space centre was used. A regularised, iterative L1 norm minimisation was then performed assuring image sparsity in the wavelet domain and data consistency [23,24,25]. Each time frame was reconstructed separately and no temporal correlations were exploited in the reconstruction. Reconstruction was CPU-parallelised on the standard reconstruction workstation of the CMR systems (32 GB RAM, Intel Xeon E5–1620 CPU).

Data analysis

Image reconstruction and concomitant field correction were performed on-line [26] while eddy current induced background phase correction [27], velocity un-aliasing and all post-processing steps were performed off-line using GTFlow (version 3.1, Gyrotools LLC, Zurich, Switzerland).

Healthy subjects

For each subject, six 2D planes were manually placed in the S2 scans. All planes intersected the aorta orthogonally to the main flow direction: one at the aortic root (Root), one at the ascending aorta (AAo), two at the aortic arch (Arch1, Arch2) and two at the descending aorta, the latter at the level of the lung-liver interface (DAo1, DAo2) (see Fig. 2). For each of these planes, contours were manually adjusted for each time point in order to compensate for aortic in-plane motion and deformation throughout the cardiac cycle. Subsequently, the contours were transferred to the CS accelerated scans. All contours were re-adjusted in the in-plane dimensions to compensate for possible bulk motion between acquisitions while their shape was kept unchanged. NF and PF were extracted from all contours whereas PV was solely extracted from the contour Root.

To assess inter- and intraobserver variability, NF and PF were additionally extracted from AAo and DAo1 in the S2 and CS6 scans from contours drawn by a second blinded observer and re-drawn by observer 1.

Finally, a quantitative particle trace analysis was performed to test for internal data consistency based on the conservation-of-mass principle. For this purpose, the percentage of particles emitted from contour AAo that reached contour DAo1 was measured. In average 398 particles were released during five time-points in early systole (time interval: 31.5 ms). Arriving particles were counted during one R-R interval.

Patients

In patients with aortopathies, a qualitative analysis of blood flow patterns was performed by particle trace visualisation and compared to previously described flow patterns in the respective diseases.

Statistical analysis

In the following, all values are given as mean ± standard deviation (SD). Bland-Altman analyses of NF, PF and PV were performed to evaluate the variability between the S2 and the CS accelerated acquisitions. Similarly, Bland-Altman analyses were performed for inter- and intraobserver variability analysis. 95% confidence intervals (CI) were calculated for mean differences (MD) and for limits of agreements (as defined by MD ± 1.96 * SD). To allow a more sensitive measure of potential temporal deviations between flow curves, the time-accumulated flow error, averaged over all contours and subjects, was calculated [10]:

$$ {E}_R=\frac{1}{n_S}\sum \limits_{S=1}^{n_S}\left(\frac{1}{n_C}\sum \limits_{C=1}^{n_C}\left(\frac{\sum \limits_{t=0}^{t_{n_T}}\left|{Q_{\mathrm{S}2}}_{S,C}(t)-{Q_{\mathrm{CS}\ \mathrm{R}}}_{S,C}(t)\right|}{\sum \limits_{t=0}^{t_{n_T}}{Q_{\mathrm{S}2}}_{S,C}(t)}\right)\right) $$

where nS, nC and nT correspond to the number of subjects, contours and time frames, respectively. QS2S, C(t) and QCS RS, C(t) denote the flow rates through the contour C of subject S at time-point t in the S2 and CS factor R measurements (R = 4, 6, 8, 10), respectively.

Results

All healthy subjects and patient data were successfully acquired and reconstructed on-line. None of the healthy subject scans needed to be aborted nor repeated. Reconstruction times were approximately 9 min for a CS6 sequence.

Healthy subjects

In 9 subjects, the examination could be extended by the additional S6 acquisition. The acquisition times including navigator efficiency, averaged over all subjects of group 1, were 15:09 (S2), 6:48 (CS4), 5:02 (CS6), 3:33 (CS8) and 2:52 min (CS10). With respect to the S2 acquisition, this corresponds to a 55, 67, 77 and 81% scan time reduction for CS4, CS6, CS8 and CS10, respectively (for absolute values, see Table 3). The mean respiratory gating efficiency was 51 ± 10%.

Table 3 Acquisition time, net flow, peak flow and peak velocity averaged over all healthy subjects of group 1 for different SENSE and CS acceleration factors

Exemplary magnitude and phase-contrast images as well as net flow curves show decreasing quality with increasing acceleration factor (Figs. 1 and 2). More caudal contours show stronger deviations (Fig. 2).

Fig. 1
figure1

Magnitude and phase-contrast images. Sagittal slices through the aorta of an exemplary dataset in peak systole. The different acceleration techniques and factors are shown in columns and labelled correspondingly. Magnitude (top row) and phase-contrast (bottom row) images with velocity encoding in foot-head direction (FH) are shown. Increased image artefacts can be observed for higher acceleration factors

Fig. 2
figure2

Net flow curves. Exemplary net flow curves extracted from datasets acquired with different acceleration factors. Each subfigure represents a different contour within the aorta. The arrows depict stronger deviations between acceleration factors. Main differences between the flow curves can be observed during peak systole. The three-dimensional figure of the aorta in the center shows the five contours where haemodynamic parameters were evaluated

These findings are confirmed by quantitative comparisons of NF, PF and PV as summarised in Table 3. Absolute discrepancies between S2 and CS8 averaged over the contours are 2.3 ml/cycle for NF 5.5 ml/s for PF and 8 cm/s for PV. Increasing deviations of NF and PF can be observed for CS10 (5.3 ml/cycle and 20.2 ml/s). A detailed analysis of MD, limits of agreement and CI is shown in the Bland-Altman plots in Fig. 3 and Additional file 1. All values are summarised in Table 4, revealing no significant differences in measured flow parameters up to a CS acceleration rate of 6 compared to S2. From an acceleration rate of 8, the NF shows a statistically significant underestimation of 2.2 ± 7.8 ml/cycle for CS8 and of 5.4 ± 8.0 ml/cycle for CS10. A significant underestimation of the PF was found for CS10 (14.4 ± 28.3 ml/s). In contrast to NF and PF, Bland-Altman analysis of PV tends towards negative biases for all acceleration factors (on average − 4%), which only for CS8 exceeds the CI of the MD. The SD of the MD show a continuous increase from CS4 to CS10 for PF (from 20.4 ml/s to 28.3 ml/s) and PV (from 8.1 cm/s to 12.2 cm/s).

Fig. 3
figure3

Bland-Altman analysis. Bland-Altman plots comparing (a) net flows in six contours and (b) peak velocities in the aortic root of S2 and compressed sensing (CS) accelerated scans. Data from different contours are indicated by different markers. The mean differences, standard deviations and their confidence intervals (red shaded areas) were calculated from the data points of all contours. The red data points denote the data of the volunteer whose net flow curves are shown in Fig. 2. The mean differences in net flow (a) between S2 and CS increase with increasing CS acceleration factor, while the standard deviations of the mean differences remain almost constant. The peak velocities (b) tend towards negative biases for all CS acceleration factors, while only exceeding the confidence interval of the mean difference for CS8. The standard deviation of the mean difference increases with increasing acceleration factor

Table 4 Summary of results of Bland-Altman analyses

In addition to the clinically relevant parameters NF, PF and PV, the accumulated flow error ER constantly increases with increasing acceleration factor (Fig. 4a). Particle trace analysis demonstrates an increased loss of particles emitted from the ascending aorta reaching the descending aorta for higher CS factors while the main difference is seen between the S2 and the CS4 scans (Fig. 4b).

Fig. 4
figure4

Advanced flow parameters. a Accumulated flow error ER as a function of CS acceleration factor. b Exemplary particle trace visualisation including the contours from which particles were ejected (ascending aorta, AAo) and arriving particles were counted (descending aorta, DAo1). c Percentage of particles ejected from the contour AAo reaching the contour DAo1 as a function of the acceleration factor

The Bland-Altman analysis of scan-rescan reproducibility demonstrates deviations for NF, PF and PV below 2.5% for all acceleration factors, while the deviations are lowest for S2 and increase with increasing acceleration factors CS6 and CS8 (Table 4, Additional file 2).

Inter- and intraobserver variabilities are similar for S2 and CS6 scans. Overall MD and SD of extracted flow parameters are under 1.6 and 4.2%, respectively (Table 4, Additional file 3).

Patients

The additional time needed for the 4D flow CMR acquisitions within the clinical examination protocol was well tolerated by all patients. The blood flow patterns of all patients are illustrated in Fig. 5 for selected time points in the cardiac cycle. Time-resolved movies are available as Additional files 4, 5 and 6. Patients #1 and #2, both with bicuspid aortic valves and ascending aortic dilatations, show pronounced right-handed helical flow patterns in the AAo; a smaller, slower helix is superimposed by a bigger, faster and less-twisted helix. In case of patient #1, helical flow in the brachiocephalic trunk and a small global vortex at the inner curvature of the arch can also be observed. In patient #3 with Marfan syndrome a vortical flow pattern can be observed at the inner curvature of the transition from the arch to the proximal DAo as well as helical flow along the DAo.

Fig. 5
figure5

Pathline screenshots. Particle traces ejected from the aorta show aberrant flow patterns in patients #1 and #2 with bicuspid aortic valves (a and b) and patient #3 with Marfan syndrome (c). Patients #1 and #2 show pronounced right-handed helical flow patterns in the AAo; a smaller, slower helix is superimposed by a larger, faster and less-twisted helix. In patient #3 a helical flow along the DAo can be observed. For the corresponding movies please refer to the Additional files 4, 5 and 6

Additional file 4: Pathline movie – Patient #1 (Bicuspid aortic valve). View of the aortic flow in patient #1 with bicuspid aortic valve. (MP4 16887 kb)

Additional file 5: Pathline movie – Patient #2 (Bicuspid aortic valve). View of the aortic flow in patient #2 with bicuspid aortic valve. (MP4 14391 kb)

Additional file 6: Pathline movie – Patient #3 (Marfan syndrome). View of the aortic flow in patient #3 with Marfan syndrome. (MP4 13540 kb)

Discussion

This study evaluated the use of a combined CS and parallel imaging technique for the acceleration of 4D flow CMR in the thoracic aorta. A thorough quantitative analysis performed in healthy subjects demonstrates its feasibility for the assessment of direct haemodynamic parameters while exemplary patient acquisitions demonstrate its potential application in clinical CMR protocols. The used CS and SENSE techniques are based on product implementations from Philips Healthcare and are therefore readily available for the Philips CMR community.

Reconstruction times of approximately 9 min can remain limiting in some cases, but a concurrent reconstruction and further acquisition of subsequent CMR scans is possible. An off-line reconstruction on dedicated hardware is also possible and could increase reconstruction speed, but remains debatable in a clinical environment due to raw data transfer speed, laborious PACS-integration and potential hardware approval issues. Finally, an acceleration of on-line reconstruction times could also be achieved by making use of higher computational power devices, e.g. graphical processor units.

Although quality of the magnitude images visually suffers from increased artefacts for higher CS factors, all flow curves extracted from the CS 4D flow data compare well with the S2 reference dataset. Deviations of the resulting flow curves are in line with deviations of k-t accelerated 4D flow CMR measurements compared to standard accelerated 4D flow or 2D PC-CMR acquisitions [7, 8, 10].

Quantitative values such as NF, PF and PV show overall small deviations when compared to S2 scans, while a trend showing increased MD for higher CS factors can be observed. Up to a CS acceleration rate of 6, Bland-Altman analyses show no significant differences in measured flow parameters. In contrast to NF and PF, Bland-Altman analysis of PV tends towards negative biases for all acceleration factors (on average − 4%), which only for CS8 exceeds the CI of the MD. The used CS sequence slightly overestimates the PV. This effect can be ascribed to higher noise values and therewith higher maximal velocities appearing in the contours.

In contrast to other studies showing increased temporal smoothing and therewith decreased PF and PV for higher acceleration rates [8, 28,29,30], this effect is not present in this approach, as no temporal correlations are used in the CS reconstruction. The potential influence of physiological variations was minimised by the use of a randomised acquisition order. Remaining eddy current related phase offsets [31] or gradient non-linearities [32] might be additional causes of the remaining systematic deviations in NF, PF and PV. However, as all parameters such as field of view, the velocity encoding gradient, repetition and echo times were kept constant between all acquisitions, an influence of such effects can be considered minimal.

The limits of agreement for all analysed flow parameters are in-line with previous studies [8, 10].

The results for the accumulated flow error ER, which is increasing with a slightly larger slope between CS8 and CS10, are consistent with the obtained results from the Bland-Altman analyses. This error metric reaches nearly 20% and implies that the acquisition of net flow at a single time frame is very sensitive to scan acceleration. Moreover, the large standard deviations in these error metrics could be attributed to additional uncertainties including physiological differences between healthy subjects, contour size, cardiac cycle length, distance between contours and flow velocities.

The increased particle loss for higher CS factors with the main difference between the S2 and the CS4 scans, is due to the higher undersampling. The spatial distribution of undersampling artefacts induced by the CS reconstruction might be an additional effect: depending on regularisation, noise-like artefacts may be increased leading to increased particle losses. A divergence-free reconstruction constraint is expected to minimise these artefacts while assuring physical properties of flowing blood [33,34,35].

An assessment of scan-rescan reproducibility and of inter- and intraobserver variabilities were performed to better classify the findings obtained from S2 versus CS scans. The scan-rescan reproducibility was highest for S2 and decreased with increasing acceleration factors CS6 and CS8 while the effect on the different, measured haemodynamic parameters is comparable and in line with findings from previous studies [36,37,38]. Inter- and intraobserver variabilities are similar for S2 and CS6 scans. Overall MD and SD of extracted flow parameters are in line with previously described results [36]. Furthermore, they remain small when compared to the scan-rescan differences.

Based on the healthy subject data, a relatively conservative CS factor of 6 was chosen for the patient acquisitions. The net scan times of aortic 4D flow CMR acquisition were in the order of 2 min resulting in total scan times of approximately 5 min, depending on heart rate and breathing navigator efficiency. The observed flow patterns in all three patients are in good agreement with previously described results: For instance, abnormal helical flow in the AAo in patients with bicuspid aortic valves, depending on the valvular phenotype has been described [19, 20]. Similarly, a significantly increased helical flow as well as a more frequently occurrence of vortical flow in the DAo in Marfan syndrome patients has been reported [17, 18]. Data quality achieved with CS therefore allows the measurement of these distinct haemodynamic features.

Although this study investigates a limited group size of 23 subjects (20 healthy subjects, 3 patients), the results demonstrate that CS allows a 4 to 8-fold acceleration of 4D flow CMR, corresponding to another 2 to 4-fold acceleration compared to typically used parallel imaging factors (S2) and is ready for clinical use for the assessment of primary haemodynamic parameters. The benefit in scan-time can in the future be traded, dependent on the clinical question, for shorter examination times, higher spatial and/or temporal resolution or acquiring several velocity-encoded volumes (Multi-Venc 4D flow CMR)[42]. Regarding advanced higher-order metrics such as wall shear stress, pulse wave velocities or pressure gradients, we hypothesize that a CS accelerated scan can be treated as a standard 4D flow scan with noise increasing with acceleration factors. As higher-order metrics are known to rely on high spatial, temporal and/or velocity resolutions, an improved calculation of these parameters might be achieved by using CS to allow for higher resolutions at constant scan time. However, this effect was not evaluated in the current study and needs to be investigated in future studies.

One limitation of the presented work is that the healthy subject study was carried out at a 3 T CMR system whereas the patients were scanned during a clinical routine CMR examination at a 1.5 T system. Apart from the field strength, both systems are from the same generation, used the same acquisition and reconstruction software and have identical coil configurations. While 3 T provides a higher signal-to-noise-ratio (SNR), the patient scans were performed after administration of contrast agent with an increased flip angle, increasing the SNR in the patient measurement. As shown by Hess et al. [39] similar results of data acquired at 1.5 T with contrast agent and at 3 T without contrast agent can be expected.

The current evaluation was focussed on aortic 4D flow CMR. Further studies are needed before transferring CS accelerated 4D flow CMR to other vascular beds and areas of interest with other velocity ranges and/or spatio-temporal dimensions.

Finally, dedicated CS approaches for velocity and phase images including exploiting spatio-temporal data redundancy, velocity-encoding-dimension sparsity constraints or physical priors are warranted to allow improved image quality and higher acceleration techniques [11, 33,34,35, 40, 41].

Conclusions

In summary, this study demonstrates that a six- to eightfold acceleration of 4D flow CMR using CS is feasible. Acquisitions in patients with aortopathies confirm the potential to integrate the proposed method in a clinical routine setting, whereby its main benefits are scan-time savings and direct on-line reconstruction.

Availability of data and materials

The imaging protocols as well as the datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

4D flow:

Three-dimensional time-resolved phase-contrast magnetic resonance

AAo:

Ascending aorta

Arch:

Aortic arch

CI:

Confidence interval

CMR:

Cardiovascular magnetic resonance

CS:

Compressed sensing

DAo:

Descending aorta

ECG:

Electrocardiogram

FH:

Foot-head

MD:

Mean difference

NF:

Net Flow

PC:

Phase-contrast

PF:

Peak flow

PV:

Peak velocity

Root:

Aortic root

SD:

Standard deviation

SENSE:

Sensitivity encoding

SNR:

Signal-to-noise-ratio

References

  1. 1.

    Kilner PJ, Gatehouse PD, Firmin DN. Flow measurement by magnetic resonance: a unique asset worth optimising. J Cardiovasc Magn Reson. 2007;9(4):723–8.

  2. 2.

    Markl M, Kilner PJ, Ebbers T. Comprehensive 4D velocity mapping of the heart and great vessels by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2011;13:7.

  3. 3.

    Kramer CM, Barkhausen J, Flamm SD, Kim RJ, Nagel E. Society for Cardiovascular Magnetic Resonance Board of trustees task force on standardized P. standardized cardiovascular magnetic resonance (CMR) protocols 2013 update. J Cardiovasc Magn Reson. 2013;15:91.

  4. 4.

    Kramer CM, Barkhausen J, Flamm SD, Kim RJ, Nagel E. Society for Cardiovascular Magnetic Resonance Board of trustees task force on standardized P. standardized cardiovascular magnetic resonance imaging (CMR) protocols, society for cardiovascular magnetic resonance: board of trustees task force on standardized protocols. J Cardiovasc Magn Reson. 2008;10:35.

  5. 5.

    Dyverfeldt P, Ebbers T. Comparison of respiratory motion suppression techniques for 4D flow MRI. Magn Reson Med. 2017;78(5):1877–82.

  6. 6.

    Dyverfeldt P, Bissell M, Barker AJ, Bolger AF, Carlhall CJ, Ebbers T, et al. 4D flow cardiovascular magnetic resonance consensus statement. J Cardiovasc Magn Reson. 2015;17:72.

  7. 7.

    Carlsson M, Toger J, Kanski M, Bloch KM, Stahlberg F, Heiberg E, et al. Quantification and visualization of cardiovascular 4D velocity mapping accelerated with parallel imaging or k-t BLAST: head to head comparison and validation at 1.5 T and 3 T. J Cardiovasc Magn Reson. 2011;13:55.

  8. 8.

    Schnell S, Markl M, Entezari P, Mahadewia RJ, Semaan E, Stankovic Z, et al. K-t GRAPPA accelerated four-dimensional flow MRI in the aorta: effect on scan time, image quality, and quantification of flow and wall shear stress. Magn Reson Med. 2014;72(2):522–33.

  9. 9.

    Jung B, Stalder AF, Bauer S, Markl M. On the undersampling strategies to accelerate time-resolved 3D imaging using k-t-GRAPPA. Magn Reson Med. 2011;66(4):966–75.

  10. 10.

    Giese D, Wong J, Greil GF, Buehrer M, Schaeffter T, Kozerke S. Towards highly accelerated Cartesian time-resolved 3D flow cardiovascular magnetic resonance in the clinical setting. J Cardiovasc Magn Reson. 2014;16:42.

  11. 11.

    Knobloch V, Boesiger P, Kozerke S. Sparsity transform k-t principal component analysis for accelerating cine three-dimensional flow measurements. Magn Reson Med. 2013;70(1):53–63.

  12. 12.

    Valvano G, Martini N, Huber A, Santelli C, Binter C, Chiappino D, et al. Accelerating 4D flow MRI by exploiting low-rank matrix structure and hadamard sparsity. Magn Reson Med. 2017;78(4):1330–41.

  13. 13.

    Cheng JY, Hanneman K, Zhang T, Alley MT, Lai P, Tamir JI, et al. Comprehensive motion-compensated highly accelerated 4D flow MRI with ferumoxytol enhancement for pediatric congenital heart disease. J Magn Reson Imaging. 2016;43(6):1355–68.

  14. 14.

    Dyvorne H, Knight-Greenfield A, Jajamovich G, Besa C, Cui Y, Stalder A, et al. Abdominal 4D flow MR imaging in a breath hold: combination of spiral sampling and dynamic compressed sensing for highly accelerated acquisition. Radiology. 2015;275(1):245–54.

  15. 15.

    Hsiao A, Lustig M, Alley MT, Murphy M, Chan FP, Herfkens RJ, et al. Rapid pediatric cardiac assessment of flow and ventricular volume with compressed sensing parallel imaging volumetric cine phase-contrast MRI. AJR Am J Roentgenol. 2012;198(3):W250–9.

  16. 16.

    Hsiao A, Yousaf U, Alley MT, Lustig M, Chan FP, Newman B, et al. Improved quantification and mapping of anomalous pulmonary venous flow with four-dimensional phase-contrast MRI and interactive streamline rendering. J Magn Reson Imaging. 2015;42(6):1765–76.

  17. 17.

    Geiger J, Hirtler D, Gottfried K, Rahman O, Bollache E, Barker AJ, et al. Longitudinal evaluation of aortic hemodynamics in Marfan syndrome: new insights from a 4D flow cardiovascular magnetic resonance multi-year follow-up study. J Cardiovasc Magn Reson. 2017;19(1):33.

  18. 18.

    Geiger J, Markl M, Herzer L, Hirtler D, Loeffelbein F, Stiller B, et al. Aortic flow patterns in patients with Marfan syndrome assessed by flow-sensitive four-dimensional MRI. J Magn Reson Imaging. 2012;35(3):594–600.

  19. 19.

    Rodriguez-Palomares JF, Dux-Santoy L, Guala A, Kale R, Maldonado G, Teixido-Tura G, et al. Aortic flow patterns and wall shear stress maps by 4D-flow cardiovascular magnetic resonance in the assessment of aortic dilatation in bicuspid aortic valve disease. J Cardiovasc Magn Reson. 2018;20(1):28.

  20. 20.

    Hope MD, Hope TA, Meadows AK, Ordovas KG, Urbania TH, Alley MT, et al. Bicuspid aortic valve: four-dimensional MR evaluation of ascending aortic systolic flow patterns. Radiology. 2010;255(1):53–61.

  21. 21.

    Pelc NJ, Bernstein MA, Shimakawa A, Glover GH. Encoding strategies for three-direction phase-contrast MR imaging of flow. J Magn Reson Imaging. 1991;1(4):405–13.

  22. 22.

    Buehrer M, Pruessmann KP, Boesiger P, Kozerke S. Array compression for MRI with large coil arrays. Magn Reson Med. 2007;57(6):1131–9.

  23. 23.

    Otazo R, Kim D, Axel L, Sodickson DK. Combination of compressed sensing and parallel imaging for highly accelerated first-pass cardiac perfusion MRI. Magn Reson Med. 2010;64(3):767–76.

  24. 24.

    Liang D, Liu B, Wang J, Ying L. Accelerating SENSE using compressed sensing. Magn Reson Med. 2009;62(6):1574–84.

  25. 25.

    Liu F, Duan Y, Peterson BS, Kangarlu A. Compressed sensing MRI combined with SENSE in partial k-space. Phys Med Biol. 2012;57(21):N391–403.

  26. 26.

    Bernstein MA, Zhou XJ, Polzin JA, King KF, Ganin A, Pelc NJ, et al. Concomitant gradient terms in phase contrast MR: analysis and correction. Magn Reson Med. 1998;39(2):300–8.

  27. 27.

    Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM, Yoganathan AP. Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. J Magn Reson Imaging. 1993;3(3):521–30.

  28. 28.

    Bollache E, Barker AJ, Dolan RS, Carr JC, van Ooij P, Ahmadian R, 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. Magn Reson Med. 2018;79(1):195–207.

  29. 29.

    Stadlbauer A, van der Riet W, Crelier G, Salomonowitz E. Accelerated time-resolved three-dimensional MR velocity mapping of blood flow patterns in the aorta using SENSE and k-t BLAST. Eur J Radiol. 2010;75(1):e15–21.

  30. 30.

    Tao Y, Rilling G, Davies M, Marshall I. Carotid blood flow measurement accelerated by compressed sensing: validation in healthy volunteers. Magn Reson Imaging. 2013;31(9):1485–91.

  31. 31.

    Gatehouse PD, Rolf MP, Graves MJ, Hofman MB, Totman J, Werner B, et al. Flow measurement by cardiovascular magnetic resonance: a multi-Centre multi-vendor study of background phase offset errors that can compromise the accuracy of derived regurgitant or shunt flow measurements. J Cardiovasc Magn Reson. 2010;12:5.

  32. 32.

    Markl M, Bammer R, Alley MT, Elkins CJ, Draney MT, Barnett A, et al. Generalized reconstruction of phase contrast MRI: analysis and correction of the effect of gradient field distortions. Magn Reson Med. 2003;50(4):791–801.

  33. 33.

    Santelli C, Loecher M, Busch J, Wieben O, Schaeffter T, Kozerke S. Accelerating 4D flow MRI by exploiting vector field divergence regularization. Magn Reson Med. 2016;75(1):115–25.

  34. 34.

    Ong F, Uecker M, Tariq U, Hsiao A, Alley MT, Vasanawala SS, et al. Robust 4D flow denoising using divergence-free wavelet transform. Magn Reson Med. 2015;73(2):828–42.

  35. 35.

    Busch J, Giese D, Wissmann L, Kozerke S. Reconstruction of divergence-free velocity fields from cine 3D phase-contrast flow measurements. Magn Reson Med. 2013;69(1):200–10.

  36. 36.

    Markl M, Wallis W, Harloff A. Reproducibility of flow and wall shear stress analysis using flow-sensitive four-dimensional MRI. J Magn Reson Imaging. 2011;33(4):988–94.

  37. 37.

    Zaman A, Motwani M, Oliver JJ, Crelier G, Dobson LE, Higgins DM, et al. 3.0T, time-resolved, 3D flow-sensitive MR in the thoracic aorta: impact of k-t BLAST acceleration using 8- versus 32-channel coil arrays. J Magn Reson Imaging. 2015;42(2):495–504.

  38. 38.

    Wentland AL, Grist TM, Wieben O. Repeatability and internal consistency of abdominal 2D and 4D phase contrast MR flow measurements. Acad Radiol. 2013;20(6):699–704.

  39. 39.

    Hess AT, Bissell MM, Ntusi NA, Lewis AJ, Tunnicliffe EM, Greiser A, et al. Aortic 4D flow: quantification of signal-to-noise ratio as a function of field strength and contrast enhancement for 1.5T, 3T, and 7T. Magn Reson Med. 2015;73(5):1864–71.

  40. 40.

    Hutter J, Schmitt P, Saake M, Stubinger A, Grimm R, Forman C, et al. Multi-dimensional flow-preserving compressed sensing (MuFloCoS) for time-resolved velocity-encoded phase contrast MRI. IEEE Trans Med Imaging. 2015;34(2):400–14.

  41. 41.

    Hutter J, Schmitt P, Aandal G, Greiser A, Forman C, Grimm R, et al. Low-rank and sparse matrix decomposition for compressed sensing reconstruction of magnetic resonance 4D phase contrast blood flow imaging (loSDeCoS 4D-PCI). Med Image Comput Comput Assist Interv. 2013;16(Pt 1:558–65.

  42. 42.

    Moersdorf R, Treutlein M, Kroeger JB, Ruijsink B, Wong J, Maintz D, Weiss K, Bunck AC, Baeßler B, Giese D. Precision, reproducibility and applicability of an undersampled multi-venc 4D flow MRI sequence for the assessment of cardiac hemodynamics. Magn Reson Imaging. 2019;61:73–82.

Download references

Acknowledgements

Not applicable.

Funding

Not applicable.

Author information

EN assisted with the design of the study, the set-up of the healthy subject scan protocol, recruited these subjects, collected and analysed all the data and drafted the manuscript. KW contributed to the design of the study and the supervision of the work and assisted with critical revision of the manuscript. RB did the second analysis of data for assessment of interobserver variability and critically revised the manuscript. JK assisted in the collection of patient data and critically revised the manuscript. DM helped to set up the patient acquisitions and assisted with drafting the work. DG designed the study, supervised the work, assisted with drafting and critical revision of the manuscript. All authors read and approved the final manuscript.

Correspondence to Elisabeth Neuhaus.

Ethics declarations

Ethics approval and consent to participate

Ethics approval was obtained by the local ethics committee of the faculty of medicine of Cologne University for the healthy subject as well as for the patient study. For each subject written informed consent was obtained.

Consent for publication

All subjects gave written consent to participate and to publish.

Competing interests

KW is an employee of Philips Healthcare Germany. All other authors declare that they have no competing interests. DG is an employee of Siemens Healthineers since December 2018, the entirety of his activities related to this work were however performed prior to or outside of his duties at Siemens.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Bland-Altman analysis of peak flows. Bland-Altman plots comparing peak flows of S2 and CS accelerated scans in six contours (indicated by different markers). The mean differences, standard deviations and their confidence intervals (red shaded areas) were calculated from the data points of all contours. The red data points denote the data of the volunteer whose net flow curves are shown in Fig. 2. (PDF 243 kb)

Additional file 2:

Bland-Altman analysis of scan-rescan agreement. Bland-Altman plots comparing (a) net flows, (b) peak flows and (c) peak velocities of in each case two repeated measurements of S2, CS6 and CS8 accelerated scans. The mean differences, standard deviations and their confidence intervals (red shaded areas) were calculated from the data points of all contours. (PDF 414 kb)

Additional file 3:

Bland-Altman analysis of inter- and intraobserver agreement. Bland-Altman plots comparing net flows and peak flows resulting from inter- and intraobserver analysis. The mean differences, standard deviations and their confidence intervals (red shaded areas) were calculated from the data points of both contours. (PDF 370 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Neuhaus, E., Weiss, K., Bastkowski, R. et al. Accelerated aortic 4D flow cardiovascular magnetic resonance using compressed sensing: applicability, validation and clinical integration. J Cardiovasc Magn Reson 21, 65 (2019) doi:10.1186/s12968-019-0573-0

Download citation

Keywords

  • 4D flow CMR
  • Aorta
  • Compressed sensing
  • SENSE
  • Phase-contrast CMR