New Results on Abstract Voronoi Diagrams

Voronoi diagrams are a fundamental structure used in many areas of science. For a given set of objects, called sites, the Voronoi diagram separates the plane into regions, such that points belonging to the same region have got the same nearest site. This definition clearly depends on the type of given objects, they may be points, line segments, polygons, etc. and the distance measure used. To free oneself from these geometric notions, Klein introduced abstract Voronoi diagrams as a general construct covering many concrete Voronoi diagrams. Abstract Voronoi diagrams are based on a system of bisecting curves, one for each pair of abstract sites, separating the plane into two dominance regions, belonging to one site each. The intersection of all dominance regions belonging to one site p defines its Voronoi region. The system of bisecting curves is required to fulfill only some simple combinatorial properties, like Voronoi regions to be connected, the union of their closures cover the whole plane, and the bisecting curves are unbounded. These assumptions are enough to show that an abstract Voronoi diagram of n sites is a planar graph of complexity O(n) and can be computed in expected time O(n log n) by a randomized incremental construction. In this thesis we widen the notion of abstract Voronoi diagrams in several senses. One step is to allow disconnected Voronoi regions. We assume that in a diagram of a subset of three sites each Voronoi region may consist of at most s connected components, for a constant s, and show that the diagram can be constructed in expected time O(s2 n ∑3 ≤ j ≤ n mj / j), where mj is the expected number of connected components of a Voronoi region over all diagrams of a subset of j sites. The case that all Voronoi regions are connected is a subcase, where this algorithm performs in optimal O(n log n) time, because here s = mj =1. The next step is to additionally allow bisecting curves to be closed. We present an algorithm constructing such diagrams which runs in expected time O(s2 n log(max{s,n}) ∑2 ≤ j≤ n mj / j). This algorithm is slower by a log n-factor compared to the one for disconnected regions and unbounded bisectors. The extra time is necessary to be able to handle special phenomenons like islands, where a Voronoi region is completely surrounded by another region, something that can occur only when bisectors are closed. However, this algorithm solves many open problems and improves the running time of some existing algorithms, for example for the farthest Voronoi diagram of n simple polygons of constant complexity. Another challenge was to study higher order abstract Voronoi diagrams. In the concrete sense of an order-k Voronoi diagram points are collected in the same Voronoi region, if they have the same k nearest sites. By suitably intersecting the dominance regions this can be defined also for abstract Voronoi diagrams. The question arising is about the complexity of an order-k Voronoi diagram. There are many subsets of size k but fortunately many of them have an empty order-k region. For point sites it has already been shown that there can be at most O(k (n-k)) many regions and even though order-k regions may be disconnected when considering line segments, still the complexity of the order-k diagram remains O(k(n-k)). The proofs used to show this strongly depended on the geometry of the sites and the distance measure, and were thus not applicable for our abstract higher order Voronoi diagrams. The proofs used to show this strongly depended on the geometry of the sites and the distance measure, and were thus not applicable for our abstract higher order Voronoi diagrams. Nevertheless, we were able to come up with proofs of purely topological and combinatorial nature of Jordan curves and certain permutation sequences, and hence we could show that also the order-k abstract Voronoi diagram has complexity O(k (n-k)), assuming that bisectors are unbounded, and the order-1 regions are connected. Finally, we discuss Voronoi diagrams having the shape of a tree or forest. Aggarwal et. al. showed that if points are in convex position, then given their ordering along the convex hull, their Voronoi diagram, which is a tree, can be computed in linear time. Klein and Lingas have generalized this idea to Hamiltonian abstract Voronoi diagrams, where a curve is given, intersecting each Voronoi region with respect to any subset of sites exactly once. If the ordering of the regions along the curve is known in advance, all Voronoi regions are connected, and all bisectors are unbounded, then the abstract Voronoi diagram can be computed in linear time. This algorithm also applies to diagrams which are trees for all subsets of sites and the ordering of the unbounded regions around the diagram is known. In this thesis we go one step further and allow the diagram to be a forest for subsets of sites as long as the complete diagram is a tree. We show that also these diagrams can be computed in linear time

Computational Geometry Column 42

A compendium of thirty previously published open problems in computational geometry is presented.Comment: 7 pages; 72 reference

Asymptotic statistics of the n-sided planar Poisson-Voronoi cell. I. Exact results

