feat(mesh-filters): start meshtoGeogramMesh

Author: thewtex

packages/mesh-filters/CMakeLists.txt

@@ -68,11 +68,12 @@ FetchContent_Declare(
 
 FetchContent_MakeAvailable(pmp geogram)
 include_directories(${pmp_SOURCE_DIR}/src)
+include_directories(${geogram_SOURCE_DIR}/src/lib)
 
 foreach(pipeline geogram-conversion)
 # foreach(pipeline mesh-filters mesh-filters-sigma gaussian-kernel-radius mesh-filters-bin-shrink mesh-filters-label-image)
   add_executable(${pipeline} ${pipeline}.cxx)
-  target_link_libraries(${pipeline} PUBLIC ${ITK_LIBRARIES} pmp)
+  target_link_libraries(${pipeline} PUBLIC ${ITK_LIBRARIES} geogram)
   target_include_directories(${pipeline} PUBLIC ${CMAKE_CURRENT_SOURCE_DIR})
 endforeach()
 

packages/mesh-filters/geogram-conversion.cxx

@@ -21,6 +21,8 @@
 #include "itkOutputMesh.h"
 #include "itkPipeline.h"
 
+#include "itkmeshToGeogramMesh.h"
+
 int main( int argc, char * argv[] )
 {
   itk::wasm::Pipeline pipeline("mesh-read-write-test", "A test for reading and writing meshes", argc, argv);

packages/mesh-filters/itkmeshToGeogramMesh.h

@@ -0,0 +1,68 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  Licensed under the Apache License, Version 2.0 (the "License");
+ *  you may not use this file except in compliance with the License.
+ *  You may obtain a copy of the License at
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  Unless required by applicable law or agreed to in writing, software
+ *  distributed under the License is distributed on an "AS IS" BASIS,
+ *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ *  See the License for the specific language governing permissions and
+ *  limitations under the License.
+ *
+ *=========================================================================*/
+#ifndef itkmeshToGeogramMesh_h
+#define itkmeshToGeogramMesh_h
+
+#include <memory>
+
+#include "geogram/mesh/mesh.h"
+
+namespace itk
+{
+
+template <typename TMesh>
+auto
+meshToGeogramMesh(const TMesh * itkMesh) -> std::unique_ptr<GEO::Mesh>
+{
+  using MeshType = TMesh;
+  static constexpr unsigned int Dimension = MeshType::PointDimension;
+  static constexpr bool SinglePrecision = std::is_same<typename MeshType::CoordRepType, float>::value;
+  std::unique_ptr<GEO::Mesh> geoMesh = std::make_unique<GEO::Mesh>(Dimension, SinglePrecision);
+
+  // Copy vertices
+  auto points = itkMesh->GetPoints();
+  geoMesh->vertices.create_vertices(points->Size());
+
+  for(auto it = points->Begin(); it != points->End(); ++it)
+  {
+    const auto& point = it.Value();
+    GEO::index_t v = it.Index();
+    for (unsigned int d = 0; d < Dimension; ++d)
+    {
+      geoMesh->vertices.point(v)[d] = point[d];
+    }
+  }
+
+  // Copy faces/cells
+  auto cells = itkMesh->GetCells();
+  for(auto it = cells->Begin(); it != cells->End(); ++it)
+  {
+    const auto& cell = it.Value();
+    GEO::vector<GEO::index_t> vertices;
+    for(auto* pit = cell->PointIdsBegin(); pit != cell->PointIdsEnd(); ++pit) {
+        vertices.push_back(*pit);
+    }
+    geoMesh->facets.create_polygon(vertices);
+  }
+
+  return geoMesh;
+}
+
+} // namespace itk
+
+#endif // itkmeshToGeogramMesh_h