Towards a glaucoma risk index based on simulated hemodynamics from fundus images

Glaucoma is the leading cause of irreversible but preventable blindness in the world. Its major treatable risk factor is the intra-ocular pressure, although other biomarkers are being explored to improve the understanding of the pathophysiology of the disease. It has been recently observed that glaucoma induces changes in the ocular hemodynamics. However, its effects on the functional behavior of the retinal arterioles have not been studied yet. In this paper we propose a first approach for characterizing those changes using computational hemodynamics. The retinal blood flow is simulated using a 0D model for a steady, incompressible non Newtonian fluid in rigid domains. The simulation is performed on patient-specific arterial trees extracted from fundus images. We also propose a novel feature representation technique to comprise the outcomes of the simulation stage into a fixed length feature vector that can be used for classification studies. Our experiments on a new database of fundus images show that our approach is able to capture representative changes in the hemodynamics of glaucomatous patients. Code and data are publicly available in https://ignaciorlando.github.io.

[1]  P J Blanco,et al.  A computational approach to generate concurrent arterial networks in vascular territories , 2013, International journal for numerical methods in biomedical engineering.

[2]  Geert Molenberghs,et al.  Ocular blood flow in glaucoma – the Leuven Eye Study , 2016, Acta ophthalmologica.

[3]  Dan Liu,et al.  Image-based Blood Flow Simulation in the Retinal Circulation , 2009 .

[4]  T. Wong,et al.  Global prevalence of glaucoma and projections of glaucoma burden through 2040: a systematic review and meta-analysis. , 2014, Ophthalmology.

[5]  Calvin R. Maurer,et al.  A Linear Time Algorithm for Computing Exact Euclidean Distance Transforms of Binary Images in Arbitrary Dimensions , 2003, IEEE Trans. Pattern Anal. Mach. Intell..

[6]  Evelien Vandewalle,et al.  Lack of spontaneous venous pulsation: possible risk indicator in normal tension glaucoma? , 2013, Acta ophthalmologica.

[7]  Evelien Vandewalle,et al.  Disturbed correlation between arterial resistance and pulsatility in glaucoma patients , 2012, Acta ophthalmologica.

[8]  Uwe Himmelreich,et al.  Clinical Metabolomics and Glaucoma , 2017, Ophthalmic Research.

[9]  Martin Rumpf,et al.  A Continuous Skeletonization Method Based on Level Sets , 2002, VisSym.

[10]  P Ganesan,et al.  Analysis of retinal circulation using an image-based network model of retinal vasculature. , 2010, Microvascular research.

[11]  Charles E. Riva,et al.  Retinal Blood Flow Evaluation , 2012, Ophthalmologica.

[12]  Miguel O Bernabeu,et al.  Computational fluid dynamics assisted characterization of parafoveal hemodynamics in normal and diabetic eyes using adaptive optics scanning laser ophthalmoscopy. , 2016, Biomedical optics express.

[13]  Giovanna Guidoboni,et al.  Ocular Hemodynamics and Glaucoma: The Role of Mathematical Modeling , 2013, European journal of ophthalmology.

[14]  B. Al-Diri,et al.  Hemodynamics in the retinal vasculature during the progression of diabetic retinopathy , 2017 .

[15]  Pietro Perona,et al.  A Bayesian hierarchical model for learning natural scene categories , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[16]  Paul Mitchell,et al.  Retinal vessel diameter and open-angle glaucoma: the Blue Mountains Eye Study. , 2005, Ophthalmology.

[17]  A. Pries,et al.  Biophysical aspects of blood flow in the microvasculature. , 1996, Cardiovascular research.

[18]  Elena De Momi,et al.  Blood vessel segmentation algorithms - Review of methods, datasets and evaluation metrics , 2018, Comput. Methods Programs Biomed..

[19]  Kirk Roberts,et al.  Representation Learning for Retinal Vasculature Embeddings , 2017, FIFI/OMIA@MICCAI.