A Novel Anisotropic Fast Marching Method and its Application to Blood Flow Computation in Phase-contrast MRI

Summary

Background: Phase-contrast MRI (PC MRI) can be used to assess blood flow dynamics noninvasively inside the human body. The acquired images can be reconstructed into flow vector fields. Traditionally, streamlines can be computed based on the vector fields to visualize flow patterns and particle trajectories.

Objectives: The traditional methods may give a false impression of precision, as they do not consider the measurement uncertainty in the PC MRI images. In our prior work, we incorporated the uncertainty of the measurement into the computation of particle trajectories.

Methods: As a major part of the contribution, a novel numerical scheme for solving the anisotropic Fast Marching problem is presented. A computing time comparison to state-of-theart methods is conducted on artificial tensor fields. A visual comparison of healthy to pathological blood flow patterns is given.

Results: The comparison shows that the novel anisotropic Fast Marching solver outperforms previous schemes in terms of computing time. The visual comparison of flow patterns directly visualizes large deviations of pathological flow from healthy flow.

Conclusions: The novel anisotropic Fast Marching solver efficiently resolves even strongly anisotropic path costs. The visualization method enables the user to assess the uncertainty of particle trajectories derived from PC MRI images.

Keywords

Cardiovascular abnormalities - vascular diseases - computer-assisted diagnosis - anisotropic fast marching - uncertainty visualization

• References

• 1 Sethian JA. Level Set Methods and Fast Marching Methods Evolving Interfaces in Computational Ge-ometry, Fluid Mechanics, Computer Vision, and Materials Science. Cambridge University Press;; 1999
• 2 Deschamps T, Cohen LD. Fast extraction of minimal paths in 3D images and application to virtual endoscopy. Med Image Anal 2001; 5 (04) 281-299.
• 3 Cohen LD, Kimmel R. Global minimum for active contour models: A minimal path approach. Int J Comput Vision 1997; 24 (01) 57-78.
• 4 Benmansour F, Cohen LD. Tubular Structure Segmentation Based on Minimal Path Method and Anisotropic Enhancement. Int J Comput Vision. 2010: 1-19.
• 5 Prados E, Lenglet C, Pons JP, Wotawa N, Deriche R, Faugeras O. et al. Control Theory and Fast Marching Techniques for Brain Connectivity Mapping. In: Proceedings of CVPR’06. 2006: 1076-1083.
• 6 Konukoglu E, Sermesant M, Clatz O, Peyrat JM, Delingette H, Ayache N. A recursive anisotropic Fast Marching approach to reaction diffusion equation: application to tumor growth modeling. In: Proceedings of IPMI’07. 2007: 687-699.
• 7 Schwenke M, Hennemuth A, Fischer B, Friman O. Blood Particle Trajectories in Phase-Contrast-MRI as Minimal Paths Computed with Anisotropic Fast Marching. In: Bildverarbeitung für die Medizin. 2011: 289-293.
• 8 Handels H, Ehrhardt J. Medical image computing for computer-supported diagnostics and therapy. Advances and perspectives. Methods Inf Med 2009; 48 (01) 11-17.
• 9 Darmofal DL, Haimes R. An analysis of 3D particle path integration algorithms. J Comput Physics 1996; 123 (01) 182-195.
• 10 Friman O, Hennemuth A, Harloff A, Bock J, Markl M, Peitgen HO. Probabilistic 4D blood flow mapping. In: Proceedings of MICCAI’10. 2010: 416-423.
• 11 Schwenke M, Hennemuth A, Fischer B, Friman O. Blood Flow Computation in Phase-Contrast MRI by Minimal Paths in Anisotropic Media. In: Proceedings of MICCAI’11 2011; 01: 436-443.
• 12 Cristiani E. A Fast Marching Method for Hamilton-Jacobi Equations Modeling Monotone Front Propagations. J Sci Comput 2009; 39 (02) 189-205.
• 13 Kao CY, Osher S, Tsai YH. Fast Sweeping Methods For Static Hamilton-Jacobi Equations. SIAM J Numer Anal 2004; 42 (06) 2612-2632.
• 14 Sethian JA, Vladimirsky A. Ordered Upwind Methods for Static Hamilton-Jacobi Equations: Theory and Algorithms. SIAM J Numer Anal 2003; 41 (01) 325-363.
• 15 Jbabdi S, Bellec P, Toro R, Daunizeau J, PélégriniIssac M, Benali H. Accurate anisotropic Fast March-ing for diffusion-based geodesic tractography. Int J Biomed Imaging. 2008