如何在CGAL中更新三角剖分修改面的Face Info数据

CGAL的Constrained_Delaunay_triangulation_2没有直接提供修改面列表的接口，但可以通过以下两种可靠方式追踪面的变化，避免全量遍历：

1. 基于Triangulation_observer_2的事件回调追踪 #

这是最精准的方案，利用CGAL内置的观察者机制捕获三角剖分的修改事件：

• 继承CGAL::Triangulation_observer_2<你的CDT类型>，重载对应事件的回调函数，记录被修改的面：
  • before_face_removal(Face_handle f)：在面被删除前标记它
  • after_face_insertion(Face_handle f)：在新面插入后标记它
  • after_flip(Edge e)：边翻转后，标记涉及的两个相邻面
• 维护一个容器（比如std::unordered_set）存储修改的面，操作前后重置容器
• 操作完成后，仅遍历这个容器更新面的info数据

示例代码：

// 假设CDT是你的约束Delaunay三角剖分类型
struct FaceChangeObserver : public CGAL::Triangulation_observer_2<CDT> {
    std::unordered_set<CDT::Face_handle> modified_faces;

    void before_face_removal(CDT::Face_handle f) override {
        modified_faces.insert(f);
    }

    void after_face_insertion(CDT::Face_handle f) override {
        modified_faces.insert(f);
    }

    void after_flip(CDT::Edge e) override {
        modified_faces.insert(e.first);
        modified_faces.insert(e.first->neighbor(e.second));
    }

    void clear() { modified_faces.clear(); }
};

// 使用流程
CDT cdt;
FaceChangeObserver observer(cdt);

// 插入点前重置观测容器
observer.clear();
cdt.insert(new_point);

// 仅更新被修改的面
for (auto face : observer.modified_faces) {
    // 重新计算并设置face->info()，比如内外标识或自定义高成本数据
    face->info() = compute_face_data(face);
}

2. 利用顶点操作的返回值缩小范围 #

如果不想实现Observer，也可以通过增删点的返回值锁定影响范围：

• 插入点：insert()返回新顶点的Vertex_handle，从该顶点出发，遍历其相邻的所有面（通过vertex->face()和neighbor()方法），这些面是插入操作的直接影响面。注意这种方式可能包含少量未修改的面，但实现简单。
• 删除点：删除前先记录待删顶点的所有相邻面，删除后再遍历这些面的邻面，作为需要更新的候选。但删除操作的影响范围通常更大，Observer机制仍是更优选择。

关键注意事项 #

• 插入约束边时同样会触发面的变化，Observer机制也能捕获这些事件
• 用std::unordered_set存储修改面可以自动去重（比如边翻转会重复加入两个面）
• 如果你的面数据依赖邻面状态，可能需要对修改面的邻面做二次检查，但尽量通过事件回调减少不必要的计算

内容的提问来源于stack exchange，提问作者ImperatorS79