gdal读写图像分块处理(精华版)

//获取格式类型 char **papszMetadata = poDriver->GetMetadata();
//根据文件路径文件名，图像宽，高，波段数，数据类型，文件类型，创建新的 数据集 poDataset2=poDriver->Create(strFilePath2,nImgSizeX,nImgSizeY,nBandCount,GDT _Byte,papszMetadata); //坐标赋值 double adfGeoTransform[6]={0，1，0，0，0，1}; poDataset2->SetGeoTransform(adfGeoTransform);
int nLineSpace, //y 方向上相邻两行之间的字节偏移, 默认为 0,则行间的实际字
节
//偏移为 eBufType * nBufXSize
int nBandSpace //相邻两波段之间的字节偏移，默认为 0,则意味着波段是顺序结
构
//的，其间字节偏移为 nLineSpace * nBufYSize
. . 『中间对图像进行处理运算』 . //将图像数据写入新图像中 poDataset2->RasterIO(GF_Write, 0,0, nImgSizeX, nImgSizeY, ppafScan,pafsizex,pafsizey, GDT_Byte,nBandCount,0,0, 0,0 );
);
下面将通过实例演示其使用方法,实现的是将 7 波段图像中的第 2 3 4 波段按照
3 4 2 的顺序读入内存中,逐像素存储: int bandmap[7]; bandmap[0]=3; bandmap[1]=4; bandmap[2]=2; Scanline=(nImagSizex*8+31)/32*4; BYTE pafScan=( BYTE )CPLMalloc(sizeof(byte) *nImgSizeX*nImgSizeY*３) poDataset->RasterIO( GF_Read, 0, 0,nImgSizeX,nImgSizeY, pafScan, nImgSizeX,nImgSizeY,GDT_Byte,３,bandmap,３, Scanline*３,1 );
void * pData,
//指向目标缓冲区的指针，由用户分配
int nBufXSize, int nBufYSize,//目标块在 x 方向上和 y 方向上的大小
GDALDataType eBufType, //目标缓冲区的数据类型，原类型会自动转换为目标
类型
int nBandCount,
//要处理的波段数
//时，这种坐标表示方法将不适用。
2．2 将图像数据按照要求读入内存
图像的读写是通过 RasterIO()实现的。RasterIO()功能十分强大,它可以把图像上指
定大小的矩形象素块以缩放的形式按指定的数据类型输出或输入到用户指定大
小的缓冲区中。
原型：
CPLErr GDALDataset::RasterIO(
CString strFilePath; StrFilePath=’d:/rsdata/2005_234.tif’; GDALDataSet *poDataset; //GDAL 数据集 GDALAllRegister(); poDataset = (GDALDataset *) GDALOpen(strFilePath, GA_ReadOnly );
int * panBandMap, //记录要操作的波段的索引(波段索引从 1 开始)的数组，若
为空
//则数组中存放的是前 nBandCount 个波段的索引
int nPixelSpace, //X 方向上两个相邻象素之间的字节偏移，默认据类型 eBufType 确定
delete poDataset2; delete poDriver; //图像另存完毕
CreateCopy 方式创建 jpeg 格式文件： 接上面的过程，先不 delete，（即已经完成用 create 方式先将运算完毕的图像创建 为 tiff 格式） fomat="Jpeg" poDriver = GetGDALDriverManager()->GetDriverByName(fomat); poDataset3=poDriver->CreateCopy(strFilePath3,poDataset2,1,papszMetadata,NULL,
//将原图像数据读出，进行相应处理后，写入新文件 BYTE *ppafScan= new BYTE [nImgSizeX * nImgSizeY *nBandCount]; poDataset1->RasterIO( GF_Read,0,0, nImgSizeX, nImgSizeY, ppafScan, nImgSizeX, nImgSizeY, GDT_Byte,nBandCount,0,0, 0,0 );
这样我们就打开了这个文件。通过数据集 poDataset 即可调用各功能函数，如: GetRasterCount();//获取图像波段数； GetRasterXSize();//获取图像宽度 GetRasterYSize();//获取图像高度 GetRasterBand();//获取图像某一波段 GetGeoTransform(double *p);//获取图像地理坐标信息长度为六的数组 RasterIO();//对图像数据进行缩放读和写 …… (更具体的 API 列表可以看这里)。
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()并不能十分合理的表示图像地理坐标，当影像范围很 大
以这种方式读取之后，直接可构建位图进行显示。这里可以按照自己的需要进行 其他方式读取。以上读取方式仅仅为了显示方便，如进行图像处理相关运算，则 按波段全部读出会比较方便: poDataset->RasterIO( GF_Read, 0, 0,nImgSizeX,nImgSizeY, pafScan, nImgSizeX,nImgSizeY,GDT_Byte,bandcount,0,0,0,0); 之前开辟内存改为： BYTE pafScan=new byte[nImgSizeX*nImgSizeY*bandcount];
将图像数据读入内存后，即可通过指针 pafScan 对图像进行你想要进行的操作 了。
3．另存图像 另存文件其实就是先创建一个新的文件，然后将数据写入新文件中。 使用 gdal 创建新文件有两种方式：Create()和 CreateCopy();有些文件格式支持
Create()函数（见第一页表格第三列），可以使用 Create（）直接创建此类格式的 文件，而其他不支持 Create（）函数的图像格式，需要先创建 tiff 格式文件，然 后复制生成目标格式文件。 CString strFilePath1;//输入图像文件路径名 CString strFilePath2;//另存为的 tiff 格式图像路径 CString strFilePath2; //另存为的 jpeg 格式图像路径 GDALDataset *poDataset1; //GDAL 数据集 GDALDataset *poDataset2; //待创建的 GDAL 数据集 GDALDataset *poDataset2; //待创建的 GDAL 数据集 GDALDriver *poDriver; //驱动，用于创建新的文件
2．获取图像信息、读取图像 打开文件后，下面要做的就是获取文件的相关信息保存在相应变量中，将图像数 据读入内存中，等待后续处理了。
2．1 获取基本信息 因为不同格式数据所包含的相关信息有所不同，一般情况下我们需要得到图像的 高、宽、波段数、地理坐标信息，数据类型等。Gdal 库中有相应的函数实现这 些功能。下面的代码实现获取这些信息：
gdal 读写图像分块处理（精华版） /* 版权所有者：赵文 */
一．gdal 进行数据操作 在安装好 gdal 后，即可调用 gdal 库中的函数。 （需要包含的头文件：gdal_priv.h） 1．打开数据集 使用 gdal 库进行数据（影像）操作的第一步就是打开一个数据集。对于“数据集” 这个名词大家可能不会太习惯，但是对于一般的格式来说，一个“数据集”就是一 个文件，比如一个 TIFF 文件就是一个以 tiff 为扩展名的文件。但是对于众多 RS 数据来说，一个数据集包含的绝对不仅仅是一个文件。对于很多 RS 数据，他们 把一张图像分成数个图像文件，然后放在一个文件夹中，用一些额外的文件来组 织它们之间的关系，形成一个“数据集”(有点难以理解，暂且放过)。下面我们由 给定的文件路径文件名打开一个 tiff(GeoTIFF)文件。（任何支持的格式，打开方 式都是这样）