gdal读取txt文件_GDAL读取HDF和GRB2

nager() -&gt GetDriverByName("GTIFF");

GDALDataset * OPTemp;

OPTemp = podrive -&gtCreate(outFpath,X,Y,1,Band-&gtGetRasterDataType(),NULL);

GDALRasterBand * TPband ;

TPband = OPTemp -&gtGetRasterBand(1);

TPband -&gt RasterIO( GF_Write, 0,0,X,Y,bandvalue,X,Y,Band-&gtGetRasterDataType(),0,0);

GDALClose(OPTemp);

delete []bandvalue;

bandvalue = NULL;

GDALClose(OPTemp);

break;

}

case (GDT_UInt16):

{unsigned short* bandvalue = new unsigned short[X * Y];

Band -&gt RasterIO( GF_Read, 0,0,X,Y,bandvalue,X,Y,Band-&gtGetRasterDataType(),0,0);

GDALDriver * podrive;

podrive = GetGDALDriverManager() -&gt GetDriverByName("GTIFF");

GDALDataset * OPTemp;

OPTemp = podrive -&gtCreate(outFpath,X,Y,1,Band-&gtGetRasterDataType(),NULL);

GDALRasterBand * TPband ;

TPband = OPTemp -&gtGetRasterBand(1);

TPband -&gt RasterIO( GF_Write, 0,0,X,Y,bandvalue,X,Y,Band-&gtGetRasterDataType(),0,0);

GDALClose(OPTemp);

delete []bandvalue;

bandvalue = NULL;

break;

}

default:

{break;}

}

GDALClose((GDALDriver * )DataSet);

break;

}

}

GDALClose((GDALDriverH)pDataSet);

return 0;

}

gdal 读取grb2格式

const int blockSize = 200;

template &lttypename T&gt

void BlockImgConvertProcess( int lblockXSize,int lblockYSize, int blockSize,int bandIndex,

int j,int k,GDALRasterBand* poBand,GDALRasterBand* poBand_img)

{

T *ppdata = new T[lblockXSize * lblockYSize];

poBand-&gtRasterIO( GF_Read, j * blockSize, k * blockSize, lblockXSize, lblockYSize,

ppdata, lblockXSize, lblockYSize, poBand-&gtGetRasterDataType(),

0, 0 );

poBand_img -&gt RasterIO ( GF_Write, j * blockSize, k * blockSize, lblockXSize, lblockYSize,

ppdata,lblockXSize, lblockYSize, poBand_img-&gtGetRasterDataType(), 0, 0);

delete[] ppdata;

ppdata = NULL;

}

bool GRBToImg(const char* pchTifPath, const char* pchImgPath)

{

string strTifPath = pchTifPath;

string strTifExt = strTifPath.substr(strTifPath.rfind('.') + 1, strTifPath.length());

string strImgPath = pchImgPath;

string strImgExt = strImgPath.substr(strImgPath.rfind('.') + 1, strImgPath.length());

transform(strTifExt.begin(),strTifExt.end(),strTifExt.begin(), ::toupper);

transform(strImgExt.begin(),strImgExt.end(),strImgExt.begin(), ::toupper);

if(strTifExt != "GRB2" || strImgExt != "TIF")

return false;

GDALDataset *poDataset;

//打开图像

poDataset = (GDALDataset *)GDALOpen( pchTifPath, GA_ReadOnly );

if( poDataset == NULL )

{

return false;

}

//打开图像

//获取图像信息

double adfGeoTransform[6];

double *maomao = new double[6];

long lWidth = poDataset-&gtGetRasterXSize();

long lHeight = poDataset-&gtGetRasterYSize();

long lBandCount = poDataset-&gtGetRasterCount();

const char* chWkt = poDataset-&gtGetProjectionRef();

poDataset-&gtGetGeoTransform( adfGeoTransform );

//获取图像信息

vector&ltstring&gt bandscrip;