CGAL中是否存在可移除Alpha Shape孔洞的函数？

要实现和MATLAB boundary.m 一致的单环无孔洞外边界，不需要调整Alpha值或者额外做孔洞填充，直接利用CGAL自带的几何元素分类属性过滤孔洞边即可，稳定性和效率都远高于自己实现点在多边形内判断的方案。

问题背景 #

当前基于CGAL官方2D Alpha Shapes示例编写的代码，默认输出会保留点云内部凹腔形成的孔洞，和目标效果的差异如下：

• 常规Alpha Shape输出：带内部孔洞的多环结构（效果对比图左侧）
• 目标输出：仅保留最外层闭合边界、无内部孔洞的单环结构（效果对比图右侧）

对照参考素材：CGAL官方2D Alpha Shapes示意图、Alpha Shape与目标边界的效果对比图

实现逻辑 #

CGAL生成2D Alpha Shape时，会自动给所有三角面、边标记分类属性：

1. 邻接无限外部区域的面为外层面
2. 被外层面完全包围、属于内部孔洞区域的面为孔洞面
3. 完全落在实体区域内部的面为内部面

我们只需要保留一侧邻接外层面、另一侧邻接实体面的边，所有和孔洞面相邻的边直接丢弃，最后把保留的边按端点邻接关系拼成闭合环，就是需要的无孔洞边界。

可直接复用的核心代码 #

基于官方示例代码修改，不需要引入额外依赖：

// 初始化Alpha Shape的逻辑和官方示例完全一致，不需要改动
typedef CGAL::Alpha_shape_2<DelaunayTriangulation> Alpha_shape_2;
Alpha_shape_2 as(points.begin(), points.end(), FT(0), Alpha_shape_2::GENERAL);

// 这里替换成你实际使用的alpha值，比如找最优alpha的逻辑可直接沿用官方示例写法
auto alpha_opt = *as.find_optimal_alpha(1);
as.set_alpha(alpha_opt);

std::vector<Segment_2> outer_boundary_edges;
for (auto edge_it = as.edges_begin(); edge_it != as.edges_end(); edge_it++) {
    // 先过滤掉非正则的退化边
    if (as.classify(edge_it) != Alpha_shape_2::REGULAR) continue;
    
    auto face = edge_it->first;
    int edge_idx = edge_it->second;
    auto neighbor_face = face->neighbor(edge_idx);
    
    // 边的两侧一个是外部区域、一个是Alpha Shape实体区域，判定为外边界边
    bool is_outer_edge = 
        (as.classify(face) == Alpha_shape_2::EXTERIOR && as.classify(neighbor_face) == Alpha_shape_2::INTERIOR) ||
        (as.classify(neighbor_face) == Alpha_shape_2::EXTERIOR && as.classify(face) == Alpha_shape_2::INTERIOR);
    
    if (is_outer_edge) {
        outer_boundary_edges.push_back(as.segment(edge_it));
    }
}

// 后续处理：把outer_boundary_edges里的线段按公共端点排序，拼接成闭合多边形环即可

避坑提示 #

• 不要靠增大Alpha值消除孔洞：Alpha值过大会导致外边界过度平滑收缩，和点云实际轮廓偏差过大，无法匹配MATLAB boundary 函数用收缩系数控制边界松紧的逻辑。
• 不需要自己实现孔洞填充算法：CGAL自带的分类标签已经把内外属性标记完成，自己写点在多边形内的判断不仅运行效率低，还容易在边重合、点共线的场景出现判断错误。

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