Analyzing Data with 1D Non-Linear Shapes Using Topological Methods

Total Page:16

File Type:pdf, Size:1020Kb

Analyzing Data with 1D Non-Linear Shapes Using Topological Methods Analyzing data with 1D non-linear shapes using topological methods Dissertation Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Suyi Wang, M.S. Graduate Program in Computer Science and Engineering The Ohio State University 2018 Dissertation Committee: Yusu Wang, Advisor Rephael Wenger Tamal K. Dey c Copyright by Suyi Wang 2018 Abstract Shape analysis has been applied in many applications across a broad range of domains. Among various different families of \complex" shapes, the ones with 1D non- linear topological structures (skeletons) are particularly interesting. These shapes are simple, as they can be decomposed into 1-d pieces, but still informative in representing the connectivity and other important information behind data. In this thesis, I focus on two objects from computational topology that have been useful for modeling the skeleton of data: the Reeb graph (and its variants) and the 1-(un)stable manifold from discrete Morse theory, and study their properties as well as applications to shape analysis. The two specific topological objects that I focus on have both already been widely used in practical applications. Further theoretical understanding and applications of these two objects in modeling and studying the skeleton of data will be provided. The first part of the dissertation work concerns the so-called Reeb graph and its loop-free variant, the contour tree, which can be used to provide a 1D tree summary of an input scalar field. It has been commonly used in computer graphics and visual- ization. I have investigated one problem regarding to the theoretical understanding of the contour tree, as well as developing a variant of the Reeb graph to address the issue of noise. Carr et al. has proposed an algorithm for computing the contour tree for a piecewise linear function defined on a simplicial complex domain. This algorithm ii is simple, efficient and has been widely used in practice. However, the algorithm is often applied even when the output contour tree may not exist, in which case the algorithm may not terminate or exit with only partial output. My work provides new understanding of this contour tree algorithm and characterizes the cause for such behavior. I also propose a simple variation of the contour tree (called JS-graph) to handle this situation in practice. The Reeb graph, which provides a general 1D summary of a scalar field, has been used to extract the skeleton behind data. However, when there is significant amount of ambient noise, the Reeb graph may no longer reflect the structures of the data. To handle such ambient noise, I propose and develop a concept called \gradational Reeb graph", which incrementally merges the Reeb graphs from different density levels while keeping structures from high density regions. I demonstrate that when extracting the road network from GPS samples, the \gradational Reeb graph" could capture finer structures, which traditional Reeb graph cannot. In the second part of my thesis, I will focus on the so-called 1-(un)stable manifold from Morse Theory, which is another topological object for modeling and extracting data skeleton. In this thesis, I develop a pipeline for the automatic map reconstruc- tion problem that aims to reconstruct the underlying road network, which can be considered as a hidden geometric graph from GPS trajectory samples. I also develop a Morse theory based approach that can integrate multiple maps (graphs) into a single one, as well as correct an existing map with a partial but more trustworthy map. The effectiveness of the method is demonstrated by reconstructing maps from GPS trajectories sampled in the cites of Athens, Beijing, Berlin and Chicago. This iii work proposes a new methodology with simpler pipeline for handling large maps with better performance, which also advances the state of art. The closing topic of this thesis is the extension of 1-stable manifold to a 3D appli- cation, neuron tracing, which asks to extract trees representing neurons from digital images. Here I aim to tackle the challenges such as the massive size and noisy nature of data. In this work I have developed a new algorithm that handles new large data of 3D meso-scale brain by divide and conquer topological method. Experiments have demonstrated that the proposed method has obtained quantitatively better results than those from existing algorithms in tracing single neurons, while my methods can also produce a summary of trends for more challenging and general injection data containing multiple neurons. iv Dedicated to my parents, advisor, those I love and those who love me v Acknowledgments First of all, I would sincerely express my gratitude to my advisor, Dr. Yusu Wang. She leaded me into the topological data analysis field and has tremendously supported me throughout the study. She helped me with her experience and insight- ful understanding of the field in selecting research topics, tackling critical problems and presenting the works. I am also grateful that she is patient and tolerant towards my mistakes and sometimes even thinks ahead of me to prevent potential detours. It is fortunate to have Dr. Wang as my advisor in the Ph.D. journey. These accom- plishments would not have been achieved without her advices. I would like to thank Dr. Tamal Dey and Dr. Rephael Wenger for their advices on my work and my defense. Thanks to the graduate stuents I have met in the group: Chuanjiang Luo, Lei Wang, Xiaoyin Ge, Fengtao Fan, Andrew Slatton, Alfred Rossi, Dayu Shi, Zhe Dong, Sayan Mandal, Jiayuan Wang, Tianqi Li, Dingkang Wang, Cheng Xin, Elena Farahbakhshtouli, Ryan Slechta and Minghao Tian. Thanks to all the friends I have made during the Ph.D. studies for their companionship. Finally, I would thank my parents and family, especially my uncle Yuan Ma and Yixing Yuan, for their advices and support. vi Vita 2011 . .B.S. Physics, Beijing Normal University 2015 . .M.S. Computer Science and Engineer- ing, The Ohio State University 2016-present . .Ph.D. Candidate, Computer Science and Engineering, The Ohio State University Publications Research Publications S. Wang, Y. Wang, Y. Li Efficient Map Reconstruction and Augmentation via Topo- logical Methods. In Proceedings of the 23rd SIGSPATIAL International Confer- ence on Advances in Geographic Information Systems (SIGSPATIAL '15), 25:1{25:10, 2015. S. Wang, Y. Wang, R. Wenger JS-Graph of Join and Split Trees. In Proceedings of the thirtieth annual symposium on Computational geometry (SOCG'14), 539:539{ 539:548, 2014. L. Che, Y. Xiao, S. Wang, Q. Jiang Application of K-means Clustering Analysis in Chinese Pronunciation Degree. Computer Technology and Development, Vol.21, 223-225, 2011. Fields of Study Major Field: Computer Science and Engineering vii Table of Contents Page Abstract . ii Dedication . .v Acknowledgments . vi Vita......................................... vii List of Tables . xi List of Figures . xii 1. Introduction . .1 1.1 Overview . .1 1.2 Representing Data Skeleton via Reeb graph . .4 1.2.1 Computing Contour Trees . .6 1.2.2 Handling Noise in Reeb graph . .8 1.3 Representing Data Skeleton via 1-stable manifold . 10 1.3.1 DiMorSC Implementation . 11 1.3.2 Automatic Map Reconstruction and Augmentation . 12 1.3.3 Neuron Tracing . 15 2. Preliminary . 19 2.1 Simplicial complex . 19 2.2 Critical points and 1-(un)stable manifold . 20 2.3 Reeb graph . 24 2.4 Persistence Simplification . 25 viii 3. JS-Graph of Join and Split trees . 28 3.1 Constructing a Contour Tree . 28 3.1.1 Computing Join / Split Tree . 29 3.1.2 Computing the Merge Tree . 31 3.1.3 Types of problematic output . 33 3.2 JS-graphs and its Characterization . 35 3.2.1 Definition of JS-graphs . 35 3.2.2 Characterization of JS-graphs . 37 3.3 An Efficient Algorithm to Compute A Linear-size JS-graph . 42 3.3.1 The algorithm . 43 3.3.2 Correctness of Algorithm 3 . 47 3.3.3 Time complexity analysis . 51 3.4 Proofs in Computing JS-Graph . 55 3.4.1 Missing Details from Section 3.2 . 55 3.4.2 Missing Details from Section 3.4.2 . 56 3.4.3 Characterization of Join and Split Trees which have Merge Trees . 63 3.5 Concluding Remarks . 68 4. Gradational Reeb Graph . 70 4.1 Representing data skeleton with Reeb graph . 70 4.2 Gradational Reeb Graph . 73 4.3 Experiment . 77 4.3.1 GPS data under threshold . 77 4.3.2 Merged Result . 79 5. Map Reconstruction . 82 5.1 Related Work . 82 5.2 Map Reconstruction via Morse Theory . 84 5.2.1 (Step 1) Density field construction . 86 5.2.2 (Step 2) Ridge extraction . 87 5.2.3 DiMorSC Implementation . 91 5.3 Map Integration and Augmentation . 98 5.4 Results . 102 6. Neuron Tracing . 112 6.1 Related Work . 113 6.2 Method . 115 ix 6.2.1 Pipeline . 115 6.2.2 Stitching . 119 6.2.3 Tree simplification . 121 6.3 Results . 123 6.3.1 Tracing a single neuron . 123 6.3.2 Tracing multiple larger neurons . 127 6.3.3 Simplifying the neuron bundle . 129 7. Conclusion and Future Work . 134 Bibliography . 139 x List of Tables Table Page 6.1 Running time in seconds . 126 6.2 DIADEM score . 126 xi List of Figures Figure Page 1.1 The Reeb graph of a uniform function on a surface. .5 1.2 Reconstruct road map from GPS traces . .8 2.1 Critical points in 3D: For a maxima(minima), the function values within a small ball are all lower(higher) than it. For a 2-critical point, the points with function values lower than it form a connected compo- nent while the points with function values higher than it are separated into two components.
Recommended publications
  • Surface Topology and Reeb Graph
    Sub-Topics • Compute bounding box • Compute Euler Characteristic • Estimate surface curvature • Line description for conveying surface shape • Extract skeletal representation of shapes • Morse function and surface topology--Reeb graph • Scalar field topology--Morse-Smale complex Surface Topology – Reeb Graph What is Topology? • Topology studies the connectedness of spaces • For us: how shapes/surfaces are connected What is Topology • The study of property of a shape that does not change under deformation – Rules of deformation • 1-1 and onto • Bicontinuous (continuous both ways) • Cannot tear, join, poke or seal holes – A is homeomorphic to B Why is Topology Important? • What is the boundary of an object? • Are there holes in the object? • Is the object hollow? • If the object is transformed in some way, are the changes continuous or abrupt? • Is the object bounded , or does it extend infinitely far? Why is Topology Important? • Inherent and basic properties of a shape • We want to accurately represent and preserve these properties in different applications – Surface reconstruction – Morphing – Texturing – Simplification – Compression Image source: http://www.utdallas.edu/~xxg06 1000/physicsmorphing.htm Topology-Basic Concepts • n-manifold – Set of points – Each point has a neighborhood homeomorphic to an open set of – An n-manifold is a topological space that “locally looks like” the Euclidian space Topology-Basic Concepts • Holes/genus – Genus of a surface is the maximal number of nonintersecting simple closed curves that can be cut on the surface without disconnecting it. • Boundaries Topology-Basic Concepts • Euler’s characteristic function – • #vertices, #edges, #faces – is independent of the polygonization – Specifically, 2 2 • What are c, g, and h? χ = 2 χ = 0 Topology-Basic Concepts • Orientability – Any surface has a triangulation – Orient all triangles CW or CCW – Orientability : any two triangles sharing an edge have opposite directions on that edge.
    [Show full text]
  • Approximation of Metric Spaces by Reeb Graphs: Cycle Rank of A
    Filomat xx (yyyy), zzz–zzz Published by Faculty of Sciences and Mathematics, DOI (will be added later) University of Niˇs, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Approximation of Metric Spaces by Reeb Graphs: Cycle Rank of a Reeb Graph, the Co-rank of the Fundamental Group, and Large Components of Level Sets on Riemannian Manifolds Irina Gelbukh CIC, Instituto Polit´ecnico Nacional, 07738, Mexico City, Mexico Abstract. For a connected locally path-connected topological space X and a continuous function f on it such that its Reeb graph R f is a finite topological graph, we show that the cycle rank of R f , i.e., the first Betti number b1(R f ), in computational geometry called number of loops, is bounded from above by the co-rank of the fundamental group π1(X), the condition of local path-connectedness being important since generally b1(R f ) can even exceed b1(X). We give some practical methods for calculating the co-rank of π1(X) and a closely related value, the isotropy index. We apply our bound to improve upper bounds on the distortion of the Reeb quotient map, and thus on the Gromov-Hausdorff approximation of the space by Reeb graphs, for the distance function on a compact geodesic space and for a simple Morse function on a closed Riemannian manifold. This distortion is bounded from below by what we call the Reeb width b(M) of a metric space M, which guarantees that any real-valued continuous function on M has large enough contour (connected component of a level set).
    [Show full text]
  • The Reeb Graph Edit Distance Is Universal
    The Reeb Graph Edit Distance Is Universal Ulrich Bauer Department of Mathematics, Technical University of Munich (TUM), Germany [email protected] Claudia Landi Dipartimento di Scienze e Metodi dell’Ingegneria, Università degli Studi di Modena e Reggio Emilia, Reggio Emilia, Italy [email protected] Facundo Mémoli Department of Mathematics, The Ohio State University, Columbus, OH, USA [email protected] Abstract We consider the setting of Reeb graphs of piecewise linear functions and study distances between them that are stable, meaning that functions which are similar in the supremum norm ought to have similar Reeb graphs. We define an edit distance for Reeb graphs and prove that it is stable and universal, meaning that it provides an upper bound to any other stable distance. In contrast, via a specific construction, we show that the interleaving distance and the functional distortion distance on Reeb graphs are not universal. 2012 ACM Subject Classification Theory of computation → Computational geometry; Mathematics of computing → Algebraic topology Keywords and phrases Reeb graphs, topological descriptors, edit distance, interleaving distance Digital Object Identifier 10.4230/LIPIcs.SoCG.2020.15 Funding This research has been partially supported by FAR 2014 (UniMORE), ARCES (University of Bologna), and the DFG Collaborative Research Center SFB/TRR 109 “Discretization in Geometry and Dynamics”. Acknowledgements We thank Barbara Di Fabio and Yusu Wang for valuable discussions. 1 Introduction The concept of a Reeb graph of a Morse function first appeared in [13] and has subsequently been applied to problems in shape analysis in [14, 10]. The literature on Reeb graphs in the computational geometry and computational topology is ever growing (see, e.g., [2, 3] for a discussion and references).
    [Show full text]
  • Realizable Piecewise Linear Paths of Persistence Diagrams with Reeb Graphs
    Realizable piecewise linear paths of persistence diagrams with Reeb graphs Rehab Alharbi1, Erin Wolf Chambers2, and Elizabeth Munch2 1Dept of Mathematics & Statistics, Saint Louis University 2Dept of Computer Science, Saint Louis University 3Dept of Computational Mathematics, Science & Engineering, Dept of Mathematics, Michigan State University Abstract Reeb graphs are widely used in a range of fields for the purposes of analyzing and com- paring complex spaces via a simpler combinatorial object. Further, they are closely related to extended persistence diagrams, which largely but not completely encode the information of the Reeb graph. In this paper, we investigate the effect on the persistence diagram of a particular continuous operation on Reeb graphs; namely the (truncated) smoothing operation. This con- struction arises in the context of the Reeb graph interleaving distance, but separately from that viewpoint provides a simplification of the Reeb graph which continuously shrinks small loops. We then use this characterization to initiate the study of inverse problems for Reeb graphs using smoothing by showing which paths in persistence diagram space (commonly known as vineyards) can be realized by a path in the space of Reeb graphs via these simple operations. This allows us to solve the inverse problem on a certain family of piecewise linear vineyards when fixing an initial Reeb graph. 1 Introduction Reeb graphs have become an important tool in topological data analysis for the purpose of visualiz- ing continuous functions on complex spaces, as they yield a simplified discrete structure. Originally developed in relation to Morse theory [29], these objects are used extensively for shape comparison, constructing skeletons of data sets, surface simplification, and visualization; for more details on arXiv:2107.04654v1 [cs.CG] 9 Jul 2021 these and more applications, we refer to recent surveys on the topic [5, 31].
    [Show full text]
  • Thermodynamic Tree: the Space of Admissible Paths † Alexander N
    SIAM J. APPLIED DYNAMICAL SYSTEMS c 2013 Society for Industrial and Applied Mathematics Vol. 12, No. 1, pp. 246–278 ∗ Thermodynamic Tree: The Space of Admissible Paths † Alexander N. Gorban Abstract. Is a spontaneous transition from a state x to a state y allowed by thermodynamics? Such a question arises often in chemical thermodynamics and kinetics. We ask the following more formal question: Is there a continuous path between these states, along which the conservation laws hold, the concentra- tions remain nonnegative, and the relevant thermodynamic potential G (Gibbs energy, for example) monotonically decreases? The obvious necessary condition, G(x) ≥ G(y), is not sufficient, and we construct the necessary and sufficient conditions. For example, it is impossible to overstep the equi- − librium in 1-dimensional (1D) systems (with n components and n 1 conservation laws). The system cannot come from a state x to a state y if they are on the opposite sides of the equilibrium even if G(x) >G(y). We find the general multidimensional analogue of this 1D rule and constructively solve the problem of the thermodynamically admissible transitions. We study dynamical systems, which are given in a positively invariant convex polyhedron D and have a convex Lyapunov function G. An admissible path is a continuous curve in D along which G does not increase. For x, y ∈ D, x y (x precedes y) if there exists an admissible path from x to y and x ∼ y if x y and y x. ThetreeofG in D is a quotient space D/ ∼ . We provide an algorithm for the construction of this tree.
    [Show full text]
  • Introduction to Topological Data Analysis Julien Tierny
    Introduction to Topological Data Analysis Julien Tierny To cite this version: Julien Tierny. Introduction to Topological Data Analysis. Doctoral. France. 2017. cel-01581941 HAL Id: cel-01581941 https://hal.archives-ouvertes.fr/cel-01581941 Submitted on 5 Sep 2017 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Sorbonne Universités, UPMC Univ Paris 06 Laboratoire d’Informatique de Paris 6 Julien Tierny Introduction to Topological Data Analysis Sorbonne Universités, UPMC Univ Paris 06 Laboratoire d’Informatique de Paris 6 UMR UPMC/CNRS 7606 – Tour 26 4 Place Jussieu – 75252 Paris Cedex 05 – France Notations X Topological space ¶X Boundary of a topological space M Manifold Rd Euclidean space of dimension d s, t d-simplex, face of a d-simplex v, e, t, T Vertex, edge, triangle and tetrahedron Lk(s), St(s) Link and star of a simplex Lkd(s), Std(s) d-simplices of the link and the star of a simplex K Simplicial complex T Triangulation M Piecewise linear manifold bi i-th Betti number c Euler characteristic th ai i barycentric coordinates of a point p relatively
    [Show full text]
  • Time-Varying Reeb Graphs: a Topological Framework Supporting the Analysis of Continuous Time-Varying Data
    Time-varying Reeb Graphs: A Topological Framework Supporting the Analysis of Continuous Time-varying Data by Ajith Arthur Mascarenhas A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial ful¯llment of the requirements for the degree of Doctor of Philosophy in the Department of Computer Science. Chapel Hill 2006 Approved by: Jack Snoeyink, Advisor Herbert Edelsbrunner, Reader Valerio Pascucci, Reader Dinesh Manocha, Committee Member Leonard McMillan, Committee Member ii iii c 2006 Ajith Arthur Mascarenhas ALL RIGHTS RESERVED iv v ABSTRACT AJITH ARTHUR MASCARENHAS: Time-varying Reeb Graphs: A Topological Framework Supporting the Analysis of Continuous Time-varying Data. (Under the direction of Jack Snoeyink.) I present time-varying Reeb graphs as a topological framework to support the anal- ysis of continuous time-varying data. Such data is captured in many studies, including computational uid dynamics, oceanography, medical imaging, and climate modeling, by measuring physical processes over time, or by modeling and simulating them on a computer. Analysis tools are applied to these data sets by scientists and engineers who seek to understand the underlying physical processes. A popular tool for analyzing scienti¯c datasets is level sets, which are the points in space with a ¯xed data value s. Displaying level sets allows the user to study their geometry, their topological features such as connected components, handles, and voids, and to study the evolution of these features for varying s. For static data, the Reeb graph encodes the evolution of topological features and compactly represents topological information of all level sets.
    [Show full text]
  • Introduction to Reeb Graphs and Contour Trees
    Introduction to Reeb Graphs and Contour Trees Lecture 15 Scribed by: ABHISEK KUNDU Sometimes we are interested in the topology of smooth functions as a means to analyze and visualize intrinsic properties of geometric models and scientific data. Creation and destruction of components of the level set could be an important visual technique for this purpose. Reeb graphs are obtained by contracting the connected components of the level sets to points. They could be useful visual analytical tool because they express the connectivity of level sets. DEFINITION: ( Level Set ) A level set of a real-valued function f of n variables is a set of the form: { (X1,...,Xn) | f (X1,...,Xn) = c} where c is a constant. In other words, it is the set where the function takes on a given constant value. When the number of variables is two, this is a level curve or contour line, if it is three a level surface, and for higher values of n the level set is a level hypersurface. Sublevel set: of f is a set of the form { (X1,...,Xn) | f (X1,...,Xn) ≤ c} Contour: at a given value c, contour is the connected component of level set Xc [Figure 1] Equivalence of points: two points x and y are equivalent if they belong to the same contour of Xc for some given constant c. [Figure 1] Quotient topology and space: Let us denote ~ to be an equivalence relation defined on a topological space X. e be the set of equivalence classes and let ψ : map each point x to its equivalence class DEFINITION: Quotient topology of consists of all subsets U belonging t whose preimages, ψ -1 (U), are open in X.
    [Show full text]
  • Reeb Graphs: Approximation and Persistence
    Reeb Graphs: Approximation and Persistence Tamal K. Dey∗ Yusu Wang∗ Abstract Given a continuous function f : X → IR on a topological space X, its level set f −1(a) changes continuously as the real value a changes. Consequently, the connected components in the level sets appear, disappear, split and merge. The Reeb graph of f summarizes this information into a graph structure. Previous work on Reeb graph mainly focused on its efficient computation. In this paper, we initiate the study of two important aspects of the Reeb graph which can facilitate its broader applications in shape and data analysis. The first one is the approximation of the Reeb graph of a function on a smooth compact manifold M without boundary. The approximation is computed from a set of points P sampled from M. By leveraging a relation between the Reeb graph and the so-called vertical homology group, as well as between cycles in M and in a Rips complex constructed from P , we compute the H1-homology of the Reeb graph from P . It takes O(n log n) expected time, where n is the size of the 2-skeleton of the Rips complex. As a by-product, when M is an orientable 2-manifold, we also obtain an efficient near-linear time (expected) algorithm for computing the rank of H1(M) from point data. The best known previous algorithm for this problem takes O(n3) time for point data. The second aspect concerns the definition and computation of the persistent Reeb graph homology for a sequence of Reeb graphs defined on a filtered space.
    [Show full text]
  • Time-Varying Reeb Graphs: a Topological Framework Supporting the Analysis of Continuous Time-Varying Data
    UCRL-TH-226559 Time-varying Reeb Graphs: A Topological Framework Supporting the Analysis of Continuous Time-varying Data A. Mascarenhas December 6, 2006 Disclaimer This document was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor the University of California nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or the University of California, and shall not be used for advertising or product endorsement purposes. This work was performed under the auspices of the U.S. Department of Energy by University of California, Lawrence Livermore National Laboratory under Contract W-7405-Eng-48. UCRL-TH-226559 Time-varying Reeb Graphs: A Topological Framework Supporting the Analysis of Continuous Time-varying Data by Ajith Arthur Mascarenhas A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial ful¯llment of the requirements for the degree of Doctor of Philosophy in the Department of Computer Science. Chapel Hill 2006 Approved by: Jack Snoeyink, Advisor Herbert Edelsbrunner, Reader Valerio Pascucci, Reader Dinesh Manocha, Committee Member Leonard McMillan, Committee Member ii iii c 2006 Ajith Arthur Mascarenhas ALL RIGHTS RESERVED iv v ABSTRACT AJITH ARTHUR MASCARENHAS: Time-varying Reeb Graphs: A Topological Framework Supporting the Analysis of Continuous Time-varying Data.
    [Show full text]
  • Categorification of Reeb Graphs
    Categorification of Reeb Graphs Elizabeth Munch University of Minnesota :: Institute for Mathematics and Its Applications February 5, 2014 Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 1 / 46 What's the plan? Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 2 / 46 What's the plan? Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 2 / 46 What's the plan? Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 2 / 46 Original Reeb Graph construction Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 3 / 46 Original Reeb Graph construction Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 3 / 46 Original Reeb Graph Definition Definition Given topological space X with a function f : X ! R, say x ∼ y if x and y are in the same connected component of f −1(a). The Reeb graph of the function f is the space X= ∼ with the quotient topology. Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 4 / 46 Morse Function Setting Let f be Morse with distinct critical values. critical values , Nodes Index of critical values , Type of node I Min/max , Deg 1 vertex I Saddle , Deg 3 vertex Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 5 / 46 Applications Handle removal, from Wood 3D shape matching, from et al, 2002 Hilaga et.al, 2001 Basis for homology groups, Shape skeletonization, from Dey et. al, 2013 from Biasotti et.al, 2008 Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 6 / 46 Smoothing Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 7 / 46 Smoothing Liz Munch (IMA) Reeb Graphs - SAMSI February 5, 2014 7 / 46 Other work Functional distortion distance [Bauer Ge Wang 2013] Topological simplification [Doraiswamy Natarajan 2012], [Ge et at 2011], [Pascucci et al 2007] Smoothing Goals Give a metric to compare Reeb graphs.
    [Show full text]
  • V.4 Reeb Graphs
    122 V Morse Functions V.4 Reeb Graphs The structure of a continuous function can sometimes be made explicit by plotting the evolution of the components of the level set. This leads to the concept of the Reeb graph of the function. It has applications in medical imaging and other areas of science and engineering. Iso-surface extraction. The practical motivation for studying Reeb graphs is the extraction of iso-surfaces for three-dimensional density data. In topolog- ical lingo, the density data is a continuous function, f : [0; 1]3 R, and an 1 ! iso-surface is a level set, f − (t). If f is smooth and t is a regular value then the level set is a 2-manifold, possibly with boundary. Similarly, if f is generic PL and t is not the value of a PL critical point then the level set is a 2-manifold, again possibly with boundary. Figure V.12 illustrates this fact for a generic PL function on the unit square. Assuming we enter a triangle at a boundary Figure V.12: The level set of a generic PL function on a triangulation of the unit square. The superlevel set is white and the sublevel set is shaded. point x with f(x) = t, there is a unique other boundary point y with f(y) = t where we exit the triangle. We draw the line segment from x to y as part of the level set and repeat the construction by entering the next triangle at y. The procedure is similar for a PL function on the unit cube except that we use a graph search algorithm to collect the triangular and quadrangular surface pieces we get for the tetrahedra.
    [Show full text]