add pattern from triangulation

Author: tomvanmele

CHANGELOG.md

@@ -9,6 +9,9 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ### Added
 
+* Added "Pattern from Triangulation".
+* Added `compas_triangle` to requirements.
+
 ### Changed
 
 ### Removed

commands/RV_pattern.py

@@ -1,12 +1,17 @@
 #! python3
 # venv: brg-csd
-# r: compas_rv>=0.5.0
+# r: compas_rv>=0.6.0
 
 import rhinoscriptsyntax as rs  # type: ignore
+from compas_triangle.delaunay import conforming_delaunay_triangulation
+from compas_triangle.rhino import discretise_boundary
+from compas_triangle.rhino import discretise_constraints
 
 import compas_rhino
 import compas_rhino.conversions
 import compas_rhino.objects
+from compas.geometry import NurbsCurve
+from compas.geometry import Point
 from compas_rv.datastructures import Pattern
 from compas_rv.session import RVSession
 
@@ -102,6 +107,53 @@ def RunCommand():
         pattern = Pattern.from_meshgrid(dx=DX, nx=NX, dy=DY, ny=NY)
 
     elif option == "Triangulation":
+        boundary_guids = compas_rhino.objects.select_curves("Select outer boundary.")
+        if not boundary_guids:
+            return
+
+        rs.UnselectAllObjects()
+        hole_guids = compas_rhino.objects.select_curves("Select inner boundaries.")
+
+        rs.UnselectAllObjects()
+        segments_guids = compas_rhino.objects.select_curves("Select constraint curves.")
+
+        rs.UnselectAllObjects()
+
+        target_length = rs.GetReal("Specifiy target edge length.", 1.0)
+        if not target_length:
+            return
+
+        boundary = discretise_boundary(boundary_guids, target_length)
+        holes = None
+        segments = None
+        curves = None
+
+        if hole_guids:
+            holes = discretise_constraints(hole_guids, target_length)
+
+        if segments_guids:
+            segments = discretise_constraints(segments_guids, target_length)
+            curves: list[NurbsCurve] = [NurbsCurve.from_interpolation(segment) for segment in segments]
+
+        points, triangles = conforming_delaunay_triangulation(
+            boundary,
+            polylines=segments,
+            polygons=holes,
+            area=target_length**2 / 2,
+        )
+        pattern = Pattern.from_vertices_and_faces(points, triangles)
+
+        fixed = [vertex for boundary in pattern.vertices_on_boundaries() for vertex in boundary]
+        if curves:
+            for index, point in enumerate(points):
+                for curve in curves:
+                    closest: Point = curve.closest_point(point)
+                    if closest.distance_to_point(point) < 0.1 * target_length:
+                        fixed.append(index)
+
+        pattern.smooth_area(fixed=fixed)
+
+    elif option == "Skeleton":
         raise NotImplementedError
 
     else:

compas-RV.rhproj

@@ -8,11 +8,12 @@
   "id": "a6dc4669-0e8e-40ea-8d71-b9b0f4764ec1",
   "identity": {
     "name": "COMPAS-RhinoVAULT",
-    "version": "0.5.34-beta",
+    "version": "0.6.4",
     "publisher": {
-      "email": "[email protected]",
-      "name": "Juney Lee",
-      "company": "Carnegie Mellon University"
+      "email": "[email protected]",
+      "name": "Tom Van Mele",
+      "company": "COMPAS Association",
+      "url": "https://compas.dev/"
     },
     "description": "Rhino plugin for form finding of compression-only structures using Thrust Network Analysis.",
     "copyright": "ETH Zurich - Block Research Group",
@@ -610,7 +611,6 @@
       }
     }
   ],
-  "libraries": [],
   "resources": [
     {
       "id": "0475d743-2539-4543-b579-a5d20301277c",

requirements.txt

@@ -1,5 +1,7 @@
-compas >=2.5
+compas >=2.9
+compas_dem
 compas_rui >=0.4.2
 compas_session >=0.4.5
 compas_skeleton
 compas_tna >=0.5
+compas_triangle>=1.2.1