GDAL:Ogr读写FileGDB文件(C#)

首先需要FileGDBAPI.dll支持(ogr_FileGDB.dll只支持读,不支持写),存放处理可参见https://blog.csdn.net/xzhh19921019/article/details/53419065。上代码样例(用于将一个文件进行坐标转换生成新的文件。shapefile读写类同于fileGDB,差别在于只有一个layer,无需循环,这里略去具体代码。mdb还没有找到好的方法——不想用AE):

private void ConversionToGeodeticWithFourParametersOfLgo() {
            string errMsg = "";
            foreach (string strFileName in lstFileName) {
                //try {
                DataSource originDataSource = Ogr.Open(strFileName, 0);//0 is for reading,1 is for update
                if (originDataSource == null) {
                    errMsg += "不能打开:" + strFileName + "\r\n";
                    continue;
                }
                Layer layer0 = originDataSource.GetLayerByIndex(0);//if it is a shapefile ,then it only has 1 layer.
                if (layer0 == null) {
                    errMsg += "获取文件:" + strFileName + "的第0个图层失败!" + "\r\n";
                    continue;
                }
                SpatialReference originSpatialReference = layer0.GetSpatialRef();
                string strTemp = originSpatialReference.GetAttrValue("SPHEROID", 0);
                if (strTemp == null) {
                    errMsg += strFileName + ":请注意,该文件无空间参考系,请选用参考系设置”工具将界面设定的“源坐标系”参数写入文件!" + "\r\n";
                    continue;
                }
                SpatialReference targetSpatialReference = originSpatialReference.CloneGeogCS();
                string suffix = Path.GetExtension(strFileName);
                string targetFileName;
                //originDataSource.Dispose();
                if (originSpatialReference.IsProjected() == 1) {//projected to geodetic
                    originSpatialReference.Dispose();
                    double e1Square = 2 * 1 / targetEllipsoid.fReciprocal - 1 / targetEllipsoid.fReciprocal * 1 / targetEllipsoid.fReciprocal;// first eccentricity square
                    double b = targetEllipsoid.a * (1 - 1 / targetEllipsoid.fReciprocal);
                    double e2Square = Math.Pow(targetEllipsoid.a / b, 2) - 1;// second eccentricity square
                    double[] gaussCoefficient = LidfSurveyFunctions.GaussCoefficient(1 / targetEllipsoid.fReciprocal);
                    if (suffix.Equals(".shp")) {                        
                       //略去
                    } else if (suffix.Equals(".gdb")) {
                        targetFileName = Path.Combine(targetDirectory, Path.GetFileName(strFileName));
                        if (Directory.Exists(targetFileName)) {
                            FileHelper.DeleteFolder(targetFileName);
                        }
                        if (!Directory.Exists(targetDirectory)) {
                            Directory.CreateDirectory(targetDirectory);
                        }
                        OSGeo.OGR.Driver driver = Ogr.GetDriverByName("FileGDB");
                        if (driver == null) {
                            errMsg += "不能获取:" + strFileName + "的驱动器,请检查文件!" + "\r\n";
                            continue;
                        }
                        DataSource targetDataSource = driver.CreateDataSource(targetFileName, null);
                        int la
  • 1
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 9
    评论
ESRI Shape文件是一种常用的矢量数据格式,它包含了点、线、面等要素和属性信息。为了处理这种格式的数据,我们可以使用GDAL/OGR库。本文将介绍如何使用C#语言和GDAL/OGRESRI Shape文件。 一、安装GDAL/OGRGDAL/OGR库是一个开源的地理信息系统扩展库,它可以各种GIS数据格式。我们可以从官网(https://gdal.org/download.html)下载最新的二进制版本,并进行安装。 二、创建工作空间 在C#中使用GDAL/OGR库进行ESRI Shape文件,需要先创建一个工作空间。工作空间是一个抽象的概念,它代表了一个包含多个数据集的环境。我们可以使用以下代码创建一个工作空间: ```csharp 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文件中的所有点要素: ```csharp 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文件,并向其中添加一个点要素: ```csharp 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/OGRESRI Shape文件。我们首先创建了一个工作空间,并使用它取了ESRI Shape文件中的点要素。然后,我们创建了一个新的ESRI Shape文件,并向其中添加了一个点要素。使用GDAL/OGR库可以方便地处理各种GIS数据格式,为GIS应用程序的开发提供了便利。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 9
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值