A Novel Point Inclusion Test for Convex Polygons Based on Voronoi Tessellations

The point inclusion tests for polygons, in other words the point in polygon (PIP) algorithms are fundamental tools for many scientific fields related to computational geometry and they have been studied for a long time. PIP algorithms get direct or indirect geometric definition of a polygonal entity and validate its containment of a given point. The PIP algorithms which are working directly on the geometric entities derive linear boundary definitions for the edges of the polygon. Moreover, almost all direct test methods rely on the two point form of the line equation to partition the space into half-spaces. Voronoi tessellations use an alternate approach for half-space partitioning. Instead of line equation, distance comparison between generator points is used to accomplish the same task. Voronoi tessellations consist of convex polygons which are defined between generator points. Therefore, Voronoi tessellations have become an inspiration for us to develop a new approach of PIP testing specialized for convex polygons. Essential equations to the conversion of a convex polygon to a voronoi polygon are derived along this paper. As a reference, a very standard convex PIP testing algorithm, the sign of offset, is selected for comparison. For generalization of the comparisons the ray crossing algorithm is used as another reference. All algorithms are implemented as vector and matrix operations without any branching. This enabled us to benefit from the CPU optimizations of the underlying linear algebra libraries. All algorithms are tested for three different polygon sizes and varying point batch sizes. Overall, our proposed algorithm has performed better with varying margin between 10% to 23% comparing to the reference methods.