boost::geometry地理坐标使用问题：断言错误、策略疑惑及sym_difference实现咨询

一、地理坐标系多边形的sym_difference实现 #

你定义的polygon_t基于地理坐标系（经纬度），而Boost.Geometry的拓扑运算（如sym_difference、intersection等）默认仅支持平面坐标系，直接调用会触发断言错误，推荐以下解决方案：

方案1：投影到平面坐标系计算（优先推荐） #

将地理坐标转换为平面投影坐标系（如UTM、摩尔魏德投影），在平面完成拓扑运算后再转回地理坐标，示例代码如下：

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <boost/geometry/strategies/transform/mollweide.hpp>
#include <boost/geometry/algorithms/transform.hpp>

namespace bg = boost::geometry;

// 地理坐标类型
using geo_point_t = bg::model::point<double, 2, bg::cs::geographic<bg::degree>>;
using geo_polygon_t = bg::model::polygon<geo_point_t>;

// 平面投影坐标类型
using plane_point_t = bg::model::point<double, 2, bg::cs::cartesian>;
using plane_polygon_t = bg::model::polygon<plane_point_t>;

int main() {
    // 初始化两个地理多边形
    geo_polygon_t poly1, poly2;
    // ... 填充poly1和poly2的顶点数据 ...

    // 定义摩尔魏德投影策略（可替换为其他平面投影）
    auto proj_strategy = bg::strategy::transform::mollweide<geo_point_t, plane_point_t>();

    // 地理坐标转平面坐标
    plane_polygon_t plane_poly1, plane_poly2;
    bg::transform(poly1, plane_poly1, proj_strategy);
    bg::transform(poly2, plane_poly2, proj_strategy);

    // 执行sym_difference运算
    std::vector<plane_polygon_t> plane_results;
    bg::sym_difference(plane_poly1, plane_poly2, std::back_inserter(plane_results));

    // 平面结果转回地理坐标
    std::vector<geo_polygon_t> geo_results;
    for (auto const& plane_poly : plane_results) {
        geo_polygon_t geo_poly;
        bg::transform(plane_poly, geo_poly, proj_strategy.inverse());
        geo_results.push_back(geo_poly);
    }

    return 0;
}

方案2：使用地理拓扑运算扩展（版本受限） #

部分新版本Boost.Geometry针对地理坐标系提供了有限拓扑运算支持，但需显式指定地理策略，不过该功能仍不完善，建议优先用方案1。

二、策略命名差异解惑 #

• boost::geometry::strategy::area::geographic<...>：旧版本Boost（1.70之前）的命名方式，属于单策略命名空间，仅为向后兼容保留。
• boost::geometry::strategies::area::geographic<...>：新版本Boost（1.70及之后）的官方推荐接口，采用复数形式strategies命名空间，功能与旧版本一致，都是计算地理多边形面积的椭球面/球面策略。

三、优质学习资源推荐 #

1. Boost.Geometry官方文档：重点阅读「Geographic Coordinate Systems」章节，覆盖地理坐标处理逻辑、策略使用、投影转换等核心内容。
2. 官方示例代码：Boost源码包中example/geographic目录下有大量可运行示例，涵盖面积计算、距离计算、投影转换等常见场景。
3. Boost.Geometry官方FAQ：针对地理坐标的常见误区有专门解答，能快速定位问题。

内容的提问来源于stack exchange，提问作者geo-amateur