private List <int> GetIdsOfContainedFeatureCentroids(IRecordSet recordset, IPolygon tilePolygon)
{
List <int> ids = new List <int>();
IFeatureCursor cursor = recordset.get_Cursor(true) as IFeatureCursor;
IFeature feature = cursor.NextFeature();
if (feature == null)
{
return(ids);
}
if (feature.ShapeCopy.GeometryType != esriGeometryType.esriGeometryMultipoint && feature.ShapeCopy.GeometryType != esriGeometryType.esriGeometryPolygon)
{
throw new ArgumentException("Only multipoint and polygon geometry types are supported by this method.");
}
this.tracePolygon(tilePolygon);
IRelationalOperator2 relationalOperator = tilePolygon as IRelationalOperator2;
try
{
// Envelope must contain centriod
IPoint centroid = null;
while (feature != null)
{
centroid = (feature.ShapeCopy as IArea).Centroid;
if (relationalOperator.Contains(centroid))
{
this.tracePolygonContainsPoint(tilePolygon, centroid, feature.OID);
ids.Add(feature.OID);
}
feature = cursor.NextFeature();
}
if (ids.Count > 0)
{
System.Diagnostics.Debug.WriteLine("OIDs with centroids contained: " + ids.ToDelimitedString <int>(","));
}
return(ids);
}
finally
{
this.ReleaseComObject(cursor);
cursor = null;
}
}
public static void PanToGeometry(IActiveView view, IGeometry geometry, bool setCenterAt = false)
{
IEnvelope env = new EnvelopeClass();
env = view.Extent;
IRelationalOperator2 operation = env as IRelationalOperator2;
if (setCenterAt || !operation.Contains(geometry))
{
ISegmentCollection poly = new PolygonClass();
IArea area = geometry.Envelope as IArea;
env.CenterAt(area.Centroid);
view.Extent = env;
view.Refresh();
view.ScreenDisplay.UpdateWindow();
}
}
private void addToNeighboringFeature(ref HashSet <IPolygon4>[] polys)
{
foreach (IPolygon4 geoS in polys[0])
{
//polys[0].
if (geoS.IsEmpty)
{
continue;
}
else
{
bool bTest = false;
foreach (IPolygon4 geoB in polys[1])
{
if (geoB.IsEmpty)
{
continue;
}
else
{
IRelationalOperator2 ro2 = (IRelationalOperator2)geoB;
if (!ro2.Disjoint(geoS))
{
try
{
ITopologicalOperator5 tp6 = (ITopologicalOperator5)geoB;
tp6.Union(geoS);
bTest = true;
break;
}
catch (Exception e)
{
Console.WriteLine(e.ToString());
}
}
}
}
}
}
}
public Dictionary <string, List <string> > PolygonNeighbors(
IFeatureClass featureClass,
string uid_field_name,
IQueryFilter queryFilter,
bool areaOverlap = false,
bool containsSelf = true)
{
if (featureClass.ShapeType != esriGeometryType.esriGeometryPolygon)
{
return(null);
}
int uid_idx = featureClass.FindField(uid_field_name);
// 创建输出字典
Dictionary <string, List <string> > outDict = new Dictionary <string, List <string> >();
var featureCursor = featureClass.Search(queryFilter, true);
IFeature feature = null;
AddStepProgressor(featureClass.FeatureCount(queryFilter), "Processing Neighbors . . . ");
while ((feature = featureCursor.NextFeature()) != null)
{
if (!ContinueStepProgressor(new object[] { featureCursor }))
{
return(null);
}
IRelationalOperator2 relOperator = (IRelationalOperator2)feature.Shape;
List <string> nbr_list = new List <string>();
var nbFeatureCursor = featureClass.Search(queryFilter, true);
IFeature nbFeature = null;
while ((nbFeature = nbFeatureCursor.NextFeature()) != null)
{
if (nbFeature.OID == feature.OID)
{
continue;
}
if (areaOverlap) // 分析重合的相邻情况
{
if (relOperator.IsNear(nbFeature.Shape, 0) || relOperator.Overlaps(nbFeature.Shape))
{
nbr_list.Add(nbFeature.Value[uid_idx].ToString());
}
}
else
{
if (relOperator.IsNear(nbFeature.Shape, 0))
{
nbr_list.Add(nbFeature.Value[uid_idx].ToString());
}
}
}
Marshal.ReleaseComObject(nbFeatureCursor);
GC.Collect();
GC.WaitForPendingFinalizers();
if (nbr_list.Count > 0) // 插入数据
{
string main_str = feature.Value[uid_idx].ToString();
if (containsSelf)
{
nbr_list.Insert(0, main_str);
}
outDict.Add(main_str, nbr_list);
}
}
Marshal.ReleaseComObject(featureCursor);
GC.Collect();
GC.WaitForPendingFinalizers();
return(outDict);
}
public Dictionary <string, List <string> > PolygonNeighbors(
IFeatureClass featureClass,
string uid_field_name)
{
if (featureClass.ShapeType != esriGeometryType.esriGeometryPolygon)
{
return(null);
}
int uid_idx = featureClass.FindField(uid_field_name);
// 创建输出字典
Dictionary <string, List <string> > outDict = new Dictionary <string, List <string> >();
var featureCursor = featureClass.Search(null, true);
IFeature feature = null;
AddStepProgressor(featureClass.FeatureCount(null), "Processing Neighbors . . . ");
while ((feature = featureCursor.NextFeature()) != null)
{
if (!ContinueStepProgressor(new object[] { featureCursor }))
{
return(null);
}
IRelationalOperator2 relOperator = (IRelationalOperator2)feature.Shape;
List <string> nbr_list = new List <string>();
ISpatialFilter spatialFilter = new SpatialFilterClass();
spatialFilter.Geometry = feature.Shape;
spatialFilter.SpatialRel = esriSpatialRelEnum.esriSpatialRelTouches;
var nbFeatureCursor = featureClass.Search(spatialFilter, true);
IFeature nbFeature = null;
while ((nbFeature = nbFeatureCursor.NextFeature()) != null)
{
if (nbFeature.OID == feature.OID)
{
continue;
}
nbr_list.Add(nbFeature.Value[uid_idx].ToString());
}
Marshal.ReleaseComObject(nbFeatureCursor);
GC.Collect();
GC.WaitForPendingFinalizers();
if (nbr_list.Count > 0) // 插入数据
{
string main_str = feature.Value[uid_idx].ToString();
nbr_list.Insert(0, main_str);
outDict.Add(main_str, nbr_list);
}
}
Marshal.ReleaseComObject(featureCursor);
GC.Collect();
GC.WaitForPendingFinalizers();
return(outDict);
}
public void OtherRule(string idname, string checkname, string IN_RuleType, string TABLENAME, string SUPTABLE, string inputtext, IHookHelper m_hookHelper)
{
CommonClass common = new CommonClass();
IFeatureClass IN_FeatureClass = common.GetLayerByName(TABLENAME).FeatureClass;
IGeoDataset pGeoDataset = IN_FeatureClass as IGeoDataset;
ISpatialReference spatialReference = pGeoDataset.SpatialReference;
IFeatureClass IN_Sup_FeatureClass = null;
if (SUPTABLE != null)
{
IN_Sup_FeatureClass = common.GetLayerByName(SUPTABLE).FeatureClass;
}
string ErrorFilePath = this.topoDir + "\\" + checkname + idname + ".shp";
if (IN_RuleType == "面多部件检查")
{
try
{
common.CreatShpFile(this.topoDir, spatialReference, esriGeometryType.esriGeometryPolygon, checkname + idname);
List <ErrorEntity> list = new List <ErrorEntity>();
IFeatureCursor cursor = IN_FeatureClass.Search(null, false);
int tempCount = 0;
IFeature pFeature = cursor.NextFeature();
while (pFeature != null)
{
IGeometry pGeo = pFeature.ShapeCopy;
ITopologicalOperator pTopoOperator = pGeo as ITopologicalOperator;
int iCount = 0;
if (IN_FeatureClass.ShapeType == esriGeometryType.esriGeometryPolygon)
{
iCount = (pGeo as IPolygon).ExteriorRingCount;
}
else if (IN_FeatureClass.ShapeType == esriGeometryType.esriGeometryPolyline)
{
iCount = ((pGeo as IPolyline) as IGeometryCollection).GeometryCount;
}
else if (IN_FeatureClass.ShapeType == esriGeometryType.esriGeometryMultipoint)
{
iCount = ((pGeo as IMultipoint) as IPointCollection).PointCount;
}
if (iCount > 1)
{
tempCount++;
common.GenerateSHPFile(pFeature.ShapeCopy, ErrorFilePath);
}
pFeature = cursor.NextFeature();
}
//new ErrorTable().AddErr(list, ErrType.MultiPart, idname);
if (!DicTopoError.ContainsKey(idname))
{
DicTopoError.Add(idname, tempCount);
}
else
{
DicTopoError[idname] = tempCount;
}
Marshal.ReleaseComObject(cursor);
if (tempCount == 0)
{
string shpFolder = System.IO.Path.GetDirectoryName(ErrorFilePath);
IWorkspaceFactory pWorkspaceFac = new ShapefileWorkspaceFactoryClass();
IWorkspace pWorkSpace = pWorkspaceFac.OpenFromFile(ErrorFilePath, 0);
IFeatureWorkspace pFeatureWorkSpace = pWorkSpace as IFeatureWorkspace;
if (System.IO.File.Exists(ErrorFilePath))
{
IFeatureClass pFCChecker = pFeatureWorkSpace.OpenFeatureClass(ErrorFilePath);
if (pFCChecker != null)
{
IDataset pds = pFCChecker as IDataset;
pds.Delete();
}
}
}
}
catch (Exception ex)
{
LogHelper.WriteLog(typeof(TopoChecker), ex);
}
}
else if (IN_RuleType == "面自相交检查")
{
int errorCount = 0;
try
{
common.CreatShpFile(this.topoDir, spatialReference, esriGeometryType.esriGeometryPoint, checkname + idname);
IFeatureLayer pLayer = new FeatureLayer();
pLayer.FeatureClass = IN_FeatureClass;
TopoClassChecker topo = new TopoClassChecker();
IFeatureCursor cursor;
IQueryFilter filter = new QueryFilterClass();
filter.SubFields = pLayer.FeatureClass.OIDFieldName + "," + pLayer.FeatureClass.ShapeFieldName;
cursor = pLayer.FeatureClass.Search(filter, true);
IFeature feature = null;
object missing = Type.Missing;
while ((feature = cursor.NextFeature()) != null)
{
IPoint tempPoint = null;
StringBuilder builder = new StringBuilder();
IGeometryCollection shape = feature.Shape as IGeometryCollection;
for (int i = 0; i < shape.GeometryCount; i++)
{
esriNonSimpleReasonEnum enum2;
IPointCollection newPoints = shape.get_Geometry(i) as IPointCollection;
IRing ring = newPoints as IRing;
int num2 = 0;
if (ring.IsClosed)
{
num2 = 1;
}
PolylineClass o = new PolylineClass();
o.AddPointCollection(newPoints);
o.SpatialReference = spatialReference;
ITopologicalOperator3 @operator = o;
@operator.IsKnownSimple_2 = false;
if ([email protected]_IsSimpleEx(out enum2) && (enum2 == esriNonSimpleReasonEnum.esriNonSimpleSelfIntersections))
{
List <string> list2 = new List <string>();
List <string> list3 = new List <string>();
for (int j = num2; j < newPoints.PointCount; j++)
{
IPoint point = newPoints.get_Point(j);
tempPoint = point;
string item = point.X.ToString() + "," + point.Y.ToString();
if (list2.Contains(item))
{
if (!list3.Contains(item))
{
builder.Append(";");
builder.Append(item);
list3.Add(item);
}
}
else
{
list2.Add(item);
}
}
}
Marshal.ReleaseComObject(o);
o = null;
}
if (builder.Length > 0)
{
errorCount++;
string[] strArray = builder.ToString().Substring(1).Split(new char[] { ';' });
ESRI.ArcGIS.Geometry.IPointCollection pPointCollection1 = new ESRI.ArcGIS.Geometry.MultipointClass();
foreach (string str in strArray)
{
if (!string.IsNullOrEmpty(str))
{
string[] strArray2 = str.Split(new char[] { ',' });
double pX = double.Parse(strArray2[0]);
double pY = double.Parse(strArray2[1]);
ESRI.ArcGIS.Geometry.IPoint point = new ESRI.ArcGIS.Geometry.PointClass();
point.X = pX;
point.Y = pY;
common.GenerateSHPFile(point, ErrorFilePath);
}
}
}
}
Marshal.ReleaseComObject(cursor);
if (errorCount == 0)
{
string shpFolder = System.IO.Path.GetDirectoryName(ErrorFilePath);
IWorkspaceFactory pWorkspaceFac = new ShapefileWorkspaceFactoryClass();
IWorkspace pWorkSpace = pWorkspaceFac.OpenFromFile(ErrorFilePath, 0);
IFeatureWorkspace pFeatureWorkSpace = pWorkSpace as IFeatureWorkspace;
if (System.IO.File.Exists(ErrorFilePath))
{
IFeatureClass pFCChecker = pFeatureWorkSpace.OpenFeatureClass(ErrorFilePath);
if (pFCChecker != null)
{
IDataset pds = pFCChecker as IDataset;
pds.Delete();
}
}
}
}
catch (Exception ex)
{
LogHelper.WriteLog(typeof(TopoChecker), ex);
}
if (!DicTopoError.ContainsKey(idname))
{
DicTopoError.Add(idname, errorCount);
}
else
{
DicTopoError[idname] = errorCount;
}
}
else if (IN_RuleType == "缝隙检查")
{
int errorCount = 0;
try
{
common.CreatShpFile(this.topoDir, spatialReference, esriGeometryType.esriGeometryPolygon, checkname + idname);
this.m_hookHelper = m_hookHelper;
IFeatureClass pFClass = IN_FeatureClass;
//获取空间参考
IGeometry geometryBag = new GeometryBagClass();
IGeoDataset geoDataset = pFClass as IGeoDataset;
geometryBag.SpatialReference = geoDataset.SpatialReference;
////属性过滤
IFeatureCursor featureCursor = pFClass.Search(null, false);
// 遍历游标
IFeature currentFeature = featureCursor.NextFeature();
IGeometryCollection geometryCollection = geometryBag as IGeometryCollection;
object missing = Type.Missing;
while (currentFeature != null)
{
geometryCollection.AddGeometry(currentFeature.Shape, ref missing, ref missing);
currentFeature = featureCursor.NextFeature();
}
// 合并要素
ITopologicalOperator unionedPolygon = null;
unionedPolygon = new Polygon() as ITopologicalOperator;
unionedPolygon.ConstructUnion(geometryCollection as IEnumGeometry);
Marshal.ReleaseComObject(featureCursor);
IPolygon4 pMergerPolygon = unionedPolygon as IPolygon4;
IGeometryBag pOutGeometryBag = pMergerPolygon.ExteriorRingBag; //获取外部环
IGeometryCollection pOutGmtyCollection = pOutGeometryBag as IGeometryCollection;
for (int i = 0; i < pOutGmtyCollection.GeometryCount; i++) //对外部环遍历
{
IGeometry pOutRing = pOutGmtyCollection.get_Geometry(i); //外部环
//【此处可以对外部环进行操作】
IPointCollection pOutRingCollection = pOutRing as IPointCollection;
for (int j = 0; j < pOutRingCollection.PointCount; j++)
{
IPoint pOutRingPoint = pOutRingCollection.get_Point(j);//获取外环上的点
}
IGeometryBag pInteriotGeometryBag = pMergerPolygon.get_InteriorRingBag(pOutRing as IRing); //获取内部环
IGeometryCollection pInteriorGeometryCollection = pInteriotGeometryBag as IGeometryCollection;
for (int j = 0; j < pInteriorGeometryCollection.GeometryCount; j++)
{
ISegmentCollection SegCol = pInteriorGeometryCollection.get_Geometry(j) as ISegmentCollection;
IPolygon PPolygon = new PolygonClass();
ISegmentCollection newSegCol = PPolygon as ISegmentCollection;
newSegCol.AddSegmentCollection(SegCol);
//pInteriorGeometry即为多边形的内部环
IGeometry inRing = PPolygon as IGeometry;
inRing.SpatialReference = geometryBag.SpatialReference;
IArea area = inRing as IArea;
Double getarea = System.Math.Abs(Convert.ToDouble(area.Area));
if (inputtext == null || inputtext == "" || getarea < Convert.ToDouble(inputtext))
{
Boolean flag = true;
ISpatialFilter filter = new SpatialFilterClass();
filter.Geometry = inRing;
filter.SubFields = pFClass.OIDFieldName + "," + pFClass.ShapeFieldName;
filter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;
IFeatureCursor o = pFClass.Search(filter, true);
IFeature feature = o.NextFeature();
while (feature != null && flag)
{
IPolygon4 pPolygon = feature.Shape as IPolygon4;
IGeometryBag iOutGeometryBag = pPolygon.ExteriorRingBag; //获取外部环
IGeometryCollection iOutGmtyCollection = iOutGeometryBag as IGeometryCollection;
for (int m = 0; m < iOutGmtyCollection.GeometryCount && flag; m++) //对外部环遍历
{
IGeometry outGeo = iOutGmtyCollection.get_Geometry(m);
IGeometryCollection polyGonGeo = new PolygonClass();
polyGonGeo.AddGeometry(outGeo);
IPolygon iPolygon = polyGonGeo as IPolygon;
iPolygon.SimplifyPreserveFromTo();
IRelationalOperator2 pRelationalOperator2 = iPolygon as IRelationalOperator2;
if (!pRelationalOperator2.Contains(inRing))
{
errorCount++;
common.GenerateSHPFile(inRing, ErrorFilePath);
flag = false;
}
}
feature = o.NextFeature();
}
Marshal.ReleaseComObject(o);
}
}
if (errorCount == 0)
{
string shpFolder = System.IO.Path.GetDirectoryName(ErrorFilePath);
IWorkspaceFactory pWorkspaceFac = new ShapefileWorkspaceFactoryClass();
IWorkspace pWorkSpace = pWorkspaceFac.OpenFromFile(ErrorFilePath, 0);
IFeatureWorkspace pFeatureWorkSpace = pWorkSpace as IFeatureWorkspace;
if (System.IO.File.Exists(ErrorFilePath))
{
IFeatureClass pFCChecker = pFeatureWorkSpace.OpenFeatureClass(ErrorFilePath);
if (pFCChecker != null)
{
IDataset pds = pFCChecker as IDataset;
pds.Delete();
}
}
}
}
}
catch (Exception ex)
{
LogHelper.WriteLog(typeof(TopoChecker), ex);
}
if (!DicTopoError.ContainsKey(idname))
{
DicTopoError.Add(idname, errorCount);
}
else
{
DicTopoError[idname] = errorCount;
}
}
else if (IN_RuleType == "面重叠检查")
{
IFeatureClass outIFC = null;
this.m_hookHelper = m_hookHelper;
ESRI.ArcGIS.AnalysisTools.Intersect intersect = new ESRI.ArcGIS.AnalysisTools.Intersect();
string outString = this.topoDir + "\\" + checkname + idname + ".shp";
intersect.in_features = common.GetPathByName(TABLENAME);
intersect.out_feature_class = outString;
Geoprocessor geoProcessor = new Geoprocessor();
geoProcessor.OverwriteOutput = true;
try
{
if (this.gp == null)
{
this.gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
}
IGeoProcessorResult result = (IGeoProcessorResult)this.gp.Execute(intersect, null);
if (result.Status != ESRI.ArcGIS.esriSystem.esriJobStatus.esriJobSucceeded)
{
Console.WriteLine("gp工具执行错误");
}
else
{
outIFC = this.gp.Open(result.ReturnValue) as IFeatureClass;
}
}
catch (Exception ex)
{
LogHelper.WriteLog(typeof(TopoChecker), ex);
Console.WriteLine(ex.Message);
}
if (outIFC.FeatureCount(null) == 0)
{
string shpFolder = System.IO.Path.GetDirectoryName(ErrorFilePath);
IWorkspaceFactory pWorkspaceFac = new ShapefileWorkspaceFactoryClass();
IWorkspace pWorkSpace = pWorkspaceFac.OpenFromFile(ErrorFilePath, 0);
IFeatureWorkspace pFeatureWorkSpace = pWorkSpace as IFeatureWorkspace;
if (System.IO.File.Exists(ErrorFilePath))
{
IFeatureClass pFCChecker = pFeatureWorkSpace.OpenFeatureClass(ErrorFilePath);
if (pFCChecker != null)
{
IDataset pds = pFCChecker as IDataset;
pds.Delete();
}
}
}
if (!DicTopoError.ContainsKey(idname))
{
DicTopoError.Add(idname, outIFC.FeatureCount(null));
}
else
{
DicTopoError[idname] = outIFC.FeatureCount(null);
}
}
else if (IN_RuleType == "面重叠检查（与其他图层）")
{
IFeatureClass outIFC = null;
this.m_hookHelper = m_hookHelper;
ESRI.ArcGIS.AnalysisTools.Intersect intersect = new ESRI.ArcGIS.AnalysisTools.Intersect();
string outString = this.topoDir + "\\" + checkname + idname + ".shp";
intersect.in_features = @common.GetPathByName(TABLENAME) + ";" + @common.GetPathByName(SUPTABLE);
intersect.out_feature_class = outString;
Geoprocessor geoProcessor = new Geoprocessor();
geoProcessor.OverwriteOutput = true;
try
{
if (this.gp == null)
{
this.gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
}
IGeoProcessorResult result = (IGeoProcessorResult)this.gp.Execute(intersect, null);
if (result.Status != ESRI.ArcGIS.esriSystem.esriJobStatus.esriJobSucceeded)
{
Console.WriteLine("gp工具执行错误");
}
else
{
outIFC = this.gp.Open(result.ReturnValue) as IFeatureClass;
}
}
catch (Exception ex)
{
Console.WriteLine(ex.Message);
LogHelper.WriteLog(typeof(TopoChecker), ex);
}
if (outIFC.FeatureCount(null) == 0)
{
string shpFolder = System.IO.Path.GetDirectoryName(ErrorFilePath);
IWorkspaceFactory pWorkspaceFac = new ShapefileWorkspaceFactoryClass();
IWorkspace pWorkSpace = pWorkspaceFac.OpenFromFile(ErrorFilePath, 0);
IFeatureWorkspace pFeatureWorkSpace = pWorkSpace as IFeatureWorkspace;
if (System.IO.File.Exists(ErrorFilePath))
{
IFeatureClass pFCChecker = pFeatureWorkSpace.OpenFeatureClass(ErrorFilePath);
if (pFCChecker != null)
{
IDataset pds = pFCChecker as IDataset;
pds.Delete();
}
}
}
if (!DicTopoError.ContainsKey(idname))
{
DicTopoError.Add(idname, outIFC.FeatureCount(null));
}
else
{
DicTopoError[idname] = outIFC.FeatureCount(null);
}
}
Marshal.ReleaseComObject(IN_FeatureClass);
if (SUPTABLE != null)
{
Marshal.ReleaseComObject(IN_Sup_FeatureClass);
}
}
/// <summary>
/// 缝隙检查。融合后去内环。
/// </summary>
/// <param name="pFeature"></param>
/// <returns></returns>
public List <Dictionary <string, IGeometry> > CheckFeatureGap(IFeatureClass pFClass, string inputtext)
{
List <Dictionary <string, IGeometry> > listGeo = new List <Dictionary <string, IGeometry> >();
if (pFClass == null)
{
return(null);
}
//Stopwatch watch = new Stopwatch();
//watch.Start();
//获取空间参考
IGeometry geometryBag = new GeometryBagClass();
IGeoDataset geoDataset = pFClass as IGeoDataset;
geometryBag.SpatialReference = geoDataset.SpatialReference;
////属性过滤
//ISpatialFilter queryFilter = new SpatialFilterClass();
//queryFilter.SubFields = "Shape";
IFeatureCursor featureCursor = pFClass.Search(null, false);
// 遍历游标
IFeature currentFeature = featureCursor.NextFeature();
if (currentFeature == null)
{
return(null);
}
IGeometryCollection geometryCollection = geometryBag as IGeometryCollection;
object missing = Type.Missing;
while (currentFeature != null)
{
geometryCollection.AddGeometry(currentFeature.Shape, ref missing, ref missing);
currentFeature = featureCursor.NextFeature();
}
// 合并要素
ITopologicalOperator unionedPolygon = null;
if (pFClass.ShapeType == esriGeometryType.esriGeometryPoint)
{
unionedPolygon = new Multipoint() as ITopologicalOperator;
unionedPolygon.ConstructUnion(geometryCollection as IEnumGeometry);
}
else if (pFClass.ShapeType == esriGeometryType.esriGeometryPolyline)
{
unionedPolygon = new Polyline() as ITopologicalOperator;
unionedPolygon.ConstructUnion(geometryCollection as IEnumGeometry);
}
else
{
unionedPolygon = new Polygon() as ITopologicalOperator;
unionedPolygon.ConstructUnion(geometryCollection as IEnumGeometry);
}
Marshal.ReleaseComObject(featureCursor);
IPolygon4 pMergerPolygon = unionedPolygon as IPolygon4;
IGeometryBag pOutGeometryBag = pMergerPolygon.ExteriorRingBag; //获取外部环
IGeometryCollection pOutGmtyCollection = pOutGeometryBag as IGeometryCollection;
for (int i = 0; i < pOutGmtyCollection.GeometryCount; i++) //对外部环遍历
{
IGeometry pOutRing = pOutGmtyCollection.get_Geometry(i); //外部环
//【此处可以对外部环进行操作】
IPointCollection pOutRingCollection = pOutRing as IPointCollection;
for (int j = 0; j < pOutRingCollection.PointCount; j++)
{
IPoint pOutRingPoint = pOutRingCollection.get_Point(j);//获取外环上的点
}
IGeometryBag pInteriotGeometryBag = pMergerPolygon.get_InteriorRingBag(pOutRing as IRing); //获取内部环
IGeometryCollection pInteriorGeometryCollection = pInteriotGeometryBag as IGeometryCollection;
for (int j = 0; j < pInteriorGeometryCollection.GeometryCount; j++)
{
ISegmentCollection SegCol = pInteriorGeometryCollection.get_Geometry(j) as ISegmentCollection;
IPolygon PPolygon = new PolygonClass();
ISegmentCollection newSegCol = PPolygon as ISegmentCollection;
newSegCol.AddSegmentCollection(SegCol);
//pInteriorGeometry即为多边形的内部环
IGeometry inRing = PPolygon as IGeometry;
inRing.SpatialReference = geometryBag.SpatialReference;
IArea area = inRing as IArea;
Double getarea = System.Math.Abs(Convert.ToDouble(area.Area));
if (inputtext == null || inputtext == "" || getarea < Convert.ToDouble(inputtext))
{
Boolean flag = true;
ISpatialFilter filter = new SpatialFilterClass();
filter.Geometry = inRing;
filter.SubFields = pFClass.OIDFieldName + "," + pFClass.ShapeFieldName;
filter.SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects;
IFeatureCursor o = pFClass.Search(filter, true);
IFeature feature = o.NextFeature();
while (feature != null && flag)
{
Console.WriteLine(feature.OID);
IPolygon4 pPolygon = feature.Shape as IPolygon4;
IGeometryBag iOutGeometryBag = pPolygon.ExteriorRingBag; //获取外部环
IGeometryCollection iOutGmtyCollection = iOutGeometryBag as IGeometryCollection;
for (int m = 0; m < iOutGmtyCollection.GeometryCount && flag; m++) //对外部环遍历
{
IGeometry outGeo = iOutGmtyCollection.get_Geometry(m);
IGeometryCollection polyGonGeo = new PolygonClass();
polyGonGeo.AddGeometry(outGeo);
IPolygon iPolygon = polyGonGeo as IPolygon;
iPolygon.SimplifyPreserveFromTo();
IRelationalOperator2 pRelationalOperator2 = iPolygon as IRelationalOperator2;
if (!pRelationalOperator2.Contains(inRing))
{
Dictionary <string, IGeometry> itemDic = new Dictionary <string, IGeometry>();
itemDic.Add(feature.OID.ToString(), inRing);
listGeo.Add(itemDic);
flag = false;
}
}
feature = o.NextFeature();
}
Marshal.ReleaseComObject(o);
}
}
}
if (listGeo.Count > 0)
{
return(listGeo);
}
else
{
return(null);
}
}
