inverse problems in tda focus on metric graphs
play

Inverse problems in TDA Focus on metric graphs Steve Oudot joint - PowerPoint PPT Presentation

Geometry and Topology of Data ICERM, December 2017 Inverse problems in TDA Focus on metric graphs Steve Oudot joint work with Elchanan (Isaac) Solomon (Brown University) Persistence diagrams as descriptors for data ) ( g Model i n n


  1. Geometry and Topology of Data ICERM, December 2017 Inverse problems in TDA Focus on metric graphs Steve Oudot — joint work with Elchanan (Isaac) Solomon (Brown University)

  2. Persistence diagrams as descriptors for data ) ( g Model i n n i f l e p r m e n a c s e ( ) Data Descriptor (TDA) - genericity - stability - invariance - · · · 1

  3. Persistence diagrams as descriptors for data ) ( g Model i n n i f l e p r m e n a c s e ( ) Data Descriptor (TDA) - genericity - stability autoencoders! - invariance - · · · 1

  4. Persistence diagrams as descriptors for data ) ( g Model i n n i f l e p r m e n a c s e ( ) noise Data Descriptor signal (TDA) Challenges: - signal vs noise discrimination - model inference i.e. infer a model explaining the data / descriptor - hypothesis testing - · · · 1

  5. Persistence diagrams as descriptors for data ) ( g Model i n n i f l e p r m e n a c s e ( ) noise Data Descriptor signal (TDA) Challenges: - signal vs noise discrimination - model inference model uniqueness? i.e. infer a model explaining the data / descriptor - hypothesis testing - · · · 1

  6. Persistence diagrams as descriptors for data ) ( g Model i n n i f l e p r m e n a c s e ( ) Data Descriptor (TDA) Lipschitz operator left inverse? 1

  7. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations offsets filtration barcode / diagram point cloud simplicial filtration 2

  8. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations t C t ( X, d X ) R 2 t ( X, d X ) 2

  9. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations 1 1 α ≥ π/ 2 dgm C ( P, ℓ 2 ) = { (0 , + ∞ ) } ⊔ { (0 , 1 2 ) } ⊔ { (0 , 1 2 ) } dgm R ( P, ℓ 2 ) = { (0 , + ∞ ) } ⊔ { (0 , 1) } ⊔ { (0 , 1) } ⇒ diagrams for different values of α are indistinguishable 2

  10. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations Prop: [Folklore] For any metric tree ( X, d X ) : dgm R ( X, d X ) = dgm C ( X, d X ) = { (0 , + ∞ ) } ⇒ no information on the metric X is 0 -hyperbolic ⇒ metric balls are convex ⇒ geodesic triangles are tripods 2

  11. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations • Reeb graphs ⇒ Reeb graphs are indistinguishable from their diagrams 2

  12. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations • Reeb graphs • Real-valued functions Prop: [Folklore] Given f : X → R and h : Y → X homeomorphism, dgm f ◦ h = dgm f this is too large a class of transformations, for our purposes we rather target isometries ⇒ Persistence is invariant under reparametrizations 2

  13. Lack of injectivity in general • Unions of (open) balls — ˇ Cech/Rips/Delaunay filtrations • Reeb graphs • Real-valued functions possible solutions: • richer topological invariants (e.g. persistent homotopy) • use several filter functions ( concatenation vs multipersistence) 2

  14. Persistent Homology Transform (PHT) implicit: PHT( X ) = PHT( X, F ) PHT( X )= { dgm f w | w ∈ W } ( X, d X ) (compact) PHT( X ) · · · note: here, as before, dgm f contains dgm f w F = { f w } w ∈ W (diagrams, d b ) R 3

  15. Persistent Homology Transform (PHT) implicit: PHT( X ) = PHT( X, F ) PHT( X )= { dgm f w | w ∈ W } ( X, d X ) (compact) PHT( X ) · · · note: here, as before, dgm f contains dgm f w F = { f w } w ∈ W (diagrams, d b ) R Thm: [Turner, Mukherjee, Boyer 2014] X S d − 1 Let F = {�· , w �} w ∈ S d − 1 , where d = 2 , 3 is fixed. Then, PHT is injective on the set of linear embeddings w of compact simplicial complexes in R d . Extension: [Turner et al., in progress] True for arbitrary d and semialgebraic compact sets. 3

  16. Persistent Homology Transform (PHT) implicit: PHT( X ) = PHT( X, F ) PHT( X )= { dgm f w | w ∈ W } ( X, d X ) (compact) PHT( X ) · · · note: here, as before, dgm f contains dgm f w F = { f w } w ∈ W Q: can we derive an intrinsic version? (diagrams, d b ) R Thm: [Turner, Mukherjee, Boyer 2014] X S d − 1 Let F = {�· , w �} w ∈ S d − 1 , where d = 2 , 3 is fixed. Then, PHT is injective on the set of linear embeddings w of compact simplicial complexes in R d . Extension: [Turner et al., in progress] True for arbitrary d and semialgebraic compact sets. 3

  17. PHT for intrinsic metrics Given ( X, d X ) compact, take F = { d X ( · , x ) } x ∈ X 4

  18. PHT for intrinsic metrics Given ( X, d X ) compact, take F = { d X ( · , x ) } x ∈ X Thm (local stability): [Carri` ere, O., Ovsjanikov 2015] this construction holds for general metric spaces, however it makes sense mostly for compact Let ( X, d X ) and ( Y, d Y ) be compact length spaces with positive convexity radius ( ̺ ( X ) , ̺ ( Y ) > 0 ). Let x ∈ X and y ∈ Y . 1 If d GH (( X, x ) , ( Y, y )) ≤ 20 min { ̺ ( X ) , ̺ ( Y ) } , then d b (dgm d X ( · , x ) , dgm d Y ( · , y )) ≤ 20 d GH (( X, x ) , ( Y, y )) . 4

  19. PHT for intrinsic metrics Given ( X, d X ) compact, take F = { d X ( · , x ) } x ∈ X Corollary (local stability of PHT): Let ( X, d X ) and ( Y, d Y ) be compact length spaces with positive convexity radius ( ̺ ( X ) , ̺ ( Y ) > 0 ). 1 If d GH ( X, Y ) ≤ 20 min { ̺ ( X ) , ̺ ( Y ) } , then d H ( PHT ( X ) , PHT ( Y )) ≤ 20 d GH (( X, x ) , ( Y, y )) . Thm (local stability): [Carri` ere, O., Ovsjanikov 2015] this construction holds for general metric spaces, however it makes sense mostly for compact Let ( X, d X ) and ( Y, d Y ) be compact length spaces with positive convexity radius ( ̺ ( X ) , ̺ ( Y ) > 0 ). Let x ∈ X and y ∈ Y . 1 If d GH (( X, x ) , ( Y, y )) ≤ 20 min { ̺ ( X ) , ̺ ( Y ) } , then d b (dgm d X ( · , x ) , dgm d Y ( · , y )) ≤ 20 d GH (( X, x ) , ( Y, y )) . 4

  20. PHT for intrinsic metrics Given ( X, d X ) compact, take F = { d X ( · , x ) } x ∈ X Corollary (local stability of PHT): Let ( X, d X ) and ( Y, d Y ) be compact length spaces with positive convexity radius ( ̺ ( X ) , ̺ ( Y ) > 0 ). 1 If d GH (( X ) , ( Y )) ≤ 20 min { ̺ ( X ) , ̺ ( Y ) } , then d H ( PHT ( X ) , PHT ( Y )) ≤ 20 d GH (( X, x ) , ( Y, y )) . a # X →∞ d GH ( T, X ) − → 0 b b d H ( PHT ( T ) , PHT ( X )) is bounded away from 0 a 4

  21. PHT for metric graphs from now on we will focus on compact metric graphs, which are length spaces that admit Focus: compact metric graphs (1-dimensional stratified length spaces) PHT: F = { d X ( · , x ) } x ∈ X , dgm = extended persistence diagram Thm (global stability): [Dey, Shi, Wang 2015] For any compact metric graphs X, Y , Note: this result does not extend to the class of compact length spaces. Indeed, a graph d H ( PHT ( X ) , PHT ( Y )) ≤ 18 d GH ( X, Y ) . Thm (density): [Gromov] Compact metric graphs are GH-dense among the compact length spaces. 5

  22. PHT for metric graphs from now on we will focus on compact metric graphs, which are length spaces that admit Focus: compact metric graphs (1-dimensional stratified length spaces) PHT: F = { d X ( · , x ) } x ∈ X , dgm = extended persistence diagram Thm (global stability): [Dey, Shi, Wang 2015] For any compact metric graphs X, Y , Note: this result does not extend to the class of compact length spaces. Indeed, a graph d H ( PHT ( X ) , PHT ( Y )) ≤ 18 d GH ( X, Y ) . Thm (density): [Gromov] Compact metric graphs are GH-dense among the compact length spaces. Q: injectivity of PHT on metric graphs? 5

  23. PHT for metric graphs Negative result: PHT is not injective on all compact metric graphs X Y PHT ( X ) = PHT ( Y ) while X �≃ Y 5

  24. PHT for metric graphs Negative result: PHT is not injective on all compact metric graphs X Y PHT ( X ) = PHT ( Y ) while X �≃ Y Note: Aut ( X ) is non-trivial, hence Ψ X : x �→ dgm d X ( · , x ) is not injective So, maybe there is a connection between Ψ X being injective and PHTitself being injective... this is precisely what our results show 5

  25. PHT for metric graphs Let Inj Ψ = { X compact metric graph s.t. Ψ X is injective } Thm 1: PHT is injective on Inj Ψ . Thm 2: Inj Ψ is GH-dense among the compact metric graphs. 1 1 0.9 1.1 5 1.1 5 10 ⇒ Note: Ψ X injective Aut ( X ) trivial Thus, Inj Ψ is a strict subset of the graphs with trivial automorphism group p q �⇒ 0.9 5 6 1 dgm d X ( · , p ) = dgm d X ( · , q ) 5

  26. PHT for metric graphs Let Inj Ψ = { X compact metric graph s.t. Ψ X is injective } Thm 1: PHT is injective on Inj Ψ . + Gromov’s density result Thm 2: Inj Ψ is GH-dense among the compact metric graphs. Corollary: There is a GH-dense subset of the compact length spaces on which PHT is injective. 1 1 0.9 1.1 5 1.1 5 10 ⇒ Note: Ψ X injective Aut ( X ) trivial Thus, Inj Ψ is a strict subset of the graphs with trivial automorphism group p q �⇒ 0.9 5 6 1 dgm d X ( · , p ) = dgm d X ( · , q ) 5

  27. PHT for metric graphs Let Inj Ψ = { X compact metric graph s.t. Ψ X is injective } Thm 1: PHT is injective on Inj Ψ . + Gromov’s density result Thm 2: Inj Ψ is GH-dense among the compact metric graphs. Corollary: There is a GH-dense subset of the compact length spaces on which PHT is injective. Thm 3: PHT is GH- locally injective on compact metric graphs. 5

Recommend


More recommend