Connected Components on Lie Groups and Applications to Multi-Orientation Image Analysis

We develop and analyze a new algorithm to find the connected components of a compact set I from a Lie group G endowed with a left-invariant Riemannian distance. For a given $\delta >0$ , the algorithm finds the largest cover of I such that all sets in the cover are separated by at least distance $\delta$ . We call the sets in the cover the $\delta$ -connected components of I (closely related to $\check{\text{C}}$ ech complexes of radius $\delta /2$ ). The grouping relies on an iterative procedure involving morphological dilations with Hamilton-Jacobi-Bellman kernels on G and notions of $\delta$ -thickened sets. We prove that the algorithm converges in finitely many iteration steps. We find the optimal value for $\delta$ using persistence diagrams. We also propose to use specific affinity matrices. This allows for grouping of $\delta$ -connected components based on their local proximity and alignment. Among the many different applications of the algorithm, in this article, we focus on illustrating that the method can efficiently identify (possibly overlapping) branches in complex vascular trees on retinal images. This is done by applying an orientation score transform to the images that allows us to view them as functions from $L_2\left(G\right)$ where $G=SE\left(2\right)$ , the Lie group of roto-translations. By applying our algorithm in this Lie group, we illustrate that we obtain $\delta$ -connected components that differentiate between crossing structures and that group well-aligned, nearby structures. This contrasts standard connected component algorithms in $R^2$ .