GDAL--API手册.pdf

在打开 GDAL 所支持的光栅数据之前需要注册驱动。这里的驱动是针对 GDAL 支持

的所有 数据格式。通常可以通过调用 GDALAllRegister() 函数来注册所有已知

#include "gdal_priv.h"

int main()

{

GDALDataset *poDataset;

GDALAllRegister();

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

if( poDataset == NULL )

int main()

{

GDALDatasetH hDataset;

GDALAllRegister();

hDataset = GDALOpen( pszFilename, GA_ReadOnly );

if( hDataset == NULL )

{

...;

dataset = gdal.Open( filename, GA_ReadOnly )

名）。它的具体解释由相应的驱动程序负责。 它可能是一个 URL，或者是文件

名以后后面带有许多用于控制打开方式的参数。通常建议， 不要在打开文件的

选择对话框中对文件的类型做太多的限制。

如果 GDAL 数据模型一节所描述的，一个 GDALDataset 包含了光栅数据的一系列

的波段信息。 同时它还包含元数据、一个坐标系统、投影类型、光栅的大小以

及其他许多信息。

adfGeoTransform[0] /* 左上角 x */

adfGeoTransform[1] /* 东西方向一个像素对应的距离 */

adfGeoTransform[5] /* 南北方向一个像素对应的距离 */

如果需要输出 dataset 的基本信息，可以这样：

double adfGeoTransform[6];

printf( "Driver: %s/%s\n",

poDataset->GetDriver()->GetDescription(),

poDataset->GetRasterXSize(), poDataset->GetRasterYSize(),

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] );

printf( "Driver: %s/%s\n",

GDALGetDriverShortName( hDriver ),

GDALGetDriverLongName( hDriver ) );

printf( "Size is %dx%dx%d\n",

GDALGetRasterXSize( hDataset ),

GDALGetRasterYSize( hDataset ),

GDALGetRasterCount( hDataset ) );

if( GDALGetProjectionRef( hDataset ) != NULL )

adfGeoTransform[0], adfGeoTransform[3] );

printf( "Pixel Size = (%.6f,%.6f)\n",

adfGeoTransform[1], adfGeoTransform[5] );

}

In Python:

print 'Driver: ', dataset.GetDriver().ShortName,'/', \

print 'Projection is ',dataset.GetProjection()

geotransform = dataset.GetGeoTransform()

if not geotransform is None:

print 'Origin = (',geotransform[0], ',',geotransform[3],')'

print 'Pixel Size = (',geotransform[1], ',',geotransform[5],')'

现在，我们可以通过 GDAL 获取光栅的一个波段。同样每个波段含有元数据、块

大小、 颜色表以前其他一些信息。下面的代码从 dataset 获取一个

GDALRasterBand 对象， 并且显示它的一些信息。

In C++:

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) )

printf( "Band has %d overviews.\n",

poBand->GetOverviewCount() );