C#使用GDAL读取DEM数据

读取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.       ReadRaster函数

声明：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 – 待读取数据x方向的起始位置，即x方向的偏移值

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

xSize – 待读取数据的宽度

ySize – 待读取数据的高度

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

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

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

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

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