7437367
Vines and vineyards by updating persistence in linear time
Persistent homology is the mathematical core of recent work on shape, including reconstruction, recognition, and matching. Its pertinent information is encapsulated by a pairing of the critical values of a function, visualized by points forming a diagram in the plane. The original algorithm in [10] computes the pairs from an ordering of the simplices in a triangulation and takes worst-case time cubic in the number of simplices. The main result of this paper is an algorithm that maintains the pairing in worst-case linear time per transposition in the ordering. A side-effect of the algorithm's analysis is an elementary proof of the stability of persistence diagrams [7] in the special case of piecewise-linear functions. We use the algorithm to compute 1-parameter families of diagrams which we apply to the study of protein folding trajectories.
14707519
Inequalities for the curvature of curves and surfaces
In this paper, we bound the difference between the total mean curvatures of two closed surfaces in R3 in terms of their total absolute curvatures and the Fréchet distance between the volumes they enclose. The proof relies on a combination of methods from algebraic topology and integral geometry. We also bound the difference between the lengths of two curves using the same methods.
10864618
Stability of persistence diagrams
The persistence diagram of a real-valued function on a topological space is a multiset of points in the extended plane. We prove that under mild assumptions on the function, the persistence diagram is stable: small changes in the function imply only small changes in the diagram. We apply this result to estimating the homology of sets in a metric space and to comparing and classifying geometric shapes.
17520762
Extraction and simplification of iso-surfaces in tandem
The tandem algorithm combines the marching cube algorithm for surface extraction and the edge contraction algorithm for surface simplification in lock-step to avoid the costly intermediate step of storing the entire extracted surface triangulation. Beyond this basic strategy, we introduce refinements to prevent artifacts in the resulting triangulation, first, by carefully monitoring the amount of simplification during the process and, second, by driving the simplification toward a compromise between shape approximation and mesh quality. We have implemented the algorithm and used extensive computational experiments to document the effects of various design options and to further fine-tune the algorithm.
6301471
Adaptive Simplicial Grids from Cross-Sections of Monotone Complexes
We study the maintenance of a simplicial grid or complex under changing density requirements. The proposed method works in any fixed dimension and generates grids by projecting cross-sections of a monotone simplicial complex that lives in one dimension higher than the grid. The density of the grid is adapted by locally moving the cross-section up or down along the extra dimension.
7249846
Space Searching for Intersecting Objects
Abstract Determining or counting geometric objects that intersect another geometric query object is at the core of algorithmic problems in a number of applied areas of computer science. This article presents a family of space-efficient data structures that realize sublinear query time for points, line segments, lines and polygons in the plane, and points, line segments, planes, and polyhedra in three dimensions.
17780302
Triangulating topological spaces
Given a subspace<inline-equation><f><blkbd>X⊆R</blkbd><sup>d</sup></f></inline-equation> and a finite set <inline-equation><f>S⊆<blkbd>R</blkbd><sup>d</sup></f></inline-equation>, we introduce the Delaunay simplicial complex, <inline-equation><f><sc>D</sc><inf><blkbd>X</blkbd></inf></f></inline-equation>, restricted by <inline-equation><f><blkbd>X</blkbd></f></inline-equation>. Its simplices are spanned by subsets<inline-equation><f>T⊆S</f></inline-equation> for which the common intersection of Voronoi cells meets <inline-equation><f><blkbd>X</blkbd></f></inline-equation> in a non-empty set. By the nerve theorem,<inline-equation><f><lim align="c"><op>⋃</op></lim><sc>D</sc><inf><blkbd>X</blkbd></inf></f></inline-equation> and <inline-equation><f><blkbd>X</blkbd></f></inline-equation> are homotopy equivalent if all such sets are contractible. This paper shows that <inline-equation><f><lim align="c"><op>⋃</op></lim><sc>D</sc><inf><blkbd>X</blkbd></inf></f></inline-equation> and <inline-equation><f><blkbd>X</blkbd></f></inline-equation> <?Pub Caret1>are homeomorphic if the sets can be further subdivided in a certain way so they form a regular CW complex.
12582971
Stabbing line segments
An algorithm for the geometric problem of determining a line (called a stabbing line) which intersects each ofn given line segments in the plane is presented. As a matter of fact, the algorithm computes a description of all stabbing lines. A purely geometric fact is proved which infers that this description requiresO(n) space to be specified. Our algorithm computes it inO(n logn) time which is optimal in the worst case.Using the description of the stabbing lines, we are able to decide inO(logn) time whether or not a specified line is a stabbing line. Finally, the problem of maintaining the description of all stabbing lines while inserting and deleting line segments is addressed.
30538700
Computing a Ham-Sandwich Cut in Two Dimensions
Let B be a set of n"b black points and W a set of n"w, white points in the Euclidean plane. A line h is said to bisect B (or W) if, at most, half of the points of B (or W) lie on any one side of h. A line that bisects both B and W is called a ham-sandwich cut of B and W. We give an algorithm that computes a ham-sandwich cut of B and W in 0((n"h+n"w) log (min {n"b, n"w}+ 1)) time. The algorithm is considerably simpler than the previous most efficient one which takes 0((n"b + n"w) log (n"b + n"w)) time.