如何使用PROJ C++ API实现CRS间坐标转换？

以下是和你提供的pyproj代码逻辑完全对应的C++实现，包含完整流程和关键注意事项：

核心流程对应关系 #

pyproj的Transformer.from_crs + transform流程，在C++ API中对应：

• 初始化PROJ上下文
• 通过EPSG代码创建源/目标CRS
• 创建CRS转换对象
• 执行坐标转换
• 清理资源

完整示例代码 #

#include <proj.h>
#include <iostream>

int main() {
    // 1. 创建PROJ上下文（pyproj自动处理，C++需显式初始化）
    proj_context* ctx = proj_context_create();
    if (!ctx) {
        std::cerr << "Failed to create PROJ context" << std::endl;
        return 1;
    }

    // 2. 从EPSG代码创建源CRS（对应pyproj的4326，WGS84经纬度）
    proj_obj* src_crs = proj_crs_from_epsg(ctx, 4326);
    // 目标CRS（对应pyproj的6677，日本平面直角坐标系IX系）
    proj_obj* dst_crs = proj_crs_from_epsg(ctx, 6677);
    if (!src_crs || !dst_crs) {
        std::cerr << "Failed to create CRS from EPSG code" << std::endl;
        proj_context_destroy(ctx);
        return 1;
    }

    // 3. 创建坐标转换对象（对应pyproj的Transformer.from_crs）
    PJ* transformer = proj_create_crs_to_crs(ctx, src_crs, dst_crs, nullptr);
    if (!transformer) {
        std::cerr << "Failed to create transformer" << std::endl;
        proj_destroy(src_crs);
        proj_destroy(dst_crs);
        proj_context_destroy(ctx);
        return 1;
    }

    // 4. 准备坐标：注意PROJ中EPSG4326默认顺序是 经度(x), 纬度(y)
    //    对应pyproj输入的(35.9032597, 139.9343755)是纬度在前，需要调换顺序
    double x = 139.9343755;
    double y = 35.9032597;

    // 5. 执行正向转换（对应pyproj的transform方法）
    int result = proj_trans(transformer, PJ_FWD, &x, &y, nullptr);
    if (result != 0) {
        std::cerr << "Transformation failed" << std::endl;
        proj_destroy(transformer);
        proj_destroy(src_crs);
        proj_destroy(dst_crs);
        proj_context_destroy(ctx);
        return 1;
    }

    // 输出转换结果，和pyproj的输出一致（x,y对应pyproj返回的第二、第一个值）
    std::cout << "转换结果：(" << x << ", " << y << ")" << std::endl;

    // 6. 清理资源（pyproj自动回收，C++需显式释放）
    proj_destroy(transformer);
    proj_destroy(src_crs);
    proj_destroy(dst_crs);
    proj_context_destroy(ctx);

    return 0;
}

关键注意事项 #

• 坐标顺序差异：pyproj默认接受(纬度, 经度)作为EPSG4326的输入，而PROJ C++ API中EPSG4326的标准顺序是(经度, 纬度)，转换时必须调换顺序才能得到正确结果
• 错误处理：示例中加入了基础错误检查，实际生产环境需根据PROJ的返回值做更细致的错误处理
• 编译链接：编译时需链接PROJ库，以Linux为例，编译命令类似：

  g++ -o proj_transform proj_transform.cpp -lproj -I/usr/include/proj -L/usr/lib/x86_64-linux-gnu/
  

  路径需根据你的PROJ安装位置调整

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