ADD subdivision and triangulation and types_std

Author: petrasvestartas

CMakeLists.txt

@@ -225,6 +225,7 @@ function(add_nanobind_module module_name source_file)
 endfunction()
 
 # Add new modules here
+add_nanobind_module(types_std src/types_std.cpp)
 add_nanobind_module(booleans_ext src/booleans.cpp)
 add_nanobind_module(meshing_ext src/meshing.cpp)
 add_nanobind_module(intersections_ext src/intersections.cpp)
@@ -233,4 +234,7 @@ add_nanobind_module(reconstruction_ext src/reconstruction.cpp)
 add_nanobind_module(skeletonization_ext src/skeletonization.cpp)
 add_nanobind_module(slicer_ext src/slicer.cpp)
 add_nanobind_module(straight_skeleton_2_ext src/straight_skeleton_2.cpp)
+add_nanobind_module(triangulation_ext src/triangulation.cpp)
+add_nanobind_module(subdivision_ext src/subdivision.cpp)
+
 

docs/examples/example_reconstruction_pointset_smoothing.py

@@ -17,7 +17,7 @@ def reconstruction_pointset_smoothing():
     ply_file_path = Path(__file__).parent.parent.parent / "data" / "box.ply"
     original_points = Pointcloud.from_ply(ply_file_path)
 
-    smoothed_points = pointset_smoothing(original_points, 1000, 3)
+    smoothed_points = pointset_smoothing(original_points, 100, 3)
 
     return Pointcloud(original_points), Pointcloud(smoothed_points)
 

docs/examples/example_straight_skeleton_2_interior_straight_skeleton_with_holes.py

@@ -1,4 +1,3 @@
-from compas.geometry import Point
 from compas.geometry import Polygon
 from compas_viewer import Viewer
 from line_profiler import profile

docs/examples/example_subdivision.py

@@ -0,0 +1,69 @@
+from compas.geometry import Box
+from compas.geometry import Polyhedron
+from compas.geometry import Translation
+from compas_viewer import Viewer
+from line_profiler import profile
+
+from compas_cgal.subdivision import mesh_subdivide_catmull_clark
+from compas_cgal.subdivision import mesh_subdivide_loop
+from compas_cgal.subdivision import mesh_subdivide_sqrt3
+
+
+@profile
+def main():
+    """Subdivide a mesh using different subdivision algorithms."""
+
+    # ==============================================================================
+    # Input geometry
+    # ==============================================================================
+
+    box = Box.from_diagonal(([0, 0, 0], [1, 1, 1]))
+
+    V0, F0 = box.to_vertices_and_faces(triangulated=False)
+    V1, F1 = box.to_vertices_and_faces(triangulated=False)
+    V2, F2 = box.to_vertices_and_faces(triangulated=True)
+    V3, F3 = box.to_vertices_and_faces(triangulated=True)
+
+    # ==============================================================================
+    # Compute
+    # ==============================================================================
+
+    V4, F4 = mesh_subdivide_catmull_clark((V1, F1), 3)
+    V5, F5 = mesh_subdivide_loop((V2, F2), 3)
+    V6, F6 = mesh_subdivide_sqrt3((V3, F3), 3)
+
+    S2 = Polyhedron(V0, F0)
+    S4 = Polyhedron(V4.tolist(), F4.tolist())
+    S6 = Polyhedron(V5.tolist(), F5.tolist())
+    S8 = Polyhedron(V6.tolist(), F6.tolist())
+
+    S2.transform(Translation.from_vector([1, 0, 0]))
+    S4.transform(Translation.from_vector([2, 0, 0]))
+    S6.transform(Translation.from_vector([3, 0, 0]))
+    S8.transform(Translation.from_vector([4, 0, 0]))
+
+    # ==============================================================================
+    # Visualize
+    # ==============================================================================
+
+    return S2, S4, S6, S8
+
+
+S2, S4, S6, S8 = main()
+
+
+# ==============================================================================
+# Visualize
+# ==============================================================================
+
+viewer = Viewer()
+
+viewer.renderer.camera.target = [3, 1, 0.5]
+viewer.renderer.camera.position = [5, -2, 0.5]
+
+viewer.scene.add(S2.to_mesh(), show_points=False)
+viewer.scene.add(S4.to_mesh(), show_points=False)
+viewer.scene.add(S6.to_mesh(), show_points=False)
+viewer.scene.add(S8.to_mesh(), show_points=False)
+
+viewer.show()

docs/examples/example_subdivision.rst

