A Hierarchical Approach for Regular Centroidal Voronoi Tessellations

International audienceIn this paper we consider Centroidal Voronoi Tessellations (CVTs) and study their regularity. CVTs are geometric structures that enable regular tessellations of geometric objects and are widely used in shape modeling and analysis. While several efficient iterative schemes, with defined local convergence properties, have been proposed to compute CVTs, little attention has been paid to the evaluation of the resulting cell decompositions. In this paper, we propose a regularity criterion that allows us to evaluate and compare CVTs independently of their sizes and of their cell numbers. This criterion allows us to compare CVTs on a common basis. It builds on earlier theoretical work showing that second moments of cells converge to a lower bound when optimising CVTs. In addition to proposing a regularity criterion, this paper also considers computational strategies to determine regular CVTs. We introduce a hierarchical framework that propagates regularity over decomposition levels and hence provides CVTs with provably better regularities than existing methods. We illustrate these principles with a wide range of experiments on synthetic and real models

A PDE approach to centroidal tessellations of domains

We introduce a class of systems of Hamilton-Jacobi equations that characterize critical points of functionals associated to centroidal tessellations of domains, i.e. tessellations where generators and centroids coincide, such as centroidal Voronoi tessellations and centroidal power diagrams. An appropriate version of the Lloyd algorithm, combined with a Fast Marching method on unstructured grids for the Hamilton-Jacobi equation, allows computing the solution of the system. We propose various numerical examples to illustrate the features of the technique

Microstructure modeling and crystal plasticity parameter identification for predicting the cyclic mechanical behavior of polycrystalline metals

Computational homogenization permits to capture the influence of the microstructure on the cyclic mechanical behavior of polycrystalline metals. In this work we investigate methods to compute Laguerre tessellations as computational cells of polycrystalline microstructures, propose a new method to assign crystallographic orientations to the Laguerre cells and use Bayesian optimization to find suitable parameters for the underlying micromechanical model from macroscopic experiments

Layered Fields for Natural Tessellations on Surfaces

Mimicking natural tessellation patterns is a fascinating multi-disciplinary problem. Geometric methods aiming at reproducing such partitions on surface meshes are commonly based on the Voronoi model and its variants, and are often faced with challenging issues such as metric estimation, geometric, topological complications, and most critically parallelization. In this paper, we introduce an alternate model which may be of value for resolving these issues. We drop the assumption that regions need to be separated by lines. Instead, we regard region boundaries as narrow bands and we model the partition as a set of smooth functions layered over the surface. Given an initial set of seeds or regions, the partition emerges as the solution of a time dependent set of partial differential equations describing concurrently evolving fronts on the surface. Our solution does not require geodesic estimation, elaborate numerical solvers, or complicated bookkeeping data structures. The cost per time-iteration is dominated by the multiplication and addition of two sparse matrices. Extension of our approach in a Lloyd's algorithm fashion can be easily achieved and the extraction of the dual mesh can be conveniently preformed in parallel through matrix algebra. As our approach relies mainly on basic linear algebra kernels, it lends itself to efficient implementation on modern graphics hardware.Comment: Natural tessellations, surface fields, Voronoi diagrams, Lloyd's algorith