Computational Topology in Reconstruction, Mesh Generation, and Data Analysis
|
|
- Laura Small
- 5 years ago
- Views:
Transcription
1 Computational Topology in Reconstruction, Mesh Generation, and Data Analysis Tamal K. Dey Department of Computer Science and Engineering The Ohio State University Dey (2014) Computational Topology CCCG 14 1 / 55
2 Outline Topological concepts: Dey (2014) Computational Topology CCCG 14 2 / 55
3 Outline Topological concepts: Topological spaces Dey (2014) Computational Topology CCCG 14 2 / 55
4 Outline Topological concepts: Topological spaces Maps Dey (2014) Computational Topology CCCG 14 2 / 55
5 Outline Topological concepts: Topological spaces Maps Complexes Dey (2014) Computational Topology CCCG 14 2 / 55
6 Outline Topological concepts: Topological spaces Maps Complexes Homology groups Dey (2014) Computational Topology CCCG 14 2 / 55
7 Outline Topological concepts: Topological spaces Maps Complexes Homology groups Applications: Dey (2014) Computational Topology CCCG 14 2 / 55
8 Outline Topological concepts: Topological spaces Maps Complexes Homology groups Applications: Manifold reconstruction Dey (2014) Computational Topology CCCG 14 2 / 55
9 Outline Topological concepts: Topological spaces Maps Complexes Homology groups Applications: Manifold reconstruction Delaunay mesh generation Dey (2014) Computational Topology CCCG 14 2 / 55
10 Outline Topological concepts: Topological spaces Maps Complexes Homology groups Applications: Manifold reconstruction Delaunay mesh generation Topological data analysis Dey (2014) Computational Topology CCCG 14 2 / 55
11 Topologcal spaces Topology Background A point set with open subsets closed under union and finite intersections Dey (2014) Computational Topology CCCG 14 3 / 55
12 Topologcal spaces Topology Background A point set with open subsets closed under union and finite intersections d-ball B d {x R d x 1} d-sphere S d {x R d x = 1} k-manifold: neighborhoods homeomorphic to open k-ball 2-sphere, torus, double torus are 2-manifolds Dey (2014) Computational Topology CCCG 14 3 / 55
13 Maps Topology Background Homeomorphism h : T 1 T 2 where h is continuous, bijective and has continuous inverse Isotopy : continuous deformation that maintains homeomorphism homotopy equivalence: map linked to continuous deformation only Homeomorphic No isotopy Dey (2014) Computational Topology CCCG 14 4 / 55
14 Simplicial complex Topology Background Abstract V (K): vertex set, k-simplex: (k + 1)-subset σ V (K) Complex K = {σ σ σ = σ K} Dey (2014) Computational Topology CCCG 14 5 / 55
15 Simplicial complex Topology Background Abstract V (K): vertex set, k-simplex: (k + 1)-subset σ V (K) Complex K = {σ σ σ = σ K} Geometric k-simplex: k + 1-point convex hull Complex K: (i) t K if t is a face of t K (ii) t 1, t 2 K t 1 t 2 is a face of both Dey (2014) Computational Topology CCCG 14 5 / 55
16 Simplicial complex Topology Background Abstract V (K): vertex set, k-simplex: (k + 1)-subset σ V (K) Complex K = {σ σ σ = σ K} Geometric k-simplex: k + 1-point convex hull Complex K: (i) t K if t is a face of t K (ii) t 1, t 2 K t 1 t 2 is a face of both Triangulation: K is a triangulation of a topological space T if T K Dey (2014) Computational Topology CCCG 14 5 / 55
17 Reconstruction Surface Reconstruction Dey (2014) Computational Topology CCCG 14 6 / 55
18 Sampling Sampling Sample P Σ R 3 Dey (2014) Computational Topology CCCG 14 7 / 55
19 Local Feature Size Sampling Lfs(x) is the distance to medial axis Dey (2014) Computational Topology CCCG 14 8 / 55
20 Sampling ε-sample (Amenta-Bern-Eppstein 98) Each x has a sample within εlfs(x) distance Dey (2014) Computational Topology CCCG 14 9 / 55
21 Sampling Crust and Cocone Guarantees Theorem (Crust: Amenta-Bern 1999) Any point x Σ is within O(ε)Lfs(x) distance from a point in the output. Conversely, any point of the output surface has a point x Σ within O(ε)Lfs(x) distance for ε < Dey (2014) Computational Topology CCCG / 55
22 Sampling Crust and Cocone Guarantees Theorem (Crust: Amenta-Bern 1999) Any point x Σ is within O(ε)Lfs(x) distance from a point in the output. Conversely, any point of the output surface has a point x Σ within O(ε)Lfs(x) distance for ε < Theorem (Cocone: Amenta-Choi-Dey-Leekha 2000) The output surface computed by Cocone from an ε sample is homeomorphic to the sampled surface for ε < Dey (2014) Computational Topology CCCG / 55
23 Sampling Restricted Voronoi/Delaunay Definition Restricted Voronoi: Vor P Σ : Intersection of Vor (P) with the surface/manifold Σ. Dey (2014) Computational Topology CCCG / 55
24 Sampling Restricted Voronoi/Delaunay Definition Restricted Delaunay: Del P Σ : dual of Vor P Σ Dey (2014) Computational Topology CCCG / 55
25 Sampling Topology Closed Ball property (Edelsbrunner, Shah 94) If restricted Voronoi cell is a closed ball in each dimension, then Del P Σ is homeomorphic to Σ. Dey (2014) Computational Topology CCCG / 55
26 Sampling Topology Closed Ball property (Edelsbrunner, Shah 94) If restricted Voronoi cell is a closed ball in each dimension, then Del P Σ is homeomorphic to Σ. Dey (2014) Computational Topology CCCG / 55
27 Sampling Topology Closed Ball property (Edelsbrunner, Shah 94) If restricted Voronoi cell is a closed ball in each dimension, then Del P Σ is homeomorphic to Σ. Theorem For a sufficiently small ε if P is an ε-sample of Σ, then (P, Σ) satisfies the closed ball property, and hence Del P Σ Σ. Dey (2014) Computational Topology CCCG / 55
28 Sampling Topology Closed Ball property (Edelsbrunner, Shah 94) If restricted Voronoi cell is a closed ball in each dimension, then Del P Σ is homeomorphic to Σ. Dey (2014) Computational Topology CCCG / 55
29 Sampling Topology Closed Ball property (Edelsbrunner, Shah 94) If restricted Voronoi cell is a closed ball in each dimension, then Del P Σ is homeomorphic to Σ. Dey (2014) Computational Topology CCCG / 55
30 Boundaries Input Variations Ambiguity in reconstruction Non-homeomorphic Restricted Delaunay [DLRW09] Non-orientabilty Dey (2014) Computational Topology CCCG / 55
31 Input Variations Boundaries Ambiguity in reconstruction Non-homeomorphic Restricted Delaunay [DLRW09] Non-orientabilty Theorem (D.-Li-Ramos-Wenger 2009) Given a sufficiently dense sample of a smooth compact surface Σ with boundary one can compute a Delaunay mesh isotopic to Σ. Dey (2014) Computational Topology CCCG / 55
32 Input Variations Open: Reconstructing nonsmooth surfaces Guarantee of homeomorphism is open Dey (2014) Computational Topology CCCG / 55
33 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Dey (2014) Computational Topology CCCG / 55
34 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Dey (2014) Computational Topology CCCG / 55
35 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Dey (2014) Computational Topology CCCG / 55
36 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Dey (2014) Computational Topology CCCG / 55
37 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Use (ε, δ)-sampling Dey (2014) Computational Topology CCCG / 55
38 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Use (ε, δ)-sampling Restricted Delaunay does not capture topology Slivers are arbitrarily oriented [CDR05] Del P Σ Σ no matter how dense P is. Dey (2014) Computational Topology CCCG / 55
39 High Dimensions High Dimensional PCD Curse of dimensionality (intrinsic vs. extrinsic) Reconstruction of submanifolds brings ambiguity Use (ε, δ)-sampling Restricted Delaunay does not capture topology Slivers are arbitrarily oriented [CDR05] Del P Σ Σ no matter how dense P is. Delaunay triangulation becomes harder Dey (2014) Computational Topology CCCG / 55
40 Reconstruction High Dimensions Theorem (Cheng-Dey-Ramos 2005) Given an (ε, δ)-sample P of a smooth manifold Σ R d for appropriate ε, δ > 0, there is a weight assignment of P so that Del ˆP Σ Σ which can be computed efficiently. Dey (2014) Computational Topology CCCG / 55
41 Reconstruction High Dimensions Theorem (Cheng-Dey-Ramos 2005) Given an (ε, δ)-sample P of a smooth manifold Σ R d for appropriate ε, δ > 0, there is a weight assignment of P so that Del ˆP Σ Σ which can be computed efficiently. Theorem (Chazal-Lieutier 2006) Given an ε-noisy sample P of manifold Σ R d, there exists r p ρ(σ) for each p P so that the union of balls B(p, r p ) is homotopy equivalent to Σ. Dey (2014) Computational Topology CCCG / 55
42 High Dimensions Reconstructing Compacts Dey (2014) Computational Topology CCCG / 55
43 High Dimensions Reconstructing Compacts Lfs vanishes, introduce µ-reach and define (ε, µ)-samples. Dey (2014) Computational Topology CCCG / 55
44 High Dimensions Reconstructing Compacts Lfs vanishes, introduce µ-reach and define (ε, µ)-samples. Theorem (Chazal-Cohen-S.-Lieutier 2006) Given an (ε, µ)-sample P of a compact K R d for appropriate ε, µ > 0, there is an α so that union of balls B(p, α) is homotopy equivalent to K η for arbitrarily small η. Dey (2014) Computational Topology CCCG / 55
45 Mesh generation Surface and volume mesh Dey (2014) Computational Topology CCCG / 55
46 Mesh generation Surface and volume mesh Dey (2014) Computational Topology CCCG / 55
47 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 Dey (2014) Computational Topology CCCG / 55
48 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 To mesh some domain D Dey (2014) Computational Topology CCCG / 55
49 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 To mesh some domain D 1 Initialize points P D, compute Del P Dey (2014) Computational Topology CCCG / 55
50 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 To mesh some domain D 1 Initialize points P D, compute Del P 2 If some condition is not satisfied, insert a point p D into P and repeat Dey (2014) Computational Topology CCCG / 55
51 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 To mesh some domain D 1 Initialize points P D, compute Del P 2 If some condition is not satisfied, insert a point p D into P and repeat 3 Return Del P D Dey (2014) Computational Topology CCCG / 55
52 Delaunay refinement Mesh generation Pioneered by Chew89,Ruppert92,Shewchuk98 To mesh some domain D 1 Initialize points P D, compute Del P 2 If some condition is not satisfied, insert a point p D into P and repeat 3 Return Del P D Burden is to show termination (by packing argument) Dey (2014) Computational Topology CCCG / 55
53 Mesh generation Sampling Theorem Theorem (Amenta-Bern 98, Cheng-D.-Edelsbrunner-Sullivan 01) If P Sigma is a discrete ε-sample of a smooth surface Σ, then for ε < 0.09, Del P Σ satisfies: Dey (2014) Computational Topology CCCG / 55
54 Mesh generation Sampling Theorem Theorem (Amenta-Bern 98, Cheng-D.-Edelsbrunner-Sullivan 01) If P Sigma is a discrete ε-sample of a smooth surface Σ, then for ε < 0.09, Del P Σ satisfies: It is homeomorphic to Σ Dey (2014) Computational Topology CCCG / 55
55 Mesh generation Sampling Theorem Theorem (Amenta-Bern 98, Cheng-D.-Edelsbrunner-Sullivan 01) If P Sigma is a discrete ε-sample of a smooth surface Σ, then for ε < 0.09, Del P Σ satisfies: It is homeomorphic to Σ Each triangle has normal aligning within O(ε) angle to the surface normals Dey (2014) Computational Topology CCCG / 55
56 Mesh generation Sampling Theorem Theorem (Amenta-Bern 98, Cheng-D.-Edelsbrunner-Sullivan 01) If P Sigma is a discrete ε-sample of a smooth surface Σ, then for ε < 0.09, Del P Σ satisfies: It is homeomorphic to Σ Each triangle has normal aligning within O(ε) angle to the surface normals Hausdorff distance between Σ and Del P Σ is O(ε 2 ) of LFS. Dey (2014) Computational Topology CCCG / 55
57 Mesh generation Sampling Theorem Modified Theorem (Boissonnat-Oudot 05) If P Σ is such that each Voronoi edge-surface intersection x lies within εlfs(x) from a sample, then for ε < 0.09, Del P Σ satisfies: It is homeomorphic to Σ Each triangle has normal aligning within O(ε) angle to the surface normals Hausdorff distance between Σ and Del P Σ is O(ε 2 ) of LFS. Dey (2014) Computational Topology CCCG / 55
58 Mesh generation Basic Delaunay Refinement 1 Initialize points P Σ, compute Del P 2 If some condition is not satisfied, insert a point c Σ into P and repeat 3 Return Del P Σ Dey (2014) Computational Topology CCCG / 55
59 Mesh generation Surface Delaunay Refinement 1 Initialize points P Σ, compute Del P 2 If some Voronoi edge intersects Σ at x with d(x, P) > εlfs(x), insert x in P and repeat 3 Return Del P Σ Dey (2014) Computational Topology CCCG / 55
60 Mesh generation Difficulty How to compute Lfs(x)? Can be approximated by computing approximatemedial axis needs a dense sample. Dey (2014) Computational Topology CCCG / 55
61 Mesh generation Difficulty How to compute Lfs(x)? Can be approximated by computing approximatemedial axis needs a dense sample. Dey (2014) Computational Topology CCCG / 55
62 A Solution Mesh generation Replace d(x, P) < εlfs(x) with d(x, P) < λ, an user parameter Dey (2014) Computational Topology CCCG / 55
63 Mesh generation A Solution Replace d(x, P) < εlfs(x) with d(x, P) < λ, an user parameter Topology guarantee is lost Dey (2014) Computational Topology CCCG / 55
64 Mesh generation A Solution Replace d(x, P) < εlfs(x) with d(x, P) < λ, an user parameter Topology guarantee is lost Require topological disks around vertices Dey (2014) Computational Topology CCCG / 55
65 Mesh generation A Solution Replace d(x, P) < εlfs(x) with d(x, P) < λ, an user parameter Topology guarantee is lost Require topological disks around vertices Guarantees manifolds Dey (2014) Computational Topology CCCG / 55
66 Mesh generation A Solution Replace d(x, P) < εlfs(x) with d(x, P) < λ, an user parameter Topology guarantee is lost Require topological disks around vertices Guarantees manifolds Dey (2014) Computational Topology CCCG / 55
67 Mesh generation A Solution 1 Initialize points P Σ, compute Del P 2 If some Voronoi edge intersects Σ at x with d(x, P) > εlfs(x), insert x in P and repeat 3 Return Del P Σ Dey (2014) Computational Topology CCCG / 55
68 Mesh generation A Solution 1 Initialize points P Σ, compute Del P 2 If some Voronoi edge intersects Σ at x with d(x, P) > λlfs(x), insert x in P and repeat 3 Return Del P Σ Dey (2014) Computational Topology CCCG / 55
69 Mesh generation A Solution 1 Initialize points P Σ, compute Del P 2 If some Voronoi edge intersects Σ at x with d(x, P) > λlfs(x), insert x in P and repeat 3 If restricted triangles around a vertex p do not form a topological disk, insert furthest x where a dual Voronoi edge of a triangle around p intersects Σ 4 Return Del P Σ Dey (2014) Computational Topology CCCG / 55
70 A Meshing Theorem Mesh generation Theorem Previous algorithm produces output mesh with the following guarantees: Dey (2014) Computational Topology CCCG / 55
71 Mesh generation A Meshing Theorem Theorem Previous algorithm produces output mesh with the following guarantees: 1 Output mesh is always a 2-manifold 2 If λ is sufficiently small, the output mesh satisfies topological and geometric guarantees: Dey (2014) Computational Topology CCCG / 55
72 Mesh generation A Meshing Theorem Theorem Previous algorithm produces output mesh with the following guarantees: 1 Output mesh is always a 2-manifold 2 If λ is sufficiently small, the output mesh satisfies topological and geometric guarantees: 1 It is related to Σ by an isotopy Dey (2014) Computational Topology CCCG / 55
73 A Meshing Theorem Mesh generation Theorem Previous algorithm produces output mesh with the following guarantees: 1 Output mesh is always a 2-manifold 2 If λ is sufficiently small, the output mesh satisfies topological and geometric guarantees: 1 It is related to Σ by an isotopy 2 Each triangle has normal aligning within O(λ) angle to the surface normals Dey (2014) Computational Topology CCCG / 55
74 A Meshing Theorem Mesh generation Theorem Previous algorithm produces output mesh with the following guarantees: 1 Output mesh is always a 2-manifold 2 If λ is sufficiently small, the output mesh satisfies topological and geometric guarantees: 1 It is related to Σ by an isotopy 2 Each triangle has normal aligning within O(λ) angle to the surface normals 3 Hausdorff distance between Σ and Del P Σ is O(λ 2 ) of LFS. Dey (2014) Computational Topology CCCG / 55
75 Topological Data Analysis Data Analysis by Persistent Homology Persistent homology [Edelsbrunner-Letscher-Zomorodian 00], [Zomorodian-Carlsson 02] Dey (2014) Computational Topology CCCG / 55
76 Chain Topological Data Analysis Let K be a finite simplicial complex Dey (2014) Computational Topology CCCG / 55
77 Chain Topological Data Analysis Let K be a finite simplicial complex Simplicial complex Dey (2014) Computational Topology CCCG / 55
78 Chain Topological Data Analysis Let K be a finite simplicial complex Simplicial complex Definition A p-chain in K is a formal sum of p-simplices: c = i the addition in a ring, Z, Z 2, R etc. a i σ i ; sum is Dey (2014) Computational Topology CCCG / 55
79 Chain Topological Data Analysis Let K be a finite simplicial complex 1-chain ab + bc + cd (a i Z 2 ) Definition A p-chain in K is a formal sum of p-simplices: c = i the addition in a ring, Z, Z 2, R etc. a i σ i ; sum is Dey (2014) Computational Topology CCCG / 55
80 Boundary Topological Data Analysis Definition A p-boundary p+1 c of a (p + 1)-chain c is defined as the sum of boundaries of its simplices Dey (2014) Computational Topology CCCG / 55
81 Boundary Topological Data Analysis Definition A p-boundary p+1 c of a (p + 1)-chain c is defined as the sum of boundaries of its simplices Simplicial complex Dey (2014) Computational Topology CCCG / 55
82 Boundary Topological Data Analysis Definition A p-boundary p+1 c of a (p + 1)-chain c is defined as the sum of boundaries of its simplices 2-chain bcd + bde (under Z 2 ) Dey (2014) Computational Topology CCCG / 55
83 Boundary Topological Data Analysis Definition A p-boundary p+1 c of a (p + 1)-chain c is defined as the sum of boundaries of its simplices 1-boundary bc +cd +db+bd +de +eb = bc +cd +de +eb = 2 (bcd +bde) (under Z 2 ) Dey (2014) Computational Topology CCCG / 55
84 Cycles Topological Data Analysis Definition A p-cycle is a p-chain that has an empty boundary Dey (2014) Computational Topology CCCG / 55
85 Cycles Topological Data Analysis Definition A p-cycle is a p-chain that has an empty boundary Simplicial complex Dey (2014) Computational Topology CCCG / 55
86 Cycles Topological Data Analysis Definition A p-cycle is a p-chain that has an empty boundary 1-cycle ab + bc + cd + de + ea (under Z 2 ) Dey (2014) Computational Topology CCCG / 55
87 Cycles Topological Data Analysis Definition A p-cycle is a p-chain that has an empty boundary 1-cycle ab + bc + cd + de + ea (under Z 2 ) Each p-boundary is a p-cycle: p p+1 = 0 Dey (2014) Computational Topology CCCG / 55
88 Groups Topological Data Analysis Definition The p-chain group C p (K) of K is formed by p-chains under addition Dey (2014) Computational Topology CCCG / 55
89 Groups Topological Data Analysis Definition The p-chain group C p (K) of K is formed by p-chains under addition The boundary operator p induces a homomorphism p : C p (K) C p 1 (K) Dey (2014) Computational Topology CCCG / 55
90 Groups Topological Data Analysis Definition The p-chain group C p (K) of K is formed by p-chains under addition The boundary operator p induces a homomorphism p : C p (K) C p 1 (K) Definition The p-cycle group Z p (K) of K is the kernel ker p Dey (2014) Computational Topology CCCG / 55
91 Groups Topological Data Analysis Definition The p-chain group C p (K) of K is formed by p-chains under addition The boundary operator p induces a homomorphism p : C p (K) C p 1 (K) Definition The p-cycle group Z p (K) of K is the kernel ker p Definition The p-boundary group B p (K) of K is the image im p+1 Dey (2014) Computational Topology CCCG / 55
92 Homology Topological Data Analysis Definition The p-dimensional homology group is defined as H p (K) = Z p (K)/B p (K) Dey (2014) Computational Topology CCCG / 55
93 Homology Topological Data Analysis Definition The p-dimensional homology group is defined as H p (K) = Z p (K)/B p (K) Definition Two p-chains c and c are homologous if c = c + p+1 d for some chain d Dey (2014) Computational Topology CCCG / 55
94 Homology Topological Data Analysis Definition The p-dimensional homology group is defined as H p (K) = Z p (K)/B p (K) Definition Two p-chains c and c are homologous if c = c + p+1 d for some chain d (a) trivial (null-homologous) cycle; (b), (c) nontrivial homologous cycles Dey (2014) Computational Topology CCCG / 55
95 Complexes Topological Data Analysis Let P R d be a point set Dey (2014) Computational Topology CCCG / 55
96 Topological Data Analysis Complexes Let P R d be a point set B(p, r) denotes an open d-ball centered at p with radius r Dey (2014) Computational Topology CCCG / 55
97 Topological Data Analysis Complexes Let P R d be a point set B(p, r) denotes an open d-ball centered at p with radius r Definition The Čech complex Cr (P) is a simplicial complex where a simplex σ C r (P) iff Vert(σ) P and p Vert(σ) B(p, r/2) 0 Dey (2014) Computational Topology CCCG / 55
98 Topological Data Analysis Complexes Let P R d be a point set B(p, r) denotes an open d-ball centered at p with radius r Definition The Čech complex Cr (P) is a simplicial complex where a simplex σ C r (P) iff Vert(σ) P and p Vert(σ) B(p, r/2) 0 Definition The Rips complex R r (P) is a simplicial complex where a simplex σ R r (P) iff Vert(σ) are within pairwise Euclidean distance of r Dey (2014) Computational Topology CCCG / 55
99 Topological Data Analysis Complexes Let P R d be a point set B(p, r) denotes an open d-ball centered at p with radius r Definition The Čech complex C r (P) is a simplicial complex where a simplex σ C r (P) iff Vert(σ) P and p Vert(σ) B(p, r/2) 0 Definition The Rips complex R r (P) is a simplicial complex where a simplex σ R r (P) iff Vert(σ) are within pairwise Euclidean distance of r Proposition For any finite set P R d and any r 0, C r (P) R r (P) C 2r (P) Dey (2014) Computational Topology CCCG / 55
100 Point set P Topological Data Analysis Dey (2014) Computational Topology CCCG / 55
101 Topological Data Analysis Balls B(p, r/2) for p P Dey (2014) Computational Topology CCCG / 55
102 Topological Data Analysis Čech complex C r (P) Dey (2014) Computational Topology CCCG / 55
103 Topological Data Analysis Rips complex R r (P) Dey (2014) Computational Topology CCCG / 55
104 Topological Data Analysis Topological persistence r(x) = d(x, P): distance to point cloud P Sublevel sets r 1 [0, a] are union of balls Evolution of the sublevel sets with increasing a left hole persists longer Persistent homology quantizes this idea Dey (2014) Computational Topology CCCG / 55
105 Topological Data Analysis Persistent Homology f : T R; T a = f 1 (, a], the sublevel set T a T b for a b provides inclusion map ι : T a T b Induced map ι : H p (T a ) H p (T b ) giving the sequence 0 H p (T a1 ) H p (T a2 ) H p (T an ) H p (T) Persistent homology classes: Image of f ij p : H p (T ai ) H p (T aj ) Dey (2014) Computational Topology CCCG / 55
106 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Dey (2014) Computational Topology CCCG / 55
107 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Union of balls with its nerve Čech complex Dey (2014) Computational Topology CCCG / 55
108 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Union of balls with its nerve Čech complex Dey (2014) Computational Topology CCCG / 55
109 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Union of balls with its nerve Čech complex (a) (b) (c) (d) Dey (2014) Computational Topology CCCG / 55
110 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Union of balls with its nerve Čech complex Evolution of sublevel sets becomes Filtration: (a) (b) = K 0 K 1 K n = K 0 H p (K 1 ) H p (K n ) = H p (K). (c) (d) Dey (2014) Computational Topology CCCG / 55
111 Topological Data Analysis Continuous to Discrete Replace T with a simplicial complex K := K(T) Union of balls with its nerve Čech complex Evolution of sublevel sets becomes Filtration: (a) (b) = K 0 K 1 K n = K 0 H p (K 1 ) H p (K n ) = H p (K). Birth and Death of homology classes (c) (d) Dey (2014) Computational Topology CCCG / 55
112 Bar Codes Topological Data Analysis birth-death and bar codes Dey (2014) Computational Topology CCCG / 55
113 Topological Data Analysis Persistence Diagram a bar [a, b] is represented as a point in the plane Dey (2014) Computational Topology CCCG / 55
114 Topological Data Analysis Persistence Diagram a bar [a, b] is represented as a point in the plane incorporate the diagonal in the diagram Dgm p (f ) Dey (2014) Computational Topology CCCG / 55
115 Topological Data Analysis Persistence Diagram a bar [a, b] is represented as a point in the plane incorporate the diagonal in the diagram Dgm p (f ) Dey (2014) Computational Topology CCCG / 55
116 Topological Data Analysis Persistence Diagram a bar [a, b] is represented as a point in the plane incorporate the diagonal in the diagram Dgm p (f ) Death Birth Dey (2014) Computational Topology CCCG / 55
117 Topological Data Analysis Stability of Persistence Diagram Bottleneck distance (C all bijections) d B (Dgm p (f ), Dgm p (g)) := inf sup x c(x) c C x Dgm p (f ) Theorem (Cohen-Steiner,Edelsbrunner,Harer 06) d B (Dgm p (f ), Dgm p (g)) f g Dey (2014) Computational Topology CCCG / 55
118 Back to Point Data Topological Data Analysis d T be the distance function from the space T. Dey (2014) Computational Topology CCCG / 55
119 Topological Data Analysis Back to Point Data d T be the distance function from the space T. d P be the distance function from sample P Dey (2014) Computational Topology CCCG / 55
120 Topological Data Analysis Back to Point Data Death d T be the distance function from the space T. d P be the distance function from sample P Birth Dey (2014) Computational Topology CCCG / 55
121 Topological Data Analysis Back to Point Data Death d T be the distance function from the space T. d P be the distance function from sample P d T d P d(t, P) = ε Birth Dey (2014) Computational Topology CCCG / 55
122 Topological Data Analysis Back to Point Data Death d T be the distance function from the space T. d P be the distance function from sample P d T d P d(t, P) = ε d B (Dgm p (d T ), Dgm (d P )) ε Birth Dey (2014) Computational Topology CCCG / 55
123 Topological Data Analysis Back to Point Data Death d T be the distance function from the space T. d P be the distance function from sample P d T d P d(t, P) = ε d B (Dgm p (d T ), Dgm (d P )) ε Birth Dey (2014) Computational Topology CCCG / 55
124 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Dey (2014) Computational Topology CCCG / 55
125 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Dey (2014) Computational Topology CCCG / 55
126 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Sparsified Rips complex [Sheehy 12] Dey (2014) Computational Topology CCCG / 55
127 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Sparsified Rips complex [Sheehy 12] Graph Induced Complex (GIC) [D.-Fang-Wang 13] R r (P ) Dey (2014) Computational Topology CCCG / 55
128 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Sparsified Rips complex [Sheehy 12] Graph Induced Complex (GIC) [D.-Fang-Wang 13] R r (P ) Q Dey (2014) Computational Topology CCCG / 55
129 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Sparsified Rips complex [Sheehy 12] Graph Induced Complex (GIC) [D.-Fang-Wang 13] R r (P ) R r (Q) Dey (2014) Computational Topology CCCG / 55
130 Topological Data Analysis Issues: Choice of Complexes Čech complexes are difficult to compute Most literature considers Rips complexes, but they are Huge Sparsified Rips complex [Sheehy 12] Graph Induced Complex (GIC) [D.-Fang-Wang 13] R r (P ) R r (Q) G r (P, Q) Dey (2014) Computational Topology CCCG / 55
131 Topological Data Analysis Issues: Filtration Maps Zigzag inclusions[carlsson-de Silva-Morozov 09] K 1 K 2 K 3... K n Dey (2014) Computational Topology CCCG / 55
132 Topological Data Analysis Issues: Filtration Maps Zigzag inclusions[carlsson-de Silva-Morozov 09] K 1 K 2 K 3... K n Dey (2014) Computational Topology CCCG / 55
133 Topological Data Analysis Issues: Filtration Maps Zigzag inclusions[carlsson-de Silva-Morozov 09] K 1 K 2 K 3... K n Simplicial maps instead of inclusions [D.-Fan-Wang 14] K f 1 f K 2 f 1 K2 n Kn = K Dey (2014) Computational Topology CCCG / 55
134 Topological Data Analysis Issues: Filtration Maps Zigzag inclusions[carlsson-de Silva-Morozov 09] K 1 K 2 K 3... K n Simplicial maps instead of inclusions [D.-Fan-Wang 14] K f 1 f K 2 f 1 K2 n Kn = K Efficient algorithm for Zigzag simplicial maps? Dey (2014) Computational Topology CCCG / 55
135 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Dey (2014) Computational Topology CCCG / 55
136 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Dey (2014) Computational Topology CCCG / 55
137 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Scalar field data analysis [CGOS 09] Dey (2014) Computational Topology CCCG / 55
138 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Scalar field data analysis [CGOS 09] Computing topological structures Dey (2014) Computational Topology CCCG / 55
139 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Scalar field data analysis [CGOS 09] Computing topological structures Reeb graphs [C-MEHNP 03, D.-Wang 11] Dey (2014) Computational Topology CCCG / 55
140 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Scalar field data analysis [CGOS 09] Computing topological structures Reeb graphs [C-MEHNP 03, D.-Wang 11] Homology cycles Dey (2014) Computational Topology CCCG / 55
141 Further Problems Further Issues/Problems What about stabilty when domains aren t same. Interleaving [Chazal-Oudot 08], [CCGGO08] Scalar field data analysis [CGOS 09] Computing topological structures Reeb graphs [C-MEHNP 03, D.-Wang 11] Homology cycles Dey (2014) Computational Topology CCCG / 55
142 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis Dey (2014) Computational Topology CCCG / 55
143 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis Dey (2014) Computational Topology CCCG / 55
144 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis Dey (2014) Computational Topology CCCG / 55
145 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis First solution for surfaces: Erickson-Whittlesey [SODA05] Dey (2014) Computational Topology CCCG / 55
146 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis First solution for surfaces: Erickson-Whittlesey [SODA05] General problem NP-hard: Chen-Freedman [SODA10] Dey (2014) Computational Topology CCCG / 55
147 Further Problems Optimal Homology Basis Problem Compute an optimal set of cycles forming a basis First solution for surfaces: Erickson-Whittlesey [SODA05] General problem NP-hard: Chen-Freedman [SODA10] H 1 basis for simplicial complexes: D.-Sun-Wang [SoCG10] Dey (2014) Computational Topology CCCG / 55
148 Optimal Localization Further Problems Compute an optimal cycle in a given class. Dey (2014) Computational Topology CCCG / 55
149 Optimal Localization Further Problems Compute an optimal cycle in a given class. Dey (2014) Computational Topology CCCG / 55
150 Optimal Localization Further Problems Compute an optimal cycle in a given class. Dey (2014) Computational Topology CCCG / 55
151 Optimal Localization Further Problems Compute an optimal cycle in a given class. Dey (2014) Computational Topology CCCG / 55
152 Optimal Localization Further Problems Compute an optimal cycle in a given class. Surfaces: Colin de Verdière-Lazarus [DCG05], Colin de Verdière-Erickson [SODA06], Chambers-Erickson-Nayyeri [SoCG09] Dey (2014) Computational Topology CCCG / 55
153 Optimal Localization Further Problems Compute an optimal cycle in a given class. Surfaces: Colin de Verdière-Lazarus [DCG05], Colin de Verdière-Erickson [SODA06], Chambers-Erickson-Nayyeri [SoCG09] General problem NP-hard: Chen-Freedman [SODA10] Dey (2014) Computational Topology CCCG / 55
154 Optimal Localization Further Problems Compute an optimal cycle in a given class. Surfaces: Colin de Verdière-Lazarus [DCG05], Colin de Verdière-Erickson [SODA06], Chambers-Erickson-Nayyeri [SoCG09] General problem NP-hard: Chen-Freedman [SODA10] Special cases: Dey-Hirani-Krishnamoorthy [STOC10] Dey (2014) Computational Topology CCCG / 55
155 Conclusions Further Problems Reconstructions : Dey (2014) Computational Topology CCCG / 55
156 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open Dey (2014) Computational Topology CCCG / 55
157 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Dey (2014) Computational Topology CCCG / 55
158 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : Dey (2014) Computational Topology CCCG / 55
159 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : piecewise smooth surfaces, complexes Dey (2014) Computational Topology CCCG / 55
160 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : piecewise smooth surfaces, complexes Dey (2014) Computational Topology CCCG / 55
161 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : piecewise smooth surfaces, complexes Data Analysis : Dey (2014) Computational Topology CCCG / 55
162 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : piecewise smooth surfaces, complexes Data Analysis : functions on spaces Dey (2014) Computational Topology CCCG / 55
163 Conclusions Further Problems Reconstructions : non-smooth surfaces remain open high dimensions still not satisfactory Mesh Generation : piecewise smooth surfaces, complexes Data Analysis : functions on spaces connecting to data mining, machine learning. Dey (2014) Computational Topology CCCG / 55
164 Thank Thank You Dey (2014) Computational Topology CCCG / 55
Estimating Geometry and Topology from Voronoi Diagrams
Estimating Geometry and Topology from Voronoi Diagrams Tamal K. Dey The Ohio State University Tamal Dey OSU Chicago 07 p.1/44 Voronoi diagrams Tamal Dey OSU Chicago 07 p.2/44 Voronoi diagrams Tamal Dey
More informationA Multicover Nerve for Geometric Inference. Don Sheehy INRIA Saclay, France
A Multicover Nerve for Geometric Inference Don Sheehy INRIA Saclay, France Computational geometers use topology to certify geometric constructions. Computational geometers use topology to certify geometric
More informationTopological Data Analysis - I. Afra Zomorodian Department of Computer Science Dartmouth College
Topological Data Analysis - I Afra Zomorodian Department of Computer Science Dartmouth College September 3, 2007 1 Acquisition Vision: Images (2D) GIS: Terrains (3D) Graphics: Surfaces (3D) Medicine: MRI
More informationTopological Inference via Meshing
Topological Inference via Meshing Benoit Hudson, Gary Miller, Steve Oudot and Don Sheehy SoCG 2010 Mesh Generation and Persistent Homology The Problem Input: Points in Euclidean space sampled from some
More informationA Geometric Perspective on Sparse Filtrations
CCCG 2015, Kingston, Ontario, August 10 12, 2015 A Geometric Perspective on Sparse Filtrations Nicholas J. Cavanna Mahmoodreza Jahanseir Donald R. Sheehy Abstract We present a geometric perspective on
More informationSimplicial Complexes: Second Lecture
Simplicial Complexes: Second Lecture 4 Nov, 2010 1 Overview Today we have two main goals: Prove that every continuous map between triangulable spaces can be approximated by a simplicial map. To do this,
More informationTopological Inference via Meshing
Topological Inference via Meshing Don Sheehy Theory Lunch Joint work with Benoit Hudson, Gary Miller, and Steve Oudot Computer Scientists want to know the shape of data. Clustering Principal Component
More information66 III Complexes. R p (r) }.
66 III Complexes III.4 Alpha Complexes In this section, we use a radius constraint to introduce a family of subcomplexes of the Delaunay complex. These complexes are similar to the Čech complexes but differ
More informationPersistence stability for geometric complexes
Lake Tahoe - December, 8 2012 NIPS workshop on Algebraic Topology and Machine Learning Persistence stability for geometric complexes Frédéric Chazal Joint work with Vin de Silva and Steve Oudot (+ on-going
More information4.2 Simplicial Homology Groups
4.2. SIMPLICIAL HOMOLOGY GROUPS 93 4.2 Simplicial Homology Groups 4.2.1 Simplicial Complexes Let p 0, p 1,... p k be k + 1 points in R n, with k n. We identify points in R n with the vectors that point
More informationA fast and robust algorithm to count topologically persistent holes in noisy clouds
A fast and robust algorithm to count topologically persistent holes in noisy clouds Vitaliy Kurlin Durham University Department of Mathematical Sciences, Durham, DH1 3LE, United Kingdom vitaliy.kurlin@gmail.com,
More informationHomology and Persistent Homology Bootcamp Notes 2017
Homology and Persistent Homology Bootcamp Notes Summer@ICERM 2017 Melissa McGuirl September 19, 2017 1 Preface This bootcamp is intended to be a review session for the material covered in Week 1 of Summer@ICERM
More informationLecture 1 Discrete Geometric Structures
Lecture 1 Discrete Geometric Structures Jean-Daniel Boissonnat Winter School on Computational Geometry and Topology University of Nice Sophia Antipolis January 23-27, 2017 Computational Geometry and Topology
More informationGeometric Inference. 1 Introduction. Frédéric Chazal 1 and David Cohen-Steiner 2
Geometric Inference Frédéric Chazal 1 and David Cohen-Steiner 2 1 INRIA Futurs, 4 rue Jacques Monod - 91893 Orsay Cedex France frederic.chazal@inria.fr 2 INRIA Sophia-Antipolis, 2004 route des lucioles
More informationData Analysis, Persistent homology and Computational Morse-Novikov theory
Data Analysis, Persistent homology and Computational Morse-Novikov theory Dan Burghelea Department of mathematics Ohio State University, Columbuus, OH Bowling Green, November 2013 AMN theory (computational
More informationPERSISTENT HOMOLOGY OF FINITE TOPOLOGICAL SPACES
PERSISTENT HOMOLOGY OF FINITE TOPOLOGICAL SPACES HANEY MAXWELL Abstract. We introduce homology and finite topological spaces. From the basis of that introduction, persistent homology is applied to finite
More informationA Multicover Nerve for Geometric Inference
A Multicover Nerve for Geometric Inference Donald R. Sheehy Abstract We show that filtering the barycentric decomposition of a Čech complex by the cardinality of the vertices captures precisely the topology
More informationAn introduction to Topological Data Analysis through persistent homology: Intro and geometric inference
Sophia-Antipolis, January 2016 Winter School An introduction to Topological Data Analysis through persistent homology: Intro and geometric inference Frédéric Chazal INRIA Saclay - Ile-de-France frederic.chazal@inria.fr
More informationGeometric Inference. 1 Introduction. Frédéric Chazal 1 and David Cohen-Steiner 2
Geometric Inference Frédéric Chazal 1 and David Cohen-Steiner 2 1 INRIA Futurs, 4 rue Jacques Monod - 91893 Orsay Cedex France frederic.chazal@inria.fr 2 INRIA Sophia-Antipolis, 2004 route des lucioles
More informationCorrectness. The Powercrust Algorithm for Surface Reconstruction. Correctness. Correctness. Delaunay Triangulation. Tools - Voronoi Diagram
Correctness The Powercrust Algorithm for Surface Reconstruction Nina Amenta Sunghee Choi Ravi Kolluri University of Texas at Austin Boundary of a solid Close to original surface Homeomorphic to original
More information4. Simplicial Complexes and Simplicial Homology
MATH41071/MATH61071 Algebraic topology Autumn Semester 2017 2018 4. Simplicial Complexes and Simplicial Homology Geometric simplicial complexes 4.1 Definition. A finite subset { v 0, v 1,..., v r } R n
More informationClustering algorithms and introduction to persistent homology
Foundations of Geometric Methods in Data Analysis 2017-18 Clustering algorithms and introduction to persistent homology Frédéric Chazal INRIA Saclay - Ile-de-France frederic.chazal@inria.fr Introduction
More informationManifolds. Chapter X. 44. Locally Euclidean Spaces
Chapter X Manifolds 44. Locally Euclidean Spaces 44 1. Definition of Locally Euclidean Space Let n be a non-negative integer. A topological space X is called a locally Euclidean space of dimension n if
More informationThe Medial Axis of the Union of Inner Voronoi Balls in the Plane
The Medial Axis of the Union of Inner Voronoi Balls in the Plane Joachim Giesen a, Balint Miklos b,, Mark Pauly b a Max-Planck Institut für Informatik, Saarbrücken, Germany b Applied Geometry Group, ETH
More informationNew Results on the Stability of Persistence Diagrams
New Results on the Stability of Persistence Diagrams Primoz Skraba Jozef Stefan Institute TAGS - Linking Topology to Algebraic Geometry and Statistics joint work with Katharine Turner and D. Yogeshwaran
More informationarxiv: v1 [math.st] 11 Oct 2017
An introduction to Topological Data Analysis: fundamental and practical aspects for data scientists Frédéric Chazal and Bertrand Michel October 12, 2017 arxiv:1710.04019v1 [math.st] 11 Oct 2017 Abstract
More informationGeometric Data. Goal: describe the structure of the geometry underlying the data, for interpretation or summary
Geometric Data - ma recherche s inscrit dans le contexte de l analyse exploratoire des donnees, dont l obj Input: point cloud equipped with a metric or (dis-)similarity measure data point image/patch,
More informationLecture notes for Topology MMA100
Lecture notes for Topology MMA100 J A S, S-11 1 Simplicial Complexes 1.1 Affine independence A collection of points v 0, v 1,..., v n in some Euclidean space R N are affinely independent if the (affine
More informationA THINNING ALGORITHM FOR TOPOLOGICALLY CORRECT 3D SURFACE RECONSTRUCTION
A THINNING ALGORITHM FOR TOPOLOGICALLY CORRECT 3D SURFACE RECONSTRUCTION Leonid Tcherniavski Department of Informatics University of Hamburg tcherniavski@informatik.uni-hamburg.de ABSTRACT The existing
More informationSurfaces, meshes, and topology
Surfaces from Point Samples Surfaces, meshes, and topology A surface is a 2-manifold embedded in 3- dimensional Euclidean space Such surfaces are often approximated by triangle meshes 2 1 Triangle mesh
More informationClique homological simplification problem is NP-hard
Clique homological simplification problem is NP-hard NingNing Peng Department of Mathematics, Wuhan University of Technology pnn@whut.edu.cn Foundations of Mathematics September 9, 2017 Outline 1 Homology
More informationReconstruction of Filament Structure
Reconstruction of Filament Structure Ruqi HUANG INRIA-Geometrica Joint work with Frédéric CHAZAL and Jian SUN 27/10/2014 Outline 1 Problem Statement Characterization of Dataset Formulation 2 Our Approaches
More informationMesh Generation. Jean-Daniel Boissonnat DataShape, INRIA
Mesh Generation Jean-Daniel Boissonnat DataShape, INRIA http://www-sop.inria.fr/datashape Algorithmic Geometry Mesh generation J-D. Boissonnat 1 / 27 Meshing surfaces and 3D domains visualization and graphics
More informationObserving Information: Applied Computational Topology.
Observing Information: Applied Computational Topology. Bangor University, and NUI Galway April 21, 2008 What is the geometric information that can be gleaned from a data cloud? Some ideas either already
More informationLecture 7: Jan 31, Some definitions related to Simplical Complex. 7.2 Topological Equivalence and Homeomorphism
CS 6170 Computational Topology: Topological Data Analysis University of Utah Spring 2017 School of Computing Lecture 7: Jan 31, 2017 Lecturer: Prof. Bei Wang Scribe: Avani Sharma,
More informationTopological estimation using witness complexes. Vin de Silva, Stanford University
Topological estimation using witness complexes, Acknowledgements Gunnar Carlsson (Mathematics, Stanford) principal collaborator Afra Zomorodian (CS/Robotics, Stanford) persistent homology software Josh
More informationDiffusion Maps and Topological Data Analysis
Diffusion Maps and Topological Data Analysis Melissa R. McGuirl McGuirl (Brown University) Diffusion Maps and Topological Data Analysis 1 / 19 Introduction OVERVIEW Topological Data Analysis The use of
More informationTHE BASIC THEORY OF PERSISTENT HOMOLOGY
THE BASIC THEORY OF PERSISTENT HOMOLOGY KAIRUI GLEN WANG Abstract Persistent homology has widespread applications in computer vision and image analysis This paper first motivates the use of persistent
More informationDetection and approximation of linear structures in metric spaces
Stanford - July 12, 2012 MMDS 2012 Detection and approximation of linear structures in metric spaces Frédéric Chazal Geometrica group INRIA Saclay Joint work with M. Aanjaneya, D. Chen, M. Glisse, L. Guibas,
More informationGeometric structures on 2-orbifolds
Geometric structures on 2-orbifolds Section 1: Manifolds and differentiable structures S. Choi Department of Mathematical Science KAIST, Daejeon, South Korea 2010 Fall, Lectures at KAIST S. Choi (KAIST)
More informationA Fast and Simple Surface Reconstruction Algorithm
A Fast and Simple Surface Reconstruction Algorithm Siu-Wing Cheng Jiongxin Jin Man-Kit Lau Abstract We present an algorithm for surface reconstruction from a point cloud. It runs in O(n log n) time, where
More informationTutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass
Tutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass University of California, Davis, USA http://www.cs.ucdavis.edu/~koehl/ims2017/ What is a shape? A shape is a 2-manifold with a Riemannian
More informationTopological Data Analysis
Proceedings of Symposia in Applied Mathematics Topological Data Analysis Afra Zomorodian Abstract. Scientific data is often in the form of a finite set of noisy points, sampled from an unknown space, and
More informationLecture 0: Reivew of some basic material
Lecture 0: Reivew of some basic material September 12, 2018 1 Background material on the homotopy category We begin with the topological category TOP, whose objects are topological spaces and whose morphisms
More informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
More information4. Definition: topological space, open set, topology, trivial topology, discrete topology.
Topology Summary Note to the reader. If a statement is marked with [Not proved in the lecture], then the statement was stated but not proved in the lecture. Of course, you don t need to know the proof.
More informationComplex conforming Delaunay triangulation
. RESEARCH PAPERS. SCIENCE CHINA Information Sciences June 2010 Vol. 53 No. 6: 1130 1140 doi: 10.1007/s11432-010-0097-6 Complex conforming Delaunay triangulation MENG XianHai 1, LI JiGang 1, YANG Qin 1,
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationEffective Computational Geometry for Curves and Surfaces. Chapter 7 Computational Topology: An Introduction
Effective Computational Geometry for Curves and Surfaces Chapter 7 Computational Topology: An Introduction Günter Rote and Gert Vegter We give an introduction to combinatorial topology, with an emphasis
More informationOpen Problems Column Edited by William Gasarch
Open Problems Column Edited by William Gasarch 1 This Issues Column! This issue s Open Problem Column is by Brittany Terese Fasy and Bei Wang and is on Open Problems in Computational Topology; however,
More informationAnalysis of high dimensional data via Topology. Louis Xiang. Oak Ridge National Laboratory. Oak Ridge, Tennessee
Analysis of high dimensional data via Topology Louis Xiang Oak Ridge National Laboratory Oak Ridge, Tennessee Contents Abstract iii 1 Overview 1 2 Data Set 1 3 Simplicial Complex 5 4 Computation of homology
More informationPersistent Homology and Nested Dissection
Persistent Homology and Nested Dissection Don Sheehy University of Connecticut joint work with Michael Kerber and Primoz Skraba A Topological Data Analysis Pipeline A Topological Data Analysis Pipeline
More informationApproximating Polygonal Objects by Deformable Smooth Surfaces
Approximating Polygonal Objects by Deformable Smooth Surfaces Ho-lun Cheng and Tony Tan School of Computing, National University of Singapore hcheng,tantony@comp.nus.edu.sg Abstract. We propose a method
More information35 CURVE AND SURFACE RECONSTRUCTION
35 CURVE AND SURFACE RECONSTRUCTION Tamal K. Dey INTRODUCTION The problem of reconstructing a shape from its sample appears in many scientific and engineering applications. Because of the variety in shapes
More informationDiscrete Surfaces. David Gu. Tsinghua University. Tsinghua University. 1 Mathematics Science Center
Discrete Surfaces 1 1 Mathematics Science Center Tsinghua University Tsinghua University Discrete Surface Discrete Surfaces Acquired using 3D scanner. Discrete Surfaces Our group has developed high speed
More informationLecture 5: Simplicial Complex
Lecture 5: Simplicial Complex 2-Manifolds, Simplex and Simplicial Complex Scribed by: Lei Wang First part of this lecture finishes 2-Manifolds. Rest part of this lecture talks about simplicial complex.
More informationCell-Like Maps (Lecture 5)
Cell-Like Maps (Lecture 5) September 15, 2014 In the last two lectures, we discussed the notion of a simple homotopy equivalences between finite CW complexes. A priori, the question of whether or not a
More informationComputational Geometry
Computational Geometry 600.658 Convexity A set S is convex if for any two points p, q S the line segment pq S. S p S q Not convex Convex? Convexity A set S is convex if it is the intersection of (possibly
More informationStability and Statistical properties of topological information inferred from metric data
Beijing, June 25, 2014 ICML workshop: Topological Methods for Machine Learning Stability and Statistical properties of topological information inferred from metric data Frédéric Chazal INRIA Saclay Introduction
More informationComputing the Betti Numbers of Arrangements. Saugata Basu School of Mathematics & College of Computing Georgia Institute of Technology.
1 Computing the Betti Numbers of Arrangements Saugata Basu School of Mathematics & College of Computing Georgia Institute of Technology. 2 Arrangements in Computational Geometry An arrangement in R k is
More informationQUANTIFYING HOMOLOGY CLASSES
QUANTIFYING HOMOLOGY CLASSES CHAO CHEN 1 AND DANIEL FREEDMAN 1 1 Rensselaer Polytechnic Institute, 110 8th street, Troy, NY 12180, U.S.A. E-mail address, {C. Chen, D. Freedman}: {chenc3,freedman}@cs.rpi.edu
More informationReeb Spaces and the Robustness of Preimages
Reeb Spaces and the Robustness of Preimages by Amit Patel Department of Computer Science Duke University Ph.D. Dissertation 2010 Copyright c 2010 by Amit Patel All rights reserved except the rights granted
More informationStratified Structure of Laplacian Eigenmaps Embedding
Stratified Structure of Laplacian Eigenmaps Embedding Abstract We construct a locality preserving weight matrix for Laplacian eigenmaps algorithm used in dimension reduction. Our point cloud data is sampled
More informationKurt Mehlhorn, MPI für Informatik. Curve and Surface Reconstruction p.1/25
Curve and Surface Reconstruction Kurt Mehlhorn MPI für Informatik Curve and Surface Reconstruction p.1/25 Curve Reconstruction: An Example probably, you see more than a set of points Curve and Surface
More informationQuadratic and cubic b-splines by generalizing higher-order voronoi diagrams
Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
More informationReflection groups 4. Mike Davis. May 19, Sao Paulo
Reflection groups 4 Mike Davis Sao Paulo May 19, 2014 https://people.math.osu.edu/davis.12/slides.html 1 2 Exotic fundamental gps Nonsmoothable aspherical manifolds 3 Let (W, S) be a Coxeter system. S
More informationHowever, this is not always true! For example, this fails if both A and B are closed and unbounded (find an example).
98 CHAPTER 3. PROPERTIES OF CONVEX SETS: A GLIMPSE 3.2 Separation Theorems It seems intuitively rather obvious that if A and B are two nonempty disjoint convex sets in A 2, then there is a line, H, separating
More informationAnalyse Topologique des Données
Collège de France 31 mai 2017 Analyse Topologique des Données Frédéric Chazal DataShape team INRIA Saclay - Ile-de-France frederic.chazal@inria.fr Introduction [Shape database] [Galaxies data] [Scanned
More informationPlanar Graphs. 1 Graphs and maps. 1.1 Planarity and duality
Planar Graphs In the first half of this book, we consider mostly planar graphs and their geometric representations, mostly in the plane. We start with a survey of basic results on planar graphs. This chapter
More informationThe Cyclic Cycle Complex of a Surface
The Cyclic Cycle Complex of a Surface Allen Hatcher A recent paper [BBM] by Bestvina, Bux, and Margalit contains a construction of a cell complex that gives a combinatorial model for the collection of
More informationComplexity of Delaunay Triangulation for Points on Lower-dimensional Polyhedra.
published in Proc. 18th ACM-SIAM Sympos. Discrete Algorithms (SODA) 2007, p.1106 1113. Complexity of Delaunay Triangulation for Points on Lower-dimensional Polyhedra. Nina Amenta UC Davis Dominique Attali
More informationOther Voronoi/Delaunay Structures
Other Voronoi/Delaunay Structures Overview Alpha hulls (a subset of Delaunay graph) Extension of Voronoi Diagrams Convex Hull What is it good for? The bounding region of a point set Not so good for describing
More informationComplexity of Delaunay triangulation for points on lower-dimensional polyhedra
Complexity of Delaunay triangulation for points on lower-dimensional polyhedra Nina Amenta UC Davis Dominique Attali CNRS-LIS Grenoble Olivier Devillers INRIA Sophia-Antipolis Abstract We show that the
More informationOn Computing Handle and Tunnel Loops
On Computing Handle and Tunnel Loops Tamal K. Dey Kuiyu Li Jian Sun Department of Computer Science and Engineering The Ohio State University Columbus, OH, 43210, USA Abstract Many applications seek to
More informationOn the Topology of Finite Metric Spaces
On the Topology of Finite Metric Spaces Meeting in Honor of Tom Goodwillie Dubrovnik Gunnar Carlsson, Stanford University June 27, 2014 Data has shape The shape matters Shape of Data Regression Shape of
More informationFinal Exam, F11PE Solutions, Topology, Autumn 2011
Final Exam, F11PE Solutions, Topology, Autumn 2011 Question 1 (i) Given a metric space (X, d), define what it means for a set to be open in the associated metric topology. Solution: A set U X is open if,
More informationA roadmap for the computation of persistent homology
A roadmap for the computation of persistent homology Nina Otter Mathematical Institute, University of Oxford (joint with Mason A. Porter, Ulrike Tillmann, Peter Grindrod, Heather A. Harrington) 2nd School
More informationOn the Locality of Extracting a 2-Manifold in IR 3
On the Locality of Extracting a 2-Manifold in IR 3 Daniel Dumitriu 1, Stefan Funke 2, Martin Kutz, and Nikola Milosavljević 3 1 Max-Planck-Institut für Informatik Campus E 1.4, 66123 Saarbrücken, Germany
More informationMath 734 Aug 22, Differential Geometry Fall 2002, USC
Math 734 Aug 22, 2002 1 Differential Geometry Fall 2002, USC Lecture Notes 1 1 Topological Manifolds The basic objects of study in this class are manifolds. Roughly speaking, these are objects which locally
More informationThe Gudhi Library: Simplicial Complexes and Persistent Homology
The Gudhi Library: Simplicial Complexes and Persistent Homology Clément Maria, Jean-Daniel Boissonnat, Marc Glisse, Mariette Yvinec To cite this version: Clément Maria, Jean-Daniel Boissonnat, Marc Glisse,
More informationSURFACE RECONSTRUCTION AND SIMPLIFICATION
NATIONAL AND KAPODISTRIAN UNIVERSITY OF ATHENS DEPARTMENT OF INFORMATICS AND TELECOMMUNICATIOS SECTOR OF THEORETICAL INFORMATICS SURFACE RECONSTRUCTION AND SIMPLIFICATION Alexander Agathos Master Thesis
More informationA Practical Guide to Persistent Homology
A Practical Guide to Persistent Homology Dmitriy Morozov Lawrence Berkeley National Lab A Practical Guide to Persistent Code snippets available at: http://hg.mrzv.org/dionysus-tutorial Homology (Dionysus
More informationtopological data analysis and stochastic topology yuliy baryshnikov waikiki, march 2013
topological data analysis and stochastic topology yuliy baryshnikov waikiki, march 2013 Promise of topological data analysis: extract the structure from the data. In: point clouds Out: hidden structure
More informationProvably Good Moving Least Squares
Provably Good Moving Least Squares Ravikrishna Kolluri Computer Science Division University of California at Berkeley 1 Problem Definition Given a set of samples on a closed surface build a representation
More informationPart 3. Topological Manifolds
Part 3 Topological Manifolds This part is devoted to study of the most important topological spaces, the spaces which provide a scene for most of geometric branches in mathematics such as Differential
More informationSOCIAL CHOICE AMONG COMPLEX OBJECTS:MATHEMATICAL TOOLS
SOCIAL CHOICE AMONG COMPLEX OBJECTS:MATHEMATICAL TOOLS SIMONA SETTEPANELLA Abstract. Here the reader can find some basic definitions and notations in order to better understand the model for social choise
More informationMeshing Volumes with Curved Boundaries
Engineering with Computers manuscript No. (will be inserted by the editor) Meshing Volumes with Curved Boundaries Steve Oudot 1, Laurent Rineau 2, Mariette Yvinec 3 1 INRIA, BP 93 06902 Sophia Antipolis,
More informationSurface reconstruction based on a dynamical system
EUROGRAPHICS 2002 / G. Drettakis and H.-P. Seidel (Guest Editors) Volume 21 (2002), Number 3 Surface reconstruction based on a dynamical system N.N. Abstract We present an efficient algorithm that computes
More informationSimplicial Objects and Homotopy Groups. J. Wu
Simplicial Objects and Homotopy Groups J. Wu Department of Mathematics, National University of Singapore, Singapore E-mail address: matwuj@nus.edu.sg URL: www.math.nus.edu.sg/~matwujie Partially supported
More informationSurfaces Beyond Classification
Chapter XII Surfaces Beyond Classification In most of the textbooks which present topological classification of compact surfaces the classification is the top result. However the topology of 2- manifolds
More informationSIMPLICIAL ENERGY AND SIMPLICIAL HARMONIC MAPS
SIMPLICIAL ENERGY AND SIMPLICIAL HARMONIC MAPS JOEL HASS AND PETER SCOTT Abstract. We introduce a combinatorial energy for maps of triangulated surfaces with simplicial metrics and analyze the existence
More informationShortest Non-trivial Cycles in Directed and Undirected Surface Graphs
Shortest Non-trivial Cycles in Directed and Undirected Surface Graphs Kyle Fox University of Illinois at Urbana-Champaign Surfaces 2-manifolds (with boundary) genus g: max # of disjoint simple cycles whose
More information25 HIGH-DIMENSIONAL TOPOLOGICAL DATA ANALYSIS
25 HIGH-DIMENSIONAL TOPOLOGICAL DATA ANALYSIS Frédéric Chazal INTRODUCTION Modern data often come as point clouds embedded in high-dimensional Euclidean spaces, or possibly more general metric spaces.
More informationHardness Results for Homology Localization
Hardness Results for Homology Localization Chao Chen Institute of Science and Technology Austria Vienna University of Technology, Austria Email: ChaoChen79@gmail.com Daniel Freedman Hewlett-Packard Laboratories,
More informationCollars and Intestines: Practical Conforming Delaunay Refinement
Collars and Intestines: Practical Conforming Delaunay Refinement Alexander Rand and Noel Walkington Carnegie Mellon University Summary. While several existing Delaunay refinement algorithms allow acute
More informationHOMOLOGY. Equivalently,
HOMOLOGY JOHN ROGNES 1. Homology 1.1. The Euler characteristic. The Euler characteristic of a compact triangulated surface X is defined to be the alternating sum χ(x) = V E + F where V, E and F are the
More informationTHE RISE, FALL AND RISE OF SIMPLICIAL COMPLEXES. Andrew Ranicki (Edinburgh) aar Warwick 2nd February, 2018
1 THE RISE, FALL AND RISE OF SIMPLICIAL COMPLEXES Andrew Ranicki (Edinburgh) http://www.maths.ed.ac.uk/ aar Warwick 2nd February, 2018 2 Simplicial complexes A simplicial complex is a combinatorial scheme
More informationAlpha-Beta Witness Complexes
Alpha-Beta Witness Complexes Dominique Attali 1, Herbert Edelsbrunner 2, John Harer 3, and Yuriy Mileyko 4 1 LIS-CNRS, Domaine Universitaire, BP 46, 38402 Saint Martin d Hères, France 2 Departments of
More informationDual Complexes of Cubical Subdivisions of R n
Dual Complexes of Cubical Subdivisions of R n Herbert Edelsbrunner and Michael Kerber Abstract We use a distortion to define the dual complex of a cubical subdivision of R n as an n-dimensional subcomplex
More informationVoronoi Diagram. Xiao-Ming Fu
Voronoi Diagram Xiao-Ming Fu Outlines Introduction Post Office Problem Voronoi Diagram Duality: Delaunay triangulation Centroidal Voronoi tessellations (CVT) Definition Applications Algorithms Outlines
More information