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

Polyhedron projection to xy plane with CGAL

本文关键字:xy 平面 投影 多面体 CGAL 使用      更新时间:2023-10-16

我正在努力计算多面体在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 投影。然后将生成的多边形连接在一起。

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

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

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

如果你的对象是水密的,你可以只取法线朝上的面。 您还可以计算面的连接分量,其中 normls 指向上方,然后为此类分量的边界创建多边形。

投影特征类用于将 3D 点解释为 2D 点,例如构造 2D 约束 Delaunay 三角测量,但我看不出它在这里可以派上用场。