C# Gdal读取和处理(字段,获取Feature、Geometry)SHP文件

31 篇文章 8 订阅

using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Windows.Forms;
//
using System.IO;
using OSGeo.GDAL;
using OSGeo.OGR;
using OSGeo.OSR;
using System.Collections;
 
namespace GdalReadSHP
{
    /// <summary>
    /// 定义SHP解析类
    /// </summary>
    public class ShpRead
    {
        /// 保存SHP属性字段
        public OSGeo.OGR.Driver oDerive;
        public List<string> m_FeildList;
        private Layer oLayer;
        public string sCoordiantes;
        public ShpRead()
        {
            m_FeildList = new List<string>();
            oLayer = null;
            sCoordiantes = null;
        }
 
        /// <summary>
        /// 初始化Gdal
        /// </summary>
        public void InitinalGdal()
        {
            // 为了支持中文路径
            Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES");
            // 为了使属性表字段支持中文
            Gdal.SetConfigOption("SHAPE_ENCODING", "");
            Gdal.AllRegister();
            Ogr.RegisterAll();
            
            oDerive = Ogr.GetDriverByName("ESRI Shapefile");
            if (oDerive == null)
            {
                MessageBox.Show("文件不能打开,请检查");
            }
        }
 
        /// <summary>
        /// 获取SHP文件的层
        /// </summary>
        /// <param name="sfilename"></param>
        /// <param name="oLayer"></param>
        /// <returns></returns>
        public bool GetShpLayer(string sfilename)
        {
            if (null == sfilename || sfilename.Length <= 3)
            {
                oLayer = null;
                return false;
            }
            if (oDerive == null)
            {
                MessageBox.Show("文件不能打开,请检查");
            }
            DataSource ds = oDerive.Open(sfilename, 1);
            if (null == ds)
            {
                oLayer = null;
                return false;
            }
            int iPosition = sfilename.LastIndexOf("\\");
            string sTempName = sfilename.Substring(iPosition + 1, sfilename.Length - iPosition - 4 - 1);
            oLayer = ds.GetLayerByName(sTempName);
            if (oLayer == null)
            {
                ds.Dispose();
                return false;
            }
            return true;
        }
        /// <summary>
        /// 获取所有的属性字段
        /// </summary>
        /// <returns></returns>
        public bool GetFeilds()
        {
            if (null == oLayer)
            {
                return false;
            }
            m_FeildList.Clear();
            wkbGeometryType oTempGeometryType = oLayer.GetGeomType();
            List<string> TempstringList = new List<string>();
 
            //
            FeatureDefn oDefn = oLayer.GetLayerDefn();
            int iFieldCount = oDefn.GetFieldCount();
            for (int iAttr = 0; iAttr < iFieldCount; iAttr++)
            {
                FieldDefn oField = oDefn.GetFieldDefn(iAttr);
                if (null != oField)
                {
                    m_FeildList.Add(oField.GetNameRef());
                }
            }
            return true;
        }
        /// <summary>
        ///  获取某条数据的字段内容
        /// </summary>
        /// <param name="iIndex"></param>
        /// <param name="FeildStringList"></param>
        /// <returns></returns>
        public bool GetFeildContent(int iIndex, out List<string> FeildStringList)
        {
            FeildStringList = new List<string>();
            Feature oFeature = null;
            if ((oFeature = oLayer.GetFeature(iIndex)) != null)
            {
                
                FeatureDefn oDefn = oLayer.GetLayerDefn();
                int iFieldCount = oDefn.GetFieldCount();
                // 查找字段属性
                for (int iAttr = 0; iAttr < iFieldCount; iAttr++)
                {
                    FieldDefn oField = oDefn.GetFieldDefn(iAttr);
                    string sFeildName = oField.GetNameRef();
 
                    #region 获取属性字段
                    FieldType Ftype = oFeature.GetFieldType(sFeildName);
                    switch (Ftype)
                    {
                        case FieldType.OFTString:
                            string sFValue = oFeature.GetFieldAsString(sFeildName);
                            string sTempType = "string";
                            FeildStringList.Add(sFValue);
                            break;
                        case FieldType.OFTReal:
                            double dFValue = oFeature.GetFieldAsDouble(sFeildName);
                            sTempType = "float";
                            FeildStringList.Add(dFValue.ToString());
                            break;
                        case FieldType.OFTInteger:
                            int iFValue = oFeature.GetFieldAsInteger(sFeildName);
                            sTempType = "int";
                            FeildStringList.Add(iFValue.ToString());
                            break;
                        default:
                            //sFValue = oFeature.GetFieldAsString(ChosenFeildIndex[iFeildIndex]);
                            sTempType = "string";
                            break;
                    }
                    #endregion
                }
            }
            return true;
        }
        /// <summary>
        /// 获取数据
        /// </summary>
        /// <returns></returns>
        public bool GetGeometry(int iIndex)
        {
            if (null == oLayer)
            {
                return false;
            }
            int iFeatureCout = oLayer.GetFeatureCount(0);
            Feature oFeature = null;
            oFeature = oLayer.GetFeature(iIndex);
            //  Geometry
            Geometry oGeometry = oFeature.GetGeometryRef();
            wkbGeometryType oGeometryType = oGeometry.GetGeometryType();
            switch (oGeometryType)
            {
                case wkbGeometryType.wkbPoint:
                    oGeometry.ExportToWkt(out sCoordiantes);
                    sCoordiantes = sCoordiantes.ToUpper().Replace("POINT (", "").Replace(")", "");
                    break;
                case wkbGeometryType.wkbLineString:
                case wkbGeometryType.wkbLinearRing:
                    oGeometry.ExportToWkt(out sCoordiantes);
                    sCoordiantes = sCoordiantes.ToUpper().Replace("LINESTRING (", "").Replace(")", "");
                    break;
                default:
                    break;
            }
            return false;
        }
 
    }//END class
}

