C#用GDAL/OGR库读写ESRI Shape小结[原创]

ESRI Shape文件是一种常用的矢量数据格式，它包含了点、线、面等要素和属性信息。为了处理这种格式的数据，我们可以使用GDAL/OGR库。本文将介绍如何使用C#语言和GDAL/OGR库读写ESRI Shape文件。

一、安装GDAL/OGR库

GDAL/OGR库是一个开源的地理信息系统扩展库，它可以读写各种GIS数据格式。我们可以从官网（https://gdal.org/download.html）下载最新的二进制版本，并进行安装。

二、创建工作空间

在C#中使用GDAL/OGR库进行ESRI Shape文件的读写，需要先创建一个工作空间。工作空间是一个抽象的概念，它代表了一个包含多个数据集的环境。我们可以使用以下代码创建一个工作空间：

using OSGeo.OGR;
using OSGeo.GDAL;

Gdal.AllRegister();
Ogr.RegisterAll();

string workspace = @"C:\data\shapefiles";
Driver driver = Ogr.GetDriverByName("ESRI Shapefile");
DataSource dataSource = driver.Open(workspace, 0);

这里使用了GDAL/OGR库中的`Driver`、`DataSource`和`Gdal`等类。首先，我们调用`Gdal.AllRegister()`方法和`Ogr.RegisterAll()`方法注册GDAL/OGR库中的所有驱动程序。然后，我们指定了一个路径作为工作空间，并使用`Ogr.GetDriverByName()`方法获取ESRI Shapefile驱动程序。最后，我们调用`driver.Open()`方法打开工作空间，并将返回的`DataSource`对象存储到`dataSource`变量中。

三、读取ESRI Shape文件

在创建了工作空间之后，我们可以使用`dataSource`对象读取ESRI Shape文件中的要素和属性信息。以下代码演示了如何读取ESRI Shape文件中的所有点要素：

Layer layer = dataSource.GetLayerByName("points");
Feature feature = null;

while ((feature = layer.GetNextFeature()) != null)
{
    Geometry geometry = feature.GetGeometryRef();
    double x = geometry.GetX(0);
    double y = geometry.GetY(0);
    
    Console.WriteLine($"Point ({x}, {y})");
}

这里使用了GDAL/OGR库中的`Layer`、`Feature`和`Geometry`等类。首先，我们使用`dataSource.GetLayerByName()`方法获取名为“points”的图层，并将返回的`Layer`对象存储到`layer`变量中。然后，我们使用`layer.GetNextFeature()`方法循环遍历图层中的所有要素，并将返回的`Feature`对象存储到`feature`变量中。接着，我们使用`feature.GetGeometryRef()`方法获取要素的几何体，并将返回的`Geometry`对象存储到`geometry`变量中。最后，我们使用`geometry.GetX(0)`和`geometry.GetY(0)`方法获取要素的坐标，并将其打印到控制台上。

四、写入ESRI Shape文件

在读取了ESRI Shape文件之后，我们可以使用`dataSource`对象写入新的要素和属性信息。以下代码演示了如何创建一个新的ESRI Shape文件，并向其中添加一个点要素：

Layer layer = dataSource.CreateLayer("new_points", null, wkbGeometryType.wkbPoint, null);
FieldDefn fieldDefn = new FieldDefn("name", FieldType.OFTString);
layer.CreateField(fieldDefn, 1);

Feature feature = new Feature(layer.GetLayerDefn());
Geometry geometry = new Point(121.5, 31.2);
feature.SetGeometry(geometry);
feature.SetField("name", "Shanghai");

layer.CreateFeature(feature);
dataSource.SyncToDisk();

这里使用了GDAL/OGR库中的`FieldDefn`、`Point`和`wkbGeometryType`等类。首先，我们使用`dataSource.CreateLayer()`方法创建一个名为“new_points”的图层，并将返回的`Layer`对象存储到`layer`变量中。然后，我们使用`new FieldDefn("name", FieldType.OFTString)`语句创建一个名为“name”的字符串型属性。接着，我们使用`layer.CreateField()`方法创建一个新的属性字段，并将`fieldDefn`变量和`1`参数传递给该方法。接下来，我们使用`new Point(121.5, 31.2)`语句创建一个新的点几何体，并将其存储到`geometry`变量中。然后，我们使用`new Feature(layer.GetLayerDefn())`语句创建一个新的要素，并将其存储到`feature`变量中。接着，我们使用`feature.SetGeometry(geometry)`方法设置要素的几何体。最后，我们使用`feature.SetField("name", "Shanghai")`方法设置要素的属性值，并使用`layer.CreateFeature(feature)`方法将要素添加到图层中。最后，我们使用`dataSource.SyncToDisk()`方法将修改保存到磁盘中。

总结

本文介绍了如何使用C#语言和GDAL/OGR库读写ESRI Shape文件。我们首先创建了一个工作空间，并使用它读取了ESRI Shape文件中的点要素。然后，我们创建了一个新的ESRI Shape文件，并向其中添加了一个点要素。使用GDAL/OGR库可以方便地处理各种GIS数据格式，为GIS应用程序的开发提供了便利。