3661

Validation of real-time phase contrast MRI with online compressed sensing reconstruction in phantom and patients
Tania Lala1, Lea Christierson2,3, Petter Frieberg1, Daniel Giese4,5, Nina Hakacova 2, Pia Sjöberg1, Ellen Ostenfeld1, and Johannes Töger1
1Clinical Physiology, Department of Clinical Sciences Lund, Skåne University Hospital, Lund University, Lund, Sweden, 2Department of Clinical Sciences Lund, Pediatric Heart Center, Skåne University Hospital, Lund University, Lund, Sweden, 3Department of Biomedical Engineering, Lund University, Lund, Sweden, 4Magnetic Resonance, Siemens Healthcare GmbH, Erlangen, Germany, 5Institute of Radiology, University Hospital Erlangen, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany

Synopsis

Keywords: Flow, Velocity & Flow

Real-time Phase Contrast MRI with Compressed Sensing reconstruction is a promising method to image cardiac flow in patients with arrhythmia, where cardiac-gated sequences fail. In this study, we compared it against the gated standard clinical flow MRI using a phantom model and in patients. The phantom underwent pulsatile periodic and non-periodic flow. For in vivo validation, non-arrhythmic patient data (N=10) and data from patient with atrial fibrillation (N=1) were collected from the ascending aorta. Non-periodic flow was captured both in the phantom and in vivo. Real-time MRI showed good accuracy in net flow quantification and underestimated peak flow.

Background

Phase Contrast (PC) MRI is a widely used clinical tool to quantify cardiac and vessel blood flow. However, standard cardiac-gated methods fail to capture flow during arrhythmias1. Real-Time (RT) PC MRI enables fast, free-breathing, arrhythmia-robust acquisitions, but may suffer from degraded temporal and spatial resolution compared to gated flow imaging2 or require offline reconstruction schemes3–5. By using compressed sensing (CS), it is possible to achieve high acceleration factors with maintained spatial and temporal resolution while keeping a fast, inline reconstruction6. We hypothesized that RT PC MRI combined with CS could be a clinically relevant method for flow imaging during arrhythmias.
Therefore, the aim of this study was to validate the performance of a research Real-Time PC sequence with CS reconstruction in a phantom and in non-arrhythmic patients against the standard cardiac-gated clinical flow sequence, for future use in arrhythmic patients.

Methods

Phantom and patient scans were performed using a 1.5T clinical scanner (MAGNETOM Sola, Siemens Healthcare, Erlangen, Germany). The research RT CS sequence was tested with two accelerations: 1) with acceleration factor R=8 and shared velocity encoding7 (SVE) and 2) with R=13.5. The pulse sequence parameters of the research RT CS and standard clinical sequences are summarized in Figure 1.
The phantom included a straight tube (20 mm inner diameter) in the feet-head direction to simulate the ascending aorta (Figure 2, panel A). The phantom was connected to a pulsatile pump driven by a programmable motor set to generate flow rates in five steps between 3 and 8 l/min, with pulsation frequency set to 60 beats per minute (N=5 experiments).
The experimental protocol included timer and beaker measurements followed by through-plane MRI at a plane representing the ascending aorta with the clinical and RT sequences, followed by a second timer and beaker measurement to assess experiment stability. The clinical sequence was gated by a signal generated by the motor at the start of each cycle. Timer and beaker measurements were compared to the clinical sequence for validation of the set up (Figure 2, Panels B, C, D).
Real-time net flow (in liters/minute as calculated from the real-time flow curves) was compared to timer and beaker measurements, and peak flow rate (mean peak flow rate over all cycles) was compared to the clinical sequence. To investigate performance in arrhythmias, two sequential flow curves with different stroke volumes were programmed to occur periodically. The clinical sequence was gated according to the periodicity. The resulting flow curve was compared to what RT CS captured.
The same MRI protocol was used in patients (N=10) without arrhythmias for in vivo validation and a patient (N=1) with atrial fibrillation to test feasibility in arrhythmia. Through-plane measurements were performed at the ascending aorta and net flow and peak flow rate were calculated, comparing RT to the clinical sequence. Written, informed consent was obtained from all patients.
Image analysis was performed using the software Segment8 version 3.3 (Medviso AB, Lund, Sweden). Phase background correction was performed in all PC images using static ROIs9. Linear regression and modified Bland-Altman analyses were done in GraphPad Prism version 9.4.1 (GraphPad Software, San Diego, California, USA), and used to compare RT flow with clinical MRI flow and timer-beaker flow.

Results

Images from the phantom scans are shown in Figure 3, Panel A. Timer and beaker before and after MRI differed less than 2% for all pump settings. Linear regression and modified Bland-Altman analysis of net and peak flow quantification are shown in Figure 4, Panel A for RT CS with R=8, SVE (net flow bias: 2.5±1.3%, peak flow rate bias: -8.7±1.9%) and Panel B for RT CS with R=13.5 (net flow bias -0.1%±1.2%, peak flow rate bias -7.1±2.1%). Figure 4, Panel C shows an example of an RT CS periodic flow curve in the phantom and in Panel D each cycle of this example is compared to the curve by the standard clinical sequence. Figure 4, Panel C shows how the sequences captured phantom non-periodic flow.
Figure 3, Panel B shows images from one patient. Linear regression and modified Bland Altman analysis plots are shown in Figure 5, Panel A for RT CS with R=8, SVE (net flow bias=-3.1±8.9%, peak flow bias= -17.0±5.78%) and Panel B for RT CS with R=13.5 (net flow bias= -6.1±9.4%, peak flow bias=-14.0±4.7%). Figure 5, panel C shows an example of periodic flow in one of the non-arrhythmic patients where each cycle was compared to the corresponding curve as assessed by the standard clinical sequence (Panel D). Finally, Panel E shows how arrhythmic flow was captured in the atrial fibrillation patient with the RT (R=13.5) sequence.