@@ -0,0 +1,18 @@
+Mesh Subdivision
+================
+
+This example demonstrates different mesh subdivision schemes available in COMPAS CGAL.
+
+Key Features:
+
+* Catmull-Clark subdivision for quad meshes
+* Loop subdivision for triangle meshes
+* √3 (sqrt3) subdivision for triangle meshes
+* Side-by-side comparison of different subdivision methods
+
+.. figure:: /_images/example_subdivision.png
+    :figclass: figure
+    :class: figure-img img-fluid
+
+.. literalinclude:: example_subdivision.py
+   :language: python

docs/examples/example_triangulation.py

@@ -0,0 +1,59 @@
+from compas.datastructures import Mesh
+from compas.geometry import Polygon
+from compas.geometry import Translation
+from compas_viewer import Viewer
+from line_profiler import profile
+
+from compas_cgal.triangulation import conforming_delaunay_triangulation
+from compas_cgal.triangulation import refined_delaunay_mesh
+
+
+@profile
+def main():
+    """Triangulate a mesh with holes."""
+
+    boundary = Polygon.from_sides_and_radius_xy(64, 4)
+
+    hole = Polygon.from_sides_and_radius_xy(128, 1)
+
+    hole1 = hole.transformed(Translation.from_vector([2, 0, 0]))
+    hole2 = hole.transformed(Translation.from_vector([-2, 0, 0]))
+    hole3 = hole.transformed(Translation.from_vector([0, 2, 0]))
+    hole4 = hole.transformed(Translation.from_vector([0, -2, 0]))
+
+    holes = [hole1, hole2, hole3, hole4]
+
+    V, F = conforming_delaunay_triangulation(
+        boundary,
+        holes=holes,
+    )
+
+    cdt = Mesh.from_vertices_and_faces(V, F)
+
+    V, F = refined_delaunay_mesh(
+        boundary,
+        holes=holes,
+        maxlength=0.5,
+        is_optimized=True,
+    )
+
+    rdm = Mesh.from_vertices_and_faces(V, F)
+
+    return cdt, rdm
+
+
+cdt, rdm = main()
+
+# ==============================================================================
+# Visualize
+# ==============================================================================
+
+viewer = Viewer()
+
+viewer.renderer.camera.target = [0, 0, 0]
+viewer.renderer.camera.position = [0, -1, 13]
+
+viewer.scene.add(cdt.transformed(Translation.from_vector([-5, 0, 0])))
+viewer.scene.add(rdm.transformed(Translation.from_vector([+5, 0, 0])))
+
+viewer.show()

docs/examples/example_triangulation.rst

@@ -0,0 +1,18 @@
+Mesh Triangulation
+==================
+
+This example demonstrates how to perform constrained and refined Delaunay triangulation using COMPAS CGAL.
+
+Key Features:
+
+* Creating polygonal boundary and holes
+* Conforming Delaunay triangulation
+* Refined Delaunay meshing with size constraints
+* Side-by-side visualization of different triangulation methods
+
+.. figure:: /_images/example_triangulation.png
+    :figclass: figure
+    :class: figure-img img-fluid
+
+.. literalinclude:: example_triangulation.py
+   :language: python

src/compas_cgal/intersections.py

@@ -2,6 +2,7 @@
 from compas.plugins import plugin
 
 from compas_cgal import intersections_ext
+from compas_cgal import types_std
 
 from .types import PolylinesNumpy
 from .types import VerticesFaces
@@ -48,6 +49,6 @@ def intersection_mesh_mesh(
     VB = np.asarray(VB, dtype=np.float64, order="C")
     FB = np.asarray(FB, dtype=np.int32, order="C")
 
-    pointsets: intersections_ext.VectorRowMatrixXd = intersections_ext.intersection_mesh_mesh(VA, FA, VB, FB)
+    pointsets: types_std.VectorRowMatrixXd = intersections_ext.intersection_mesh_mesh(VA, FA, VB, FB)
 
     return pointsets

src/compas_cgal/measure.py

@@ -3,6 +3,7 @@
 from compas.plugins import plugin
 
 from compas_cgal import measure_ext
+from compas_cgal import types_std
 
 from .types import VerticesFaces
 
@@ -76,6 +77,6 @@ def mesh_centroid(mesh: VerticesFaces) -> list[float]:
     V, F = mesh
     V = np.asarray(V, dtype=np.float64, order="C")
     F = np.asarray(F, dtype=np.int32, order="C")
-    vector_of_double: measure_ext.VectorDouble = measure_ext.centroid(V, F)
+    vector_of_double: types_std.VectorDouble = measure_ext.centroid(V, F)
     point = Point(*vector_of_double)
     return point

src/compas_cgal/skeletonization.py

@@ -2,6 +2,7 @@
 from compas.plugins import plugin
 
 from compas_cgal import skeletonization_ext
+from compas_cgal import types_std  # noqa: F401
 
 from .types import PolylinesNumpySkeleton
 from .types import VerticesFaces