2,100 research outputs found
Manifold Reconstruction in Arbitrary Dimensions using Witness Complexes
A preliminary version of this paper has been presented at the ACM Symp. on Computational Geometry 2007International audienceIt is a well-established fact that the witness complex is closely related to the restricted Delaunay triangulation in low dimensions. Specifically, it has been proved that the witness complex coincides with the restricted Delaunay triangulation on curves, and is still a subset of it on surfaces, under mild sampling conditions. In this paper, we prove that these results do not extend to higher-dimensional manifolds, even under strong sampling conditions such as uniform point density. On the positive side, we show how the sets of witnesses and landmarks can be enriched, so that the nice relations that exist between restricted Delaunay triangulation and witness complex hold on higher-dimensional manifolds as well. We derive from our structural results an algorithm that reconstructs manifolds of any arbitrary dimension or co-dimension at different scales. The algorithm combines a farthest-point refinement scheme with a vertex pumping strategy. It is very simple conceptually, and it does not require the input point sample to be sparse. Its running time is bounded by c(d)n2, where n is the size of the input point cloud, and c(d) is a constant depending solely (yet exponentially) on the dimension d of the ambient space. Although this running time makes our reconstruction algorithm rather theoretical, recent work has shown that a variant of our approach can be made tractable in arbitrary dimensions, by building upon the results of this paper
Towards Persistence-Based Reconstruction in Euclidean Spaces
Manifold reconstruction has been extensively studied for the last decade or
so, especially in two and three dimensions. Recently, significant improvements
were made in higher dimensions, leading to new methods to reconstruct large
classes of compact subsets of Euclidean space . However, the complexities
of these methods scale up exponentially with d, which makes them impractical in
medium or high dimensions, even for handling low-dimensional submanifolds. In
this paper, we introduce a novel approach that stands in-between classical
reconstruction and topological estimation, and whose complexity scales up with
the intrinsic dimension of the data. Specifically, when the data points are
sufficiently densely sampled from a smooth -submanifold of , our
method retrieves the homology of the submanifold in time at most ,
where is the size of the input and is a constant depending solely on
. It can also provably well handle a wide range of compact subsets of
, though with worse complexities. Along the way to proving the
correctness of our algorithm, we obtain new results on \v{C}ech, Rips, and
witness complex filtrations in Euclidean spaces
Only distances are required to reconstruct submanifolds
In this paper, we give the first algorithm that outputs a faithful
reconstruction of a submanifold of Euclidean space without maintaining or even
constructing complicated data structures such as Voronoi diagrams or Delaunay
complexes. Our algorithm uses the witness complex and relies on the stability
of power protection, a notion introduced in this paper. The complexity of the
algorithm depends exponentially on the intrinsic dimension of the manifold,
rather than the dimension of ambient space, and linearly on the dimension of
the ambient space. Another interesting feature of this work is that no explicit
coordinates of the points in the point sample is needed. The algorithm only
needs the distance matrix as input, i.e., only distance between points in the
point sample as input.Comment: Major revision, 16 figures, 47 page
A Growing Self-Organizing Network for Reconstructing Curves and Surfaces
Self-organizing networks such as Neural Gas, Growing Neural Gas and many
others have been adopted in actual applications for both dimensionality
reduction and manifold learning. Typically, in these applications, the
structure of the adapted network yields a good estimate of the topology of the
unknown subspace from where the input data points are sampled. The approach
presented here takes a different perspective, namely by assuming that the input
space is a manifold of known dimension. In return, the new type of growing
self-organizing network presented gains the ability to adapt itself in way that
may guarantee the effective and stable recovery of the exact topological
structure of the input manifold
Approximating Loops in a Shortest Homology Basis from Point Data
Inference of topological and geometric attributes of a hidden manifold from
its point data is a fundamental problem arising in many scientific studies and
engineering applications. In this paper we present an algorithm to compute a
set of loops from a point data that presumably sample a smooth manifold
. These loops approximate a {\em shortest} basis of the
one dimensional homology group over coefficients in finite field
. Previous results addressed the issue of computing the rank of
the homology groups from point data, but there is no result on approximating
the shortest basis of a manifold from its point sample. In arriving our result,
we also present a polynomial time algorithm for computing a shortest basis of
for any finite {\em simplicial complex} whose edges have
non-negative weights
Approximating Local Homology from Samples
Recently, multi-scale notions of local homology (a variant of persistent
homology) have been used to study the local structure of spaces around a given
point from a point cloud sample. Current reconstruction guarantees rely on
constructing embedded complexes which become difficult in high dimensions. We
show that the persistence diagrams used for estimating local homology, can be
approximated using families of Vietoris-Rips complexes, whose simple
constructions are robust in any dimension. To the best of our knowledge, our
results, for the first time, make applications based on local homology, such as
stratification learning, feasible in high dimensions.Comment: 23 pages, 14 figure
Constructing Intrinsic Delaunay Triangulations of Submanifolds
We describe an algorithm to construct an intrinsic Delaunay triangulation of
a smooth closed submanifold of Euclidean space. Using results established in a
companion paper on the stability of Delaunay triangulations on -generic
point sets, we establish sampling criteria which ensure that the intrinsic
Delaunay complex coincides with the restricted Delaunay complex and also with
the recently introduced tangential Delaunay complex. The algorithm generates a
point set that meets the required criteria while the tangential complex is
being constructed. In this way the computation of geodesic distances is
avoided, the runtime is only linearly dependent on the ambient dimension, and
the Delaunay complexes are guaranteed to be triangulations of the manifold
- …