ADS: Random Sampling of Occupancy Functions using Adaptive Delaunay Scaffolding

Dense random sampling and surfacing of shapes encoded via implicit occupancy functions (OFs) are critical elements of many applications. Existing methods largely provide either one or the other of random sampling or mesh surfaces: ray shooting approaches deliver random samples with no connectivity, and grid-based methods deliver mesh surfaces but their sampling is highly biased. We propose a new method which delivers both pseudo-random OF surface samples and an isosurface mesh connecting them. Our method achieves these goals while requiring an order of magnitude fewer function evaluations than prior approaches. Key to our Adaptive Delaunay Sampling (ADS) approach is a progressively computed Delaunay tetrahedralization of points in 3D space, which we use as a sampling and surfacing scaffold. Starting from an initial coarse Delaunay scaffold, we repeatedly refine crossing edges, ones whose end vertices lie on opposite sides of the surface, augmenting the scaffold with points closer and closer to the surface. Each refinement step uses the Delaunay criterion to incorporate the newly added vertices into the scaffold, introducing new crossing edges. We use the intersections of fine crossing edges with the OF surface as the output samples, and use the marching tetrahedra method to surface these samples. We subsequently use normal estimation to densify the sampling near fine features and in areas of high surface curvature. We validate ADS by sampling 150 inputs at different resolutions, and provide extensive comparisons to existing alternatives. Our experiments demonstrate significant improvement in accuracy/function evaluation count trade-off, and showcase downstream applications.