Hello GDAL（简单的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);
}