aabb

Author: tomvanmele

docs/_images/example_surface_aabb.png

@@ -0,0 +1,43 @@
+from math import radians
+from compas.geometry import Translation, Rotation
+from compas.geometry import Point, Polyline
+from compas_occ.geometry import NurbsSurface
+
+from compas_view2.app import App
+
+points = [
+    [Point(0, 0, 0), Point(1, 0, 0), Point(2, 0, 0), Point(3, 0, 0), Point(4, 0, 0)],
+    [Point(0, 1, 0), Point(1, 1, 2), Point(2, 1, 2), Point(3, 1, 0), Point(4, 1, 0)],
+    [Point(0, 2, 0), Point(1, 2, 2), Point(2, 2, 2), Point(3, 2, 0), Point(4, 2, 0)],
+    [Point(0, 3, 0), Point(1, 3, 0), Point(2, 3, 0), Point(3, 3, 0), Point(4, 3, 0)],
+]
+
+surface = NurbsSurface.from_points(points=points)
+
+T = Translation.from_vector([0, -1.5, 0])
+R = Rotation.from_axis_and_angle([0, 0, 1], radians(45))
+
+surface.transform(R * T)
+
+# ==============================================================================
+# AABB
+# ==============================================================================
+
+box = surface.aabb()
+
+# ==============================================================================
+# Visualisation
+# ==============================================================================
+
+view = App()
+
+for row in surface.points:
+    view.add(Polyline(row), show_points=True, pointsize=20, pointcolor=(1, 0, 0), linewidth=2, linecolor=(1.0, 0, 0))
+
+for col in zip(* surface.points):
+    view.add(Polyline(col), linewidth=2, linecolor=(0, 1.0, 0))
+
+view.add(surface.to_mesh(u=50))
+view.add(box, show_faces=False)
+
+view.run()

docs/examples/surface_aabb.py

@@ -0,0 +1,10 @@
+********************************************************************************
+Surface Axis-Aligned Bounding Box
+********************************************************************************
+
+.. figure:: /_images/example_surface_aabb.png
+    :figclass: figure
+    :class: figure-img img-fluid
+
+.. literalinclude:: surface_aabb.py
+    :language: python

docs/examples/surface_aabb.rst

@@ -527,8 +527,8 @@ def copy(self) -> NurbsSurface:
     def transform(self, T: Transformation) -> None:
         """Transform this surface."""
         _T = gp_Trsf()
-        _T.SetValues(* T.list)
-        self.occ_surface.Transform(T)
+        _T.SetValues(* T.list[:12])
+        self.occ_surface.Transform(_T)
 
     def transformed(self, T: Transformation) -> NurbsSurface:
         """Transform an independent copy of this surface."""
@@ -609,7 +609,7 @@ def closest_point(self, point, distance=None) -> Point:
     def aabb(self, precision: float = 0.0) -> Box:
         """Compute the axis aligned bounding box of the surface."""
         box = Bnd_Box()
-        BndLib_AddSurface_Add(GeomAdaptor_Surface(self.occ_curve), precision, box)
+        BndLib_AddSurface_Add(GeomAdaptor_Surface(self.occ_surface), precision, box)
         return Box.from_diagonal((
             Point.from_occ(box.CornerMin()),
             Point.from_occ(box.CornerMax())