从网格中读取顶点和三角形

Question

我只是花了几个小时在Visual Studio c++中的CGAL试图了解网格是如何工作的。我想要的是访问顶点和三角形的列表(顶点是double[3]的形式，三角形是int[3]的形式)。以下是我正在编写的脚本:

http://doc.cgal.org/latest/Surface_mesher/Surface_mesher_2mesh_a_3d_gray_image_8cpp-example.html

关键是-函数CGAL::output_surface_facets_to_off (out, c2t3);以。off格式输出我一个很好的文件(由MeshLab访问)，但我不能通过操纵c2t3或tr变量来做任何类似的事情。我期待的是这样的内容:

c2t3。顶点(从0到N)和c2t3。三角形(从0到M)，值为整数的三倍。我得到的是顶点列表，切面列表，单元格列表，边列表。除了在非排序的顶点列表中查找每个顶点编号，没有其他方法可以从facet中获得顶点编号。

谁能解决我的问题，指出我做错了什么?此外，CGAL的API非常…生的。挖掘源代码也是非常硬核的——我找不到output_surface函数体。

Answer 1

好的，我在Complex_2_in_triangulation_3_file_writer.h文件中找到了答案。显然，CGAL库正在创建整个地图的顶点和寻找面顶点在这个地图。代码部分:

using CGAL::Surface_mesher::number_of_facets_on_surface;
typedef typename C2t3::Triangulation Tr;
typedef typename Tr::Finite_facets_iterator Finite_facets_iterator;
typedef typename Tr::Finite_vertices_iterator Finite_vertices_iterator;
typedef typename Tr::Facet Facet;
typedef typename Tr::Edge Edge;
typedef typename Tr::Vertex_handle Vertex_handle;
std::map<Vertex_handle, int> V;
int inum = 0;
for(Finite_vertices_iterator vit = tr.finite_vertices_begin();
vit != tr.finite_vertices_end();
++vit)
{
    V[vit] = inum++;
}
for( Finite_facets_iterator fit = tr.finite_facets_begin();
fit != tr.finite_facets_end(); ++fit)
{
const typename Tr::Cell_handle cell = fit->first;
const int& index = fit->second;
    if (cell->is_facet_on_surface(index)==true)
    {
    const int index1 = V[cell->vertex(tr.vertex_triple_index(index, 0))];
    const int index2 = V[cell->vertex(tr.vertex_triple_index(index, 1))];
    const int index3 = V[cell->vertex(tr.vertex_triple_index(index, 2))];
    }
}