Towards persistence-based reconstruction in euclidean spaces

Frédéric Chazal, Steve Yann Oudot
2008 Proceedings of the twenty-fourth annual symposium on Computational geometry - SCG '08  
Manifold reconstruction has been extensively studied for the last decade or so, especially in two and three dimensions. Recent advances in higher dimensions have led to new methods to reconstruct large classes of compact subsets of R d . However, the complexities of these methods scale up exponentially with d, making them impractical in medium or high dimensions, even on data sets of low intrinsic dimensionality. In this paper, we introduce a novel approach that stands in-between classical
more » ... struction and topological estimation, and whose complexity scales up with the intrinsic dimension of the data. Our algorithm combines two paradigms: greedy refinement, and topological persistence. Given a point cloud in R d , we build a set of landmarks iteratively, while maintaining a nested pair of abstract complexes, whose images in R d lie close to the data, and whose persistent homology eventually coincides with the homology of the underlying shape. When the data points are densely sampled from a smooth m-submanifold X of R d , our method retrieves the homology of X in time at most c(m)n 5 , where n is the size of the input and c(m) is a constant depending solely on m. To prove the correctness of our algorithm, we investigate onČech, Rips, and witness complex filtrations in Euclidean spaces. More precisely, we show how previous results on unions of balls can be transposed toČech filtrations, and from there to Rips and witness complex filtrations. Finally, investigating further on witness complexes, we quantify a conjecture of Carlsson and de Silva, which states that witness complex filtrations should have cleaner persistence barcodes thanČech or Rips filtrations, at least on smooth submanifolds of Euclidean spaces.
doi:10.1145/1377676.1377719 dblp:conf/compgeom/ChazalO08 fatcat:vge7xs7oprdythqpvvwad3blla