rhino surface kinks, closest point on surface or boundaries and mapping/remapping objects on surface

Author: Robin Oval

src/compas_rhino/geometry/surface.py

@@ -6,8 +6,13 @@
 import compas_rhino
 
 from compas_rhino.geometry import RhinoGeometry
+from compas_rhino.geometry import RhinoCurve
+
+from compas_rhino.utilities import delete_objects
 
 from compas.geometry import subtract_vectors
+from compas.geometry import angle_vectors
+
 
 try:
     import rhinoscriptsyntax as rs
@@ -229,6 +234,33 @@ def borders(self, type=1):
         curves = rs.ExplodeCurves(border, delete_input=True)
         return curves
 
+    def kinks(self, threshold=1e-3):
+        """Return the XYZ coordinates of kinks, i.e. tangency discontinuities, along the surface's boundaries.
+
+        Returns
+        -------
+        list
+            The list of XYZ coordinates of surface boundary kinks.
+
+        """
+
+        kinks = []
+        borders = self.borders(type = 0)
+        
+        for border in borders:
+            border = RhinoCurve(border)
+            extremities = map(lambda x: rs.EvaluateCurve(border.guid, rs.CurveParameter(border.guid, x)), [0., 1.])
+        
+            if border.is_closed():
+                start_tgt, end_tgt = border.tangents(extremities)
+                if angle_vectors(start_tgt, end_tgt) > threshold:
+                    kinks += extremities 
+        
+            else:
+                kinks += extremities
+
+        return list(set(kinks))
+
     def project_point(self, point, direction=(0, 0, 1)):
         projections = rs.ProjectPointToSurface(point, self.guid, direction)
         if not projections:
@@ -242,13 +274,6 @@ def project_points(self, points, direction=(0, 0, 1), include_none=True):
         projections[:] = [self.closest_point(point) if not point else point for point in projections]
         return map(list, projections)
 
-    def closest_point(self, point, maxdist=None):
-        point = self.geometry.ClosestPoint(Point3d(*point))
-        return list(point)
-
-    def closest_points(self, points, maxdist=None):
-        return [self.closest_point(point) for point in points]
-
     def pull_point(self, point):
         pass
 
@@ -278,7 +303,143 @@ def pull_mesh(self, mesh, fixed=None, d=1.0):
     def pull_meshes(self, meshes):
         pass
 
+    def closest_point(self, xyz):
+        """Return the XYZ coordinates of the closest point on the surface from input XYZ-coordinates.
+    
+        Parameters
+        ----------
+        xyz : list
+            XYZ coordinates.
+
+        Returns
+        -------
+        list
+            The XYZ coordinates of the closest point on the surface.
+
+        """
+
+        return rs.EvaluateSurface(self.guid, *rs.SurfaceClosestPoint(self.guid, xyz))
+
+    def closest_points(self, points):
+        return [self.closest_point(point) for point in points]
+
+    def closest_point_on_boundaries(self, xyz):
+        """Return the XYZ coordinates of the closest point on the boundaries of the surface from input XYZ-coordinates.
+    
+        Parameters
+        ----------
+        xyz : list
+            XYZ coordinates.
+
+        Returns
+        -------
+        list
+            The XYZ coordinates of the closest point on the boundaries of the surface.
+
+        """
+
+        borders = self.borders(type = 0)
+        proj_dist = {tuple(proj_xyz): distance_point_point(xyz, proj_xyz) for proj_xyz in [RhinoCurve(border).closest_point(xyz) for border in borders]}
+        delete_objects(borders)
+        return min(proj_dist, key = proj_dist.get)
+
+    def closest_points_on_boundaries(self, points):
+        return [self.closest_point_on_boundaries(point) for point in points]
+
+    # --------------------------------------------------------------------------
+    # mapping
+    # --------------------------------------------------------------------------
+
+    def map_uv0(self, xyz):
+        """Return the UV(0) point from the mapping of a XYZ point based on the UV parameterisation of the surface.
+
+        Parameters
+        ----------
+        xyz : list
+            XYZ coordinates.
+
+        Returns
+        -------
+        list
+            The UV(0) coordinates of the mapped point.
+
+        """
+
+        return rs.SurfaceClosestPoint(self.guid, xyz) + (0.,)
+
+    def remap_xyz_point(self, uv):
+        """Return the XYZ point from the re-mapping of a UV point based on the UV parameterisation of the surface.
+
+        Parameters
+        ----------
+        uv : list
+            UV(0) coordinates.
+
+        Returns
+        -------
+        list
+            The XYZ coordinates of the re-mapped point.
+
+        """
+
+        return tuple(rs.EvaluateSurface(self.guid, *uv))
+
+    def remap_xyz_line(self, line):
+        """Return the XYZ points from the re-mapping of a UV line based on the UV parameterisation of the surface.
+
+        Parameters
+        ----------
+        uv : list
+            List of UV(0) coordinates.
+
+        Returns
+        -------
+        list
+            The list of XYZ coordinates of the re-mapped line.
+
+        """
+
+        return (self.remap_xyz_point(line[0][:2]), self.remap_xyz_point(line[1][:2]))
+
+    def remap_xyz_polyline(self, polyline):
+        """Return the XYZ points from the re-mapping of a UV polyline based on the UV parameterisation of the surface.
+
+        Parameters
+        ----------
+        uv : list
+            List of UV(0) coordinates.
+
+        Returns
+        -------
+        list
+            The list of XYZ coordinates of the re-mapped polyline.
+
+        """
+
+        return [self.remap_xyz_point(vertex[:2]) for vertex in polyline]
+
+    def remap_xyz_mesh(self, mesh, cls=None):
+        """Return the mesh from the re-mapping of a UV mesh based on the UV parameterisation of the surface.
+
+        Parameters
+        ----------
+        mesh : Mesh
+            A mesh.
+
+        Returns
+        -------
+        Mesh, cls
+            The re-mapped mesh.
+
+        """
+
+        if cls is None:
+            cls = type(mesh)
 
+        vertices, faces = mesh.to_vertices_and_faces()
+        vertices = [self.remap_xyz_point(uv0[:2]) for uv0 in vertices]
+        return cls.from_vertices_and_faces(vertices, faces)
+ 
 # ==============================================================================
 # Main
 # ==============================================================================