CGAL 6.1 - Kinetic Space Partition
|
This CGAL component implements the kinetic space partition proposed by Bauchet et. al [1]. It takes as input a set of non-coplanar convex polygons and partitions the bounding box of the input into polyhedra, where each polyhedron and its facets are convex. Each facet of the partition is part of the input polygons or an extension of them.
The kinetic partition homogeneously expands each input polygon along its support plane until collisions occur between them. At each collision, their expansion may stop, they may restrict the propagation along an intersection line between two support planes or they may continue beyond the intersection line. The polygons are split at the beginning of the kinetic simulation and at each collision along the intersection line to be intersection-free.
Whether a polygon is expanded beyond an intersection with another polygon, depends on the user specified k
parameter. Choosing k = 1
will cause the expansion of polygons to locally stop at the first intersection with another polygon and choosing k
equal to the number of input polygons will lead to a full expansion of each polygon to the bounding box.
The first step of the method creates a plane arrangement between the support planes of the input polygons. The computation of a plane arrangement has \(O(n^3)\). To speed up the computation the input data can be decomposed into separate kinetic partitions using an octree. The decomposition limits the expansion of an input polygon to octree nodes it initially intersects. The usage of an octree significantly speeds up the creation of the kinetic partition. However, it adds facets to the partition which originate from the octree and do not belong to an input polygon. The plane arrangement contains all points and lines of the intersections between the support planes and the bounding box. For each support plane, all intersections with the bounding box and other support planes are given by lines, edges and vertices of the arrangement. The kinetic partition created in the second step is a subset of faces of the arrangement depending on the k
parameter.
The kinetic partition for a chosen k
is obtained by propagating each polygon within its support plane. As intersections with other polygons can only occur at the known edges in the plane arrangement, the 3D collision problem can be solved as separate 2D polygon edge collisions.
The algorithm has five parameters:
k
: unsigned intk
, the maximum number of intersections that can occur for a polygon before its expansion stops. The initial intersections of the original input polygons are not considered. Thus increasing the k
leads to a higher complexity of the partitioning, i.e., a higher number of facets and a higher number of volumes. For a certain k
the partition can be considered to be complete and an increase in k
will not further increase the complexity of the partition. A typical choice of k
is in the range of 1 to 3.reorient_bbox
: booleanreorient_bbox
to true aligns the x-axis of the bounding box with the direction of the largest variation in horizontal direction of the input data while maintaining the z-axis.bbox_dilation_ratio
: FTmax_octree_node_size
: unsigned intmax_octree_depth
: unsigned intThe kinetic partition can be accessed as a LinearCellComplex
via CGAL::Kinetic_space_partition_3::get_linear_cell_complex()
.
The following example reads a set of polygons from a file and creates a kinetic partition. Increasing the k
parameter to 2 or 3 leads to a more detailed kinetic partition.
File Kinetic_space_partition/kinetic_partition.cpp
This package is an implementation of Bauchet et. al [1] with an octree replacing the grid subdivision. A proof of concept of the kinetic partition was developed by Simon Giraudot and Dmitry Anisimov.