GDAL中GDALDataType中值与其在C++中数据类型对应

GDAL中的GDALDataType是一个枚举型，其中的值为：

• GDT_Unknown : 未知数据类型
• GDT_Byte : 8bit正整型 (C++中对应unsigned char)
• GDT_UInt16 : 16bit正整型 (C++中对应 unsigned short)
• GDT_Int16 : 16bit整型 (C++中对应 short 或 short int)
• GDT_UInt32 : 32bit 正整型 (C++中对应unsigned long)
• GDT_Int32 : 32bit整型 (C++中对应int 或 long 或 long int)
• GDT_Float32 : 32bit 浮点型 (C++中对应float)
• GDT_Float64 : 64bit 浮点型 (C++中对应double)
• GDT_CInt16 : 16bit复整型 (C++中对应Complex Int16)
• GDT_CInt32 : 32bit复整型 (C++中对应Complex Int32)
• GDT_CFloat32 : 32bit复浮点型 (C++中对应Complex Float32)
• GDT_CFloat64 : 64bit复浮点型 (C++中对应Complex Float64)

今天在使用GDAL读取tif格式DEM/DSM数据时，发现读取到的高程值总是不对。在网上搜索，是因为定义的数据类型不对导致，在获取高程前，需要获取被读取的数据存储类型

GDALDataType gdal_data_type = poDataset->GetRasterBand(1)->GetRasterDataType();

然后根据不同数据类型，定义*pafScanblock1

    switch (gdal_data_type)
    {
    case GDT_Byte:
    {
        unsigned char *pafScanblock1;
        pafScanblock1 = (unsigned char *)CPLMalloc(sizeof(unsigned char)*(1)*(1));//建议能小则小否则会造成内存不足的情况
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_UInt16:
    {
        unsigned short *pafScanblock1;
        pafScanblock1 = (unsigned short *)CPLMalloc(sizeof(unsigned short)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_Int16:
    {
        short int *pafScanblock1;
        pafScanblock1 = (short int *)CPLMalloc(sizeof(short int)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_UInt32:
    {
        unsigned long *pafScanblock1;
        pafScanblock1 = (unsigned long *)CPLMalloc(sizeof(unsigned long)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_Int32:
    {
        long *pafScanblock1;
        pafScanblock1 = (long *)CPLMalloc(sizeof(long)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_Float32:
    {
        float *pafScanblock1;
        pafScanblock1 = (float *)CPLMalloc(sizeof(float)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    case GDT_Float64:
    {
        double *pafScanblock1;
        pafScanblock1 = (double *)CPLMalloc(sizeof(double)*(1)*(1));
        poBand1->RasterIO(GF_Read, dx, dy, 1, 1, pafScanblock1, 1, 1, GDALDataType(poBand1->GetRasterDataType()), 0, 0);
        cout << Xgeo << "  " << Ygeo << "  " << *pafScanblock1 << endl;
        break;
    }

    default: break;

    }