2269

4D flow MRI vs. Doppler Echocardiography: Data Integration Framework for Direct Comparison of Intra-cardiac Flow Dynamics
Jeesoo Lee1,2, Nadia El Hangouche2, Alex J Barker3, James D Thomas2, and Michael Markl1,4
1Radiology, Northwestern University, Chicago, IL, United States, 2Cardiology, Northwestern Memorial Hospital, Chicago, IL, United States, 3Radiology, University of Colorado Denver, Denver, CO, United States, 4Biomedical Engineering, Northwestern University, Chicago, IL, United States

Synopsis

We developed a framework to integrate 4D flow MRI with color Doppler echocardiography (CDE) for direct localized inter-modal velocity comparison. The framework combines CDE and 4D flow MRI data by using cine-MRI and a two-step image registration; 1) registration of CDE to cine-MRI and 2) registration of cine-MRI to 4D flow MRI. As a result, a voxel-wise inter-modal velocity difference map overlaid to a tissue image was constructed. Six healthy subjects with a same-day CDE and 4D flow MRI were utilized and a strong correlation with moderate agreement was found in the LVOT.

Introduction

4D flow MRI is an emerging technique capable of visualizing complex 3D cardiovascular flow. Previous studies have demonstrated that 4D flow MRI can assess 3D hemodynamic parameters (e.g. wall shear stress distribution1, 3D flow patterns2, and viscous energy loss3) associated with cardiovascular diseases. In order to demonstrate clinical utility, comparisons must be made with Doppler echocardiography, the reference standard for clinical cardiovascular flow assessment. Previous comparison studies have focused on peak velocity4-8 and shown 4D flow MRI can measure comparable velocity when encoding velocity was appropriately set prior to a scan4. However, those studies did not match the measurement location, and thus, are not fully investigating inter-modality agreement. This limitation can be tackled by using 2D color Doppler echocardiography (CDE) which measures the planar velocity field. The goal of this study was to develop a framework for matching measurement locations of ‘same-day’ 4D flow MRI and CDE to perform voxel-wise velocity comparison.

Methods

Six healthy volunteers underwent free-breathing 4D flow MRI (spatial resolution=3.4-3.8 x 2.4-2.5 x 2.4-2.8 mm3, temporal resolution=38.4-40.6 ms, venc=150 cm/s, 1.5T Aera: Siemens Healthcare, Erlangen, Germany) and transthoracic echocardiography on the same day. Baseline characteristics of the subjects are listed in Table 1. 4D flow MRI datasets were preprocessed by using MATLAB (R2018b, Mathworks Inc., Natick, USA) for phase offset correction, noise masking, anti-aliasing as described previously9. MRI velocity datasets were segmented in 3D, including the aorta and left ventricular outflow tract (LVOT) using a previously described machine learning approach10. Three-chamber (3C) CDE datasets with a clear appearance of velocity in the LVOT were exported from an ultrasound scanner and reconstructed with MATLAB. A workflow for registering CDE image pixels to 4D flow MRI voxels is illustrated in Figure 1. First, 3C CDE tissue image was manually registered with 3C standard cine-MRI by matching aortic valve leaflet hinge points and a posterior mitral annulus junction point seen in both images (Figure 1A). Second, the cine image was registered with 4D flow MRI by using rigid transformation combined with optimization algorithm in MATLAB. This step was necessary since cine-MRI was acquired at end inspiration. The left ventricle with the aorta in the cine image was segmented and the mask data was used to perform registration with a 3D aorta segmentation mask. A generalized pattern search algorithm11 was implemented to automatically find a location of the global maximum cross-correlation coefficient between the masks (Figure 1B). Then, CDE image pixels were transferred to 4D flow MRI space by using the image transformation matrices obtained after the previous steps. Lastly, 4D flow MRI velocity vectors were interpolated to CDE pixel locations and projected to ultrasound beamlines (Figure 1C) and compared with CDE velocity data. Peak systole was chosen as a reference cardiac phase for comparison. A velocity difference map was constructed in the data overlapping region. Linear regression and Bland-Altman analysis was performed for the velocities in the LVOT.

Results

Qualitative inspection of CDE image registered to cine image showed a reasonable overlap of anatomical features such as mitral valve leaflets and aorta wall (Figure 2). Cine to 4D flow MRI registration improved Dice score between the cine segmentation masks and the 3D aorta mask cross-sections (initially 0.49 ± 0.30, after 0.80 ± 0.04, p = 0.035, Figure 3). Figure 4A shows a reformatted 4D flow MRI velocity (VMRI) and a CDE velocity (VCDE) map from a subject who had the best overall velocity resemblance. Figure 4B shows voxel-wise velocity difference maps from all subjects. Large discrepant voxels (|VMRI-VCDE| > 0.6 m/s) were mostly located in the downstream of the aortic valve. There was a strong correlation of velocities in the LVOT (r = 0.79, p < 0.001) with a moderate agreement (bias -13.4% and limits of agreement ± 41.4%).

Discussion

