《影像数据格式访问开源包GDAL.ppt》由会员分享，可在线阅读，更多相关《影像数据格式访问开源包GDAL.ppt（16页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、影像数据格式访问开源包GDAL Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望GDAL简介官网：http:/www.gdal.org/GDAL是一个开源的栅格空间数据和影像数据的读写、转换器。目前支持100多种栅格格式（http:/www.gdal.org/formats_list.html）C/C+版GDAL数据模型核心数据Dataset（数据集）Raster Band（波段）Metadata（元数据）坐标系统Coordinate System（坐标系统）Aff
2、ine GeoTransform（仿射变换）GCP（地面控制点）Color Table（颜色表）Overviews（金字塔层）让格式得到支持编译时打开格式开关有些是默认的编译的，有些是需修改后重新编译的参看http:/www.gdal.org/formats_list.html程序中获取格式的Driver注册格式的Driver注册所有格式注册所有格式GDALAllRegister();注意：必须在编译时打开格式开关。通常在程序中调用其他GDAL函数之前，必须首先调用此函数反动作GDALDestroyDriverManager();打开数据集GDALOpenGDALDatasetH GDALOp
3、en(const char*pszFilename,GDALAccess eAccess)GDALAccess:GA_ReadOnly：只读 GA_Update ：读写可自动识别数据格式。成功返回句柄，失败返回NULL。GDALClose关闭数据集获取基本信息GDALGetRasterXSize获取列数GDALGetRasterYSize获取行数GDALGetRasterCount获取波段数GDALGetGeoTransform获取仿射变换参数GDALGetGCPs获取控制点获取波段GDALGetRasterBandHGDALRasterBandH GDALGetRasterBand (GDA
4、LDatasetH hDS,int nBandId)nBandId：从1起算成功返回波段句柄，失败返回NULL获取波段基本信息GDALGetRasterBandXSizeGDALGetRasterBandYSizeGDALGetRasterDataTypeGDALGetRasterColorTableGDALGetDefaultRAT按块(block)读写栅格数据按block读写可以获得很高的性能GDALGetBlockSize：获取一个标准block的大小（横向和纵向像素数）ReadBlock：读一个blockWriteBlock：写一个block计算block数目：GDALGetBlock
5、Size(hBand,&nXBlockSize,&nYBlockSize);nXBlocks=(GDALRasterBandGetXSize(hBand)+nXBlockSize-1)/nXBlockSize;nYBlocks=(GDALRasterBandGetYSize(hBand)+nYBlockSize-1)/nYBlockSize;GDALReadBlock (GDALRasterBandH hBand,int nXOff,int nYOff,void*pData )必须预先分配足够的内存给pData。CPLErr GDALWriteBlock (GDALRasterBandH hB
6、and,int nXOff,int nYOff,void*pData )读写任意区域栅格数据CPLErr GDALRasterIO (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进行类型转换。像
7、素坐标和地理坐标仿射转换获取转换系数：doublea6;GDALGetGeoTransform(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+a4*c+a5*
8、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。