简介ARCGIS检查锐角工具是一份面向ARCMAP桌面端用户的多边形几何质量检查功能包定位在空间数据编辑与成果质检环节。工具通过设定角度阈值自动识别多边形中的尖锐角帮助测绘、规划及GIS数据生产人员提前发现可能影响分析精度或后续地理处理结果的几何缺陷。包体非常精简共3个文件其中锐角检查工具.tbx用于在ArcToolbox中直接调用jianjiao_check.py为可定制参数的Python脚本使用说明.pdf则对工具原理、参数设置与操作流程作了详细演示兼顾操作与二次开发。资源大小仅521KB轻量实用已有5436人学习下载。借助该功能包用户可以快速掌握锐角检查的完整思路并结合手动编辑工具对问题顶点进行修正提升空间数据入库前的质量管控效率。 在国土调查、不动产确权、规划数据整理这些活儿里泡久了的兄弟大概率都有过这种经历从CAD转过来的地块或者早年手勾的边界好好的一张面放大一看边角上全是尖尖的锐角有的角度小到只有几度。这种图斑平时看着没事可一旦轮到入库质检、拓扑检查、面积统计轮番上阵它就变成了一串串报错逼着你回头一点一点清理。我刚做数据整理那几年第一次批量处理农村土地确权的几千个疑似图斑就是被这批锐角磨到没脾气才开始满世界找检查锐角的插件最后万不得已决定自己拿ArcPy写一个。这篇就把这个“ARCGIS检查锐角工具”从原理到落地完整拆一遍代码直接拿走就能用。1. 锐角从哪来不改会出什么事1.1 三类最常见的“产锐角”场景锐角不是凭空长出来的大部分情况下都能找到源头。我经手过的数据里锐角几乎都来自下面三个场景。第一个是CAD转GIS的遗留问题。很多地块边线在CAD里是PLINE绘制的节点特别密转成SHP之后这些冗余折点全部保留下来两条原本该是平滑衔接的边中间硬生生多出一串小折线角度自然就失控了。尤其是地形图转地类图斑、规划红线转入库数据这种情况占比非常大。第二个是人工勾绘时的手抖。用编辑器沿屏幕描边界的时候连续点击鼠标稍微犹豫一下就可能留下Z字形的折线放大看角度极小。这类问题在做外业调绘数字化时特别常见勾的人自己都没注意到。第三个是几何操作后的残边。面要素裁剪、合并、擦除、拓扑修复之后边界处经常产生细小的“碎屑”折角。比如两个面做Union之后公共边处会出现很短的线段这些短边和相邻边之间就容易形成锐角。1.2 锐角不改会出哪些事很多人觉得锐角只是“难看”不影响面积所以懒得处理。但实际业务里它带来的麻烦远不止观感问题。首先是质检报错。国土空间调查、确权登记、年度变更调查这些任务都有一套严格的数据质检规则锐角是常见的检查项之一。我记得有一年做省级汇交质检软件一次性报了几千个尖锐角问题涉及图斑几百个只能一个个退回处理工期直接被拖了两天。其次是拓扑与入库问题。锐角往往伴随着短边而短边在拓扑检查里很容易触发“不能有悬挂边”“不能有伪节点”这类规则直接导致数据入库失败。尤其是一些基于ArcGIS Server发布的数据服务对图层几何质量要求更高有锐角的数据发上去前端编辑经常出现奇怪的吸附行为。再一个是制图和空间分析的失真。制图时锐角处标注压盖严重晕线和填充符号在尖角位置显示得乱七八糟。空间分析时如果你要做基于角度阈值的缓冲区、视线分析或者要素折弯分析锐角会影响计算结果的合理性。2. 检查锐角之前先把原理和方案想明白2.1 向量夹角公式锐角检查唯一的数学依据锐角检查说白了就是算“角”。一个面要素在数学上是由封闭折线围成的每个折点的左右两侧各有一条边这两条边从折点出发分别连向前一个点和后一个点。我们现在要算的就是这两个相邻边在折点处的夹角。计算夹角的公式是向量点积的标准写法cosθ (v1 · v2) / (|v1| × |v2|)其中v1是“当前点→前一个点”的向量v2是“当前点→后一个点”的向量θ就是我们要的夹角。这个公式在初中数学里就有但真正在ArcGIS里落地时会有几个细节要注意。第一个细节是闭合环。面要素的边界是闭合的所以第一个点的前一个点其实是最后一个点最后一个点的后一个点其实是第一个点。用Python处理的时候我会用取模运算 (i-1) % n 和 (i1) % n 来拿前后点这样就不用单独处理首尾。第二个细节是重复点。如果数据里存在相邻两个点完全重合的情况那其中一个向量长度为0点积公式就失效了。这种情况我一般在计算前先过滤掉重复点或者在计算时加一个判断返回空值跳过。第三个细节是角度的取值范围。点积公式算出来的是0到180度之间的夹角我们关心的是“小于某个阈值”的锐角所以直接在结果上做小于比较就行不用考虑大于180度的情况。2.2 阈值定多少没有万能值但有推荐基线每次有人问我“阈值该设成几度”我都会反问一句你这批数据是干什么用的因为不同业务对锐角的要求差别挺大。我根据自己的使用经验整理了下面这个参考表业务场景推荐阈值说明国土调查/变更调查5°~10°图斑形状允许较自然的尖角阈值太大会误伤真实边界不动产确权/宗地10°~15°宗地边界强调规整锐角往往代表勾绘错误规划入库/报批15°~20°规划地块形状普遍规整角度太小会被判定不合格一般数据预处理10°折中选项适合不确定业务标准时先跑一轮摸底需要注意的是阈值设置得越小查出来的数量越少但漏检的风险越高设置得越大误报越多后期人工复核成本越高。建议第一次先用10度跑一遍看看数量和数据特征再根据结果调整。2.3 方案对比为什么推荐ArcPy脚本在动手之前我先对比过几种可行的方案。用字段计算器可以算吗理论上可以在“计算几何”里能拿到每个点坐标但你要在字段里维护多个点的坐标、再做向量运算、还要处理闭合环写出来的表达式又长又难维护性能也差。我试过一次就放弃了。用ArcGIS自带的拓扑规则可以吗不行。ArcGIS的拓扑规则集合里并没有专门针对“最小夹角”的规则你要做锐角检查只能借助其他间接方式去推断。很多人搜遍了工具箱也没找到对应功能最后才去找插件。用第三方插件/Add-In呢网上确实流传过一些“检查尖锐角插件”但版本兼容是个大问题10.x的addin拿到Pro里根本装不上而且来路不明的工具在数据生产环境里本身就是风险。所以我最终选择自己写ArcPy脚本——可控、可复现、改参数很方便而且是纯Python迁移成本低。这也是我想在这篇里重点分享的方案。3. 手把手用ArcPy实现一个可复用的锐角检查工具3.1 脚本思路与核心代码我的核心想法很简单把所有面要素的每个折点都遍历一遍如果某个折点的角度小于阈值就把这个点的位置、所属要素编号、折点编号、角度值全部记录下来输出成一个点要素。这样你既能快速定位到问题图斑又能直接在结果里看到每个角的具体度数。下面是完整脚本我用的是ArcGIS的arcpy.da游标性能和稳定性都比老版游标好得多。# -*- coding: utf-8 -*- import arcpy import math def calc_angle(px, py, qx, qy, rx, ry): 计算向量qp与qr的夹角返回角度制数值若存在零向量返回None ax, ay px - qx, py - qy # 当前点 - 前一点 bx, by rx - qx, ry - qy # 当前点 - 后一点 len_a math.hypot(ax, ay) len_b math.hypot(bx, by) if len_a 0 or len_b 0: return None cos_theta (ax * bx ay * by) / (len_a * len_b) cos_theta max(-1.0, min(1.0, cos_theta)) # 防止浮点误差越界 return math.degrees(math.acos(cos_theta)) def check_sharp_corners(in_fc, out_fc, threshold10.0): 检查面要素锐角将小于阈值的折点输出为点要素 desc arcpy.Describe(in_fc) sr desc.spatialReference out_path, out_name arcpy.os.path.split(out_fc) arcpy.CreateFeatureclass_management(out_path, out_name, POINT, spatial_referencesr) arcpy.AddField_management(out_fc, SourceOID, LONG) arcpy.AddField_management(out_fc, PartIdx, LONG) arcpy.AddField_management(out_fc, VtxIdx, LONG) arcpy.AddField_management(out_fc, Angle, DOUBLE) fields [OID, SHAPE] out_fields [SHAPE, SourceOID, PartIdx, VtxIdx, Angle] with arcpy.da.SearchCursor(in_fc, fields) as s_cur, \ arcpy.da.InsertCursor(out_fc, out_fields) as i_cur: for oid, shape in s_cur: if shape is None: continue part_idx 0 for part in shape: pts [(pnt.X, pnt.Y) for pnt in part] # 去掉闭合环末尾的重复点 while len(pts) 1 and pts[0] pts[-1]: pts.pop() n len(pts) if n 3: part_idx 1 continue for i in range(n): px, py pts[(i - 1) % n] qx, qy pts[i] rx, ry pts[(i 1) % n] ang calc_angle(px, py, qx, qy, rx, ry) if ang is not None and ang threshold: i_cur.insertRow(( arcpy.PointGeometry(arcpy.Point(qx, qy), sr), oid, part_idx, i, ang )) part_idx 1 if __name__ __main__: in_fc rD:\data\polygons.shp # 改成你的面数据路径 out_fc rD:\data\sharp_corners.shp # 输出点数据路径 threshold 10.0 # 角度阈值 check_sharp_corners(in_fc, out_fc, threshold) arcpy.AddMessage(检查完成结果已输出到 out_fc)这里有几个实现细节值得说一下。用for part in shape可以直接同时处理多部件要素和带洞要素。多部件面的每个外环或内环都是一个part脚本会逐part遍历。带洞要素的洞也是part计算逻辑一致不用担心漏掉。闭合环的处理上我先判断pts[0] pts[-1]相等就把最后一个点去掉这样后续用取模运算的时候不会把同一个点算两遍。这个细节不处理的话闭合点的角度会被错误地计算在“第一个点”和“最后一个点”上。零向量防御是必须的。数据里经常有相邻点坐标完全相同的情况不判断的话division by zero异常会直接让脚本中断。我的做法是返回None跳过这个点不影响整体执行。3.2 在ArcMap和ArcGIS Pro里跑起来脚本写好后最快的运行方式就是在ArcMap或者ArcGIS Pro的Python窗口里直接粘贴执行。打开ArcMap在“地理处理 Python”里打开Python窗口把if __name__之前的函数定义贴进去再单独执行带上实际路径的那几行调用代码。ArcGIS Pro也是一样在“分析 Python”里打开Python窗格。有一点要注意ArcMap 10.x用的Python 2.7ArcGIS Pro用的是Python 3.x我上面的代码在两个环境里都是兼容的只要你用的是ArcGIS 10.4及以上版本arcpy.da模块都能用。如果你要长期复用我推荐封装成脚本工具。做法是在ArcToolbox里新建一个脚本工具箱添加脚本参数设置四个输入面要素、输出点要素、角度阈值。这样以后在工具箱里点两下就能跑不用每次都打开Python窗口改路径。我自己的工具箱里面还加了个“阈值默认值10”的提示避免同事用的时候不清不楚。3.3 结果怎么用定位、筛选与修复跑完之后输出的是一个点要素每个点代表一个小于阈值的锐角。我习惯把它直接叠加到原始数据上用Select By Location或者属性查询就能快速找到问题图斑。输出字段里SourceOID对应原始要素的FID方便你回原表里定位PartIdx是部件编号VtxIdx是折点编号这两个字段在人工编辑时特别有用——你打开编辑工具跳到指定部件和折点眨眼就能找到问题位置。修复锐角最常用的有三种方式。第一种是手动编辑用“编辑折点”工具把问题折点拖动或者删除适合数量少的时候。第二种是使用“消除”工具Eliminate把短边合并到相邻边适合由短边引起的锐角。第三种是“简化面”Simplify Polygon通过设置容差重新概括边界适合数据量大、批量处理的情况。但要注意简化面是个“动刀”的操作容差太大会明显改变面积所以操作前一定要备份原始数据。我的习惯是先跑脚本摸清锐角分布再用简化面处理处理完再跑一遍脚本直到结果数量降到可接受范围。4. 实操中的高频问题与排查技巧4.1 坐标系没设对角度结果就废了这是最坑、也最容易忽略的问题。如果数据是地理坐标系GCS也就是经纬度存储的直接按经纬度坐标计算向量夹角会得到完全失真的结果。原因很简单经度方向和纬度方向在地球表面上并不“等长”同一个角度值在不同纬度下的实际视觉大小完全不同。所以正确做法是检查锐角之前把数据投影到合适的投影坐标系再动手。比如数据在某个城市就用当地的高斯-克吕格投影或者UTM分带投影。投影之后角度的几何意义才成立。运行脚本前我建议先看一眼数据的属性确保Spatial Reference一栏是投影坐标系别拿经纬度数据硬跑。4.2 多部件、带洞、弧段要素分别怎么处理我的脚本已经兼容多部件和带洞要素但弧段真曲线情况不一样。ArcGIS里有些要素的边界并不是折线而是贝塞尔曲线或者圆弧比如从CAD转过来的部分数据。这种情况下直接遍历折点会把弧段当成一条直线边处理锐角自然查不出来。怎么判断你的数据里有没有弧段可以用“要素折点转点”工具看点的密度或者用arcpy查看几何属性里的curve字段。如果确实有弧段先对要素类使用“加密”Densify工具把曲线转成密集折点再跑锐角检查。加密时容差建议设置为0.5米或1米太小数据量暴涨太大查不出细节。4.3 数据量特别大运行卡死怎么办我最早跑20万多个面的检查时也遇到过卡顿后来优化了两点基本解决了。一是全程使用da.SearchCursor和da.InsertCursor这两个新游标比老版ArcObjects接口快很多是首选的性能方案。二是不要每查到一个角就立刻写入要素类先攒在内存的列表里最后一次性插入。脚本里我直接用了InsertCursor逐行插入但如果你发现很慢可以把结果先存在Python列表最后统一写库。4.4 修完一批又冒出一批这是正常现象很多人修复完锐角后发现简化面工具又制造了新的锐角于是怀疑工具不稳。这其实是几何概括的必然结果——当你删除一个折点时相邻两个边被拉直有可能形成新的夹角。所以锐角处理是个“检查—修复—再检查”的循环过程不用指望一遍跑完。我的建议是设一个目标值比如跑完剩余锐角少于10个就算合格然后人工收尾。最后说点个人体会。锐角检查这件事查出来只是第一步真正花时间的是修复和复核。但如果没有一个趁手的工具光定位问题就能让人崩溃。我现在做入库质检标准动作是先跑这个脚本出锐角分布再跑拓扑检查最后按问题清单逐条过。脚本封装成工具箱里的按钮之后团队里的新人也能直接上手。建议有条件的话你也把它做成固定工具留在自己的工具箱里以后什么数据进来先扫一遍能省下不少晚上加班的功夫。本文还有配套的精品资源点击获取