add volume measurement

Author: tomvanmele

CHANGELOG.md

@@ -0,0 +1,17 @@
+# Changelog
+
+All notable changes to this project will be documented in this file.
+
+The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
+and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
+
+
+## Unreleased
+
+### Added
+
+* Volume calculation of closed triangle mesh.
+
+### Changed
+
+### Removed

scripts/test_volume.py

@@ -0,0 +1,13 @@
+from compas.geometry import Box, Frame
+from compas.datastructures import Mesh
+
+from compas_cgal.measure import volume
+
+box = Box(Frame.worldXY(), 1, 1, 1)
+
+mesh = Mesh.from_shape(box)
+mesh.quads_to_triangles()
+
+V = volume(mesh.to_vertices_and_faces())
+
+print(V)

setup.py

@@ -59,6 +59,7 @@ def get_library_dirs():
             'src/booleans.cpp',
             'src/slicer.cpp',
             'src/intersections.cpp',
+            'src/measure.cpp',
         ]),
         include_dirs=[
             './include',
@@ -166,6 +167,7 @@ def build_extensions(self):
         "Programming Language :: Python :: 3",
         "Programming Language :: Python :: 3.6",
         "Programming Language :: Python :: 3.7",
+        "Programming Language :: Python :: 3.8",
         "Programming Language :: Python :: Implementation :: CPython",
     ],
     keywords=[],

src/compas_cgal.cpp

@@ -7,6 +7,7 @@ void init_meshing(py::module&);
 void init_booleans(py::module&);
 void init_slicer(py::module&);
 void init_intersections(py::module&);
+void init_measure(py::module&);
 
 PYBIND11_MODULE(_cgal, m) {
     m.doc() = "";
@@ -19,4 +20,5 @@ PYBIND11_MODULE(_cgal, m) {
     init_booleans(m);
     init_slicer(m);
     init_intersections(m);
+    init_measure(m);
 }

src/compas_cgal/__init__.py

@@ -68,7 +68,7 @@
     'compas_cgal.booleans',
     'compas_cgal.intersections',
     'compas_cgal.meshing',
-    'compas_cgal.slicer'
+    'compas_cgal.slicer',
 ]
 
 __all__ = ["HOME", "DATA", "DOCS", "TEMP"]

src/compas_cgal/measure.py

@@ -0,0 +1,29 @@
+import numpy as np
+from compas_cgal._cgal import measure
+from compas.plugins import plugin
+
+
+__all__ = [
+    'volume'
+]
+
+
+@plugin(category='trimesh', pluggable_name='trimesh_volume')
+def volume(mesh):
+    """Compute the volume of a closed triangle mesh.
+
+    Parameters
+    ----------
+    mesh : tuple of vertices and faces
+        The mesh.
+
+    Returns
+    -------
+    float
+        The volume of the mesh.
+
+    """
+    V, F = mesh
+    V = np.asarray(V, dtype=np.float64)
+    F = np.asarray(F, dtype=np.int32)
+    return measure.volume(V, F)

src/measure.cpp

@@ -0,0 +1,29 @@
+#include "measure.h"
+
+#include <CGAL/Polygon_mesh_processing/measure.h>
+
+namespace PMP = CGAL::Polygon_mesh_processing;
+namespace py = pybind11;
+
+double
+pmp_volume(
+    Eigen::Ref<const compas::RowMatrixXd> & V,
+    Eigen::Ref<const compas::RowMatrixXi> & F)
+{
+    Mesh mesh = compas::mesh_from_vertices_and_faces(V, F);
+
+    double volume = PMP::volume(mesh);
+
+    return volume;
+};
+
+void init_measure(py::module & m) {
+    py::module submodule = m.def_submodule("measure");
+
+    submodule.def(
+        "volume",
+        &pmp_volume,
+        py::arg("V").noconvert(),
+        py::arg("F").noconvert()
+    );
+};

src/measure.h

@@ -0,0 +1,11 @@
+#ifndef COMPAS_MEASURE_H
+#define COMPAS_MEASURE_H
+
+#include <compas.h>
+
+double
+pmp_volume(
+    Eigen::Ref<const compas::RowMatrixXd> & V,
+    Eigen::Ref<const compas::RowMatrixXi> & F);
+
+#endif /* COMPAS_MEASURE_H */