Hello GDAL（简单的GDAL入门程序）

怼烎@ 2022-07-12 10:49 186阅读 0赞

Windows下的开发环境搭建参考博文：[GDAL开发环境搭建（VS2010 C++版）][GDAL_VS2010 C]

下面是一个很简单的读取GeoTiff格式的图像，然后输出其元数据信息的示例。  
其中比较重要的是使用GDAL的六参数模型，计算图像的左上角和右下角信息。

#include <iostream>
    #include <iomanip>
    #include "gdal_priv.h"
    
    using std::cout;
    
    /* @brief 计算图像行列号在给定坐标系下对应的地理坐标 @param x 行号 @param y 列号 @param coords 返回的地理坐标 @param transform 变换的六参数 */
    void toGeoCoord(int x, int y, double* coords, double* transform)
    {
        coords[0] = transform[0] + x * transform[1] + y * transform[2];
        coords[1] = transform[3] + x * transform[4] + y * transform[5];
    }
    
    /* @brief 读取数据并输出元数据信息 @param fileName 数据文件路径 @return 数据读取失败返回1，读取成功返回0 */
    int readGeoRaster(const char* fileName)
    {
        GDALAllRegister();  // 注册所有支持的数据格式驱动
        GDALDataset* dataset = static_cast<GDALDataset*>(GDALOpen(fileName, GA_ReadOnly));  // 以只读模式打开给定的数据文件
        if (!dataset)
        {
            cout << "影像读取失败：（" << fileName << "）！" << '\n';
            return EXIT_FAILURE;
        }
    
        int nSizeX = dataset->GetRasterXSize(); // 影像的宽度（像元数目）
        int nSizeY = dataset->GetRasterYSize(); // 影像的高度（像元数目）
        int nBandCount = dataset->GetRasterCount(); // 影像波段数
    
        cout << "图像大小：(" << nSizeX << ", " << nSizeY << ")\n";
        cout << "波段数：" << nBandCount << "\n";
    
        double adfGeoTransform[6];  // 存储图像的六参数信息
        dataset->GetGeoTransform(adfGeoTransform);  // 获取数据的六参数信息
    
        double adfULCoord[2];   // 图像左上角坐标
        double adfLRCoord[2];   // 图像右下角坐标
        // 获取左上角和右下角的地理坐标
        toGeoCoord(0, 0, adfULCoord, adfGeoTransform);
        toGeoCoord(nSizeX - 1, nSizeY - 1, adfLRCoord, adfGeoTransform);
    
        cout << std::fixed; // 对于double类型的数据可以使其正常输入，阻止默认的科学计数法输出，
        cout << "左上角坐标：(" << adfULCoord[0] << ", " << adfULCoord[1] << ")\n";
        cout << "右下角坐标：(" << adfLRCoord[0] << ", " << adfLRCoord[1] << ")\n";
    
        const char* pszProj = dataset->GetProjectionRef(); // 获得WKT形式的投影信息
        cout << "投影信息：" << pszProj << '\n';
    
        GDALClose(dataset); // 关闭数据集
    
        return EXIT_SUCCESS;
    }
    
    int main()
    {
        const char* pszFileName = R"(F:\2016\Data\Landsat\LT51240392010131BKT01\LT51240392010131BKT01_B1.TIF)";
        return readGeoRaster(pszFileName);
    }

[GDAL_VS2010 C]: http://blog.csdn.net/theonegis/article/details/25590269