如何访问 CGAL 3D 三角剖分中的构面？

Question

我正在使用CGAL 来计算一组点的 3D 三角剖分：

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Delaunay_triangulation_3<K>   CGALTriangulation;
typedef CGALTriangulation::Point            Point;

// Construction from a list of points
std::list<Point> points;
points.push_front(Point(0, 0, 0));
points.push_front(Point(2, 0, 0));
points.push_front(Point(0, 2, 0));
points.push_front(Point(2, 2, 0));
points.push_front(Point(1, 1, 1));

// Perform triangulation
CGALTriangulation T(points.begin(), points.end());

访问三角形（面）

我需要在 Unity 中创建一个网格，所以我使用 CGAL，因为它有很多算法可以处理这个复杂的问题。问题是在 API 中很难找到一种方法来访问构成三角剖分的不同三角形（以及它们的顶点）我还没有找到一种方法：（强>

注意请注意，仅访问顶点对我来说是不够的：

for (CGALTriangulation::Finite_vertices_iterator it = T.finite_vertices_begin(); 
    it != T.finite_vertices_end(); 
    it++) 
{
    CGALTriangulation::Triangulation_data_structure::Vertex v = *it;
    // Do something with the vertex
}

因为我没有得到任何关于每个顶点属于哪个面（三角形）的信息。而三角形正是我所需要的！

如何访问三角剖分的三角形（面）？如何从每个 facet 中取出顶点？

Answer 1

我不知道你能达到什么目标。 3D Delaunay 三角剖分是将点的凸包分解为四面体。无论如何，如果您想访问三角剖分的各个方面，请使用 Finite_facets_iterator。

类似：

for (CGALTriangulation::Finite_facets_iterator it = T.finite_facets_begin(); 
    it != T.finite_facets_end(); 
    it++) 
{
    std::pair<CGALTriangulation::Cell_handle, int> facet = *it;
    CGALTriangulation::Vertex_handle v1 = facet.first->vertex( (facet.second+1)%4 );
    CGALTriangulation::Vertex_handle v2 = facet.first->vertex( (facet.second+2)%4 );
    CGALTriangulation::Vertex_handle v3 = facet.first->vertex( (facet.second+3)%4 );
}

如果您对表面网格感兴趣，您可能希望查看重构算法，例如 Poisson surface reconstruction 或 Advancing Front Reconstruction。