arcgis面重叠检查.docx

using System;using System.Collections.Generic;using System.Linq;using System.Text;using System.Data;using ESRI.ArcGIS.Geometry;using ESRI.ArcGIS.DataSourcesGDB;using ESRI.ArcGIS.Geodatabase;namespace PolygonCoveredPolygon public class PolygonCoveredPolygon /* 功能函数：面重叠检查 /* 参 数：源图层、目标图层 /* 返 回 值：检查结果（源图层名、源要素OID、目标图层名、目标要素的OID） public DataTable DataCheck(IFeatureClass orgFeaCls, IFeatureClass desFeaCls, out Exception outError) DataTable dt = null; outError = null; if (orgFeaCls.ShapeType != esriGeometryType.esriGeometryPolygon) return null; /源图层不是面 if (desFeaCls.ShapeType != esriGeometryType.esriGeometryPolygon) return null; /目标图层不是面 dt = new DataTable(); dt.Columns.Add(源图层名); dt.Columns.Add(源要素OID); dt.Columns.Add(目标图层); dt.Columns.Add(目标要素OID); /查找源要素 IFeatureCursor pOrgCursor = orgFeaCls.Search(null, false); IFeature pOrgFea = pOrgCursor.NextFeature();/源要素 /查找目标要素 IFeatureCursor pDesCursor = desFeaCls.Search(null, false); IFeature pDesFea = pDesCursor.NextFeature(); /目标要素 List OID = new List(); /要保存目标的OID while (pDesFea != null) if(!OID.Contains(pDesFea.OID) OID.Add(pDesFea.OID); pDesFea = pDesCursor.NextFeature(); if(OID=null|OID.Count=0) return null; while (pOrgFea!= null) /源要素不为空 IGeometry pOrgFeaGeo = pOrgFea.Shape; /源要素的几何形状 for (int i = 0; i OID.Count; i+) IFeature mDesFea=desFeaCls.GetFeature(OIDi); if(mDesFea=null) continue; IRelationalOperator pRelation = mDesFea.Shape as IRelationalOperator; if (pRelation.Overlaps(pOrgFeaGeo) /判断 源要素和目标要素 是否重叠 如果面重叠 DataRow dr = dt.NewRow(); dr源图层名 = pOrgFea.Class.AliasName; dr源要素OID = pOrgFea.OID; dr目标图层 = mDesFea.Class.AliasName; dr目标要素OID = mDesFea.OID; dt.Rows.Add(dr); pOrgFea = pOrgCursor.NextFeature(); return dt; using System;using System.Collections.Generic;using System.Linq;using System.Text;using System.Data;using ESRI.ArcGIS.Geometry;using ESRI.ArcGIS.DataSourcesGDB;using ESRI.ArcGIS.Geodatabase; namespace LineCoveredCheck public class LineCoveredCheckCls /* 功能函数：重复性检查 /* 参 数：源图层、目标图层 /* 返 回 值：检查结果（源图层名、源要素OID、目标图层名、目标要素OID） public DataTable DataCheck(IFeatureClass orgFeaCls, IFeatureClass desFeaCls, out Exception outError) DataTable dt = null; outError = null; dt = new DataTable(); dt.Columns.Add(源图层名); dt.Columns.Add(源要素OID); dt.Columns.Add(目标图层名); dt.Columns.Add(目标要素OID); /查找源要素 IFeatureCursor pOrgCursor = orgFeaCls.Search(null, false); IFeature pOrgFea = pOrgCursor.NextFeature(); /源要素 /查找目标要素 IFeatureCursor pDesCursor = desFeaCls.Search(null, false); IFeature pDesFea = pDesCursor.NextFeature(); /源要素 List OID = new List(); /保存目标要素的OID while (pDesFea != null) if (!OID.Contains(pDesFea.OID) OID.Add(pDesFea.OID); pDesFea = pDesCursor.NextFeature(); if (OID = null | OID.Count = 0) return null; while (pOrgFea != null) /源要素不为空 IGeometry pOrgFeaGeo = pOrgFea.Shape; /源要素的几何形状 for (int i = 0; i OID.Count; i+) /循环目标要素 IFeature mDesFea = desFeaCls.GetFeature(OIDi);/ 获取目标要素 if (mDesFea = null) continue; IRelationalOperator pRelation = mDesFea.Shape as IRelationalOperator; /几何关系 if (pRelation.Equals(pOrgFeaGeo) /判断空间位置是否重合 if (pOrgFea.OID = OIDi) continue;/如果是同一个要素就跳出 if (BeSameProperties1(pOrgFea, mDesFea) DataRow dr = dt.NewRow(); dr源图层名 = pOrgFea.Class.AliasName; /显示错误信息 dr源要素OID = pOrgFea.OID; dr目标图层名 = mDesFea.Class.AliasName; dr目标要素OID = mDesFea.OID; dt.Rows.Add(dr); pOrgFea = pOrgCursor.NextFeature(); return dt; private bool BeSameProperties1(IFeature pOrgFea, IFeature mDesFea) bool b = true; string fieldNames = ; /字段字符串返回值 for (int m = 0; m mDesFea.Fields.FieldCount; m+) IField pField = mDesFea.Fields.get_Field(m); /字段 if (pField = null) continue; if (pField.Type = esriFieldType.esriFieldTypeOID) continue;/OID 字段排除 if (pField.Type = esriFieldType.esriFieldTypeGeometry) continue;/Geometry 排除 if (pField.Name.ToUpper().StartsWith(SHAPE) continue; string pfieldName = pField.Name; fieldNames += pfieldName + ,; fieldNames = fieldNames.Substring(0, fieldNames.Length - 1);/所有字段的字符串 string fieldsArr = fieldNames.Split(new char , );/需要检查的字段数组 int orgindex = -1; int desindex = -1; for (int j = 0; j fieldsArr.Length; j+) orgindex = pOrgFea.Fields.FindField(fieldsArrj); desindex = mDesFea.Fields.FindField(fieldsArrj); if (pOrgFea.get_Value(orgindex).ToString() != mDesFea.get_Value(desindex).ToString() return false; /属性不一致; return b; public DataTable DataCheck1(IFeatureClass pFeaCls, out Exception pError) DataTable dt = null; pError = null; dt = new DataTable(); dt.Columns.Add(源图层名); dt.Columns.Add(源要素OID); dt.Columns.Add(目标图层名); dt.Columns.Add(目标要素OID); IFeatureCursor pCursor = pFeaCls.Search(null, false); IFeature pFea = pCursor.NextFeature(); List LstFea = new List(); while (pFea != null) if (!LstFea.Contains(pFea) LstFea.Add(pFea); pFea = pCursor.NextFeature(); foreach (IFeature pOrgFea in LstFea) foreach (IFeature pDesFea in LstFea) if (pOrgFea.OID = pDesFea.OID) continue; /判断空间位置是否重合，若重合，则进行线面的检查 IGeometry pGeo=pOrgFea.Shape; IRelationalOperator pRelation = pDesFea.Shape as IRelationalOperator; /几何关系 if (pRelation.Equals(pGeo) /检查属性是否一致 if (BeSameProperties(pOrgFea, pDesFea) DataRow dr = dt.NewRow(); dr源图层名 = pOrgFea.Class.AliasName; /显示错误信息 dr源要素OID = pOrgFea.OID; dr目标图层名 = pDesFea.Class.AliasName; dr目标要素OID = pDesFea.OID; dt.Rows.Add(dr); return dt; private bool BeSameProperties(IFeature pOrgFea, IFeature pDesFea) bool beSucced = true; for (int i = 0; i pOrgFea.Fields.FieldCount; i+) if (pOrgFea.Fields.get_Field(i).Type = esriFieldType.esriFieldTypeOID) continue; if (pOrgFea.Fields.get_F