如何使用boost::geometry::line_interpolate在椭球上沿大圆非最短路径生成导航点

你当前用boost::geometry::line_interpolate结合椭球策略计算地理点间的最短路径插值点，但需要在特定场景下生成沿地球另一侧绕行的长路径（比如(20, 20)与(30, 20)之间，走跨越350度经度的路径而非10度的最短路径）。现有代码基于官方示例实现，默认生成最短路径的插值点。

要解决这个问题，关键在于理解Boost.Geometry的line_interpolate对于地理坐标系下的线段，默认会选择最短的大圆弧（弧长≤180度对应的路径）进行插值。要生成绕地球另一侧的长路径，我们需要构造一个能让Boost识别为“长弧”的线段，具体步骤如下：

修改思路 #

Boost.Geometry判断地理线段的路径时，会选择两点间经度差绝对值更小的那个弧。所以只要让两点的经度差绝对值大于180度，它就会自动选择长路径的大圆弧。我们可以通过调整其中一个端点的经度（使用等价的、超出常规范围的经度值）来实现这一点。

修改后的完整代码 #

#include <iostream>
#include <cmath>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>

// 辅助函数：将经度归一化到[-180, 180]的常规范围（可选）
double normalize_longitude(double lon) {
    lon = fmod(lon, 360.0);
    if (lon > 180.0) {
        lon -= 360.0;
    } else if (lon < -180.0) {
        lon += 360.0;
    }
    return lon;
}

int main() {
    typedef boost::geometry::model::d2::point_xy<double, boost::geometry::cs::geographic<boost::geometry::degree>> Point_Type;
    using Segment_Type = boost::geometry::model::segment<Point_Type>;
    using Multipoint_Type = boost::geometry::model::multi_point<Point_Type>;

    boost::geometry::srs::spheroid<double> spheroid(6378137.0, 6356752.3142451793);
    boost::geometry::strategy::line_interpolate::geographic<boost::geometry::strategy::vincenty> str(spheroid);

    // 调整终点经度为-330（等价于30度，但与起点20度的经度差为350度，触发长路径计算）
    Segment_Type const start_end_points { {20, 20}, {-330, 20} }; 
    double distance { 50000 }; // 每50km生成一个插值点

    Multipoint_Type mp;
    boost::geometry::line_interpolate(start_end_points, distance, mp, str);

    // 输出归一化后的插值点（如果不需要归一化，直接输出原始坐标即可）
    std::cout << "Long path interpolation points (normalized longitude):\n";
    for (auto const& p : mp) {
        double lon = normalize_longitude(p.x());
        double lat = p.y();
        std::cout << "POINT(" << lon << " " << lat << ")\n";
    }

    return 0;
}

代码说明 #

• 端点调整：把终点(30,20)的经度改为-330（30-360=-330），这样起点(20,20)和终点的经度差绝对值为350度，超过180度，Boost会自动选择长路径的大圆弧进行插值。
• 经度归一化：normalize_longitude函数将插值点的经度转换到[-180,180]的常规范围，方便后续使用。如果不需要这个处理，可以直接输出原始的插值点坐标。

这样修改后，line_interpolate就会生成你需要的沿地球另一侧绕行的长路径插值点了。

内容的提问来源于stack exchange，提问作者J'e