使用 CGAL 将多面体投影到 xy 平面

Polyhedron projection to xy plane with CGAL

我正在努力计算多面体在 xy 平面上的投影。目前我使用以下代码:

vector<Polygon_2> ii;
vector<Polygon_with_holes_2> oi;
for (Facet_iterator s = polyhedron.facets_begin();
                            s != polyhedron.facets_end(); ++s) {
    Halfedge_facet_circulator h = s->facet_begin(), he(h);

    Polygon_2 polygon;
    do {
        Point p1 = h->vertex()->point();
        polygon.insert(polygon.vertices_end(), Point_2(p1.x(), p1.y()));
    } while (++h != he);
    if (polygon.orientation() == CGAL::NEGATIVE)
        polygon.reverse_orientation();
    ii.push_back(polygon);
}
CGAL::join(ii.begin(), ii.end(), std::back_inserter(oi));

这将遍历多面体表面,并为每个表面手动执行 2d 投影。然后将生成的多边形连接在一起。

但是,从手册中我觉得这不是使用库来完成此类任务的预期方式。 class Project_traits_xy_3 提示,有一种常规方法可以使用 CGAL 实现投影。但是,我找不到合适的文档或示例。

谁能给我指出正确的方向?我觉得,这应该是一个优雅的标准任务。

您正在寻找的是轮廓的计算,而我们在 CGAL 中没有将其作为高级函数。 您的代码完成了这项工作。

如果你的物体是防水的,你可以只取法线向上的面。您还可以计算法线向上的连接面组件,然后为这些组件的边界创建多边形。

投影特征class用于将3D点解释为2D点,例如构造2D Constrained Delaunay三角剖分,但我看不出它能在这里派上用场。