This is a clean re-implementation of work I did years ago at a visual effects company. This library implements the three core CSG operations (intersection, union, and difference) on triangulated meshes. Unlike typical CSG implementations, the meshes to be operated on do not need to be solid (i.e. watertight) or even manifold.
All major pieces of functionality are working in limited tests. Only 'difference' operations are currently supported, with or without watertight "capping." Intersection and union operations are variations of the existing code, but I want to do more debugging and optimizations before implementation. No support for preserving normals or UVs. Poor handling of edge and corner cases. Output is a triangulated mesh, with the "triangle sliver" problem that comes along with that.
Requires Eigen.
In my previous implementation, I used a propietary retriangulation algorithm to break up the cut triangles. However, for this public implementation I am using Jonathan Shewchuk's excellent Triangle library. This code is free for personal and research use but not for commercial use, and does not fall under the license that the rest of this code uses.
If you wish to create the Python bindings or the Blender integration, SWIG also must be installed.
In the top-level csg directory:
mkdir build
cd build
cmake ..
make
After making the project, as described above, run make install to install the DSOs into your local Python 3.6 site packages. Then, in Blender, go to 'User Preferences' and select 'Install Addon From File.' Install the csg.py file, found in the blender directory and enable it (by clicking the empty square next to the name). To use the CSG tool, select two triangulated meshes: the clay first, then the knife. Execue the CSG operation (space bar, then type 'CSG') and two triangulated meshes will be created: one for the portion of the clay mesh above the knife, one for the portion of the clay mesh below the knife.
Essentially, the CSG algorithm finds all intersections of the triangles in mesh "A" and mesh "B", retriangulates the intersected faces to preserve these new edges, and then seperates the cut mesh (what we called the "clay") into two new meshes.
In more detail:
- Create AABB trees of the two meshes (the "clay" and the "knife").
- Using the AABB trees, do a broadphase intersection step to find pairs of bounding boxes that may contain intersecting triangles
- Determine exact triangle-triangle intersections, yeilding potentially new points that describe new edges.
- Retriangulate the cut triangles with any added points and edges from step 3.
- Categorize the cut triangles into two new surfaces, based on whether they are above or below the face that cut them. We then flood-fill the membership among the uncut triangles.
- Generate the output meshes by combining the four mesh fragment results (clay above the knife, clay below the knife, knife above the clay, knife below the clay) in various ways, depending on the operation we desire.