Conclusion

Real-time PC MRI with online compressed sensing reconstruction can quantify net flow accurately but underestimates peak flow measurements. This is in line with previous studies using offline reconstructions2,3. Moreover, RT CS was feasible in a patient with arrhythmia, where gated flow imaging would fail.

Acknowledgements

No acknowledgement found.

References

1.Gatehouse PD, Keegan J, Crowe LA, et al. Applications of phase-contrast flow and velocity imaging in cardiovascular MRI. Eur Radiol. 2005;15(10):2172-2184.

2.Traber J, Wurche L, Dieringer MA, et al. Real-time phase contrast magnetic resonance imaging for assessment of haemodynamics: from phantom to patients. Eur Radiol. 2016;26(4):986-996.

3.Haji‐Valizadeh H, Guo R, Kucukseymen S, et al. Highly accelerated free‐breathing real‐time phase contrast cardiovascular MRI via complex‐difference deep learning. Magn Reson Med. 2021;86(2):804-819.

4.Sun A, Zhao B, Li Y, He Q, Li R, Yuan C. Real-time phase-contrast flow cardiovascular magnetic resonance with low-rank modeling and parallel imaging. J Cardiovasc Magn Reson. 2017;19(1):19.

5.Haji‐Valizadeh H, Feng L, Ma LE, et al. Highly accelerated, real‐time phase‐contrast MRI using radial k ‐space sampling and GROG‐GRASP reconstruction: a feasibility study in pediatric patients with congenital heart disease. NMR in Biomedicine. 2020;33(5). doi:10.1002/nbm.4240

6.Xiong Fei, et al. Highly Accelerated Free-Breathing Real-Time 2D Flow Imaging using Compressed Sensing and Shared Velocity Encoding. SCMR 25th annual scientific sessions (2022)

7.Lin HY, Bender JA, Ding Y, et al. Shared velocity encoding: A method to improve the temporal resolution of phase-contrast velocity measurements. Magnetic Resonance Medicine. 2012;68(3):703-710.

8.Heiberg E, Sjögren J, Ugander M, Carlsson M, Engblom H, Arheden H. Design and validation of Segment - freely available software for cardiovascular image analysis. BMC Med Imaging. 2010;10(1):1.

9.Lankhaar JW, Hofman MBM, Marcus JT, Zwanenburg JJM, Faes TJC, Vonk-Noordegraaf A. Correction of phase offset errors in main pulmonary artery flow quantification. J Magn Reson Imaging. 2005;22(1):73-79.

Figures

Typical pulse sequence parameters for the standard clinical and real-time flow sequences

Experimental phantom set-up: A) Sagittal view of the phantom. The black line specifies the imaging plane of interest (ascending aorta plane). Arrows indicate the flow direction B) MRI images of the coronal view of the phantom at the plane of interest with the standard clinical flow sequence C) pulsatile flow curve in the phantom and D) Standard clinical MRI compared to timer and beaker measurements. Linear regression: solid line is the identity line and dashed line is the regression line, and modified Bland Altman: solid line shows bias and dashed lines are the 95% limits of agreement.

Example images of the phantom (A) and in-vivo (B) at peak systole with the standard clinical sequence (top row), RT, CS R=8, SVE (middle row) and RT, CS, R=13.5 (bottom row). The gray colormap on the bottom indicates that the pixels in the velocity encoded images range between [-200 cm/s, 200 cm/s].

Flow analysis in pulsatile phantom. RT CS compared to timer and beaker for net flow and to standard clinical MRI for peak flow quantification. Linear regression and modified Bland-Altman plots for A) RT CS, R=8, SVE and B) RT CS, R=13.5. Linear regression: line of identity is solid, and the regression line is dashed. Bland-Altman: bias shown by the solid line and the 95% limits of agreement by the dashed. C) RT CS in periodic pulsatile flow, D) individual cycles from example C compared to standard clinical MRI, E) RT CS and standard clinical MRI in non-periodic flow.

Flow analysis in patients. RT CS compared to standard clinical MRI for net flow and peak flow quantification. Linear regression and modified Bland-Altman plots for A) RT CS, R=8, SVE and B) RT CS, R=13.5. Linear regression: line of identity is solid, and the regression line is dashed. Bland-Altman: bias shown by the solid line and the 95% limits of agreement by the dashed. C) RT CS in a non-arrhythmic patient, D) individual cycles from example C compared to standard clinical MRI, E) RT CS, R=13.5 in atrial fibrillation (double arrows emphasize different RR intervals).

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)
3661
DOI: https://doi.org/10.58530/2023/3661