要在 C#读取 GeoPackage 文件,可以使用 GDAL 库。GDAL 是一个开源的地理空间数据处理库,它支持多种格式的地理空间数据,包括 GeoPackage。 以下是使用 GDALC#读取 GeoPackage 的示例代码: ```csharp using OSGeo.OGR; using OSGeo.GDAL; // 注册所有的驱动程序 Gdal.AllRegister(); // 打开数据源 DataSource ds = Ogr.Open("path/to/your/geopackage.gpkg", 0); // 获取第一个图层 Layer layer = ds.GetLayerByIndex(0); // 获取图层的元数据 FeatureDefn defn = layer.GetLayerDefn(); // 输出所有字段名 for (int i = 0; i < defn.GetFieldCount(); i++) { FieldDefn fieldDefn = defn.GetFieldDefn(i); Console.WriteLine(fieldDefn.GetName()); } // 逐个读取要素 Feature feature; while ((feature = layer.GetNextFeature()) != null) { // 输出要素的属性值 for (int i = 0; i < defn.GetFieldCount(); i++) { FieldDefn fieldDefn = defn.GetFieldDefn(i); Console.WriteLine(feature.GetFieldAsString(i)); } // 获取要素的几何形状 Geometry geom = feature.GetGeometryRef(); if (geom != null) { // 在这里进行几何形状的处理 } feature.Dispose(); } // 关闭数据源 ds.Dispose(); ``` 以上代码使用 OSGeo.OGR 和 OSGeo.GDAL 命名空间中的类来读取 GeoPackage 文件。首先需要调用 Gdal.AllRegister() 方法来注册所有的驱动程序,然后通过 Ogr.Open() 方法打开 GeoPackage 文件,并通过 GetLayerByIndex() 方法获取第一个图层。接着可以通过 GetLayerDefn() 方法获取图层的元数据,包括字段名和数据类型等信息。使用 GetNextFeature() 方法逐个读取要素,可以通过 GetFieldAsString() 方法获取要素的属性值,通过 GetGeometryRef() 方法获取要素的几何形状。最后需要记得关闭数据源。 需要注意的是,在运行以上代码之前,需要将 GDAL 的 C/C++ 库和 C# 绑定库(gdal_csharp.dll 和 osr_csharp.dll)添加到项目中,并将它们的依赖项(gdalXXX.dll 和 osrXXX.dll)复制到输出目录。可以从 GDAL 的官方网站(http://www.gdal.org/)下载最新版本的 GDAL 库。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值