ADD normals, boundary, and grouping.

Author: pv

docs/examples/example_mapping.py

@@ -1,20 +1,23 @@
 from pathlib import Path
 
 from compas.colors import Color
+from compas.geometry import Polyline
 from compas.datastructures import Mesh
 from compas_viewer import Viewer
 from compas_viewer.config import Config
 from tessagon.adaptors.list_adaptor import ListAdaptor
-from tessagon.types.hex_tessagon import HexTessagon
+from tessagon.types.zig_zag_tessagon import ZigZagTessagon
 
 from compas_libigl.mapping import map_mesh
 from compas_libigl.parametrisation import trimesh_lsc_mapping
 
+from compas import json_dump
+
 # ==============================================================================
 # Input geometry: 3D Mesh
 # ==============================================================================
 
-mesh = Mesh.from_obj(Path(__file__).parent.parent.parent / "data" / "minimal_surface.obj")
+mesh = Mesh.from_off(Path(__file__).parent.parent.parent / "data" / "beetle.off")
 
 for vertex in mesh.vertices():
     x, y, z = mesh.vertex_attributes(vertex, "xyz")  # type: ignore
@@ -32,24 +35,52 @@
     "function": lambda u, v: [u * 1, v * 1, 0],
     "u_range": [-0.255, 1.33],
     "v_range": [-0.34, 1.33],
-    "u_num": 16,
-    "v_num": 10,
+    "u_num": 20,
+    "v_num": 20,
     "u_cyclic": False,
     "v_cyclic": False,
     "adaptor_class": ListAdaptor,
 }
-tessagon = HexTessagon(**options)
+tessagon = ZigZagTessagon(**options)
 tessagon_mesh = tessagon.create_mesh()
 pv = tessagon_mesh["vert_list"]
 pf = tessagon_mesh["face_list"]
 
 # ==============================================================================
 # Mapping: 3D Mesh, 2D Pattern, UV
+# mv - mapped vertices
+# mf - mapped faces
+# mn - mapped normals
+# mg - mapped boundaries (True or False)
+# mb - mapped groups for polygons with holes (list of int)
 # ==============================================================================
-
-mv, mf = map_mesh((v, f), (pv, pf))
+mv, mf, mn, mb, mg = map_mesh((v, f), (pv, pf))
 mesh_mapped = Mesh.from_vertices_and_faces(mv, mf)
 
+
+# ==============================================================================
+# Offset mesh by normals, normals are interpolated from the original mesh.
+# ==============================================================================
+mesh_mapped_offset = mesh_mapped.copy()
+for i in range(mesh_mapped.number_of_vertices()):   
+    mesh_mapped_offset.vertex_attributes(i, "xyz", mesh_mapped.vertex_attributes(i, "xyz") - mn[i]*0.001)
+
+# ==============================================================================
+# Get Boundary Polylines
+# ==============================================================================
+boundaries = []
+for i in range(len(mb)):
+    if not mb[i]:
+        continue
+    points = []
+    for j in range(len(mf[i])):
+        id = mf[i][j]
+        points.append(mesh_mapped.vertex_attributes(id, "xyz") + mn[id]*0.002)
+    points.append(points[0])
+    polyline = Polyline(points)
+    boundaries.append(polyline)
+
+
 # ==============================================================================
 # Viewer
 # ==============================================================================
@@ -59,15 +90,20 @@
 config.camera.position = [5, 2, 1.5]
 
 viewer = Viewer(config=config)
-viewer.scene.add(mesh, name="mesh", show_faces=False, linecolor=Color.grey(), opacity=0.2)
+# viewer.scene.add(mesh, name="mesh", show_faces=False, linecolor=Color.grey(), opacity=0.2)
 viewer.scene.add(Mesh.from_vertices_and_faces(pv, pf), name="pattern2d")
-viewer.scene.add(mesh_mapped, name="mesh_mapped", facecolor=Color.red())
+viewer.scene.add(mesh_mapped, name="mesh_mapped", facecolor=Color.pink())
+viewer.scene.add(mesh_mapped_offset, name="mesh_mapped", facecolor=Color.blue())
+for boundary in boundaries:
+    viewer.scene.add(boundary, name="boundary", linecolor=Color.yellow(), linewidth=3)
 
 # To see where the pattern is mapped:
 uv = trimesh_lsc_mapping((v, f))
 mesh_flattened = mesh.copy()
 for i in range(mesh.number_of_vertices()):
     mesh_flattened.vertex_attributes(i, "xyz", [uv[i][0], uv[i][1], 0])
 
