使用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进行栅格数据读取示例

推荐阅读

广东有哪些比较好的中职学校

怎么看电脑的详细配置显卡，电脑怎么查看配置显卡

2023济南哪些技能证书能申请企业职工技能提升补贴？

夜书所见书的意思夜书所见书的意思简述

mysql安装后没有mysql数据库安装的mysql没有服务

浅析ADAMoracle多重优势特征

尼康85mm移轴镜头尼康85mm

爱普生4100清零软件

吃猪血可以清肺吗

西门子冰箱冷冻吸不紧常见故障,以及处理方式？

秋梨膏怎么熬如何在家熬制秋梨膏

一般纳税人信用等级怎么查纳税信用评价信息在哪里查询，纳税信用评价结果报告

九寨沟出现蓝冰的原因

折耳猫养殖注意事项_折耳猫养殖注意事项有哪些

佳能单反档位使用技巧佳能单反常用档

怎么把ftp服务器的图片在网页中显示？？不是本地显,ftp服务器实现图片缩略图的解决办法

mysql数据库视图的优缺点视图mysql数据库

什么是响应式布局响应式布局的解释

ppt图表制作如何添加数据分析,PPT如何添加图表

简单IOC的实现