如何从CGAL的General_polygon_set_2中解析出单个曲线

嘿，我完全理解你刚从MATLAB转C++接触CGAL的困惑——毕竟CGAL的模板和迭代器体系确实和MATLAB的数组思维差异挺大的。不过别担心，你已经走对路了，接下来咱们一步步把每个曲线拆出来：

核心思路回顾 #

你已经正确提取了Polygon_with_holes_2对象，接下来的关键是：outer_boundary()和每个孔洞返回的General_polygon_2本质是曲线的容器，每个元素是CGAL的Curve_2类型（在你的场景下就是线段或圆弧）。我们只需要遍历这些容器，再区分曲线类型即可。

第一步：定义必要的类型别名（避免冗长代码） #

首先确保你已经定义了匹配的Kernel和多边形类型（根据你用的CGAL布尔运算配置调整）：

// 根据你的实际Kernel调整，这里用带精确构造和平方根的版本
typedef CGAL::Exact_predicates_exact_constructions_kernel_with_sqrt Kernel;
typedef Kernel::Point_2 Point;
typedef Kernel::Segment_2 Segment;
typedef Kernel::Circle_2 Circle;
typedef CGAL::Circular_arc_2<Kernel> CircularArc;
typedef CGAL::General_polygon_2<Kernel> GeneralPolygon;
typedef CGAL::Polygon_with_holes_2<Kernel> PolygonWithHoles;
typedef CGAL::General_polygon_set_2<Kernel> GeneralPolygonSet;

第二步：遍历外边界的所有曲线 #

你已经拿到了minkOuter（GeneralPolygon类型），直接用范围for循环或者迭代器遍历每个曲线，然后判断类型并提取参数：

// 遍历外边界的每条曲线
std::cout << "=== 外边界曲线 ===" << std::endl;
for (const auto& curve : minkOuter) {
    // 判断是否为线段
    if (CGAL::is_segment(curve)) {
        Segment seg = CGAL::to_segment(curve);
        Point start = seg.source();
        Point end = seg.target();
        
        // 这里执行你的自定义操作，比如打印、存储坐标等
        std::cout << "线段: (" << start.x() << ", " << start.y() << ") --> (" 
                  << end.x() << ", " << end.y() << ")" << std::endl;
    }
    // 判断是否为圆弧
    else if (CGAL::is_circular_arc(curve)) {
        CircularArc arc = CGAL::to_circular_arc(curve);
        Circle supporting_circle = arc.supporting_circle();
        Point center = supporting_circle.center();
        // 注意：CGAL的Circle默认返回平方半径，如需实际半径用CGAL::sqrt()转换
        Kernel::FT radius_sq = supporting_circle.squared_radius();
        bool is_ccw = arc.is_counterclockwise_oriented();
        Point arc_start = arc.source();
        Point arc_end = arc.target();
        
        // 自定义操作
        std::cout << "圆弧: 圆心(" << center.x() << ", " << center.y() << "), 半径" 
                  << CGAL::sqrt(radius_sq) << ", 方向: " << (is_ccw ? "逆时针" : "顺时针") 
                  << ", 起止点: (" << arc_start.x() << ", " << arc_start.y() << ") --> (" 
                  << arc_end.x() << ", " << arc_end.y() << ")" << std::endl;
    }
}

第三步：遍历所有孔洞的曲线 #

用holes_begin()和holes_end()迭代器遍历每个孔洞，每个孔洞也是一个GeneralPolygon，遍历逻辑和外边界完全一致：

// 遍历所有孔洞
std::cout << "\n=== 孔洞曲线 ===" << std::endl;
int hole_index = 0;
for (auto hole_it = mink.holes_begin(); hole_it != mink.holes_end(); ++hole_it) {
    const GeneralPolygon& hole = *hole_it;
    std::cout << "\n孔洞 " << ++hole_index << "（包含 " << hole.size() << " 条曲线）:" << std::endl;
    
    for (const auto& curve : hole) {
        // 重复外边界的判断逻辑，这里省略重复代码，直接复用上面的判断块即可
        if (CGAL::is_segment(curve)) {
            Segment seg = CGAL::to_segment(curve);
            // 自定义操作...
        } else if (CGAL::is_circular_arc(curve)) {
            CircularArc arc = CGAL::to_circular_arc(curve);
            // 自定义操作...
        }
    }
}

关键注意事项 #

1. Kernel的匹配：确保你用的Kernel支持圆弧类型（比如带Circular_kernel的变种），否则无法正确识别圆弧。
2. 精确数值类型：CGAL的FT是精确数类型，如果你需要转换为double用于后续计算，可以用CGAL::to_double()函数，但注意可能丢失精度。
3. 曲线方向：孔洞的曲线默认是顺时针方向（和外边界的逆时针相反），这是CGAL的多边形约定，如果你需要统一方向可以调用reverse()方法反转。

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