GDAL实现投影转换功能

使用GDAL库来实现坐标投影转换是一个非常好的选择，因为GDAL提供了强大的地理空间数据处理功能。以下是一个使用Qt和C++结合GDAL实现的完整代码示例，满足你的需求，包括设置坐标系转换、自定义投影、文件支持和转换参数等功能。

环境准备
确保安装以下依赖项：

GDAL - 用于处理地理空间数据的库
Qt - 用于界面开发
C++编译器 - 用于编译代码
安装GDAL
可以使用以下命令安装GDAL（以Ubuntu为例）：
sudo apt-get install gdal-bin libgdal-dev
Qt工程设置
在你的CMakeLists.txt或pro文件中添加GDAL库的链接：

cmake
复制代码
find_package(GDAL REQUIRED)
target_link_libraries(your_project_name GDAL)
C++代码示例
以下是一个实现投影转换功能的完整代码示例：

#include <QCoreApplication>
#include <QString>
#include <iostream>
#include <gdal_priv.h>
#include <ogr_spatialref.h>

// 投影转换类
class ProjectionTransformer {
public:
    // 构造函数，初始化源和目标投影
    ProjectionTransformer(const std::string& srcProj, const std::string& dstProj) {
        // 初始化GDAL
        GDALAllRegister();

        // 创建源投影
        srcSRS = new OGRSpatialReference();
        srcSRS->importFromProj4(srcProj.c_str());

        // 创建目标投影
        dstSRS = new OGRSpatialReference();
        dstSRS->importFromProj4(dstProj.c_str());

        // 创建转换对象
        transform = OGRCreateCoordinateTransformation(srcSRS, dstSRS);
        if (transform == nullptr) {
            throw std::runtime_error("无法创建坐标转换对象");
        }
    }

    // 析构函数，释放资源
    ~ProjectionTransformer() {
        delete srcSRS;
        delete dstSRS;
        if (transform != nullptr) {
            OCTDestroyCoordinateTransformation(transform);
        }
    }

    // 坐标转换函数
    bool transformPoint(double& x, double& y, double& z) {
        return transform->Transform(1, &x, &y, &z);
    }

    // 文件转换函数
    bool transformFile(const std::string& inputFile, const std::string& outputFile) {
        GDALDataset* dataset = (GDALDataset*)GDALOpen(inputFile.c_str(), GA_ReadOnly);
        if (!dataset) {
            std::cerr << "无法打开文件: " << inputFile << std::endl;
            return false;
        }

        // 设置投影
        dataset->SetProjection(srcSRS->ExportToWkt());
        GDALDataset* outDataset = (GDALDataset*)GDALCreateCopy(
            GDALGetDriverByName("GTiff"), outputFile.c_str(), dataset, FALSE, nullptr, nullptr, nullptr);

        GDALClose(dataset);
        GDALClose(outDataset);
        return true;
    }

private:
    OGRSpatialReference* srcSRS;
    OGRSpatialReference* dstSRS;
    OGRCoordinateTransformation* transform;
};

int main(int argc, char *argv[]) {
    QCoreApplication app(argc, argv);

    // 示例：自定义投影转换参数
    const std::string srcProj = "+proj=longlat +datum=WGS84 +no_defs";
    const std::string dstProj = "+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84";

    try {
        ProjectionTransformer transformer(srcProj, dstProj);

        // 示例：单点坐标转换
        double x = 116.397, y = 39.907, z = 0;
        if (transformer.transformPoint(x, y, z)) {
            std::cout << "转换后的坐标：" << x << ", " << y << ", " << z << std::endl;
        } else {
            std::cerr << "坐标转换失败" << std::endl;
        }

        // 示例：文件转换
        if (transformer.transformFile("input.tif", "output.tif")) {
            std::cout << "文件转换成功" << std::endl;
        } else {
            std::cerr << "文件转换失败" << std::endl;
        }

    } catch (const std::exception& ex) {
        std::cerr << "错误: " << ex.what() << std::endl;
    }

    return app.exec();
}

代码说明
构造函数：ProjectionTransformer类的构造函数接收源和目标投影参数的字符串（Proj.4格式），并创建相应的投影系统和转换对象。
transformPoint：将单个点（x, y, z）从源坐标系转换到目标坐标系。
transformFile：将整个文件的投影从源坐标系转换到目标坐标系。
异常处理：在无法创建转换对象时抛出异常，以便捕获和处理错误。
注意事项
坐标参数：代码中的srcProj和dstProj为WGS84和墨卡托投影的Proj.4格式，可以根据需求更改。
投影文件：GDAL支持多种文件格式，GDALCreateCopy用于将源文件复制为新的投影。
多线程：如果需要多线程，可以将ProjectionTransformer类的实例化和调用移到线程内。