gdal库rasteri函数

GDAL (Geospatial Data Abstraction Library) 是一个开源的地理空间数据处理库，其中包含了很多处理栅格数据的函数，其中一个常用的函数就是 `RasterIO`。

`RasterIO` 函数用于读取或写入栅格数据，其语法如下：

CPLErr RasterIO( GDALRWFlag eRWFlag,
                 int nXOff, int nYOff, int nXSize, int nYSize,
                 void * pData, int nBufXSize, int nBufYSize,
                 GDALDataType eBufType,
                 int nBandCount, int *panBandMap,
                 int nPixelSpace, int nLineSpace, int nBandSpace );

参数说明：

- `eRWFlag`: 读取或写入标志，可选值为 `GF_Read` 或 `GF_Write`。
- `nXOff`、`nYOff`: 要读取或写入的栅格数据的左上角像素的 X、Y 坐标。
- `nXSize`、`nYSize`: 要读取或写入的栅格数据的宽度和高度（单位为像素）。
- `pData`: 存储读取或写入数据的缓冲区，必须是一个指向足够大的内存块的指针。
- `nBufXSize`、`nBufYSize`: 缓冲区的宽度和高度（单位为像素）。
- `eBufType`: 缓冲区的数据类型，可选值包括：`GDT_Byte`、`GDT_Int16`、`GDT_UInt16`、`GDT_Int32`、`GDT_UInt32`、`GDT_Float32`、`GDT_Float64`、`GDT_CInt16`、`GDT_CInt32`、`GDT_CFloat32` 和 `GDT_CFloat64`。
- `nBandCount`、`panBandMap`: 要读取或写入的波段数和波段映射数组。
- `nPixelSpace`、`nLineSpace`、`nBandSpace`: 存储像素、行和波段之间的字节间距（单位为字节）。

示例代码：

#include "gdal.h"

int main()
{
    GDALAllRegister();

    GDALDataset *poDS = (GDALDataset*) GDALOpen("test.tif", GA_ReadOnly);
    if (poDS == NULL)
    {
        // 打开数据集失败
        return 1;
    }

    int nXSize = poDS->GetRasterXSize();
    int nYSize = poDS->GetRasterYSize();
    int nBands = poDS->GetRasterCount();

    int nBufXSize = nXSize;
    int nBufYSize = nYSize;
    void *pData = new unsigned char[nXSize * nYSize * nBands];

    int nPixelSpace = GDALGetDataTypeSize(poDS->GetRasterBand(1)->GetRasterDataType()) / 8;
    int nLineSpace = nXSize * nPixelSpace;
    int nBandSpace = nXSize * nYSize * nPixelSpace;

    CPLErr eErr = poDS->RasterIO(GF_Read, 0, 0, nXSize, nYSize, pData, nBufXSize, nBufYSize, GDT_Byte, nBands, NULL, nPixelSpace, nLineSpace, nBandSpace);

    if (eErr != CE_None)
    {
        // 读取数据失败
        delete[] pData;
        GDALClose(poDS);
        return 1;
    }

    // 处理读取到的数据

    delete[] pData;
    GDALClose(poDS);

    return 0;
}

该示例代码用 `RasterIO` 函数读取了一个名为 "test.tif" 的栅格数据，并将数据存储在一个缓冲区中。要注意的是，在使用 `RasterIO` 函数之前，需要调用 `GDALAllRegister` 函数进行初始化，并且需要在最后调用 `GDALClose` 函数释放资源。