如何使用C++编程语言与GDAL库进行空间参考坐标系统（SCS）转换，并提供完整的示例代码？

在C++中，GDAL（Geospatial Data Abstraction Library）是一个广泛使用的库，用于处理地理空间数据。如果你想要进行空间参考坐标系统的转换，你可以使用GDAL提供的函数，如`OGRCoordinateTransformation`或`OSRTransform`。下面是一个简单的示例，演示了如何使用GDAL库从一个坐标系转换到另一个坐标系：

首先，你需要包含必要的GDAL头文件，并链接GDAL库。假设你已经安装了GDAL并设置好环境变量。

#include <gdal/gdal.h>
#include <gdal/ogr_api.h>

然后，你需要创建一个OGRCoordinateTransformation对象，它代表两个坐标系之间的转换。这通常通过`OGRCreateCoordinateTransformation`函数完成：

// 创建源坐标系对象（例如EPSG:4326）
const char* src_srs = "EPSG:4326";
OGRSpatialReference* src_srs_obj = new OGRSpatialReference();
if (src_srs_obj->importFromWkt(src_srs) != OGRERR_NONE) {
    // 处理错误...
}

// 创建目标坐标系对象（例如EPSG:3857）
const char* dst_srs = "EPSG:3857";
OGRSpatialReference* dst_srs_obj = new OGRSpatialReference();
if (dst_srs_obj->importFromWkt(dst_srs) != OGRERR_NONE) {
    // 处理错误...
}

// 创建转换对象
OGRCoordinateTransformation* transform = OGRCreateCoordinateTransformation(src_srs_obj, dst_srs_obj);
if (!transform) {
    // 处理错误...
}

接下来，你可以使用这个转换对象将坐标从一种体系转换到另一种体系：

// 假设你有一个经纬度点（lon, lat）
double lon = -74.006017;
double lat = 40.712776;

// 转换
double x, y;
if (!transform->Transform(&x, &y, &lon, &lat)) {
    // 处理转换失败...
}

std::cout << "转换后的经度:" << x << ", 纬度:" << y << std::endl;

最后，记得在使用完转换对象后释放内存：

delete transform;
delete src_srs_obj;
delete dst_srs_obj;