A roadmap for the computation of persistent homology
|
|
- Julia Boone
- 6 years ago
- Views:
Transcription
1 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 and Conference on Topological Data Analysis Querétaro, 11 December 2015 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
2 What is out there? Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
3 Perseus javaplex Dionysus phom TDA GAP persistence PHAT DIPHA GUDHI Kenzo module The persistence landscape toolbox SimpPers jholes HoPeS RIVET Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
4 A brief (biased) history of PH software PH algorithm H. Edelsbrunner, D. Letscher, G. Carlsson and A. Zomorodian 2005 Plex V. de Silva, P. Perry, L. Kettner, A. Zomorodian 2011 javaplex A. Tausz, M. Vejdemo-Johansson, H. Adams 2012 Perseus 2013 PHAT V. Nanda M. Kerber, J. Reininghaus, U. Bauer, H. Wagner 2014 DIPHA 2014 GUDHI 2016 RIVET Nina Otter (Oxford) M. Kerber, J. Reininghaus, U. Bauer C. Maria, J.D. Boissonnat, M. Glisse, M. Yvinec M. Lesnick, M. Wright Roadmap for PH computation Quere taro, 11 December / 45
5 Pipeline for PH computation Data (1) Filtered complex (2) (3) Barcodes Interpretation Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
6 Step 1: from data to a filtered complex Point cloud data Point cloud S Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
7 Step 1: from data to a filtered complex Point cloud data Point cloud S Filtered simplicial complex Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
8 Čech complex C ɛ (S) = { σ S x σ B x(ɛ) } Nerve theorem: C ɛ (S) x S B x(ɛ) Complexity: 2 O(n) Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
9 Vietoris-Rips complex VR ɛ (S) = {σ S x, y σ : d(x, y) 2ɛ} Nerve theorem does not hold Complexity: 2 O(n) For all ɛ 0 we have C ɛ (S) VR ɛ (S) C 2ɛ (S). Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
10 α complex α ɛ (S) = { σ S x σ B x(ɛ) V x } Nerve theorem: α ɛ (S) x S B x(ɛ) Complexity: n O( d/2 ) (n points in R d ) For all ɛ 0 we have α ɛ (S) C ɛ (S) VR ɛ (S). Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
11 Optimisations and speed-ups (A) Sparse filtrations: construct sparse filtered complexes which approximate homology of a given complex. Witness complex [de Silva, Carlsson 2004] Graph induced complex [Dey, Fang, Wang 2014] Linear size approximation of VR [Sheehy 2013] Linear size approximation of Čech [Kerber, Sharathkumar 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
12 Optimisations and speed-ups (A) Sparse filtrations: construct sparse filtered complexes which approximate homology of a given complex. Witness complex [de Silva, Carlsson 2004] Graph induced complex [Dey, Fang, Wang 2014] Linear size approximation of VR [Sheehy 2013] Linear size approximation of Čech [Kerber, Sharathkumar 2013] Graph induced complex on point data, Dey, Fan, Wang 2014 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
13 Optimisations and speed-ups (A) Sparse filtrations: construct sparse filtered complexes which approximate homology of a given complex. Witness complex [de Silva, Carlsson 2004] Graph induced complex [Dey, Fang, Wang 2014] linear size approximation of VR [Sheehy 2013] linear size approximation of Čech [Kerber, Sharathkumar 2013] (B) Heuristics: given a filtered simplicial complex, reduce number of simplices without changing homology. Tidy set [Zomorodian 2010] Morse reductions [Mischaikow, Nanda 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
14 Optimisations summary Complex K Size of K (A) Approximation (B) Reduction Čech O(2 n ) linear size approx. tidy set Morse red. Vietoris Rips (VR) O(2 n ) linear size approx. tidy set (lazy) witness Morse red. GIC α Data structures: n O( d/2 ) (n points in R d ) Simplex tree [Boissonnat, Maria 2012] - tidy set Morse red. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
15 Implementations Complex K (A) Approximation (B) Reduction Čech Dionysus linear size approx.? tidy set not open source Morse red. Perseus Vietoris Rips (VR) Dionysus,Perseus, PHAT DIPHA, javaplex, GUDHI α Dionysus, HoPeS (1-D) linear size approx.? (lazy) witness javaplex GIC? - tidy set Morse red. tidy set Morse red. Data structure: Simplex tree GUDHI Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
16 Other data types 1 Undirected (weighted) networks Weight rank clique filtration [Petri et al 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
17 Other data types 1 Undirected (weighted) networks Weight rank clique filtration [Petri et al 2013] 2 Image data: cubical complexes Lower star filtration [Wagner et al 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
18 Other data types: implementations 1 Undirected (weighted) networks Weight rank clique filtration jholes 2 Image data: cubical complexes DIPHA, Perseus Lower star filtration DIPHA Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
19 Step 2: from a filtered complex to barcodes Standard algorithm Step 1: Order the simplices so that the total order is compatible with the filtration. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
20 Step 2: from a filtered complex to barcodes Standard algorithm Step 1: Order the simplices so that the total order is compatible with the filtration. Step 2: Construct the boundary matrix. Let n be the number of simplices. The boundary matrix B is an n n-matrix defined by: { 1 σ i σ j and dim(σ i ) = dim(σ j ) 1 B(i, j) = 0 otherwise Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
21 Step 2: from a filtered complex to barcodes Standard algorithm Step 1: Order the simplices so that the total order is compatible with the filtration. Step 2: Construct the boundary matrix. Let n be the number of simplices. The boundary matrix B is an n n-matrix defined by: { 1 σ i σ j and dim(σ i ) = dim(σ j ) 1 B(i, j) = 0 otherwise Step 3: Reduce the matrix using column additions from left to right. Step 4: Read the endpoints of the persistence intervals from the matrix to get the barcode. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
22 Standard algorithm Step 1. Given a filtered simplicial complex = K 0 K 1 K m = K put an order on its simplices such that: A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
23 Standard algorithm Step 1. Given a filtered simplicial complex = K 0 K 1 K m = K put an order on its simplices such that: A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. Example: K 1 K 2 K 3 K 4 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
24 Standard algorithm Put an order on the simplices such that: σ 1 A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. σ 2 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
25 Standard algorithm Put an order on the simplices such that: A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. σ 2 σ 1 σ 4 σ 1 σ 2 σ 3 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
26 Standard algorithm Put an order on the simplices such that: A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. σ 2 σ 1 σ 2 σ 2 σ 4 σ 1 σ 3 σ 4 σ5 σ 1 σ 6 σ 3 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
27 Standard algorithm Put an order on the simplices such that: A face of a simplex precedes the simplex. A simplex in K i precedes simplices in K \ K i. σ 2 σ 1 σ 4 σ 6 σ 2 σ 2 σ 7 σ 1 σ 3 σ 4 σ5 σ 1 σ 3 σ 1 σ 2 σ 4 σ5 σ 6 σ 3 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
28 Standard algorithm Step 2. We obtain the boundary matrix: For j = 1,..., 7 define low(j) = i if i is the position of the lowest 1 in column j. If the column is zero leave low(j) undefined. e.g. low(4) = 2 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
29 Standard algorithm Step 3. Reduce the boundary matrix. Let n be the number of simplices. Algorithm: for j = 1 to n do while there exists i < j with low(i) = low(j) do add column i to column j end while end for Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
30 Standard algorithm Step 3. Reduce the boundary matrix. Algorithm: for j = 1 to n do while there exists i < j with low(i) = low(j) do add column i to column j end while end for Add column 5 to column 6: Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
31 Standard algorithm Step 3. Reduce the boundary matrix. Algorithm: for j = 1 to n do while there exists i < j with low(i) = low(j) do add column i to column j end while end for Add column 4 to column 6: Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
32 Standard algorithm Step 4. Read the persistence pairs. If low(j) = i then σ j is negative and paired with the positive σ i. If low(j) is undefined then σ j is positive. If there exists k such that low(k) = j then σ j is paired with the negative simplex σ k. If no such k exists σ j is unpaired σ 1 positive, unpaired σ 2 positive, paired with σ 4 σ 3 positive, paired with σ 5 σ 4 negative, paired with σ 2 σ 5 negative, paired with σ 3 σ 6 positive, paired with σ 7 σ 7 negative, paired with σ 6 Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
33 Standard algorithm σ 2 σ 1 σ 4 σ 6 σ 2 σ 2 σ 2 σ 7 σ 1 σ 3 σ 4 σ5 σ 1 σ 3 σ 4 σ5 σ 1 σ 6 σ 3 σ 1 positive, unpaired interval [1, ) in H 0. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
34 Standard algorithm σ 2 σ 1 σ 4 σ 6 σ 2 σ 2 σ 2 σ 7 σ 1 σ 3 σ 4 σ5 σ 1 σ 3 σ 4 σ5 σ 1 σ 6 σ 3 σ 1 positive, unpaired interval [1, ) in H 0. σ 2 positive, paired with σ 4 no interval, since σ 2 and σ 4 enter at the same time in the filtration Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
35 Standard algorithm σ 2 σ 1 σ 4 σ 6 σ 2 σ 2 σ 2 σ 7 σ 1 σ 3 σ 4 σ5 σ 1 σ 3 σ 4 σ5 σ 1 σ 6 σ 3 σ 1 positive, unpaired interval [1, ) in H 0. σ 2 positive, paired with σ 4 no interval, since σ 2 and σ 4 enter at the same time in the filtration σ 3 positive, paired with σ 5 interval [2, 3) in H 0. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
36 Standard algorithm σ 2 σ 1 σ 4 σ 6 σ 2 σ 2 σ 2 σ 7 σ 1 σ 3 σ 4 σ5 σ 1 σ 3 σ 4 σ5 σ 1 σ 6 σ 3 σ 1 positive, unpaired interval [1, ) in H 0. σ 2 positive, paired with σ 4 no interval, since σ 2 and σ 4 enter at the same time in the filtration σ 3 positive, paired with σ 5 interval [2, 3) in H 0. σ 6 positive, paired with σ 7 interval [3, 4) in H 1. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
37 Twist algorithm Note: If column j is reduced and low(j) = i > 0 then simplex σ i is positive and reducing column i will result in setting column i to zero. σ i is a codimension 1 face of σ j Algorithm: Reduce columns from right to left If column j is reduced and low(j) = i > 0 then set column i to zero. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
38 Twist algorithm Previous example: Column 7 is reduced and low(7) = 6, so set column 6 to zero. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
39 Twist algorithm Previous example: Column 7 is reduced and low(7) = 6, so set column 6 to zero Now we are done. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
40 Spectral sequence algorithm Let k j denote the number of simplices in subcomplex K j. Let B j denote the columns numbered k j to k j. Let B i denote the rows numbered k i to k i. Idea: Reduce the matrix in phases: in each phase r, reduce columns in B j by adding columns in the blocks from B j r+1 to B j. Optimisation: The reduction in each block and each phase is independent, and can be executed in parallel. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
41 Spectral sequence algorithm Phase r = 1: Phase r = 2: B B j B k B B j B k B B B j B j B k B k Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
42 Spectral sequence algorithm Phase r = 3: Phase r = k: B B j B k B B j B k B B B j B j B k B k Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
43 Heuristic speed-ups Sequential: Dual algorithm [de Silva, Morozov, Vejdemo-Johansson 2011] Twist algorithm [Chen, Kerber 2011] Parallel: Spectral sequence algorithm [Edelsbrunner, Harer 2008] Chunk algorithm [Bauer, Kerber, Reininghaus 2013] Distributed algorithm [Bauer, Kerber, Reininghaus 2014] Blow-up complex [Morozov, Lewis 2015] Data structure: Bit tree pivot column [Bauer, Kerber, Reininghaus, Wagner 2013] Compressed annotation matrix [Boissonnat, Dey, Maria 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
44 Speed-ups: implementations Sequential: Dual algorithm Dionysus, javaplex, GUDHI, PHAT, DIPHA Twist algorithm PHAT,DIPHA Parallel: Spectral sequence algorithm PHAT Chunk algorithm PHAT Distributed algorithm DIPHA Blow-up complex open source? Data structure: Bit tree pivot column PHAT, DIPHA Compressed annotation matrix GUDHI Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
45 Step 3: interpretation (A) Comparing diagrams: p-th Wasserstein distance W p [σ](x, Y ) := inf φ: X Y [ ] 1/p σ(x, φ(x)) p x X for 1 p < and σ a metric on the plane. W [σ](x, Y ) := Example: W [L ] is the bottleneck distance. inf sup σ(x, φ(x) φ: X Y x X Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
46 (B) Statistics: The space (D, W p [σ]) (D, W p [L ]) is a Polish space. Fréchet expectations exist but are not unique nor continuous. [Mileyko et al 2011] Same is true for (D, W 2 [L 2 ]). Algorithm for estimation of Fréchet mean [Turner et al 2014] Fréchet mean for continuously varying diagrams as a probability distribution for (D, W 2 [L 2 ]). [Munch et al 2015] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
47 (C) Alternative: map persistence diagrams to simpler space Persistence landscapes [Bubenik 2014] Space of algebraic functions [Adcock et al 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
48 Implementations Dionysus and javaplex: bottleneck distance TDA package: bottleneck and Wasserstein distances, persistence landscapes, Kernel distances,... Persistence Landscape Toolbox: persistence landscapes (means, norms,... ), T -student test,... Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
49 Summary 1 of 2 Software Precomp. filt. Compl. Parallel PH algorithms Interpr. javaplex VR, LW, W, CW Perseus VR Dionysus α, VR, Čech PHAT shared DIPHA VR, lower star distr. GUDHI VR standard, dual zig-zag Morse reductions + standard standard, dual zig-zag standard, dual, twist chunk, spectral seq. dual, twist distributed multifield dual bottleneck bottleneck Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
50 Summary 2 of 2 Software Precomp. filt. Compl. Parallel PH algorithms Interpr. jholes WRCF shared standard (javaplex) HoPeS α (only 1-D) TDA VR dual (GUDHI, Dionysus, PHAT) Wasserstein, landscapes, bootstraps, etc. phom VR,LW standard, dual Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
51 Benchmarking We study the libraries javaplex, Perseus, Dionysus, DIPHA and GUDHI Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
52 Benchmarking We study the libraries javaplex, Perseus, Dionysus, DIPHA and GUDHI... by using complexes and algorithms implemented by the largest number of libraries: Vietoris-Rips complex and standard and dual algorithms Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
53 Methods We study the software from three different points of view: 1 Performance measured in CPU and real time 2 Memory usage 3 Maximum size of simplicial complex allowed by the software Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
54 Data Synthetic data Points sampled from Klein bottle Random Vietoris Rips complexes Real world data Genomic sequence of HIV virus [Chan et al, 2013] 3D scans of Stanford dragon C. Elegans neuronal network [Petri et al, 2013] Human genome network [Petri et al, 2013] Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
55 Machines Cluster 1 : 1728 (180*16) cores of 2.0GHz RAM : 64 GiB x80 nodes, 128 GiB x 4 nodes 1 Advanced Research Computing (ARC), Oxford Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
56 Results: performance Data set C. elegans Klein HIV Dragon 1 Dragon 2 Size of complex javaplex st Dionysus st DIPHA st Perseus Dionysus d DIPHA d GUDHI Table: Performance measured in wall-time and CPU seconds for the computation of PH with the Vietoris Rips complex. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
57 Results: memory usage Data set C. elegans Klein HIV Dragon 1 Dragon 2 Size of complex javaplex st < 5 < 15 > 120 > 120 > 120 Dionysus st DIPHA st Perseus Dionysus d DIPHA d GUDHI Table: Memory usage in GB. For javaplex, we indicate the value of the maximum heap size that was sufficient to perform the computation. For DIPHA, we indicate the maximum memory used by a single core. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
58 Results: maximal size of simplicial complex Software javaplex Dionysus DIPHA Perseus GUDHI st st d st d st d max. size Table: Maximal size of simplicial complex supported by the software. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
59 Conclusion Future directions: Creation of a computational topology library. User-friendly interfaces. Uniformization of input type across different implementations. Focus on step (1) and (3) of pipeline: implementation and construction of simpler simplicial complexes and interpretation of results. Multi-parameter PH. Nina Otter (Oxford) Roadmap for PH computation Querétaro, 11 December / 45
Topological Persistence
Topological Persistence Topological Persistence (in a nutshell) X topological space f : X f persistence X Dg f signature: persistence diagram encodes the topological structure of the pair (X, f) 1 Topological
More informationarxiv: v7 [math.at] 12 Sep 2017
A Roadmap for the Computation of Persistent Homology arxiv:1506.08903v7 [math.at] 12 Sep 2017 Nina Otter (otter@maths.ox.ac.uk) Mason A. Porter (mason@math.ucla.edu) Ulrike Tillmann (tillmann@maths.ox.ac.uk)
More informationWhat will be new in GUDHI library version 2.0.0
What will be new in GUDHI library version 2.0.0 Jean-Daniel Boissonnat, Paweł Dłotko, Marc Glisse, François Godi, Clément Jamin, Siargey Kachanovich, Clément Maria, Vincent Rouvreau and David Salinas DataShape,
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 informationPSB Workshop Big Island of Hawaii Jan.4-8, 2015
PSB Workshop Big Island of Hawaii Jan.4-8, 2015 Part I Persistent homology & applications Part II TDA & applications Motivating problem from microbiology Algebraic Topology 2 Examples of applications
More informationJPlex. CSRI Workshop on Combinatorial Algebraic Topology August 29, 2009 Henry Adams
JPlex CSRI Workshop on Combinatorial Algebraic Topology August 29, 2009 Henry Adams Plex: Vin de Silva and Patrick Perry (2000-2006) JPlex: Harlan Sexton and Mikael Vejdemo- Johansson (2008-present) JPlex
More informationJPlex Software Demonstration. AMS Short Course on Computational Topology New Orleans Jan 4, 2011 Henry Adams Stanford University
JPlex Software Demonstration AMS Short Course on Computational Topology New Orleans Jan 4, 2011 Henry Adams Stanford University What does JPlex do? Input: a filtered simplicial complex or finite metric
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 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 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 informationarxiv: v1 [math.at] 1 Nov 2018
Persistent-Homology-based Machine Learning and its Applications A Survey Persistent-Homology-based Machine Learning and its Applications A Survey arxiv:1811.00252v1 [math.at] 1 Nov 2018 Chi Seng Pun School
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 informationParallel & scalable zig-zag persistent homology
Parallel & scalable zig-zag persistent homology Primoz Skraba Artificial Intelligence Laboratory Jozef Stefan Institute Ljubljana, Slovenia primoz.skraba@ijs.sl Mikael Vejdemo-Johansson School of Computer
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 informationTopological Inference for Modern Data Analysis
Topological Inference for Modern Data Analysis An Introduction to Persistent Homology Giancarlo Sanchez A project presented for the degree of Master s of Science in Mathematics Department of Mathematics
More informationComputational Topology in Reconstruction, Mesh Generation, and Data Analysis
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
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 informationA persistence landscapes toolbox for topological statistics
A persistence landscapes toolbox for topological statistics Peter Bubenik, Dlotko Pawel To cite this version: Peter Bubenik, Dlotko Pawel. A persistence landscapes toolbox for topological statistics. Journal
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 informationStatistical properties of topological information inferred from data
Saclay, July 8, 2014 DataSense Research day Statistical properties of topological information inferred from data Frédéric Chazal INRIA Saclay Joint works with B.T. Fasy (CMU), M. Glisse (INRIA), C. Labruère
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 informationTOPOLOGICAL DATA ANALYSIS
TOPOLOGICAL DATA ANALYSIS BARCODES Ghrist, Barcodes: The persistent topology of data Topaz, Ziegelmeier, and Halverson 2015: Topological Data Analysis of Biological Aggregation Models 1 Questions in data
More informationTopological Data Analysis
Mikael Vejdemo-Johansson School of Computer Science University of St Andrews Scotland 22 November 2011 M Vejdemo-Johansson (St Andrews) TDA 22 November 2011 1 / 32 e Shape of Data Outline 1 e Shape of
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 informationSimBa: An Efficient Tool for Approximating Rips-filtration Persistence via Simplicial Batch-collapse
SimBa: An Efficient Tool for Approximating Rips-filtration Persistence via Simplicial Batch-collapse Tamal K. Dey 1, Dayu Shi 2, and Yusu Wang 3 1 Dept. of Computer Science and Engineering, and Dept. of
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 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 informationAn Analysis of Spaces of Range Image Small Patches
Send Orders for Reprints to reprints@benthamscience.ae The Open Cybernetics & Systemics Journal, 2015, 9, 275-279 275 An Analysis of Spaces of Range Image Small Patches Open Access Qingli Yin 1,* and Wen
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 informationWestern TDA Learning Seminar. June 7, 2018
The The Western TDA Learning Seminar Department of Mathematics Western University June 7, 2018 The The is an important tool used in TDA for data visualization. Input point cloud; filter function; covering
More informationPersistent homology in TDA
Barcelona, June, 2016 Persistent homology in TDA Frédéric Chazal and Bertrand Michel Persistent homology R X topological space f : X R X persistence Persistence diagram Nested spaces A general mathematical
More informationPackage TDA. December 9, 2017
Type Package Package TDA December 9, 2017 Title Statistical Tools for Topological Data Analysis Version 1.6 Date 2017-12-09 Author Brittany T. Fasy, Jisu Kim, Fabrizio Lecci, Clement Maria, Vincent Rouvreau.
More informationRecent Advances and Trends in Applied Algebraic Topology
Recent Advances and Trends in Applied Algebraic Topology Mikael Vejdemo-Johansson School of Computer Science University of St Andrews Scotland July 8, 2012 M Vejdemo-Johansson (St Andrews) Trends in applied
More informationTowards topological analysis of high-dimensional feature spaces.
Towards topological analysis of high-dimensional feature spaces. Hubert Wagner and Pawe l D lotko January 7, 2014 Abstract In this paper we present ideas from computational topology, applicable in analysis
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 informationGoodies in Statistic and ML. DataShape, Inria Saclay
Goodies in Statistic and ML. DataShape, Inria Saclay Gudhi is in Statistic and ML. DataShape, Inria Saclay Where did I start playing with statistics? Analysis of time varying patterns from dynamical systems,
More informationIntroduction to the R package TDA
Introduction to the R package TDA Brittany T. Fasy Computer Science Department,Tulane University Jisu Kim arxiv:1411.1830v2 [cs.ms] 29 Jan 2015 Abstract Department of Statistics, Carnegie Mellon University
More informationMetrics on diagrams and persistent homology
Background Categorical ph Relative ph More structure Department of Mathematics Cleveland State University p.bubenik@csuohio.edu http://academic.csuohio.edu/bubenik_p/ July 18, 2013 joint work with Vin
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 informationarxiv: v1 [cs.it] 10 May 2016
Separating Topological Noise from Features using Persistent Entropy Nieves Atienza 1, Rocio Gonzalez-Diaz 1, and Matteo Rucco 2 arxiv:1605.02885v1 [cs.it] 10 May 2016 1 Applied Math Department, School
More informationTopological Classification of Data Sets without an Explicit Metric
Topological Classification of Data Sets without an Explicit Metric Tim Harrington, Andrew Tausz and Guillaume Troianowski December 10, 2008 A contemporary problem in data analysis is understanding the
More informationStatistical Methods in Topological Data Analysis for Complex, High-Dimensional Data
Libraries Conference on Applied Statistics in Agriculture 2015-27th Annual Conference Proceedings Statistical Methods in Topological Data Analysis for Complex, High-Dimensional Data Patrick S. Medina Purdue
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 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 informationPersistent Homology and Nested Dissection
Persistent Homology and Nested Dissection Michael Kerber Donald R. Sheehy Primoz Skraba Abstract Nested dissection exploits the underlying topology to do matrix reductions while persistent homology exploits
More informationPersistent homology: a step-by-step introduction for newcomers
STAG: Smart Tools and Apps in computer Graphics (2016) Giovanni Pintore and Filippo Stanco (Editors) Persistent homology: a step-by-step introduction for newcomers U. Fugacci 1, S. Scaramuccia 2, F. Iuricich
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 informationTopology and the Analysis of High-Dimensional Data
Topology and the Analysis of High-Dimensional Data Workshop on Algorithms for Modern Massive Data Sets June 23, 2006 Stanford Gunnar Carlsson Department of Mathematics Stanford University Stanford, California
More informationPersistent Homology of Directed Networks
Persistent Homology of Directed Networks Samir Chowdhury and Facundo Mémoli Abstract While persistent homology has been successfully used to provide topological summaries of point cloud data, the question
More informationRandom Simplicial Complexes
Random Simplicial Complexes Duke University CAT-School 2015 Oxford 9/9/2015 Part II Random Geometric Complexes Contents Probabilistic Ingredients Random Geometric Graphs Definitions Random Geometric Complexes
More informationTopological Data Analysis
Topological Data Analysis Deepak Choudhary(11234) and Samarth Bansal(11630) April 25, 2014 Contents 1 Introduction 2 2 Barcodes 2 2.1 Simplical Complexes.................................... 2 2.1.1 Representation
More informationJPlex with BeanShell Tutorial
JPlex with BeanShell Tutorial Henry Adams henry.adams@colostate.edu December 26, 2011 NOTE: As of 2011, JPlex is being phased out in favor of the software package Javaplex, which is also written in Java
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 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 informationCOMPUTATIONAL TOPOLOGY TO MONITOR HUMAN OCCUPANCY. Institute of Information Science and Technologies - CNR, Via G. Moruzzi 1, Pisa, ITALY
COMPUTATIONAL TOPOLOGY TO MONITOR HUMAN OCCUPANCY Paolo Barsocchi, Pietro Cassará, Daniela Giorgi, Davide Moroni, Maria Antonietta Pascali Institute of Information Science and Technologies - CNR, Via G.
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 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 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 informationPersistent Homology for Defect Detection in Non-Destructive Evaluation of Materials
More Info at Open Access Database www.ndt.net/?id=18663 Persistent Homology for Defect Detection in Non-Destructive Evaluation of Materials Jose F. CUENCA 1, Armin ISKE 1 1 Department of Mathematics, University
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 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 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 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 informationPersistent Homology in Complex Network Analysis
Persistent Homology Summer School - Rabat Persistent Homology in Complex Network Analysis Ulderico Fugacci Kaiserslautern University of Technology Department of Computer Science July 7, 2017 Anything has
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 informationAn algorithm for calculating top-dimensional bounding chains
An algorithm for calculating top-dimensional bounding chains J. Frederico Carvalho Corresp., 1, Mikael Vejdemo-Johansson 2, Danica Kragic 1, Florian T. Pokorny 1 1 CAS/RPL, KTH, Royal Institute of Technology,
More informationApplied Topology. Instructor: Sara Kališnik Verovšek. Office Hours: Room 304, Tu/Th 4-5 pm or by appointment.
Applied Topology Instructor: Sara Kališnik Verovšek Office Hours: Room 304, Tu/Th 4-5 pm or by appointment. Textbooks We will not follow a single textbook for the entire course. For point-set topology,
More informationOn the Local Behavior of Spaces of Natural Images
On the Local Behavior of Spaces of Natural Images Gunnar Carlsson Tigran Ishkhanov Vin de Silva Afra Zomorodian Abstract In this study we concentrate on qualitative topological analysis of the local behavior
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 informationarxiv: v2 [math.at] 25 Oct 2012
COMPUTING HOMOLOGY AND PERSISTENT HOMOLOGY USING ITERATED MORSE DECOMPOSITION arxiv:0.49v [math.at] 5 Oct 0 PAWE L D LOTKO AND HUBERT WAGNER Abstract. In this paper we present a new approach to computing
More informationarxiv: v2 [stat.ml] 22 Feb 2018
Bartosz Zielin ski Mateusz Juda Matthias Zeppelzauer 2 Clusters within clusters.38.38.3.3.8.8.6.4.2 - -.3.34.46.6.4... Persistence diagram. - -.3.34.46 Persistence diagram.2.. Spa al Quanza on (PI) BoW
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 informationTutorial on the R package TDA
Tutorial on the R package TDA Jisu Kim Brittany T. Fasy, Jisu Kim, Fabrizio Lecci, Clément Maria, Vincent Rouvreau Abstract This tutorial gives an introduction to the R package TDA, which provides some
More informationComputational Topology in Neuroscience
Computational Topology in Neuroscience Bernadette Stolz Lincoln College University of Oxford A dissertation submitted for the degree of Master of Science in Mathematical Modelling & Scientific Computing
More informationThe Topological Data Analysis of Time Series Failure Data in Software Evolution
The Topological Data Analysis of Time Series Failure Data in Software Evolution Joao Pita Costa Institute Jozef Stefan joao.pitacosta@ijs.si Tihana Galinac Grbac Faculty of Engineering University of Rijeka
More informationHomology of Simplicial Complexes
Homology of Simplicial Complexes Math, David Perkinson Introduction. This is an introduction to the homology of simplicial complexes suitable for a first course in linear algebra. It uses little more than
More informationTopology and Data. Gunnar Carlsson Department of Mathematics, Stanford University Stanford, California October 2, 2008
Topology and Data Gunnar Carlsson Department of Mathematics, Stanford University Stanford, California 94305 October 2, 2008 1 Introduction An important feature of modern science and engineering is that
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 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 informationPackage TDA. April 27, 2017
Type Package Package TDA April 27, 2017 Title Statistical Tools for Topological Data Analysis Version 1.5.1 Date 2017-04-27 Author Brittany T. Fasy, Jisu Kim, Fabrizio Lecci, Clement Maria, Vincent Rouvreau.
More informationBARCODES: THE PERSISTENT TOPOLOGY OF DATA
BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 45, Number 1, January 2008, Pages 61 75 S 0273-0979(07)01191-3 Article electronically published on October 26, 2007 BARCODES: THE PERSISTENT
More informationTopology and fmri Data
Topology and fmri Data Adam Jaeger Statistical and Applied Mathematical Sciences Institute & Duke University May 5, 2016 fmri and Classical Methodology Most statistical analyses examine data with a variable
More informationS 1 S 3 S 5 S 7. Homotopy type of VR(S 1,r) as scale r increases
Research Statement, 2016 Henry Adams (adams@math.colostate.edu) I am interested in computational topology and geometry, combinatorial topology, and topology applied to data analysis and to sensor networks.
More informationHomological Reconstruction and Simplification in R3
Homological Reconstruction and Simplification in R3 Dominique Attali, Ulrich Bauer, Olivier Devillers, Marc Glisse, André Lieutier To cite this version: Dominique Attali, Ulrich Bauer, Olivier Devillers,
More informationOn the Nonlinear Statistics of Range Image Patches
SIAM J. IMAGING SCIENCES Vol. 2, No. 1, pp. 11 117 c 29 Society for Industrial and Applied Mathematics On the Nonlinear Statistics of Range Image Patches Henry Adams and Gunnar Carlsson Abstract. In [A.
More informationStochastic Multiresolution Persistent Homology Kernel
Stochastic Multiresolution Persistent Homology Kernel Xiaojin Zhu, Ara Vartanian, Manish Bansal, Duy Nguyen, Luke Brandl Department of Computer Sciences, University of Wisconsin Madison 1210 W. Dayton
More informationPersistence Terrace for Topological Inference of Point Cloud Data
Persistence Terrace for Topological Inference of Point Cloud Data Chul Moon 1, Noah Giansiracusa 2, and Nicole A. Lazar 1 arxiv:1705.02037v2 [stat.me] 23 Dec 2017 1 University of Georgia 2 Swarthmore College
More informationRandom Simplicial Complexes
Random Simplicial Complexes Duke University CAT-School 2015 Oxford 8/9/2015 Part I Random Combinatorial Complexes Contents Introduction The Erdős Rényi Random Graph The Random d-complex The Random Clique
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 informationarxiv: v2 [cs.cg] 8 Feb 2010 Abstract
Topological De-Noising: Strengthening the Topological Signal Jennifer Kloke and Gunnar Carlsson arxiv:0910.5947v2 [cs.cg] 8 Feb 2010 Abstract February 8, 2010 Topological methods, including persistent
More informationAn output-sensitive algorithm for persistent homology
An output-sensitive algorithm for persistent homology Chao Chen Institute of Science and Technology (IST) Austria, Klosterneuburg, Austria Michael Kerber Institute of Science and Technology (IST) Austria,
More informationBuilding a coverage hole-free communication tree
Building a coverage hole-free communication tree Anaïs Vergne, Laurent Decreusefond, and Philippe Martins LTCI, Télécom ParisTech, Université Paris-Saclay, 75013, Paris, France arxiv:1802.08442v1 [cs.ni]
More informationTopological Perspectives On Stratification Learning
Topological Perspectives On Stratification Learning Bei Wang School of Computing Scientific Computing and Imaging Institute (SCI) University of Utah www.sci.utah.edu/~beiwang August 9, 2018 Talk Overview
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 informationA Stable Multi-Scale Kernel for Topological Machine Learning
A Stable Multi-Scale Kernel for Topological Machine Learning Jan Reininghaus, Stefan Huber IST Austria Ulrich Bauer IST Austria, TU München Roland Kwitt University of Salzburg, Austria Abstract Topological
More information3. Persistence. CBMS Lecture Series, Macalester College, June Vin de Silva Pomona College
Vin de Silva Pomona College CBMS Lecture Series, Macalester College, June 2017 o o The homology of a planar region Chain complex We have constructed a diagram of vector spaces and linear maps @ 0 o 1 C
More informationStatistical topological data analysis using. Persistence landscapes
Motivation Persistence landscape Statistical topological data analysis using persistence landscapes Cleveland State University February 3, 2014 funded by AFOSR Motivation Persistence landscape Topological
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 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 informationHomology cycle bases from acyclic matchings
Homology cycle bases from acyclic matchings Dmitry Feichtner-Kozlov University of Bremen Kyoto Workshop, January 019 What is... Applied Topology? studying global features of shapes applications in other
More information