CGAL :: Surface_mesh_parameterization: записать выключенные вершины в исходном порядке - PullRequest
Новая
       39

CGAL :: Surface_mesh_parameterization: записать выключенные вершины в исходном порядке

0 голосов
/ 06 февраля 2019

Я пытаюсь изменить CGAL-4.13/examples/Surface_mesh_parameterization/lscm.cpp так, чтобы порядок вершин в результирующем отключенном файле был таким же, как и во входном файле.

Обновление: пример

Взять файлinput.off со следующим простым содержанием: 

OFF
4 2 0
-0.9310345 0.4333333 0 
-1 0.4333333 0 
-0.9310345 0.5 0 
-1 0.5 0 
3 1 0 2
3 2 3 1

Когда я вызываю стандарт lscm из CGAL с 

/path/to/CGAL-4.13-build/examples/Surface_mesh_parameterization/lscm input.off

, я получаю coords.off, содержащий 

OFF
4 2 0
-1 0.5 0
-0.931034 0.5 0
-0.931034 0.433333 0
-1 0.433333 0
3 3 2 1
3 1 0 3

и uvmap.off с 

OFF
4 2 0
-0.0166567 0.982769 0
1 1 0
1.01666 0.0172311 0
0 0 0
3 3 2 1
3 1 0 3

Файлы coords.off и uvmap.off содержат вершины и их пары параметров в одном и том же порядке (который отличается от того, что вinput.off).Вместо этого я бы хотел, чтобы параметры в uvmap.off были в порядке, соответствующем input.off.В частности, я хочу, чтобы uvmap.off выглядел так: 

OFF
4 2 0
1.01666 0.0172311 0
0 0 0
1 1 0
-0.0166567 0.982769 0
3 1 0 2
3 2 3 1

По сути, это делает coords.off избыточным, поскольку я могу использовать input.off в его роли.

Усилие врешение

Насколько я понимаю, этого можно достичь, вызвав output_uvmap_to_off(...) с 6 параметрами вместо 4 (обе версии можно найти в CGAL-4.13/include/CGAL/Surface_mesh_parameterization/IO/File_off.h).Поскольку одним из этих параметров является VertexIndexMap, я, вероятно, должен также использовать 

CGAL::Surface_mesh_parameterization::LSCM_parameterizer_3< TriangleMesh_, BorderParameterizer_, SolverTraits_ >::parameterize(...)

вместо 

CGAL::Surface_mesh_parameterization::parameterize(...)

, использованного в примере.

Вотминимальный (не очень рабочий) пример.Он получен из lscm.cpp, но я выбросил много вещей, чтобы остаться кратким. 

#include <CGAL/Simple_cartesian.h>
#include <CGAL/Surface_mesh.h>
#include <CGAL/boost/graph/Seam_mesh.h>
#include <CGAL/Surface_mesh_parameterization/IO/File_off.h>
#include <CGAL/Surface_mesh_parameterization/parameterize.h>
#include <CGAL/Surface_mesh_parameterization/Two_vertices_parameterizer_3.h>
#include <CGAL/Surface_mesh_parameterization/LSCM_parameterizer_3.h>
#include <CGAL/Polygon_mesh_processing/measure.h>
#include <boost/foreach.hpp>
#include <iostream>
#include <fstream>
#include <sstream>
#include <string>

typedef CGAL::Simple_cartesian<double>      Kernel;
typedef Kernel::Point_2                     Point_2;
typedef Kernel::Point_3                     Point_3;
typedef CGAL::Surface_mesh<Kernel::Point_3> SurfaceMesh;
typedef boost::graph_traits<SurfaceMesh>::edge_descriptor SM_edge_descriptor;
typedef boost::graph_traits<SurfaceMesh>::halfedge_descriptor SM_halfedge_descriptor;
typedef boost::graph_traits<SurfaceMesh>::vertex_descriptor SM_vertex_descriptor;
typedef SurfaceMesh::Property_map<SM_halfedge_descriptor, Point_2> UV_pmap;
typedef SurfaceMesh::Property_map<SM_edge_descriptor, bool> Seam_edge_pmap;
typedef SurfaceMesh::Property_map<SM_vertex_descriptor, bool> Seam_vertex_pmap;
typedef CGAL::Seam_mesh<SurfaceMesh, Seam_edge_pmap, Seam_vertex_pmap> Mesh;
typedef boost::graph_traits<Mesh>::vertex_descriptor vertex_descriptor;
typedef boost::graph_traits<Mesh>::halfedge_descriptor halfedge_descriptor;
typedef boost::graph_traits<Mesh>::face_descriptor face_descriptor;

namespace SMP = CGAL::Surface_mesh_parameterization;

int main(int argc, char** argv)
{
  std::ifstream in_mesh((argc>1) ? argv[1] : "data/lion.off");
  if(!in_mesh){
    std::cerr << "Error: problem loading the input data" << std::endl;
    return EXIT_FAILURE;
  }

  SurfaceMesh sm;
  in_mesh >> sm;

  Seam_edge_pmap seam_edge_pm = sm.add_property_map<SM_edge_descriptor, bool>("e:on_seam", false).first;
  Seam_vertex_pmap seam_vertex_pm = sm.add_property_map<SM_vertex_descriptor, bool>("v:on_seam", false).first;

  Mesh mesh(sm, seam_edge_pm, seam_vertex_pm);
  UV_pmap uv_pm = sm.add_property_map<SM_halfedge_descriptor, Point_2>("h:uv").first;

  halfedge_descriptor bhd = CGAL::Polygon_mesh_processing::longest_border(mesh, CGAL::Polygon_mesh_processing::parameters::all_default()).first;

  typedef SMP::Two_vertices_parameterizer_3<Mesh>                Border_parameterizer;
  typedef SMP::LSCM_parameterizer_3<Mesh, Border_parameterizer>  Parameterizer;

  // Here's where the big changes start.
  SurfaceMesh::Property_map<SM_halfedge_descriptor, int> vimap = sm.add_property_map<SM_halfedge_descriptor, int>("h:vi").first;
  SurfaceMesh::Property_map<SM_halfedge_descriptor, bool> vpmap = sm.add_property_map<SM_halfedge_descriptor, bool>("h:vp").first;
  Parameterizer parameterizer;
  parameterizer.parameterize(mesh, bhd, uv_pm, vimap, vpmap);

  const char* uvmap_file = "uvmap.off";
  std::ofstream uvmap_out(uvmap_file);
  SMP::IO::output_uvmap_to_off(mesh,sm.vertices(),sm.faces(),uv_pm,vimap,uvmap_out);

  return EXIT_SUCCESS;
}

Это не компилируется, жалуясь на требуемое преобразование в строке 131 в File_off.h.

Актуальные вопросы

  • Правильно ли инициализирована vimap?
  • Является ли это разумным направлением к моей цели написания вершин в том же порядке?
  • Еслида, как мне передать правильные аргументы output_uvmap_to_off(...)?Например, он запрашивает VertexContainer, а я предоставляю Vertex_range (отсюда и ошибка компиляции, я полагаю).Должен ли я просто собрать вершины, как предложено здесь или здесь или есть более элегантный способ?
  • Если нет, каков правильный курс действий?

1 Ответ

0 голосов
/ 27 февраля 2019

Следующая проблема отвечает на этот вопрос.

...