C#使用GDAL读取DEM数据

雨点打透心脏的1/2处 2022-08-14 00:49 274阅读 0赞

**读取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.