基于ArcGIS10.0和Oracle10g的空间数据管理平台六（C#开发）-空间数据操作类

本系统的一个重点功能就是管理空间数据，所以很多空间数据的操作，而且这些操作可能随时都可能用到，所以单独封装为一个独立的空间数据操作类（地图操作类）。封装的函数都是一些最基本的操作功能，具体实现基于ArcGIS平台的空间数据管理还需要很多复杂的功能和逻辑，后面介绍空间数据管理和编辑的时候会详细介绍。

1.首先我们需要介绍的是连接空间数据库的代码函数，具体实现如下：

/// <summary> /// 连接空间数据库并得到一个SDE工作空间的实例 /// </summary> /// <returns></returns> public static IFeatureWorkspace GetFeatrueWorkspace() { //定义一个属性 IPropertySet pPS = new PropertySetClass(); //设置数据库服务器名，服务器所在的IP地址 pPS.SetProperty("SERVER", ConfigurationSettings.AppSettings["SERVER"]); //propset.SetProperty("SERVER", "10.51.6.233"); //设置SDE的端口，这是安装时指定的，默认安装时"port:5151" pPS.SetProperty("INSTANCE", ConfigurationSettings.AppSettings["INSTANCE"]); //SDE的用户名 pPS.SetProperty("USER", ConfigurationSettings.AppSettings["USER"]); //密码 pPS.SetProperty("PASSWORD", ConfigurationSettings.AppSettings["PASSWORD"]); //设置数据库的名字,只有SQL Server Informix 数据库才需要设置 pPS.SetProperty("DATABASE", ConfigurationSettings.AppSettings["DATABASE"]); //SDE的版本,在这为默认版本 pPS.SetProperty("VERSION", ConfigurationSettings.AppSettings["VERSION"]); //定义一个工作空间,并实例化为SDE的工作空间 IWorkspaceFactory pWF = new SdeWorkspaceFactoryClass(); //打开SDE工作空间,并转化为地物工作空间 IFeatureWorkspace pFW = (IFeatureWorkspace)pWF.Open(pPS, 0); return pFW; }

上面这个函数是采用SDE连接方式，也可以采用直接连接的方式，下面一个函数通过一个if分支语句来判断可以采用SDE也可以采用直接连接的方式实现，具体请看代码及详细的注释：