+json_dump([Mesh.from_vertices_and_faces(pv, pf), mesh_mapped, mesh_mapped_offset, boundaries, mesh_flattened], "mesh_flattened.json")
+
 viewer.scene.add(mesh_flattened, name="mesh_flattened")
 viewer.show()

docs/examples/example_mapping_all_patterns.py renamed to docs/examples/example_mapping_patterns.py

@@ -7,11 +7,13 @@
 
 from compas_libigl.mapping import map_pattern_to_mesh
 
+
+
 # ==============================================================================
 # Input geometry: 3D Mesh
 # ==============================================================================
 
-mesh = Mesh.from_off(Path(__file__).parent.parent.parent / "data" / "beetle.off")
+mesh = Mesh.from_obj(Path(__file__).parent.parent.parent / "data" / "minimal_surface.obj")
 
 for vertex in mesh.vertices():
     x, y, z = mesh.vertex_attributes(vertex, "xyz")  # type: ignore
@@ -20,7 +22,8 @@
 # ==============================================================================
 # Mapping: 3D Mesh, 2D Pattern, UV
 # ==============================================================================
-mesh_mapped = map_pattern_to_mesh("ZigZag", mesh, clip_boundaries=True, tolerance=1e-6, pattern_u=10, pattern_v=10)
+
+mesh_mapped = map_pattern_to_mesh("ZigZag", mesh, clip_boundaries=True, tolerance=1e-6, pattern_u=16, pattern_v=16)
 
 # ==============================================================================
 # Viewer
@@ -32,8 +35,7 @@
 
 viewer = Viewer(config=config)
 
-viewer.scene.add(mesh, name="mesh", show_faces=False, linecolor=Color.grey(), opacity=0.2)
+# viewer.scene.add(mesh, name="mesh", show_faces=False, linecolor=Color.grey(), opacity=0.2)
 viewer.scene.add(mesh_mapped, name="mesh_mapped", facecolor=Color.red())
 
-
 viewer.show()

mesh_flattened.json

@@ -21,6 +21,8 @@
 from tessagon.types.zig_zag_tessagon import ZigZagTessagon
 
 from compas_libigl import _mapping
+from compas_libigl._types_std import VectorBool  # noqa: F401
+from compas_libigl._types_std import VectorInt  # noqa: F401
 from compas_libigl._types_std import VectorVectorInt  # noqa: F401
 
 
@@ -50,7 +52,7 @@ def map_mesh(target_mesh, pattern_mesh, clip_boundaries=True, tolerance=1e-6):
     Returns
     -------
     tuple
-        A tuple containing (vertices, faces) of the mapped pattern mesh.
+        A tuple containing: mesh vertices, normals, grouping of polygons (to form holes), boundary flags.
     """
     # Unpack mesh tuples
     v, f = target_mesh
@@ -62,10 +64,12 @@ def map_mesh(target_mesh, pattern_mesh, clip_boundaries=True, tolerance=1e-6):
     pattern_v_numpy = np.array(pv, dtype=np.float64)
 
     # Perform the mapping
-    pattern_v_numpy_copy, pattern_f_numpy_cleaned = _mapping.map_mesh_with_automatic_parameterization(v_numpy, f_numpy, pattern_v_numpy, pf, clip_boundaries, tolerance)
+    pv_numpy_copy, pf_numpy_cleaned, p_normals, pattern_is_boundary, pattern_groups = _mapping.map_mesh_with_automatic_parameterization(
+        v_numpy, f_numpy, pattern_v_numpy, pf, clip_boundaries, tolerance
+    )
 
     # Return the result as a tuple
-    return pattern_v_numpy_copy, pattern_f_numpy_cleaned
+    return pv_numpy_copy, pf_numpy_cleaned, p_normals, pattern_is_boundary, pattern_groups
 
 
 def map_pattern_to_mesh(name, mesh, clip_boundaries=True, tolerance=1e-6, pattern_u=16, pattern_v=16):
@@ -151,6 +155,6 @@ def map_pattern_to_mesh(name, mesh, clip_boundaries=True, tolerance=1e-6, patter
     pf = tessagon_mesh["face_list"]
 
     v, f = mesh.to_vertices_and_faces()
-    mv, mf = map_mesh((v, f), (pv, pf), clip_boundaries=clip_boundaries, tolerance=tolerance)
+    mapped_vertices, mapped_faces, mapped_normals, mapped_is_boundary, mapped_groups = map_mesh((v, f), (pv, pf), clip_boundaries=clip_boundaries, tolerance=tolerance)
 
-    return Mesh.from_vertices_and_faces(mv, mf)
+    return Mesh.from_vertices_and_faces(mapped_vertices, mapped_faces)