Statistical Shape Analysis of Shape Graphs with Applications to Retinal Blood-Vessel Networks

This paper provides theoretical and computational developments in statistical shape analysis of shape graphs, and demonstrates them using analysis of complex data from retinal blood-vessel (RBV) networks. The shape graphs are represented by a set of nodes and edges (planar articulated curves) connecting some of these nodes. The goals are to utilize shapes of edges and connectivities and locations of nodes to: (1) characterize full shapes, (2) quantify shape differences, and (3) model statistical variability. We develop a mathematical representation, elastic Riemannian shape metrics, and associated tools for such statistical analysis. Specifically, we derive tools for shape graph registration, geodesics, summaries, and shape modeling. Geodesics are convenient for visualizing optimal deformations, and PCA helps in dimension reduction and statistical modeling. One key challenge here is comparisons of shape graphs with vastly different complexities (in number of nodes and edges). This paper introduces a novel multi-scale representation of shape graphs to handle this challenge. Using the notions of (1) ``effective resistance" to cluster nodes and (2) elastic shape averaging of edge curves, one can reduce shape graph complexity while maintaining overall structures. This way, we can compare shape graphs by bringing them to similar complexity. We demonstrate these ideas on Retinal Blood Vessel (RBV) networks taken from the STARE and DRIVE databases.