Validation of compressed sensing accelerated 2D flow MRI in the common carotid arteries
Eva S. Peper1, Wouter V. Potters1, Bram F. Coolen1, Henk A. Marquering1, Gustav J. Strijkers1, Pim van Ooij1, and Aart J. Nederveen1

1Radiology, Academic Medical Center (AMC), Amsterdam, Netherlands

Synopsis

Flow MRI of the carotid arteries has emerged as an important imaging field during the last decade and has been shown valuable for the assessment of hemodynamics in the context of atherosclerosis. In this study a 2D flow acquisition was accelerated using random undersampling and compressed sensing reconstruction. For validation of the reconstructed flow a phantom experiment was performed at different acceleration factors followed by an in vivo scan of the carotid arteries.

Target audience:

MR scientists, Clinical researchers interested in accelerated flow MRI

Purpose:

Flow MRI of the carotid arteries has emerged as an important imaging field during the last decade and has been shown valuable for the assessment of hemodynamics in the context of atherosclerosis1. Applications of volumetric flow measurements in the carotids are nowadays still hampered by excessive scan durations and limited temporal resolutions. We hypothesize that flow measurements can be accelerated by using random undersampling of retrospectively triggered acquisitions in combination with a compressed sensing reconstruction and that validations of accelerated flow measurements can be performed using a pulsatile flow loop setup.

Methods:

Two different experiments were conducted on a 3T MRI scanner (Ingenia Philips Healthcare, NL, software release 5.1.8). First a phantom validation was performed using a pulsatile flow loop setup (LifeTec, Eindhoven, NL). The phantom mimics a human carotid artery and consists of two tubes in a rigid 3D printed block surrounded by agar (Figure 1). An air pressure driven pump created a pulsating water flow between 0.3-6.0 ml/s. The tubes were connected so that flow is occurring in opposite directions in the imaging plane. A conventional head coil was used for signal acquisition. In a second experiment flow measurements were performed in a healthy volunteer using a 8 channel surface receive neck coil. A transversal slice in the neck through left and right common carotid artery was chosen as a measurement plane. For both setups a 2D phase contrast (PC) sequence with a foot-head flow encoding direction was used (VENC 100 cm/s). Two different spatial resolutions were chosen for phantom and volunteer scans (Table 1). 24 cardiac phases were acquired using retrospective PPU gating for the volunteer scan and a peripheral pulse signal for the phantom, leading to a temporal resolution of 40 ms. In all experiments undersampling was performed by acquiring a subset of k-lines in phase encoding direction. The scanner software was adapted to acquire a predefined number of k-lines. K-space was fully sampled at the center within 20% of the matrix size. All remaining k-lines were undersampled randomly (Figure 2). The resulting acceleration factors were 1.5, 1.7, 2.0, 2.2, 2.5, 3.0 and 3.3. Raw data was imported to MATLAB using ReconFrame (Gyrotools, Zurich, CH). Compressed sensing (CS) reconstruction was performed using the BART toolbox2. Sensitivity maps were estimated from the data3. An l1 wavelet regularization was performed4 using a regularization parameter of 0.01. After reconstruction a linear phase offset correction was performed on the phase difference maps using static regions as a reference. To quantify flow in the regions of interest masking was done by thresholding magnitude images. The segmentation was performed on every dataset and was used for all cardiac phases.

Results:

Figure 1 shows representative velocity images as obtained in the phantom and in the volunteer. Flow curves for accelerated acquisitions did not show substantial deviations from the fully sampled measurement (Figure 2). The deviation from the fully sampled dataset was on average -0.12+-0.13 ml/s for a higher resolution (0.64x0.64 mm2) and 0.07+-0.23 ml/s for a lower resolution (0.8x0.8 mm2). No dependency between acceleration factor and the accuracy of the flow measurement with respect to the fully sampled dataset could be observed. The in vivo experiment shows successful application of CS acceleration for 2D flow in the carotid arteries at two different resolutions (Figure 3). As shown in Figure 4 no clear trend could be observed between the amount of acceleration and the deviation of the fully sampled flow curve. The deviation from the fully sampled dataset in vivo was -0.20+-0.50 ml/s for a higher resolution (0.8x0.8 mm2) and 0.03+-0.50 ml/s for a lower resolution (1x1 mm2).

Discussion:

The variation between the different acceleration factors in vivo can only partially be explained by the undersampled acquisition. Here also physiological variations play a role in combination with the limited temporal resolution of the measurement. The phantom validation convincingly illustrates the reliability of the proposed CS undersampling strategy for flow measurements.

Conclusion:

Flow measurements can be accelerated by using random undersampling and compressed sensing reconstruction. 2D flow can be accelerated 3-fold using compressed sensing without losing accuracy in flow estimation. Future application of this strategy will focus on 4D flow measurements using the same flow loop setup for validation purposes.

Acknowledgements

No acknowledgement found.

References

1. Potters W.V., Marquering H.A., VanBavel E., Nederveen A. J., Measuring Wall Shear Stress Using Velocity-Encoded MRI., Curr Cardiovasc Imaging Rep 2014, 7:9257

2. Takaya N., Yuan C., Chu B. et al., Association Between Carotid Plaque Characteristics and Subsequent Ischemic Cerebrovascular Events: A Prospective Assessment With MRI − Initial Results., Stroke 2006; 37(3):818–823.

3. Uecker M., Ong F., Tamir J.I., Bahri D., Virtue P., Cheng J.Y., Zhang T., and Lustig M., Berkeley Advanced Reconstruction Toolbox, Annual Meeting ISMRM, Toronto 2015, In Proc. Intl. Soc. Mag. Reson. Med. 23:2486

4. Lustig M., Donoho D. and Pauly J.M., Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med 2007, 58: 1182-1195.

5. Uecker M., Lai P., Murphy M.J., Virtue P., Elad M., Pauly J.M., Vasanawala S.S.,Lustig M.. ESPIRiT - An Eigenvalue Approach to Autocalibrating Parallel MRI: Where SENSE meets GRAPPA. Magn Reson Med 2014; 71: 990-1001.

Figures

Table 1: Scan parameters for the phantom and volunteer scans using a 2D phase contrast sequence.

Figure 1: Randomly undersampled k-space of matrix size 148x148 with an acceleration factor of 1.7 (a). Velocity map of the phantom (b) and the carotids (c) using the same scaling. The red tube shows inverted velocities, due to a returning flow.

Figure 2: Flow curves of the phantom tubes at 0.64x0.64 mm2 (A) and 0.8x0.8 mm2 (B) in plane resolution. Red curves show the returning flow. Their absolute values were inverted for better comparison. Images C and D show the differences in absolute flow between the accelerated scans and the fully sampled acquisitions.

Figure 3: Flow curves of the in vivo experiment at 0.8x0.8 mm2 (A) and 1x1 mm2 (B) in plane resolution for left and right carotid artery. Each flow curve was derived from images acquired with a different acceleration factor.

Figure 4: Mean flow over the entire heart cycle for all acceleration factors in the in vivo experiment.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
2703