影像数据格式访问开源包GDAL

《影像数据格式访问开源包GDAL》由会员分享，可在线阅读，更多相关《影像数据格式访问开源包GDAL（16页珍藏版）》请在装配图网上搜索。

1、GDALGDAL简介官网：http:/www.gdal.org/GDAL是一个开源的栅格空间数据和影像数据的读写、转换器。目前支持100多种栅格格式（http:/www.gdal.org/formats_list.html）C/C+版GDAL数据模型核心数据Dataset（数据集）Raster Band（波段）Metadata（元数据）坐标系统Coordinate System（坐标系统）Affine GeoTransform （仿射变换）GCP（地面控制点）Color Table（颜色表）Overviews（金字塔层）让格式得到支持编译时打开格式开关有些是默认的编译的，有些是需修改后重新编译

2、的参看http:/www.gdal.org/formats_list.html程序中获取格式的Driver注册格式的Driver注册所有格式注册所有格式GDALAllRegister();注意：必须在编译时打开格式开关。通常在程序中调用其他GDAL函数之前，必须首先调用此函数反动作GDALDestroyDriverManager();打开数据集GDALOpenGDALDatasetH GDALOpen (const char *pszFilename, GDALAccess eAccess) GDALAccess:GA_ReadOnly ：只读 GA_Update ：读写可自动识别数据格式。成

3、功返回句柄，失败返回NULL。GDALClose关闭数据集获取基本信息GDALGetRasterXSize获取列数GDALGetRasterYSize获取行数GDALGetRasterCount获取波段数GDALGetGeoTransform获取仿射变换参数GDALGetGCPs获取控制点获取波段GDALGetRasterBandHGDALRasterBandH GDALGetRasterBand ( GDALDatasetH hDS, int nBandId) nBandId：从1起算成功返回波段句柄，失败返回NULL获取波段基本信息GDALGetRasterBandXSizeGDALGet

4、RasterBandYSizeGDALGetRasterDataTypeGDALGetRasterColorTableGDALGetDefaultRAT按块(block)读写栅格数据按block读写可以获得很高的性能GDALGetBlockSize：获取一个标准block的大小（横向和纵向像素数）ReadBlock：读一个blockWriteBlock：写一个block计算block数目： GDALGetBlockSize(hBand, &nXBlockSize, &nYBlockSize ); nXBlocks = (GDALRasterBandGetXSize(hBand) + nXBlo

5、ckSize - 1) / nXBlockSize; nYBlocks = (GDALRasterBandGetYSize(hBand) + nYBlockSize - 1) / nYBlockSize;GDALReadBlock ( GDALRasterBandH hBand, int nXOff, int nYOff, void * pData ) 必须预先分配足够的内存给pData。CPLErr GDALWriteBlock ( GDALRasterBandH hBand, int nXOff, int nYOff, void * pData ) 读写任意区域栅格数据CPLErr GDA

6、LRasterIO ( GDALRasterBandH hBand, GDALRWFlag eRWFlag, int nXOff, int nYOff, int nXSize, int nYSize, void * pData, int nBufXSize, int nBufYSize, GDALDataType eBufType, int nPixelSpace, int nLineSpace ) 注意：性能可能不高。如果是读数据，必须分配足够的内存给pData。本函数会自动根据eBufType进行类型转换。像素坐标和地理坐标仿射转换获取转换系数：doublea6;GDALGetGeoTra

7、nsform( hDataset, a );转换函数void World2Pixel(double a6, double x, double y, int* r, int* c)*r = (int)(a4*x - a1*y - a4*a0 + a1*a3)/(a4*a2 - a1*a5);*c = (int)(a5*x - a2*y - a5*a0 + a2*a3)/(a5*a1 - a2*a4);void Pixel2World(double a6, int r, int c, double* x, double* y) *x = a0 + a1 * c + a2 * r; *y = a3

8、+ a4 * c + a5 * r;获取栅格总范围将（0，0）和（XSize, YSize）传递给Pixel2World函数即可得到对角坐标。获取格式驱动GDALDriverH GDALGetDriverByName (const char * pszName)pszName可从源代码或http:/www.gdal.org/formats_list.html的Code项获得。创建新的栅格文件GDALDatasetH GDALCreate ( GDALDriverH hDriver, const char * pszFilename, int nXSize, int nYSize, int nBands, GDALDataType eBandType, char * papszOptions ) 若成功返回数据集句柄，失败则返回NULL。papszOptions与具体的格式有关，在很多情况下可为NULL。