Kinetic Stable Delaunay Graphs

TL;DR

This paper introduces a stable Delaunay graph for moving points that is less volatile than the full triangulation, maintains many useful properties, and can be efficiently maintained with kinetic data structures.

Contribution

It defines the stable Delaunay graph, proves its properties, and presents kinetic data structures for its efficient maintenance during point motion.

Findings

Contains at least one third of the Delaunay edges.

Includes the $eta$-skeleton for $eta=1+ ext{Omega}( ext{alpha}^2)$.

Remains stable under small changes in point positions and norm approximations.

Abstract

We consider the problem of maintaining the Euclidean Delaunay triangulation $\DT$ of a set $P$ of $n$ moving points in the plane, along algebraic trajectories of constant description complexity. Since the best known upper bound on the number of topological changes in the full $\DT$ is nearly cubic, we seek to maintain a suitable portion of it that is less volatile yet retains many useful properties. We introduce the notion of a stable Delaunay graph, which is a dynamic subgraph of the Delaunay triangulation. The stable Delaunay graph (a) is easy to define, (b) experiences only a nearly quadratic number of discrete changes, (c) is robust under small changes of the norm, and (d) possesses certain useful properties. The stable Delaunay graph ($\SDG$ in short) is defined in terms of a parameter $\alpha>0$, and consists of Delaunay edges $pq$ for which the angles at which $p$ and $q$ see their Voronoi edge $e_{pq}$ are at least $\alpha$. We show that (i) $\SDG$ always contains at least roughly one third of the Delaunay edges; (ii) it contains the $\beta$-skeleton of $P$, for $\beta=1+\Omega(\alpha^2)$; (iii) it is stable, in the sense that its edges survive for long periods of time, as long as the orientations of the segments connecting (nearby) points of $P$ do not change by much; and (iv) stable Delaunay edges remain stable (with an appropriate redefinition of stability) if we replace the Euclidean norm by any sufficiently close norm. In particular, we can approximate the Euclidean norm by a polygonal norm (namely, a regular $k$-gon, with $k=\Theta(1/\alpha)$), and keep track of a Euclidean $\SDG$ by maintaining the full Delaunay triangulation of $P$ under the polygonal norm. We describe two kinetic data structures for maintaining $\SDG$. Both structures use $O^*(n)$ storage and process $O^*(n^2)$ events during the motion, each in $O^*(1)$ time.