GDAL对遥感图像的读取过程

GDAL可以读取多种遥感图像格式，以下是读取遥感图像的基本步骤：

1. 打开遥感图像文件：使用GDAL的Open函数打开遥感图像文件，例如GeoTIFF格式的遥感图像。

2. 获取遥感图像信息：使用GDAL的GetMetadata函数获取遥感图像的元数据信息，例如图像大小、坐标系、波段数等。

3. 读取遥感图像数据：使用GDAL的RasterIO函数读取遥感图像的像元值，可以按照指定的位置和大小读取遥感图像数据。

4. 关闭遥感图像文件：使用GDAL的Close函数关闭遥感图像文件。

以下是一个读取GeoTIFF格式遥感图像的示例代码：

#include "gdal_priv.h"
#include "cpl_conv.h" // for CPLMalloc()

int main()
{
    GDALDataset  *poDataset;
    GDALAllRegister();
    poDataset = (GDALDataset *) GDALOpen("example.tif", GA_ReadOnly );
    if( poDataset == NULL )
    {
        printf( "Open failed.\n" );
        exit( 1 );
    }
    int nBands = poDataset->GetRasterCount();
    int nXSize = poDataset->GetRasterXSize();
    int nYSize = poDataset->GetRasterYSize();
    double adfGeoTransform[6];
    if( poDataset->GetGeoTransform( adfGeoTransform ) == CE_None )
    {
        printf( "Origin = (%.6f,%.6f)\n",
                adfGeoTransform[0], adfGeoTransform[3] );
        printf( "Pixel Size = (%.6f,%.6f)\n",
                adfGeoTransform[1], adfGeoTransform[5] );
    }
    const char* pszProjection = poDataset->GetProjectionRef();
    printf("Projection = %s\n", pszProjection );
    float* pafScanline;
    pafScanline = (float*) CPLMalloc(sizeof(float)*nXSize);
    for(int i=0;i<nBands;i++)
    {
        GDALRasterBand  *poBand;
        poBand = poDataset->GetRasterBand(i+1);
        poBand->RasterIO( GF_Read, 0, 0, nXSize, nYSize,
                          pafScanline, nXSize, nYSize,
                          GDT_Float32, 0, 0 );
    }
    CPLFree(pafScanline);
    GDALClose( (GDALDatasetH) poDataset );
    return 0;
}

该示例代码中，首先使用GDALOpen函数打开名为"example.tif"的GeoTIFF格式遥感图像文件。然后使用GDAL的GetRasterCount、GetRasterXSize、GetRasterYSize、GetGeoTransform、GetProjectionRef等函数获取遥感图像的信息。最后使用RasterIO函数读取遥感图像数据，这里只读取了一个波段的数据。读取完数据后，使用GDALClose函数关闭遥感图像文件。