The framework proposed in this study enabled the comparison of 4D flow MRI with Doppler echocardiography at the same anatomical location. The constructed velocity difference map was beneficial to assess inter-modal agreement with respect to anatomical locations. We confirmed that both modalities provide comparable velocities in the LVOT where velocities are relatively uniform and well-aligned with each other. On the contrary, two modalities seem to disagree when the flow is disturbed by anatomical structures as seen by the presence of large discrepant voxels near and behind the aortic valve. However, this study is limited by the number of subjects (n=6) and the absence of patients. Future work will focus on a larger cohort including patients with aortic valve diseases.

Conclusion

We developed a novel framework that integrates 4D flow MRI and CDE and investigated voxel-wise velocity inter-modal agreement in LVOT. The framework is able to generate a velocity difference map superimposed onto a tissue image and this data can provide useful information to better understand the regional reliability of both modalities.

Acknowledgements

No acknowledgement found.

References

1. Van Ooij P, Potters WV, Nederveen AJ, Allen BD, Collins J, Carr J, Malaisrie SC, Markl M and Barker AJ. A methodology to detect abnormal relative wall shear stress on the full surface of the thoracic aorta using four‐dimensional flow MRI. Magnetic resonance in medicine. 2015;73:1216-1227.

2. Mahadevia R, Barker AJ, Schnell S, Entezari P, Kansal P, Fedak PW, Malaisrie SC, McCarthy P, Collins J and Carr J. Bicuspid aortic cusp fusion morphology alters aortic three-dimensional outflow patterns, wall shear stress, and expression of aortopathy. Circulation. 2014;129:673-682.

3. Barker AJ, van Ooij P, Bandi K, Garcia J, Albaghdadi M, McCarthy P, Bonow RO, Carr J, Collins J and Malaisrie SC. Viscous energy loss in the presence of abnormal aortic flow. Magnetic resonance in medicine. 2014;72:620-628.

4. Gabbour M, Schnell S, Jarvis K, Robinson JD, Markl M and Rigsby CK. 4-D flow magnetic resonance imaging: blood flow quantification compared to 2-D phase-contrast magnetic resonance imaging and Doppler echocardiography. Pediatric radiology. 2015;45:804-813.

5. Harloff A, Albrecht F, Spreer J, Stalder A, Bock J, Frydrychowicz A, Schöllhorn J, Hetzel A, Schumacher M and Hennig J. 3D blood flow characteristics in the carotid artery bifurcation assessed by flow‐sensitive 4D MRI at 3T. Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine. 2009;61:65-74.

6. Huh H, Kinno M, Thomas JD, Markl M and Barker AJ. Estimation of aortic valve effective orifice area: a same day comparison between Doppler echocardiography and 4D flow MRI. International Society for Magnetic Resonance in Medicine. 2018.

7. Markl M, Lee DC, Furiasse N, Carr M, Foucar C, Ng J, Carr J and Goldberger JJ. Left atrial and left atrial appendage 4D blood flow dynamics in atrial fibrillation. Circulation: Cardiovascular Imaging. 2016;9:e004984.

8. Rose MJ, Jarvis K, Chowdhary V, Barker AJ, Allen BD, Robinson JD, Markl M, Rigsby CK and Schnell S. Efficient method for volumetric assessment of peak blood flow velocity using 4D flow MRI. Journal of Magnetic Resonance Imaging. 2016;44:1673-1682.

9. Walker PG, Cranney GB, Scheidegger MB, Waseleski G, Pohost GM and Yoganathan AP. Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. Journal of Magnetic Resonance Imaging. 1993;3:521-530.

10. Berhane H, Scott M, Robinson JD, Rigsby CK and Markl M. 3D U-Net for Automated Segmentation of the Thoracic Aorta in 4D-Flow derived 3D PC-MRA. International Society for Magnetic Resonance in Medicine. 2019.

11. Audet C and Dennis Jr JE. Analysis of generalized pattern searches. SIAM Journal on optimization. 2002;13:889-903.

Figures

Table 1. Subject characteristics (n = 6)

Figure 1. 4D flow MRI and color Doppler echocardiography (CDE) data integration workflow. (A) A CDE tissue image was manually registered to a 3C cine MRI by matching anatomical landmarks (white arrows). (B) The cine image was registered with 4D flow MRI by matching the segmentation masks. (C) 4D flow MRI velocity vectors were interpolated to CDE pixels and projected to ultrasound beamlines for comparison.

Figure 2. 2D CDE to 2D cine MRI image registration results. The orange-colored regions indicate CDE background tissue and the white lines indicate the CDE field-of-view.

Figure 3. 2D cine MRI to 4D flow MRI registration results. The cyan colored region indicates the portion of a cine segmentation mask within a 4D flow MRI measurement volume. The pink-colored region indicates the cross-section of the 3D aorta mask at the cine image plane.

Figure 4. Voxel-wise velocity comparison results. (A1) Reformatted 4D flow MRI velocity, VMRI, map and (A2) original color Doppler echocardiography (CDE) velocity, VCDE map from the subject #5. (B) Peak systole velocity difference (VMRI−VCDE) map from all subjects. Black dotted line box represents the left ventricular outflow tract (LVOT) region. (C1) Linear regression plot with the identity curve (dotted line) and (C2) Bland-Altman analysis results for LVOT velocity comparison. All subjects are included in the analysis.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
2269