ArcGIS消除叠加操作生成的碎图斑或狭长面

在ArcGIS中叠加操作会带来碎图斑或狭长面，而这些碎图斑或狭长面又会影响后续的判断。所以就需要消除。

例如：

具体操作有以下方式：

1、最简单的的方法就是使用ArcGIS提供的“消除工具”。在"数据管理工具"——“制图综合”——“消除”。

2、使用Python脚本

# -*-coding:gbk-*-import arcpy
from arcpy import env
import time# 融合指定条件的图斑try:source_gdb_path = "F:/GIS测试数据/测试.gdb"  # 原始图层工作空间in_features = "POLYGON"   # 输入图层expression = "SHAPE_AREA <50 or SHAPE_Area/ SHAPE_Length<0.2"  # 融合限制条件env.workspace = source_gdb_pathtime_begin = time.time()print "消除指定条件下的要素{0}".format(expression)tempLayer = "block_layer"# SearchCursorfcount = len([feature[0] for feature in arcpy.da.SearchCursor(in_features, "SHAPE_AREA", expression)])while (fcount > 0):print('fcount:{0}  Time : {1} s'.format(fcount, time.time() - time_begin))tempLayer = "{0}_b".format(tempLayer)# MakeFeatureLayerarcpy.MakeFeatureLayer_management(in_features, tempLayer)in_features = "{0}_E".format(in_features)if arcpy.Exists(in_features):arcpy.Delete_management(in_features)# SelectLayerByAttribute to define feature to be eliminatedarcpy.SelectLayerByAttribute_management(tempLayer, "NEW_SELECTION", expression)# Eliminatearcpy.Eliminate_management(tempLayer, in_features, "AREA")fcount2 = len([feature[0] for feature in arcpy.da.SearchCursor(in_features, "SHAPE_AREA", expression)])if fcount2 == fcount:breakelse:fcount = fcount2print "Total Time {0}".format(time.time() - time_begin)
except Exception, e:print e.message.decode("utf_8")

3、AE+C#实现

/// <summary>
/// 消除碎图斑
/// </summary>
/// <param name="pFeatureClass">要素类</param>
/// <param name="whereClause">消除条件："SHAPE_AREA <50 or SHAPE_AREA/SHAPE_Length<0.2"</param>
public static void Eliminate(IFeatureClass pFeatureClass, string whereClause)
{Dictionary<int, IGeometry> dicOidsGeo = new Dictionary<int, IGeometry>();ISpatialFilter pSpatialFilter = new SpatialFilterClass { WhereClause = whereClause, SubFields = $"SHAPE,{pFeatureClass.OIDFieldName}" };IFeatureCursor pFeatureCursor = pFeatureClass.Search(pSpatialFilter, true);IFeature pFeature = null;while ((pFeature = pFeatureCursor.NextFeature()) != null){(pFeature.ShapeCopy as ITopologicalOperator).Simplify();dicOidsGeo[pFeature.OID] = pFeature.ShapeCopy;}System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatureCursor);int count = dicOidsGeo.Count;while (count > 0){Console.WriteLine($"个数：{count}，时间：{DateTime.Now}");int count2 = ExcudeEliminate(dicOidsGeo, pFeatureClass);if (count2 == count)break;elsecount = count2;}
}/// <summary>
/// 消除碎图斑，融到相邻面积最大面
/// </summary>
/// <param name="dicOidsGeo">几何集合</param>
/// <param name="pFeatureClass">源要素类</param>
/// <returns>未消除的个数</returns>
public static int ExcudeEliminate(Dictionary<int, IGeometry> dicOidsGeo, IFeatureClass pFeatureClass)
{Func<IFeature, IGeometry, IGeometry> func = (pFeature1, ge) =>{ITopologicalOperator pTopolo = pFeature1.ShapeCopy as ITopologicalOperator;IGeometry geo = pTopolo.Union(ge);(geo as ITopologicalOperator).Simplify();return geo;};List<int> lisOids = new List<int>();  //记录已经融合的要素IDIFeatureCursor pFeatureCursor = null;IFeature pFeature = null;ISpatialFilter pSpatialFilter = new SpatialFilterClass { SpatialRel = esriSpatialRelEnum.esriSpatialRelIntersects, GeometryField = pFeatureClass.ShapeFieldName, PostfixClause = "Order by SHAPE_AREA desc" };foreach (var ge in dicOidsGeo){pSpatialFilter.Geometry = ge.Value;pFeatureCursor = pFeatureClass.Update(pSpatialFilter, false);pFeature = pFeatureCursor.NextFeature();if (pFeature != null){pFeature.Shape = func(pFeature, ge.Value);pFeatureCursor.UpdateFeature(pFeature);lisOids.Add(ge.Key);}}System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatureCursor);if (lisOids.Count == 0)return dicOidsGeo.Count;var str = string.Empty;//删除lisOids.ForEach(x =>{str += $",{x}";dicOidsGeo.Remove(x);});str = str.Substring(1);IQueryFilter pQueryFilter = new QueryFilterClass();pQueryFilter.WhereClause = $"{pFeatureClass.OIDFieldName} in ({str})";  //str超过1000时需要分组pFeatureCursor = pFeatureClass.Update(pQueryFilter, false);while((pFeature= pFeatureCursor.NextFeature())!= null){pFeatureCursor.DeleteFeature();}System.Runtime.InteropServices.Marshal.ReleaseComObject(pFeatureCursor);return dicOidsGeo.Count;
}

