Question

我使用Polyhedron_3作为表面。我扭曲表面并确保质量我想翻转边缘以避免坏三角形。到目前为止，我的代码看起来像：

std::vector<std::pair<PlaneMeshAPI::Polyhedron::Halfedge_handle, double> > vEdgeToFlip;
for (PlaneMeshAPI::Polyhedron::Edge_iterator e = P.edges_begin(); e != P.edges_end(); ++e)
{
    // Edge_iterator so that we consider only one of the 2 possible halfedges
    bool bFlippable = true;
    if (e->is_border_edge()) bFlippable = false;
    if (bFlippable && e->facet()->marked() == -1) bFlippable = false;
    if (bFlippable && e->facet()->marked() != e->opposite()->facet()->marked()) bFlippable = false;
    // Marked() returns an int, I want to flip edges between two triangles of the same component

    if (bFlippable)
    {
        PlaneMeshAPI::Polyhedron::Facet_iterator f1, f2;
        PlaneMeshAPI::Polyhedron::Halfedge_handle heh = e;
        double lowestBef = lowestAngle(e->facet(), e->opposite()->facet()); // returns the lowest angle of the two adjacent triangles
        vEdgeToFlip.push_back(std::make_pair(e, lowestBef));
    }
}

for (int i = 0; i < vEdgeToFlip.size(); ++i)
{
    PlaneMeshAPI::Polyhedron::Halfedge_handle e = vEdgeToFlip[i].first;
    e = P.flip_edge(e);
    double lowestNow = lowestAngle(e->facet(), e->opposite()->facet());
    if (lowestNow < vEdgeToFlip[i].second)
        P.flip_edge(e);
}

代码运行正常但是当我运行P.is_valid(true)时出现此错误信息：

halfedge 7504 previous pointer integrity corrupted. summe border halfedges (2*nb) = 0 end of CGAL::HalfedgeDS_const_decorator<HDS>::is_valid(): structure is NOT VALID . counting halfedges failed. end of CGAL::Polyhedron_3<...>::is_valid(): structure is NOT VALID.

flip_edge上的文档很少见。我不知道是否需要翻转两个半边，如果它在迭代器中打破了某些东西（这样一旦我翻了一个，所有其他的都不能被翻转）。

Answer 1

我们终于找到了为什么边缘翻转会导致表面破裂。在翻转小平面e = P.flip_edge(e);之前，您必须确保它不会产生奇点：

   // Do not flip if this would create two triangle with the same vertices
    if (e->next()->opposite()->next()->opposite() == e->opposite()->next()->next()) continue;
    if (e->opposite()->next()->opposite()->next()->opposite() == e->next()->next()) continue;

    // Do not flip if it would create an edge linking a vertex with itself
    if (e->next()->vertex() == e->opposite()->next()->vertex()) continue;

CGAL Polyhedron_3 flip_edge函数打破了表面

1 个答案: