使用PROJ库实现CGCS2000坐标转WGS84坐标方法

PROJ库（Projection Library）是一个广泛使用的开源库，用于地理坐标系统的转换。要将CGCS2000（China Geocentric Coordinate System 2000）坐标系转换为WGS84（World Geodetic System 1984），你需要利用PROJ库中的特定投影参数和函数。以下是一般的步骤：

1. **安装和引入PROJ库**：首先确保已经安装了PROJ库。在C++项目中，你可以使用vcpkg或直接链接预编译的静态/动态库。

#include "proj.h" // 或者 #include <proj.h> 如果使用的是Linux环境

2. **创建一个`PJ_CONTEXT`结构体**：这用于管理PROJ操作。

PJ_CONTEXT* ctx = pj_create_context();
pj_init(ctx);

3. **选择投影**：设置从CGCS2000到WGS84的投影。CGCS2000通常用UTM或其他大地坐标系表示，而WGS84通常使用经纬度。你需要找到正确的UTM带号对应的参数，因为两者可能通过UTM进行间接转换。例如，如果CGCS2000在UTM65N，先转到WGS84的对应UTM带，再从UTM到纬度经度。

PJ *src_cgs2000, *dst_wgs84;
pj_str_t projargs;

// 设置源坐标系统（假设为UTM65N）
if (pj_set_errno(0)) {
    src_cgs2000 = pj_create("utm", &projargs); // 用实际的UTM带号替换65N
}

// 设置目标坐标系统
pj_strcpy(projargs, "+init=epsg:4326"); // WGS84的EPSG码
dst_wgs84 = pj_create("longlat", &projargs);

4. **进行坐标转换**：创建一个临时的投影，然后执行转换。

PJ *tmp_proj;
pj_transform tmp, ctx;

tmp.proj = dst_wgs84;
pj_transform_init_ctx(tmp, ctx);

// 假设src_point 是包含CGCS2000坐标的结构体
pj_transform(tmp, src_point.lon, src_point.lat, &dst_point.lng, &dst_point.lat, NULL);
pj_free_transform(&tmp);
pj_destroy(dst_wgs84);
pj_destroy(src_cgs2000);
pj_free_context(ctx);

5. **处理错误**：记得检查每个调用后的`pj_get_errno()`值，确保没有错误发生。