test_curvature_filter_openmesh.cpp

Go to the documentation of this file.

// Copyright (c) 2012-2019 University of Lyon and CNRS (France).
// All rights reserved.
//
// This file is part of MEPP2; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 3 of
// the License, or (at your option) any later version.
//
// This file is provided AS IS with NO WARRANTY OF ANY KIND, INCLUDING THE
// WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
#include "FEVV/Wrappings/Geometry_traits_openmesh.h" // Always in first position for MVS 4996 warnings
 
#include <fstream>
#include <string>
 
 
#define CGAL_USE_OM_POINTS
#include <CGAL/boost/graph/graph_traits_PolyMesh_ArrayKernelT.h>
#include <CGAL/boost/graph/properties_PolyMesh_ArrayKernelT.h>
#include "FEVV/Wrappings/Geometry_traits_openmesh.h"
 
 
#include "FEVV/Filters/Generic/Manifold/Curvature/curvature.hpp" // for calculate_curvature()
#include "FEVV/Filters/Generic/calculate_face_normals.hpp"
 
#include <OpenMesh/Core/IO/MeshIO.hh> // for OpenMesh::IO::read_mesh())
 
#include "FEVV/Wrappings/properties_openmesh.h" // for datastructure default property map
 
struct MyTraitsDouble : public OpenMesh::DefaultTraits
{
  typedef OpenMesh::Vec3d Point; // use double-values points
  typedef OpenMesh::Vec3d Normal;
};
 
typedef OpenMesh::PolyMesh_ArrayKernelT< MyTraitsDouble >
    OMMesh; // use MyTraitsDouble to remove a warning
 
//------------------------------------------------------------------------------
 
// code from Visualization/PluginFilters/curvature/CurvaturePlugin.h
void
curvature(OMMesh *mesh)
{
  typedef OMMesh HalfedgeGraph;
 
  std::cout << "Asking to Curvature mesh ! " << std::endl;
 
  auto pm = get(boost::vertex_point, *mesh);
 
  // ---
 
  //using Vector = typename FEVV::Geometry_traits< HalfedgeGraph >::Vector;
 
  // Face normal map
  // this one is a standard property map
  auto f_nm = make_property_map(FEVV::face_normal, *mesh);
 
  FEVV::Filters::calculate_face_normals(*mesh, pm, f_nm);
 
  // ---
 
  // Vertex curvature map
  // this one is a NON-standard property map!
  auto v_cm =
      FEVV::make_vertex_property_map< HalfedgeGraph,
                                      FEVV::Filters::v_Curv< HalfedgeGraph > >(
          *mesh);
 
  double min_nrm_min_curvature, max_nrm_min_curvature, min_nrm_max_curvature,
      max_nrm_max_curvature;
 
  bool value_is_geod = false;
  double value_radius = 0.001;
  FEVV::Filters::calculate_curvature( // B) call the filter corresponding to
                                      // your operation
      *mesh,
      v_cm,
      pm,
      f_nm,
      value_is_geod,
      value_radius,
      min_nrm_min_curvature,
      max_nrm_min_curvature,
      min_nrm_max_curvature,
      max_nrm_max_curvature);
 
  // ---
 
  std::cout << "Curvature mesh, isGeod: " << value_is_geod
            << " - radius: " << value_radius << "." << std::endl;
}
 
//------------------------------------------------------------------------------
 
int
main(int narg, char **argv)
{
  if(narg < 2)
  {
    std::cout << "Usage: " << argv[0] << " mesh_filename" << std::endl;
    exit(EXIT_FAILURE);
  }
  std::string filename(argv[1]);
 
 
  OMMesh m;
  if(!OpenMesh::IO::read_mesh(m, filename))
  {
    std::cout << "Failed to read file." << std::endl;
    return -1;
  }
 
 
  curvature(&m);
 
 
  std::cout << "Done." << std::endl;
 
  return 0;
}