Anisotropic Surface Remeshing without Obtuse Angles

We present a novel anisotropic surface remeshing method that can efficiently eliminate obtuse angles. Unlike previous work that can only suppress obtuse angles with expensive resampling and Lloyd-type iterations, our method relies on a simple yet efficient connectivity and geometry refinement, which can not only remove all the obtuse angles, but also preserves the original mesh connectivity as much as possible. Our method can be directly used as a post-processing step for anisotropic meshes generated from existing algorithms to improve mesh quality. We evaluate our method by testing on a variety of meshes with different geometry and topology, and comparing with representative prior work. The results demonstrate the effectiveness and efficiency of our approach

Surface segmentation for improved remeshing

Many remeshing techniques sample the input surface in a meaningful way and then triangulate the samples to produce an output triangulated mesh. One class of methods samples in a parametrization of the surface. Another class samples directly on the surface. These latter methods must have sufficient density of samples to achieve outputs that are homeomorphic to the input. In many datasets samples must be very dense even in some nearly planar regions due to small local feature size. We present an isotropic remeshing algorithm called κCVT that achieves topological correctness while sampling sparsely in all flat regions, regardless of local feature size. This is accomplished by segmenting the surface, remeshing the segmented subsurfaces individually and then stitching them back together. We show that κCVT produces quality meshes using fewer triangles than other methods. The output quality meshes are both homeomorphic and geometrically close to the input surface.postprin

Toward mixed-element meshing based on restricted Voronoi diagrams

In this paper we propose a method to generate mixed-element meshes (tetrahedra, triangular prisms, square pyramids) for B-Rep models. The vertices, edges, facets, and cells of the final volumetric mesh are determined from the combinatorial analysis of the intersections between the model components and the Voronoi diagram of sites distributed to sample the model. Inside the volumetric regions, Delaunay tetrahedra dual of the Voronoi diagram are built. Where the intersections of the Voronoi cells with the model surfaces have a unique connected component, tetrahedra are modified to fit the input triangulated surfaces. Where these intersections are more complicated, a correspondence between the elements of the Voronoi diagram and the elements of the mixedelement mesh is used to build the final volumetric mesh. The method which was motivated by meshing challenges encountered in geological modeling is demonstrated on several 3D synthetic models of subsurface rock volumes

Isotropic Surface Remeshing

International audienceThis paper proposes a new method for isotropic remeshing of tri- angulated surface meshes. Given a triangulated surface mesh to be resampled and a user-specified density function defined over it, we first distribute the desired number of samples by generalizing error diffusion, commonly used in image halftoning, to work directly on mesh triangles and feature edges. We then use the resulting sam- pling as an initial configuration for building a weighted centroidal Voronoi tessellation in a conformal parameter space, where the specified density function is used for weighting. We finally create the mesh by lifting the corresponding constrained Delaunay trian- gulation from parameter space. A precise control over the sampling is obtained through a flexible design of the density function, the latter being possibly low-pass filtered to obtain a smoother grada- tion. We demonstrate the versatility of our approach through vari- ous remeshing examples

A feature extracting and meshing approach for sheet-like structures in rocks

Meshing rock samples with sheet-like structures based their CT scanned volumetric images, is a crucial component for both visualization and numerical simulation. In rocks, fractures and veins commonly exist in the form of sheet-like objects (e.g. thin layers and distinct flat shapes), which are much smaller than the rock mass dimensions. The representations of such objects require high-resolution 3D images with a huge dataset, which are difficult and even impossible to visualize or analyze by numerical methods. Therefore, we develop a microscopic image based meshing approach to extract major sheet-like structures and then preserve their major geometric features at the macroscale. This is achieved by the following four major steps: (1) extracting major objects through extending, separation and recovering operations based on the CT scanned data/microscopic images; (2) simplifying and constructing a simplified centroidal Voronoi diagram on the extracted structures; (3) generating triangular meshes to represent the structure; (4) generating volume tetrahedron meshes constrained with the above surface mesh as the internal surfaces. Moreover, a shape similarity approach is proposed to measure and evaluate how similar the generated mesh models to the original rock samples. It is applied as criteria for further mesh generation to better describe the rock features with fewer elements. Finally, a practical CT scanned rock is taken as an application example to demonstrate the usefulness and capability of the proposed approach

Lp Centroidal Voronoi Tesselation and its applications

