写一段不同坐标系经纬度转换的代码

在地理信息系统（GIS）中，将地球上的经纬度从一种坐标系转换到另一种常用坐标系是非常常见的需求。例如，从WGS84（世界大地坐标系统1984版本）转换到UTM（ Universal Transverse Mercator）或者墨卡托投影。这里提供一个简单的示例，使用Python中的`pyproj`库来进行这种转换，因为这是一个跨语言支持的库，适用于C++程序员。

首先，确保你已经安装了`pyproj`库。如果没有，可以通过pip安装：

pip install pyproj

然后你可以使用以下C++代码片段来进行转换，这里假设你已经有`proj4cpp`库，它是`pyproj`的C++接口：

#include "proj.h"
#include <string>

// 假设你有一个包含经度、纬度的结构体或类
struct GeoCoordinates {
    double latitude;
    double longitude;
};

GeoCoordinates wgs84ToUTM(double wgs84Lon, double wgs84Lat) {
    projPJ pj = proj_create Proj4Def("+init=EPSG:4326"); // WGS84 (lat/lon)
    projPJ utmPj = proj_create Proj4Def("+proj=utm +zone=00 +datum=WGS84 +units=m +no_defs"); // UTM (米)

    // 创建一个临时坐标点
    geo_point_t point;
    proj_to_meter(pj, &point, wgs84Lon, wgs84Lat);

    // 使用临时坐标点进行转换
    doubleutm utmCoord[2];
    proj_transform(utmPj, pj, utmCoord[0], utmCoord[1]);

    // 释放资源
    proj_destroy(pj);
    proj_destroy(utmPj);

    // 返回新的UTM坐标
    GeoCoordinates result;
    result.latitude = utmCoord[1]; // 注意UTM的Y轴是纬度向下
    result.longitude = utmCoord[0];
    return result;
}

// 如果你需要从UTM回转到WGS84，只需交换pj定义即可
GeoCoordinates utmToWGS84(double utmLon, double utmLat) {
    projPJ pj = proj_create Proj4Def("+proj=utm +zone=00 +datum=WGS84 +units=m +no_defs"); // UTM
    projPJ wgs84Pj = proj_create Proj4Def("+init=EPSG:4326"); // WGS84

    // ...其他步骤类似
}