public static IWorkspace FindWsByDefault() { //定义一个属性 IPropertySet Propset = new PropertySetClass(); if (true) // 采用SDE连接 { //设置数据库服务器名，服务器所在的IP地址 Propset.SetProperty("SERVER", "192.168.1.111"); //设置SDE的端口，这是安装时指定的，默认安装时"port:5151" Propset.SetProperty("INSTANCE", "port:5151"); //SDE的用户名 Propset.SetProperty("USER", "sde"); //密码 Propset.SetProperty("PASSWORD", "sde"); //设置数据库的名字,只有SQL Server Informix 数据库才需要设置 //Propset.SetProperty("DATABASE", "sde"); //SDE的版本,在这为默认版本 Propset.SetProperty("VERSION", "SDE.DEFAULT"); } else // 直接连接 { //设置数据库服务器名,如果是本机可以用"sde:sqlserver:." Propset.SetProperty("INSTANCE", "sde:sqlserver:zhpzh"); //SDE的用户名 Propset.SetProperty("USER", "sa"); //密码 Propset.SetProperty("PASSWORD", "sa"); //设置数据库的名字,只有SQL Server Informix 数据库才需要设置 Propset.SetProperty("DATABASE", "sde"); //SDE的版本,在这为默认版本 Propset.SetProperty("VERSION", "SDE.DEFAULT"); } //定义一个工作空间,并实例化为SDE的工作空间 IWorkspaceFactory factory = new SdeWorkspaceFactoryClass(); //打开SDE工作空间,并转化为地物工作空间 IWorkspace workspace = factory.Open(Propset, 0); return workspace; }

代码注释里面提到的工作空间就相当于一个数据库，里面包含了各种可以操作的图层（相当于数据库中的表）。

2.在一个指定的工作空间中查询一个要素类

/// <summary> /// 查找指定要素 /// </summary> /// <param name="ws"></param> /// <param name="className"></param> /// <param name="dsName"></param> /// <returns></returns> public static IFeatureClass FindClassByName(IWorkspace ws, string className, string dsName) { IEnumDataset enumDs; if (dsName != "") { enumDs = ws.get_Datasets(esriDatasetType.esriDTFeatureDataset); IFeatureDataset featureDs = enumDs.Next() as IFeatureDataset; while (featureDs != null) { string strTemp = featureDs.Name.Substring(featureDs.Name.IndexOf('.') + 1); if (strTemp == dsName) { return GetFcFromDataset(featureDs, className); } featureDs = enumDs.Next() as IFeatureDataset; } } else { enumDs = ws.get_Datasets(esriDatasetType.esriDTFeatureClass); return GetFcFromEnumDataset(enumDs,className); } return null; }

然后利用这个函数实现了另外一个函数：通过要素类名和数据集名在指定的工作空间中寻找要素类

/// <summary> /// 通过要素类名和数据集名在指定的工作空间中寻找要素类 /// </summary> /// <param name="className"></param> /// <param name="datasetName"></param> /// <returns></returns> public static IFeatureClass FindClassByName(string className, string datasetName) { IWorkspace ws = FindWsByDefault(); IFeatureClass featClass = FindClassByName(ws, className, datasetName); return featClass; }

3.在数据集中查找要素类，一个数据集可能包含多个要素类，数据集就是把一些具有共同特征（共同的数据结构类型，空间参考相同等）的数据（要素类）划分成一个一个的集合。

/// <summary> /// 在数据集中查找要素类 /// </summary> /// <param name="pFD"></param> /// <param name="className"></param> /// <returns></returns> public static IFeatureClass GetFcFromDataset(IFeatureDataset pFD, string className) { IFeatureClass pFC; IFeatureClassContainer pFCC = pFD as IFeatureClassContainer; for (int i = 0; i < pFCC.ClassCount; i++) { pFC = pFCC.get_Class(i); string strTemp = pFC.AliasName.Substring(pFC.AliasName.IndexOf('.') + 1); if (strTemp == className) { return pFC; } } return null; }

4.将一个要素类从一个工作空间转移到另外一个工作空间：例如文件系统目录（Shape格式文件组成）、mdb或其他格式，文件系统目录下的所有Shape文件（一个等同于表结构的数据文件，内容等同于表里面的数据）构成一个工作空间。

/// <summary> /// 将一个要素类从一个工作空间转移到另外一个工作空间 /// 注意目标工作空间不能有改要素类，必须先清除 /// </summary> /// <param name="sourceWorkspace">源工作空间</param> /// <param name="targetWorkspace">目标工作空间</param> /// <param name="nameOfSourceFeatureClass">源要素类名</param> /// <param name="nameOfTargetFeatureClass">目标要素类名</param> public static void ConvertFeatureClass(IWorkspace sourceWorkspace, IWorkspace targetWorkspace, string nameOfSourceFeatureClass, string nameOfTargetFeatureClass, int sc_id) { //create source workspace name //创建源工作空间名称 IDataset sourceWorkspaceDataset = (IDataset)sourceWorkspace; IWorkspaceName sourceWorkspaceName = (IWorkspaceName)sourceWorkspaceDataset.FullName; //create source dataset name //创建源数据集名称 IFeatureClassName sourceFeatureClassName = new FeatureClassNameClass(); IDatasetName sourceDatasetName = (IDatasetName)sourceFeatureClassName; sourceDatasetName.WorkspaceName = sourceWorkspaceName; sourceDatasetName.Name = nameOfSourceFeatureClass; //create target workspace name //创建目标工作空间名称 IDataset targetWorkspaceDataset = (IDataset)targetWorkspace; IWorkspaceName targetWorkspaceName = (IWorkspaceName)targetWorkspaceDataset.FullName; //create target dataset name //创建目标数据集名称 IFeatureClassName targetFeatureClassName = new FeatureClassNameClass(); IDatasetName targetDatasetName = (IDatasetName)targetFeatureClassName; targetDatasetName.WorkspaceName = targetWorkspaceName; targetDatasetName.Name = nameOfTargetFeatureClass; //Open input Featureclass to get field definitions. //打开输入的要素类以得到字段定义 IName sourceName = (IName)sourceFeatureClassName; //打开源要素类 IFeatureClass sourceFeatureClass = (IFeatureClass)sourceName.Open(); //Validate the field names because you are converting between different workspace types. //验证字段名称，因为你正在不同类型的工作空间之间进行数据转换 IFieldChecker fieldChecker = new FieldCheckerClass(); IFields targetFeatureClassFields; IFields sourceFeatureClassFields = sourceFeatureClass.Fields; IEnumFieldError enumFieldError; // Most importantly set the input and validate workspaces! //最重要的设置输入和验证工作空间 fieldChecker.InputWorkspace = sourceWorkspace; fieldChecker.ValidateWorkspace = targetWorkspace; fieldChecker.Validate(sourceFeatureClassFields, out enumFieldError, out targetFeatureClassFields); // Loop through the output fields to find the geomerty field //遍历所有输出字段找到几何字段 IField geometryField; for (int i = 0; i < targetFeatureClassFields.FieldCount; i++) { if (targetFeatureClassFields.get_Field(i).Type == esriFieldType.esriFieldTypeGeometry) { geometryField = targetFeatureClassFields.get_Field(i); // Get the geometry field's geometry defenition //得到几何字段的几何定义 IGeometryDef geometryDef = geometryField.GeometryDef; //Give the geometry definition a spatial index grid count and grid size //赋予几何定义一个空间索引格网数目和格网大小值 IGeometryDefEdit targetFCGeoDefEdit = (IGeometryDefEdit)geometryDef; targetFCGeoDefEdit.GridCount_2 = 1; targetFCGeoDefEdit.set_GridSize(0, 0); //Allow ArcGIS to determine a valid grid size for the data loaded //允许ArcGIS为数据加载确定一个有效的格网大小 if (sc_id == 0) { targetFCGeoDefEdit.SpatialReference_2 = geometryField.GeometryDef.SpatialReference; } else { ISpatialReferenceFactory2 pSRF2 = new SpatialReferenceEnvironmentClass(); IGeographicCoordinateSystem pGCS = pSRF2.CreateGeographicCoordinateSystem(sc_id); ISpatialReference pSR = (ISpatialReference)pGCS; targetFCGeoDefEdit.SpatialReference_2 = pSR; } //geometryField.GeometryDef.SpatialReference; // we want to convert all of the features //转换要素类中所有的要素 IQueryFilter queryFilter = new QueryFilterClass(); queryFilter.WhereClause = ""; // Load the feature class //加载要素类 IFeatureDataConverter fctofc = new FeatureDataConverterClass(); IEnumInvalidObject enumErrors = fctofc.ConvertFeatureClass(sourceFeatureClassName, queryFilter, null, targetFeatureClassName, geometryDef, targetFeatureClassFields, "", 1000, 0); break; } } }

5.将mdb或Shape文件中要素类转换追加到sde数据库已有的要素类中，追加的mdb和Shape格式必须同SDE（空间数据库插件）数据库中表结构的数据格式相同

/// <summary> /// 将mdb或Shape文件中要素类转换追加到sde数据库已有的要素类中 /// </summary> /// <param name="sourceWorkspace"></param> /// <param name="targetWorkspace"></param> /// <param name="nameOfSourceFeatureClass"></param> /// <param name="nameOfTargetFeatureClass"></param> public static void AddFeatureClassToSDE(IWorkspace sourceWorkspace, IWorkspace targetWorkspace, string nameOfSourceFeatureClass, string nameOfTargetFeatureClass) { IFeatureWorkspace pSourceFW = sourceWorkspace as IFeatureWorkspace; IFeatureClass pSourceFc = pSourceFW.OpenFeatureClass(nameOfSourceFeatureClass); IFeatureWorkspace pTargetFW = targetWorkspace as IFeatureWorkspace; IFeatureClass pTargetFc = pTargetFW.OpenFeatureClass(nameOfTargetFeatureClass); mdb2fc(pSourceFc, pTargetFc); }

函数通过各自的工作空间和要素类名称等到具体需要转换的要素类（包括源要素类：mdb或shape格式中的表结构和目标要素类：SDE中的表结构），然后调用mdb2fc完成具体的转换工作，实现代码如下：

/// <summary> /// 追加要素到SDE的featureclass中 /// </summary> /// <param name="pSourceFc">源要素类</param> /// <param name="pTargetFc">目的要素类</param> private static void mdb2fc(IFeatureClass pSourceFc, IFeatureClass pTargetFc) { IFeatureCursor pFeaCursor = pSourceFc.Search(null, false); IFeature pFeature = pFeaCursor.NextFeature(); IField pField = new FieldClass(); int iIndex = 0; while (pFeature != null) { IFeature tempFeature = pTargetFc.CreateFeature(); tempFeature.Shape = pFeature.Shape; try { //添加字段值 for (int j = 0; j < pSourceFc.Fields.FieldCount; j++) { pField = pSourceFc.Fields.get_Field(j); iIndex = tempFeature.Fields.FindField(pField.Name); if (pField.Editable && iIndex != -1) { tempFeature.set_Value(iIndex, pFeature.get_Value(j)); } } tempFeature.Store(); } catch (System.Exception ex) { MessageBox.Show(ex.Message); return; } pFeature = pFeaCursor.NextFeature(); } }

上面根据源要素类的每一个字段的每一个值依次遍历并且插入SDE的目的要素类（表）中。

6.在图层中搜索要素

/// <summary> /// 在图层中搜索要素 /// </summary> /// <param name="pFeatLyr">需要搜索的图层</param> /// <param name="pGeo">搜索的集合条件</param> /// <returns>包含结果集的游标对象</returns> public static IFeatureCursor SearchFeatures(IFeatureLayer pFeatLyr, IGeometry pGeo) { try { ISpatialFilter spatialFilter = new SpatialFilterClass(); spatialFilter.Geometry = pGeo; spatialFilter.GeometryField = pFeatLyr.FeatureClass.ShapeFieldName; spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects; IFeatureCursor featureCursor = pFeatLyr.FeatureClass.Search(spatialFilter, false); return pFeatLyr.FeatureClass.Search(spatialFilter, false); } catch (Exception ex) { MessageBox.Show(ex.Message); return null; } }

7.从选择集中创建要素集

/// <summary> /// 从选择集中创建要素集 /// </summary> /// <param name="pFeatureLayer">来源图层</param> /// <param name="pGeometry">选择范围</param> /// <returns>FeatureClass</returns> public static IFeatureClass CreateSelectionLayer(IFeatureLayer pFL, IGeometry pGeometry) { try { IFeatureLayerDefinition pFLD = pFL as IFeatureLayerDefinition; ISpatialFilter pSF = new SpatialFilterClass(); pSF.Geometry = pGeometry; pSF.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects; IQueryFilter pQFilter = pSF; IFeatureSelection pFeatureSel = pFL as IFeatureSelection; pFeatureSel.SelectFeatures(pSF, esriSelectionResultEnum.esriSelectionResultNew, false); IFeatureLayer pSelFL = pFLD.CreateSelectionLayer("Selected Layer", true, null, null); return pSelFL.FeatureClass; } catch (Exception e) { throw (new Exception(e.Message)); } }

8.得到一个要素类和数据集的空间参考系

/// <summary> /// 得到一个要素类的空间参考系 /// </summary> /// <param name="pFeatureClass">要素类</param> /// <returns>空间参考系</returns> public static ISpatialReference getSpatialReference(IFeatureClass pFeatureClass) { IGeoDataset pGeoDataset = pFeatureClass as IGeoDataset; return pGeoDataset.SpatialReference; } /// <summary> /// 得到一个数据集的空间参考系 /// </summary> /// <param name="pFeatureDataset">数据集</param> /// <returns>空间参考系</returns> public static ISpatialReference getSpatialReference(IFeatureDataset pFeatureDataset) { IGeoDataset pGeoDataset = pFeatureDataset as IGeoDataset; return pGeoDataset.SpatialReference; }

9.改变要素类的空间参考

/// <summary> /// 改变要素类的空间参考 /// </summary> /// <param name="pFeatureClass">要素类</param> /// <param name="pGeoType">空间参考类型</param> public static void ChangeFeatureClassRef(IFeatureClass pFeatureClass, int gcsType) { IGeoDataset pGeoDataset = pFeatureClass as IGeoDataset; //QI 到IGeoDatasetSchemaEdit IGeoDatasetSchemaEdit pGeoDatasetSchemaEdit = pGeoDataset as IGeoDatasetSchemaEdit; if (pGeoDatasetSchemaEdit.CanAlterSpatialReference == true) { //创建SpatialReferenceEnvironmentClass对象 ISpatialReferenceFactory2 pSpaRefFactory = new SpatialReferenceEnvironmentClass(); //创建地理坐标系对象 IGeographicCoordinateSystem pNewGeoSys = pSpaRefFactory.CreateGeographicCoordinateSystem(gcsType); pGeoDatasetSchemaEdit.AlterSpatialReference(pNewGeoSys); } }

10.根据各种信息创建要素类

/// <summary> /// 创建一个要素类 /// </summary> /// <param name="workspace">创建要素类的工作空间</param> /// <param name="featureDataset">数据集</param> /// <param name="featureClassName">要素类名称</param> /// <param name="fields">字段集合</param> /// <param name="CLSID"></param> /// <param name="CLSEXT"></param> /// <param name="strConfigKeyword"></param> /// <returns></returns> public static IFeatureClass CreateFeatureClass(IWorkspace2 workspace, IFeatureDataset featureDataset, String featureClassName, IFields fields, UID CLSID, UID CLSEXT, String strConfigKeyword) { // name was not passed in if (featureClassName == "") { return null; } IFeatureClass featureClass; IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspace; //feature class with that name already exists if (workspace.get_NameExists(esriDatasetType.esriDTFeatureClass, featureClassName)) { featureClass = featureWorkspace.OpenFeatureClass(featureClassName); return featureClass; } // assign the class id value if not assigned if (CLSID == null) { CLSID = new UIDClass(); CLSID.Value = "esriGeoDatabase.Feature"; } IObjectClassDescription objectClassDescription = new FeatureClassDescriptionClass(); // if a fields collection is not passed in then supply our own if (fields == null) { return null; } System.String strShapeField = ""; // locate the shape field for (int j = 0; j < fields.FieldCount; j++) { if (fields.get_Field(j).Type == esriFieldType.esriFieldTypeGeometry) { strShapeField = fields.get_Field(j).Name; break; } } // Use IFieldChecker to create a validated fields collection. IFieldChecker fieldChecker = new FieldCheckerClass(); IEnumFieldError enumFieldError = null; IFields validatedFields = null; fieldChecker.ValidateWorkspace = (IWorkspace)workspace; fieldChecker.Validate(fields, out enumFieldError, out validatedFields); // The enumFieldError enumerator can be inspected at this point to determine // which fields were modified during validation. // finally create and return the feature class if (featureDataset == null)// if no feature dataset passed in, create at the workspace level { featureClass = featureWorkspace.CreateFeatureClass(featureClassName, validatedFields, CLSID, CLSEXT, esriFeatureType.esriFTSimple, strShapeField, strConfigKeyword); } else { featureClass = featureDataset.CreateFeatureClass(featureClassName, validatedFields, CLSID, CLSEXT, esriFeatureType.esriFTSimple, strShapeField, strConfigKeyword); } return featureClass; }

11.导入一个要素类到目标工作空间中

/// <summary> /// 导入一个要素类到目标工作空间中 /// </summary> /// <param name="pFW">目标工作空间</param> /// <param name="pFD">目标工作空间中的数据集</param> /// <param name="strFCName">导入的目标工作空间的要素类名称</param> /// <param name="pFC">被导入的要素类</param> public static void ImportFeatureClassToSDE(IFeatureWorkspace pFW, IFeatureDataset pFD, string strFCName, IFeatureClass pSourceFC) { if (pSourceFC == null) { MessageBox.Show("被导入的要素类为空了!"); return; } IWorkspace2 pW2 = pFW as IWorkspace2; IFeatureClass pTargetFC = null; pTargetFC = CreateFeatureClass(pW2, pFD, strFCName, pSourceFC.Fields, null, null, ""); if (pTargetFC != null) { mdb2fc(pSourceFC, pTargetFC); } else { MessageBox.Show("打开或创建目标要素类失败!"); } }

12.编辑空间参照和设置IFeatureDataset空间参照

/// <summary> /// 编辑空间参照 /// </summary> /// <param name="inputSpr"></param> /// <returns></returns> public static ISpatialReference EditSpatialReferenceDlg(ISpatialReference inputSpr) { ESRI.ArcGIS.CatalogUI.ISpatialReferenceDialog2 spatialReferenceDialog2 = new ESRI.ArcGIS.CatalogUI.SpatialReferenceDialogClass(); ISpatialReference spf = spatialReferenceDialog2.DoModalEdit(inputSpr, false, false, false, false, false, false, false, 0); return spf; } /// <summary> /// 设置IFeatureDataset空间参照 /// </summary> /// <param name="ds"></param> /// <returns></returns> public static void SetFdsSpatialReference(IFeatureClass ds) { IGeoDataset gds = ds as IGeoDataset; // Cast the geodataset to the IGeoDatasetSchemaEdit interface and set the Spatial Reference. IGeoDatasetSchemaEdit geoDatasetSchemaEdit = (IGeoDatasetSchemaEdit)gds; if (geoDatasetSchemaEdit.CanAlterSpatialReference) { geoDatasetSchemaEdit.AlterSpatialReference(EditSpatialReferenceDlg(gds.SpatialReference)); } }

这里实现所有的函数都是在类MapOperation中，要看懂上面的所有代码必须要有一些ArcObejct编程的基础知识，当然还包括我前一篇博文介绍的有关于ArcGIS的相关知识，为了给大家提供一些ArcObject的编程知识，特别是用C#讲解的，我上传一本电子书，这本书我觉得很不错，讲到了很多关于ArcGIS二次开发的各个方面。电子书的下载地址：

http://download.csdn.net/detail/wanweiaiaqiang/3824454；或在下载页搜索：ArcGIS Engine+C#实例开发教程。

凌晨1:50，睡觉了。昨天晚上我们整个项目组在公司加班到晚上10点以后哦！