使用GDAL进行栅格数据读取示例使用GDAL进行栅格数据读取示例

当然，使用GDAL可以对其支持的数据进行操作，不局限在栅格数据。但目前主要解决的问题是栅格数据，更具体一点是针对ENVI的标准数据格式。
以下程序仅作为一个简单的应用示例，主要参考了http://www.gdal.org/gdal_tutorial.html中的内容。使用时需要用到一下头文件：
cpl_config.h
cpl_conv.h
cpl_error.h
cpl_minixml.h
cpl_port.h
cpl_string.h
cpl_vsi.h
gdal.h
gdal_frmts.h
gdal_priv.h
gdal_version.h
另外需要连接gdal.lib（静态链接）或者gdal_i.lib（动态链接）。程序没有特别的功能，仅作示例而已。
gdal_sample.cpp
///
#include "./gdal/gdal_priv.h"
int main(int argc, char **argv)
{
if(argc != 2)
{
printf( "Usage: gdal_sample /n" );
return 0;
}
char *pszFilename;
pszFilename = argv[1];

// Opening the File
GDALDataset*poDataset;
GDALAllRegister();
poDataset = (GDALDataset *) GDALOpen( pszFilename, GA_ReadOnly );
if( poDataset != NULL )
{
// Getting Dataset Information
doubleadfGeoTransform[6];
printf( "Driver: %s/%s/n",
poDataset->GetDriver()->GetDescription(),
poDataset->GetDriver()->GetMetadataItem( GDAL_DMD_LONGNAME ) );
printf( "Size is %dx%dx%d/n",
poDataset->GetRasterXSize(), poDataset->GetRasterYSize(),
poDataset->GetRasterCount() );
if( poDataset->GetProjectionRef()!= NULL )
printf( "Projection is `%s'/n", poDataset->GetProjectionRef() );
if( poDataset->GetGeoTransform( adfGeoTransform ) == CE_None )
{
printf( "Origin = (%.6f,%.6f)/n",
adfGeoTransform[0], adfGeoTransform[3] );
printf( "Pixel Size = (%.6f,%.6f)/n",
adfGeoTransform[1], adfGeoTransform[5] );
}
// Fetching a Raster Band
GDALRasterBand*poBand;
intnBlockXSize, nBlockYSize;
intbGotMin, bGotMax;
doubleadfMinMax[2];

poBand = poDataset->GetRasterBand( 1 );
poBand->GetBlockSize( &nBlockXSize, &nBlockYSize );
printf( "Block=%dx%d Type=%s, ColorInterp=%s/n",
nBlockXSize, nBlockYSize,
GDALGetDataTypeName(poBand->GetRasterDataType()),
GDALGetColorInterpretationName(
poBand->GetColorInterpretation()) );
adfMinMax[0] = poBand->GetMinimum( &bGotMin );
adfMinMax[1] = poBand->GetMaximum( &bGotMax );
if( ! (bGotMin && bGotMax) )
GDALComputeRasterMinMax((GDALRasterBandH)poBand, TRUE, adfMinMax);
printf( "Min=%.3fd, Max=%.3f/n", adfMinMax[0], adfMinMax[1] );

if( poBand->GetOverviewCount() > 0 )
printf( "Band has %d overviews./n", poBand->GetOverviewCount() );
if( poBand->GetColorTable() != NULL )
printf( "Band has a color table with %d entries./n",
poBand->GetColorTable()->GetColorEntryCount() );
// Reading Raster Data
float *pafScanline;
intnXSize = poBand->GetXSize();
pafScanline = (float *) CPLMalloc(sizeof(float)*nXSize);
poBand->RasterIO( GF_Read, 0, 0, nXSize, 1,
pafScanline, nXSize, 1, GDT_Float32,
0, 0 );
// Closing the Dataset
GDALClose(poDataset);
}
【使用GDAL进行栅格数据读取示例】return 0;
}

使用GDAL进行栅格数据读取示例

推荐阅读

洛阳牡丹怎么养，洛阳牡丹花如何养

虎头茉莉怎么施肥虎头茉莉怎么养

战歌竞技场亡灵战阵容怎么组战歌竞技场亡灵战阵容分享

小学毕业寄语50句话小学毕业寄语50句

脂肪型肥胖的危害有哪些引起肥胖的三类脂肪

光火神殿为什么在秩序之战闭门不出？

我最近经常吃辣条，会口臭，是辣条造成的吗？

胡忠信名字打分102分

翠玉是玉还是翡翠翠玉属于翡翠吗

精油皂和洗面奶哪个洗脸好

吃菠萝蜜会发胖吗

gmc是什么车gmc是哪个公司的车

淘宝不认tm商标淘宝如何申请商标入驻，入驻淘宝品牌库显示近似商标

空调外机没有水滴下来怎么回事美的空调制冷室外机不滴水怎么回事

抖音极速版怎么合拍不了抖音极速版怎么合拍

信息流推荐系统智能交付解决方案探索

佳能废墨垫清零教程视频下载安装

到山西做客，主家很少用青花汾招待，这4款好酒才是宴请的首选

欧洲的牧场和食品库之称的国家是

520和521代表什么意思 520和521区别是什么