java gdal dataset.readraster_Gdal浅析（一）读取栅格数据到Dataset

gdal读取栅格数据的流程分析

今天，利用gdal写了一个demo用来读取栅格数据、并显示其波段、色调等信息。

在此分享一下，我对gdal读取栅格数据的流程的浅析。

1、使用gdal的第一步肯定是先要注册。使用静态函数gdal.AllRegister();

2、gdal将所有的栅格数据都是读取到一个它自己定义的dataset中，所以OSGeo.GDAL.Dataset ds=new OSGeo.GDAL.Dataset();

3、使用open函数 ds = gdal.Open(inputFileName, gdalconst.GA_ReadOnly);

4、至此，全路径为inputFileName的栅格文件已经读到ds中，下面介绍OSGeo.GDAL.Dataset类

总体来说，OSGeo.GDAL.Dataset实际上就是一系列相关联的波段的组合。

其共有成员如下：

virtual  ~GDALDataset ()

int  GetRasterXSize (void)

int  GetRasterYSize (void)

int  GetRasterCount (void)

GDALRasterBand *  GetRasterBand (int)

virtual void  FlushCache (void)

virtual const char *  GetProjectionRef (void)

virtual CPLErr  SetProjection (const char *)

virtual CPLErr  GetGeoTransform (double *)

virtual CPLErr  SetGeoTransform (double *)

virtual CPLErr  AddBand (GDALDataType eType, char **papszOptions=NULL)

virtual void *  GetInternalHandle (const char *)

virtual GDALDriver *  GetDriver (void)

virtual int  GetGCPCount ()

virtual const char *  GetGCPProjection ()

virtual const GDAL_GCP *  GetGCPs ()

virtual CPLErr  SetGCPs (int nGCPCount, const GDAL_GCP *pasGCPList, const char *pszGCPProjection)

virtual CPLErr  AdviseRead (int nXOff, int nYOff, int nXSize, int nYSize, int nBufXSize, int nBufYSize, GDALDataType eDT, int nBandCount, int *panBandList, char **papszOptions)

CPLErr  RasterIO (GDALRWFlag, int, int, int, int, void *, int, int, GDALDataType, int, int *, int, int, int)

int  Reference ()

int  Dereference ()

GDALAccess  GetAccess ()

int  GetShared ()

void  MarkAsShared ()

CPLErr  BuildOverviews (const char *, int, int *, int, int *, GDALProgressFunc, void *)

上述代码是C语言的，有比较烦躁的指针。不过只要大家仔细思考，还是可以很轻松的将其转换成C#语言的。

5、其中我们要使用的、或者说比较常用的方法有以下几个。

获取驱动，例如对于tif文件的驱动就是GeoTiff。

virtual GDALDriver *  GetDriver (void)

获取投影方式，在.NET中，返回值是string类型的投影描述语句

virtual const char *  GetProjectionRef (void)

获取波段数目，这个值决定了下面的分别获取每一个波段的循环次数

int  GetRasterCount (void)

for (int iBand = 1; iBand <= ds.RasterCount; iBand++)

{

Band band = ds.GetRasterBand(iBand);//获取制定波段

}

int  GetRasterXSize (void)

int  GetRasterYSize (void)

获取波段的size

最后一个：

CPLErr  BuildOverviews (const char *, int, int *, int, int *, GDALProgressFunc, void *)

这个是在ds上构建Overviews，我猜想是建类似金字塔的抽取，或称之为视图。目前还在摸索，希望高手指教！

6、这次就先研究到Band，下次从Band入手，继续摸索。

Demo的效果图如下：

--