International audienceThis paper introduces Lp -Centroidal Voronoi Tessellation (Lp -CVT), a generalization of CVT that minimizes a higher-order moment of the coordinates on the Voronoi cells. This generalization allows for aligning the axes of the Voronoi cells with a predefined background tensor field (anisotropy). Lp -CVT is computed by a quasi-Newton optimization framework, based on closed-form derivations of the objective function and its gradient. The derivations are given for both surface meshing (Ω is a triangulated mesh with per-facet anisotropy) and volume meshing (Ω is the interior of a closed triangulated mesh with a 3D anisotropy field). Applications to anisotropic, quad-dominant surface remeshing and to hex-dominant volume meshing are presented. Unlike previous work, Lp -CVT captures sharp features and intersections without requiring any pre-tagging

A unified framework for isotropic meshing based on narrow-band Euclidean distance transformation

In this paper, we propose a simple-yet-effective method for isotropic meshing relying on Euclidean distance transformation based centroidal Voronoi tessellation (CVT). Our approach improves the performance and robustness of computing CVT on curved domains while simultaneously providing high-quality output meshes. While conventional extrinsic methods compute CVTs in the entire volume bounded by the input model, we restrict the computation to a 3D shell of user-controlled thickness. Taking voxels which contain surface samples as sites, we compute the exact Euclidean distance transform on the GPU. Our algorithm is parallel and memory-efficient, and can construct the shell space for resolutions up to 20483 at interactive speed. The 3D centroidal Voronoi tessellation and restricted Voronoi diagrams are also computed efficiently on the GPU. Since the shell space can bridge holes and gaps smaller than a certain tolerance, and tolerate non-manifold edges and degenerate triangles, our algorithm can handle models with such defects, which typically cause conventional remeshing methods to fail. Our method can process implicit surfaces, polyhedral surfaces, and point clouds in a unified framework. Computational results show that our GPU-based isotropic meshing algorithm produces results comparable to state-of- the-art techniques, but is significantly faster than conventional CPU-based implementations.MOE (Min. of Education, S’pore)Published versio

Anisotropic geometry-conforming d-simplicial meshing via isometric embeddings

We develop a dimension-independent, Delaunay-based anisotropic mesh generation algorithm suitable for integration with adaptive numerical solvers. As such, the mesh produced by our algorithm conforms to an anisotropic metric prescribed by the solver as well as the domain geometry, given as a piecewise smooth complex. Motivated by the work of Lévy and Dassi [10-12,20], we use a discrete manifold embedding algorithm to transform the anisotropic problem to a uniform one. This work differs from previous approaches in several ways. First, the embedding algorithm is driven by a Riemannian metric field instead of the Gauss map, lending itself to general anisotropic mesh generation problems. Second we describe our method for computing restricted Voronoi diagrams in a dimension-independent manner which is used to compute constrained centroidal Voronoi tessellations. In particular, we compute restricted Voronoi simplices using exact arithmetic and use data structures based on convex polytope theory. Finally, since adaptive solvers require geometry-conforming meshes, we offer a Steiner vertex insertion algorithm for ensuring the extracted dual Delaunay triangulation is homeomorphic to the input geometries. The two major contributions of this paper are: a method for isometrically embedding arbitrary mesh-metric pairs in higher dimensional Euclidean spaces and a dimension-independent vertex insertion algorithm for producing geometry-conforming Delaunay meshes. The former is demonstrated on a two-dimensional anisotropic problem whereas the latter is demonstrated on both 3d and 4d problems. Keywords: Anisotropic mesh generation; metric; Nash embedding theorem; isometric; geometry-conforming; restricted Voronoi diagram; constrained centroidal Voronoi tessellation; Steiner vertices; dimension-independen

Differentiable Surface Triangulation

Triangle meshes remain the most popular data representation for surface geometry. This ubiquitous representation is essentially a hybrid one that decouples continuous vertex locations from the discrete topological triangulation. Unfortunately, the combinatorial nature of the triangulation prevents taking derivatives over the space of possible meshings of any given surface. As a result, to date, mesh processing and optimization techniques have been unable to truly take advantage of modular gradient descent components of modern optimization frameworks. In this work, we present a differentiable surface triangulation that enables optimization for any per-vertex or per-face differentiable objective function over the space of underlying surface triangulations. Our method builds on the result that any 2D triangulation can be achieved by a suitably perturbed weighted Delaunay triangulation. We translate this result into a computational algorithm by proposing a soft relaxation of the classical weighted Delaunay triangulation and optimizing over vertex weights and vertex locations. We extend the algorithm to 3D by decomposing shapes into developable sets and differentiably meshing each set with suitable boundary constraints. We demonstrate the efficacy of our method on various planar and surface meshes on a range of difficult-to-optimize objective functions. Our code can be found online: https://github.com/mrakotosaon/diff-surface-triangulation