ArcGIS消除叠加操作生成的碎图斑或狭长面相关推荐

Arcgis中消除子流域划分时出现的零碎图斑或狭长面（Eliminate）
在ArcGIS中叠加操作会带来碎图斑或狭长面,而这些碎图斑或狭长面又会影响后续的判断.所以就需要消除. 例如: 具体操作有以下方式: 1.最简单的的方法就是使用ArcGIS提供的"消除工具& ...
Arcgis 熟练和操作
Arcgis 熟练和操作 arcgis操作实验注: 对于arcgis使用有一段时间了,前期是对于矢量,栅格数据的处理有了一个初步的了解,中期是在摸索中继续多层次的使用数据,以及逐步延伸到了空间分析, ...
Arcgis利用dem数据生成等高线
Arcgis利用dem数据生成等高线问题描述如果你在网上查,如何利用dem生成等高线,那么所有教程基本都是一样的,告诉你在工作箱中找到等值线工具,填入参数即可生成等高线.贴几个网上的链接教程: a ...
Python自动生成代码 - 通过tkinter图形化操作生成代码框架
Python自动生成代码 - 通过tkinter图形化操作生成代码框架背景脚本代码 Demo_CodeGenerator.py display.py FileHandler.py: 脚本运行结果: ...
使用ArcGIS为tif影像生成tfw文件
使用ArcGIS为tif影像生成tfw文件有时候我们拿到的tif影像文件的地理坐标信息直接存储在tif文件头中,没有配套的tfw文件.然而,某些软件要求必须有tfw文件,所以就会产生这一需求:将ti ...
AS下怎样录制APP操作生成Gif动态图
本文主要讲解Android下怎样录制APP操作生成gif动态图. Android App开发完成后,我们希望录制个gif做个展示.视频也可以展示,但是它没有gif轻量简单省流量.所以这篇文章就是讲解下 ...
【ARCGIS矢量化操作相关】
这里写自定义目录标题 ARCGIS矢量化操作相关 ARCGIS矢量化操作相关 ArcScan一直是灰色在自定义的扩展模块查看,打开以上窗口,查看是否勾选了ArcScan 确定勾选了后发现还是灰色的再 ...
poi实现操作生成word表格和操作word中的图表数据
poi实现操作生成word表格和操作word中的图表数据写入word表格 1.1首先要将写入的数据整理好,将其封装在一个类中: 1.2 其次将数据写入word,写成表格 2 在word表格中动态修改 ...
关于Arcgis(geoscene) 的objectid生成
objectid是arcgis自己定义的生成序列,不能随便填写,不然会和arcmap.arcgis pro .geoscene pro 导入的数据冲突 1.oracle获取Objectid序列的方法 ...

ArcGIS消除叠加操作生成的碎图斑或狭长面

具体操作有以下方式：

ArcGIS消除叠加操作生成的碎图斑或狭长面相关推荐

最新文章

热门文章