Manifold Reconstruction Jean-Daniel Boissonnat Geometrica, INRIA http://www-sop.inria.fr/geometrica Winter School, University of Nice Sophia Antipolis January 26-30, 2015 Winter School 4 Manifold Reconstruction Sophia Antipolis 1 / 33
Geometric data analysis Images, text, speech, neural signals, GPS traces,... Geometrisation : Data = points + distances between points Hypothesis : Data lie close to a structure of “small” intrinsic dimension Problem : Infer the structure from the data Winter School 4 Manifold Reconstruction Sophia Antipolis 2 / 33
Submanifolds of R d A compact subset M ⊂ R d is a submanifold without boundary of (intrinsic) dimension k < d , if any p ∈ M has an open (topological) k -ball as a neighborhood in M R N M W φ R m U Intuitively, a submanifold of dimension k is a subset of R d that looks locally like an open set of an affine space of dimension k A curve a 1 -dimensional submanifold A surface is a 2 -dimensional submanifold More generally, manifolds are defined in an intrinsic way, Winter School 4 Manifold Reconstruction Sophia Antipolis 3 / 33 d
Triangulation of a submanifold We call triangulation of a submanifold M ⊂ R d a (geometric) simplicial complex ˆ M such that ˆ M is embedded in R d its vertices are on M it is homeomorphic to M Submanifold reconstruction The problem is to construct a triangulation ˆ M of some unknown submanifold M given a finite set of points P ⊂ M Winter School 4 Manifold Reconstruction Sophia Antipolis 4 / 33
Triangulation of a submanifold We call triangulation of a submanifold M ⊂ R d a (geometric) simplicial complex ˆ M such that ˆ M is embedded in R d its vertices are on M it is homeomorphic to M Submanifold reconstruction The problem is to construct a triangulation ˆ M of some unknown submanifold M given a finite set of points P ⊂ M Winter School 4 Manifold Reconstruction Sophia Antipolis 4 / 33
Issues in high-dimensional geometry Dimensionality severely restricts our intuition and ability to visualize data ⇒ need for automated and provably correct methods methods Complexity of data structures and algorithms rapidly grow as the dimensionality increases ⇒ no subdivision of the ambient space is affordable ⇒ data structures and algorithms should be sensitive to the intrinsic dimension (usually unknown) of the data Inherent defects : sparsity, noise, outliers Winter School 4 Manifold Reconstruction Sophia Antipolis 5 / 33
Issues in high-dimensional geometry Dimensionality severely restricts our intuition and ability to visualize data ⇒ need for automated and provably correct methods methods Complexity of data structures and algorithms rapidly grow as the dimensionality increases ⇒ no subdivision of the ambient space is affordable ⇒ data structures and algorithms should be sensitive to the intrinsic dimension (usually unknown) of the data Inherent defects : sparsity, noise, outliers Winter School 4 Manifold Reconstruction Sophia Antipolis 5 / 33
Issues in high-dimensional geometry Dimensionality severely restricts our intuition and ability to visualize data ⇒ need for automated and provably correct methods methods Complexity of data structures and algorithms rapidly grow as the dimensionality increases ⇒ no subdivision of the ambient space is affordable ⇒ data structures and algorithms should be sensitive to the intrinsic dimension (usually unknown) of the data Inherent defects : sparsity, noise, outliers Winter School 4 Manifold Reconstruction Sophia Antipolis 5 / 33
Looking for small and faithful simplicial complexes Need to compromise Size of the complex ◮ can we have dim ˆ M = dim M ? Efficiency of the construction algorithms and of the representations ◮ can we avoid the exponential dependence on d ? ◮ can we minimize the number of simplices ? Quality of the approximation ◮ Homotopy type & homology (RIPS complex, persistence) ◮ Homeomorphism (Delaunay-type complexes) Winter School 4 Manifold Reconstruction Sophia Antipolis 6 / 33
Sampling and distance functions [Niyogi et al.], [Chazal et al.] d K : x → inf p ∈ K � x − p � Distance to a compact K : Stability If the data points C are close (Hausdorff) to the geometric structure K , the topology and the geometry of the offsets K r = d − 1 ([ 0 , r ]) and C r = d − 1 ([ 0 , r ]) are close Winter School 4 Manifold Reconstruction Sophia Antipolis 7 / 33
Distance functions and triangulations ˇ Nerve theorem (Leray) The nerve of the balls (Cech complex) and the union of balls have the same homotopy type Winter School 4 Manifold Reconstruction Sophia Antipolis 8 / 33
Questions + The homotopy type of a compact set X can be computed from the ˘ Cech complex of a sample of X + The same is true for the α -complex – The ˘ Cech and the α -complexes are huge ( O ( n d ) and O ( n ⌈ d / 2 ⌉ ) ) and very difficult to compute – Both complexes are not in general homeomorphic to X (i.e. not a triangulation of X ) – The ˘ Cech complex cannot be realized in general in the same space as X Winter School 4 Manifold Reconstruction Sophia Antipolis 9 / 33
˘ Cech and Rips complexes The Rips complex is easier to compute but still very big, and less precise in approximating the topology α Winter School 4 Manifold Reconstruction Sophia Antipolis 10 / 33
An example where no offset has the right topology ! Persistent homology at rescue ! Winter School 4 Manifold Reconstruction Sophia Antipolis 11 / 33
The curses of Delaunay triangulations in higher dimensions Their complexity depends exponentially on the ambient dimension. Robustness issues become very tricky Higher dimensional Delaunay triangulations are not thick even if the vertices are well-spaced The restricted Delaunay triangulation is no longer a good approximation of the manifold even under strong sampling conditions (for d > 2 ) Winter School 4 Manifold Reconstruction Sophia Antipolis 12 / 33
3D Delaunay Triangulations are not thick even if the vertices are well-spaced Each square face can be circumscribed by an empty sphere This remains true if the grid points are slightly perturbed therefore creating thin simplices Winter School 4 Manifold Reconstruction Sophia Antipolis 13 / 33
Badly-shaped simplices Badly-shaped simplices lead to bad geometric approximations which in turn may lead to topological defects in Del |M ( P ) [Oudot] see also [Cairns], [Whitehead], [Munkres], [Whitney] Winter School 4 Manifold Reconstruction Sophia Antipolis 14 / 33
Tangent space approximation M is a smooth k -dimensional manifold ( k > 2 ) embedded in R d Bad news [Oudot 2005] The Delaunay triangulation restricted to M may be a bad approximation of the manifold even if the sample is dense t = ∆ + δ / 2 t = ∆ c w v p p 0 c 0 u t t y y z z x x Winter School 4 Manifold Reconstruction Sophia Antipolis 15 / 33
Whitney’s angle bound and tangent space approximation Lemma [Whitney 1957] If σ is a j -simplex whose vertices all lie within a distance η from a hyperplane H ⊂ R d , then sin ∠ ( aff ( σ ) , H ) ≤ 2 j η D ( σ ) Corollary If σ is a j -simplex, j ≤ k , vert ( σ ) ⊂ M , ∆( σ ) ≤ δ rch ( M ) δ ∀ p ∈ σ, sin ∠ ( aff ( σ ) , T p ) ≤ Θ( σ ) ∆( σ ) 2 ( η ≤ 2 rch ( M ) by the Chord Lemma) Winter School 4 Manifold Reconstruction Sophia Antipolis 16 / 33
The assumptions M is a differentiable submanifold of positive reach of R d The dimension k of M is small P is an ε -net of M , i.e. ∀ x ∈ M , ∃ p ∈ P , � x − p � ≤ ε rch ( M ) ◮ ◮ ∀ p , q ∈ P , � p − q � ≥ ¯ η ε ε is small enough Winter School 4 Manifold Reconstruction Sophia Antipolis 17 / 33
The tangential Delaunay complex [B. & Ghosh 2010] M T p p Construct the star of p ∈ P in the Delaunay triangulation Del Tp ( P ) 1 of P restricted to T p Del T M ( P ) = � p ∈ P star ( p ) 2 Winter School 4 Manifold Reconstruction Sophia Antipolis 18 / 33
+ Del T M ( P ) ⊂ Del ( P ) + star ( p ) , Del T p ( P ) and therefore Del T M ( P ) can be computed without computing Del ( P ) – Del T M ( P ) is not necessarily a triangulated manifold Winter School 4 Manifold Reconstruction Sophia Antipolis 19 / 33
Construction of Del T p ( P ) Given a d -flat H ⊂ R , Vor ( P ) ∩ H is a weighted Voronoi diagram in H p ′ j � x − p i � 2 ≤ � x − p j � 2 x i � 2 + � p i − p ′ i � 2 ≤ � x − p ′ j � 2 + � p j − p ′ p ′ � x − p ′ j � 2 ⇔ i p j H p i Corollary: construction of Del T p ψ p ( p i ) = ( p ′ i , −� p i − p ′ i � 2 ) (weighted point) project P onto T p which requires O ( Dn ) time 1 construct star ( ψ p ( p i )) in Del ( ψ p ( p i )) ⊂ T p i 2 star ( p i ) ≈ star ( ψ p ( p i )) (isomorphic ) 3 Winter School 4 Manifold Reconstruction Sophia Antipolis 20 / 33
Recommend
More recommend
