如何让PROJ C++识别网格目录 实现EGM96与WGS84坐标转换

核心错误点如下：

• CRS维度不匹配：EPSG:4326是纯2D经纬度坐标参考系，不带高程维度。pyproj做了自动兼容，会将其隐式升级为带椭球高的3D坐标系EPSG:4979，但PROJ C API不会做自动兼容，导致2D转带高程的EPSG:5773时转换链创建失败。
• 未判空转换对象：你没有对proj_create_crs_to_crs的返回值做非空校验，当转换链创建失败时，后续proj_trans调用自然返回inf无效值。
• 搜索路径设置逻辑冗余：你不需要手动拼接默认数据路径，proj_context_set_search_paths会自动保留默认搜索路径，额外添加自定义路径即可。
• 代码不符合C++标准：C++标准要求main函数返回值为int，你写的void main属于非标准写法，部分编译器下会出现异常。
• 错误输出异常导致内容被吞：转换失败时proj_errno_string可能返回带控制字符的字符串，触发控制台输出截断，所以看不到有效日志。

修复后的可运行代码如下：

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

int main(){
    // 创建PROJ上下文
    PJ_CONTEXT* proj_context = proj_context_create();
    // 禁用网络下载，仅使用本地栅格文件
    proj_context_set_enable_network(proj_context, 0);
    // 替换为你存放us_nga_egm96_15.tif的目录绝对路径
    const char * grid_path =  "/absolute_directory_to/proj";
    // 仅添加自定义栅格搜索路径，默认数据路径会自动保留
    proj_context_set_search_paths(proj_context, 1, &grid_path);

    // 检查路径设置错误
    int err = proj_context_errno(proj_context);
    if (err != 0) {
        std::cout << "路径设置错误: " << proj_errno_string(err) << std::endl;
        proj_context_destroy(proj_context);
        return 1;
    }

    // 源CRS用带3D维度的EPSG:4979（WGS84 经纬度+椭球高），匹配高程转换需求
    PJ* transformer = proj_create_crs_to_crs (proj_context,
                                        "EPSG:4979",
                                        "EPSG:5773",
                                        NULL);
    // 校验转换链是否创建成功
    if (transformer == nullptr) {
        err = proj_context_errno(proj_context);
        std::cout << "转换链创建失败: " << proj_errno_string(err) << std::endl;
        proj_context_destroy(proj_context);
        return 1;
    }

    // 执行逆转换：从EPSG:5773（EGM96高程）转EPSG:4979（WGS84椭球高）
    // 坐标顺序遵循EPSG原生顺序：纬度、经度、高程
    PJ_COORD input = proj_coord(43.70012234, -79.41629234, 100, 0);
    PJ_COORD output = proj_trans(transformer, PJ_INV, input);

    err = proj_errno(transformer);
    if (err != 0) {
        std::cout << "转换失败: " << proj_errno_string(err) << std::endl;
        proj_destroy(transformer);
        proj_context_destroy(proj_context);
        return 1;
    }

    // 输出结果和Python代码运行结果完全一致
    std::cout << output.v[0] << "," << output.v[1] << "," << output.v[2] << std::endl;

    // 释放资源
    proj_destroy(transformer);
    proj_context_destroy(proj_context);
    return 0; 
}

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