We achieve a detailed understanding of the $n$-sided planar Poisson-Voronoi cell in the limit of large $n$. Let ${p}\_n$ be the probability for a cell to have $n$ sides. We construct the asymptotic expansion of $\log {p}\_n$ up to terms that vanish as $n\to\infty$. We obtain the statistics of the lengths of the perimeter segments and of the angles between adjoining segments: to leading order as $n\to\infty$, and after appropriate scaling, these become independent random variables whose laws we determine; and to next order in $1/n$ they have nontrivial long range correlations whose expressions we provide. The $n$-sided cell tends towards a circle of radius (n/4\pi\lambda)^{\half}, where $\lambda$ is the cell density; hence Lewis' law for the average area $A\_n$ of the $n$-sided cell behaves as $A\_n \simeq cn/\lambda$ with $c=1/4$. For $n\to\infty$ the cell perimeter, expressed as a function $R(\phi)$ of the polar angle $\phi$, satisfies $d^2 R/d\phi^2 = F(\phi)$, where $F$ is known Gaussian noise; we deduce from it the probability law for the perimeter's long wavelength deviations from circularity. Many other quantities related to the asymptotic cell shape become accessible to calculation.Comment: 54 pages, 3 figure

Spatial Aggregation: Theory and Applications

Visual thinking plays an important role in scientific reasoning. Based on the research in automating diverse reasoning tasks about dynamical systems, nonlinear controllers, kinematic mechanisms, and fluid motion, we have identified a style of visual thinking, imagistic reasoning. Imagistic reasoning organizes computations around image-like, analogue representations so that perceptual and symbolic operations can be brought to bear to infer structure and behavior. Programs incorporating imagistic reasoning have been shown to perform at an expert level in domains that defy current analytic or numerical methods. We have developed a computational paradigm, spatial aggregation, to unify the description of a class of imagistic problem solvers. A program written in this paradigm has the following properties. It takes a continuous field and optional objective functions as input, and produces high-level descriptions of structure, behavior, or control actions. It computes a multi-layer of intermediate representations, called spatial aggregates, by forming equivalence classes and adjacency relations. It employs a small set of generic operators such as aggregation, classification, and localization to perform bidirectional mapping between the information-rich field and successively more abstract spatial aggregates. It uses a data structure, the neighborhood graph, as a common interface to modularize computations. To illustrate our theory, we describe the computational structure of three implemented problem solvers -- KAM, MAPS, and HIPAIR --- in terms of the spatial aggregation generic operators by mixing and matching a library of commonly used routines.Comment: See http://www.jair.org/ for any accompanying file

Multifractal analysis of electronic states on random Voronoi-Delaunay lattices

We consider the transport of non-interacting electrons on two- and three-dimensional random Voronoi-Delaunay lattices. It was recently shown that these topologically disordered lattices feature strong disorder anticorrelations between the coordination numbers that qualitatively change the properties of continuous and first-order phase transitions. To determine whether or not these unusual features also influence Anderson localization, we study the electronic wave functions by multifractal analysis and finite-size scaling. We observe only localized states for all energies in the two-dimensional system. In three dimensions, we find two Anderson transitions between localized and extended states very close to the band edges. The critical exponent of the localization length is about 1.6. All these results agree with the usual orthogonal universality class. Additional generic energetic randomness introduced via random potentials does not lead to qualitative changes but allows us to obtain a phase diagram by varying the strength of these potentials

Correcting curvature-density effects in the Hamilton-Jacobi skeleton

The Hainilton-Jacobi approach has proven to be a powerful and elegant method for extracting the skeleton of two-dimensional (2-D) shapes. The approach is based on the observation that the normalized flux associated with the inward evolution of the object boundary at nonskeletal points tends to zero as the size of the integration area tends to zero, while the flux is negative at the locations of skeletal points. Nonetheless, the error in calculating the flux on the image lattice is both limited by the pixel resolution and also proportional to the curvature of the boundary evolution front and, hence, unbounded near endpoints. This makes the exact location of endpoints difficult and renders the performance of the skeleton extraction algorithm dependent on a threshold parameter. This problem can be overcome by using interpolation techniques to calculate the flux with subpixel precision. However, here, we develop a method for 2-D skeleton extraction that circumvents the problem by eliminating the curvature contribution to the error. This is done by taking into account variations of density due to boundary curvature. This yields a skeletonization algorithm that gives both better localization and less susceptibility to boundary noise and parameter choice than the Hamilton-Jacobi method