C# GDAL 读取DEM栅格数据(Dataset.ReadRaster())

准备文件

编译好的gdal核心库gdal180.dll以及C#封装库gdal_wrap.dll、gdal_csharp.dll

 

 读取DEM数据步骤及说明

1.      将gdal180.dll、gdal_wrap.dll、 gdal_csharp.dll拷贝到程序的生成目录，并在项目里添加对gdal_csharp.dll库的引用。

2.      在要使用gdal的文件头部加上如下命名空间的声明：

using OSGeo.GDAL;

3.      读取代码如下：

        

Gdal.AllRegister();        // 初始化gdal库

           Dataset ds = Gdal.Open(path, Access.GA_ReadOnly);     // 读取数据到gdal的dataset中

           int XSize = ds.RasterXSize;        

           int YSize = ds.RasterYSize;        // 获取栅格数据的长和宽

           int count = ds.RasterCount;        // 获取栅格数据的点的数量

           Band demband = ds.GetRasterBand(1); // 获取第一个band         

           double[] gt = new double[6];          

           ds.GetGeoTransform(gt);            // 获取屏幕坐标转换到实际地理坐标的参数

           int Xpixel = 100;

           int Yline = 100;

           double  Xgeo = gt[0] + Xpixel*gt[1] + Yline*gt[2];

           double  Ygeo = gt[3] + Xpixel*gt[4] + Yline*gt[5];  // 计算像素坐标为(100, 100)的点的实际地理坐标

           double nodatavalue;

           int hasval;

           demband.GetNoDataValue(out nodatavalue, out hasval);   // 获取没有数据的点的值

           double[] databuf = new double[XSize * YSize]; 

           demband.ReadRaster(0, 0, ds.RasterXSize, ds.RasterYSize, databuf, ds.RasterXSize, ds.RasterYSize, 0, 0); // 读取数据到缓冲区中

4.      主要函数说明。

1． Open函数

声明：

public static Dataset Open(string utf8_path, Access eAccess);

        该函数用于从指定路径读取栅格数据。该函数有两个参数，第一个参数是栅格数据的路径，可以为本地或网络路径。第二个参数是打开数据的模式，GA_ReadOnly表示以只读方式打开，GA_Update表示以读写方式打开。该函数返回一个DataSet对象，gdal用dataset对象为所有栅格数据提供统一的访问接口。Dataset中包含1个或多个band，用于存放实际数据。Dataset还包含了栅格数据的大小、nodatavalue等信息。

2.      GetRasterBand函数

声明：

public Band GetRasterBand(int nBand);

该函数用于从dataset中取得一个band。该函数有一个参数，用于指定要取的band的序号。该函数返回一个Band对象，Band是实际存储栅格数据的对象。栅格数据的一个点上可以存储多个值，一个栅格数据可以有1个或多个Band。一个Band可以理解为栅格数据的一组值。比如说DEM数据的每个点只有一个高程值，因此只有一个Band，而一幅bmp图像每个像素点有RGB三种颜色的值，因而有3个Band，使用GDAL可以方便的单独操作每个Band。

3.      GetGeoTransform函数

声明：

public void GetGeoTransform(double[] argout);

该函数用于取得从像素坐标转换到地理坐标的6个参数。该函数有一个数组参数， 用于输出转换系数。实际应用中，首先调用该函数得到6个系统，然后按如下式子进行转换：

       Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2)

      Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5)

其中Xgeo, Ygeo是地理坐标， Xpixel, Yline是像素坐标。GT是系数数组。如果影像正上为北方向，则GT2和GT4参数为0，而GT1是象元宽，GT5是象元高， （GT0，GT3）点位置是影像的左上角。 注意，上面所说的点/线坐标系是从左上角(0,0)点到右下角，也就是坐标轴从 左到右增长，从上到下增长的坐标系（即影象的行列从左下角开始计算）。

4.      GetGeoTransform函数

声明：

public CPLErr ReadRaster(int xOff, int yOff, int xSize, int ySize, double[] buffer, int buf_xSize, int buf_ySize, int pixelSpace, int lineSpace);

该函数用于从Band中读取栅格数据。该函数有9个参数，可以对如何读取数据进行精确描述。该函数返回一个CPLErr类型，用于指示读取数据是否成功。GDAL C#封装库提供了多个重载用于读取不同类型的数据，即buffer参数的类型可以是double, byte, int, float, short或者一个指针。

 

xOff和yOff是指偏移量，即从影像的左上角起始坐标（xOff,yOff）开始读取数据。

xSize和ySize是指读取图像数据的行列数，即宽度和高度，单位都是像素。

Buffer是图像数据缓存。

buf_xSize和buf_ySize是缓存区的大小，它们须与buffer申请的大小保持一致，通过这两个参数可以控制缩放，如果它们小于xSize和ySize就是将原图缩小，反之如果它们大于xSize和ySize就是将原图放大。

pixelSpace和lineSpace一般默认取0即可。

参数说明：

xOff – 待读取数据x方向的起始位置，即x方向的偏移值

yOff – 待读取数据y方向的起始位置，即y方向的偏移值

xSize – 待读取数据的宽度

ySize – 待读取数据的高度

buffer – 缓冲区用于存放读取出来的数据

buf_xSize – 缓冲区x方向的大小，一般等于xSize

buf_ySize – 缓冲区y方向的大小，一般等于ySize

pixelSpace – 两个像素横向之间的间隔大小，一般取默认值0，即buffer的类型的大小。

lineSpace – 两行之间第一个数据之间的间隔大小，一般取默认值0，即buffer类型大小*buf_xSize.