gdal读写图像分块处理（精华版）.doc

gdal读写图像分块处理（精华版） /* 版权所有者：赵文??*/? ????????????????????????????????????????????????????????????? 一．gdal进行数据操作在安装好gdal后，即可调用gdal库中的函数。（需要包含的头文件：gdal_priv.h）1．打开数据集使用gdal库进行数据（影像）操作的第一步就是打开一个数据集。对于“数据集”这个名词大家可能不会太习惯，但是对于一般的格式来说，一个“数据集”就是一个文件，比如一个TIFF文件就是一个以tiff为扩展名的文件。但是对于众多RS数据来说，一个数据集包含的绝对不仅仅是一个文件。对于很多RS数据，他们把一张图像分成数个图像文件，然后放在一个文件夹中，用一些额外的文件来组织它们之间的关系，形成一个“数据集”(有点难以理解，暂且放过)。下面我们由给定的文件路径文件名打开一个tiff(GeoTIFF)文件。（任何支持的格式，打开方式都是这样）CString strFilePath;StrFilePath=’d:/rsdata/2005_234.tif’;GDALDataSet *poDataset; //GDAL数据集GDALAllRegister();poDataset = (GDALDataset *) GDALOpen(strFilePath, GA_ReadOnly );这样我们就打开了这个文件。通过数据集poDataset即可调用各功能函数，如:GetRasterCount();//获取图像波段数；GetRasterXSize();//获取图像宽度GetRasterYSize();//获取图像高度GetRasterBand();//获取图像某一波段GetGeoTransform(double *p);//获取图像地理坐标信息长度为六的数组RasterIO();//对图像数据进行缩放读和写……(更具体的API列表可以看这里)。 2．获取图像信息、读取图像打开文件后，下面要做的就是获取文件的相关信息保存在相应变量中，将图像数据读入内存中，等待后续处理了。2．1 获取基本信息因为不同格式数据所包含的相关信息有所不同，一般情况下我们需要得到图像的高、宽、波段数、地理坐标信息，数据类型等。Gdal库中有相应的函数实现这些功能。下面的代码实现获取这些信息：int????nBandCount=poDataset-GetRasterCount();int??nImgSizeX=poDataset-GetRasterXSize();int???? nImgSizeY=poDataset-GetRasterYSize();double????????adfGeoTransform[6];poDataset-GetGeoTransform( adfGeoTransform );//如果图像不含地理坐标信息，默认返回值是：（0、1、0、0、0、1），表中第四列表示了带//有地理坐标信息的数据格式// adfGeoTransform[0]是左上角像元的东坐标；// adfGeoTransform[3]是左上角像元的北坐标；// adfGeoTransform[1]是像元宽度；// adfGeoTransform[5]是像元高度；//图像其他点坐标计算公式：//Xp = adfGeoTransform [0] + P*adfGeoTransform [1]+L*adfGeoTransform [2];//Yp = adfGeoTransform [3] + P*adfGeoTransform [4] + L*adfGeoTransform [5];//注：用GetGeoTransform()并不能十分合理的表示图像地理坐标，当影像范围很大//时，这种坐标表示方法将不适用。2．2 将图像数据按照要求读入内存图像的读写是通过RasterIO()实现的。RasterIO()功能十分强大,它可以把图像上指定大小的矩形象素块以缩放的形式按指定的数据类型输出或输入到用户指定大小的缓冲区中。原型：CPLErr GDALDataset::RasterIO(GDALRWFlag eRWFlag,??????//读写标记如果为GF_Read，则是将影像内容写入内存，如果//为GF_Write，则是将内存中内容写入文件。int nXOff, int nYOff,????//相对于图像左上角顶点(从零开始)的行列偏移量 int nXSize, int nYSize,?????? //要读写的块在x方向的象素个数和y方向的象素列数void * pData,??