读取和创建光栅图（TIF）

GDALAllRegister();
OGRRegisterAll();
//读取光栅数据图（tif）
GDALDriver* pDriver=GetGDALDriverManager()->GetDriverByName("GTiff");
GDALDataset* dataset = (GDALDataset *)GDALOpen(m_sTifFile.c_str(),GA_ReadOnly);
if(!dataset)
{
string strError(CPLGetLastErrorMsg());
cout<<strError<<endl;
}
GDALRasterBand* pBand = dataset->GetRasterBand(1);
GDALDataType datatype = pBand->GetRasterDataType();
long *pBuf = new long[200*100];
dataset->RasterIO(GF_Read,200, 100, 200, 100,pBuf, 200, 100, GDT_Int32, 1, 0,0,0,0);

//存储光栅数据（tif）
GDALDriver* pDriverw=GetGDALDriverManager()->GetDriverByName("GTiff");
GDALDataset *pDataSetw=pDriverw->Create(m_sNewTif.c_str(),200,300,1,GDT_Int32,NULL);
GDALRasterBand* wband = pDataSetw->GetRasterBand(1);
wband->RasterIO(GF_Write,200, 100, 200, 100,pBuf, 200, 100, GDT_Int32, 0, 0);
GDALClose((GDALDatasetH)pDataSetw);