简介面向GIS从业者与ArcGIS初学者的编号辅助工具针对mdb、gdb、shp三类常见地理数据格式解决批量要素缺少顺序标识、手工编号低效易错的问题。工具调用Esri ArcGIS API组件可读取并更新数据源属性字段选定数据源、编号字段与起始值即可自动追加序列号也支持按图层或指定要素定向处理适用于数据整理、制图与空间规划等场景。压缩包共74个文件、约7.1MB以32个dll动态链接库和32个配套xml文档为主体另有4个config配置、4个exe程序与2个manifest清单。dll覆盖Carto、Geodatabase、Geometry、Controls、Display等模块xml为各组件接口说明便于对照理解调用关系。已有3992人学习下载。对需要反复处理要素编号的用户工具可直接落地使用省去自行编写脚本的环节组件与文档的对应关系亦可作为了解ArcGIS二次开发接口的参考。1. 一条燃气管网要编八千个号我为什么不再手动拉序号要素一多编号这件事就从体力活变成玄学。某公司做管线普查八千多条管段要在属性表里按巡检路线从头编一遍手动拉序号拉了三天改一个图斑就得全表重来。ARCGIS 自带的字段计算器能填 OBJECTID可那个号跟着要素创建顺序走删几个补几个就断了交出去被打回来两次。这套 ARCGIS 编号小工具干的就是把这件事固定下来给定排序键和分组键稳定地吐出一串连续、可复现的编号顺序编号、分组编号、沿线编号三种口径都能覆盖。适合手里有几千到几十万条要素、又要反复重编的测绘、规划、管网、林业调查类工作——尤其是那种「下周还要再交一版」的项目。2. 编号工具的三根支柱排序键、分组键、输出字段类型写编号脚本之前先把三件事定死按什么排、按什么分组、号写成什么类型。这三件事定错了代码写得再漂亮也是白跑。我见过太多人一上来就CalculateField填!OBJECTID!跑完看着挺顺一删要素全乱回头重来。2.1 OBJECTID 为什么不能直接当编号用OBJECTIDshapefile 里叫 FID是要素类内部维护的行标识不是业务字段。它的行为有几个坑地理数据库里从 1 开始shapefile 里从 0 开始删掉要素后不重排新增要素可能复用被删掉的号导出、投影、合并、裁剪每一步都可能重新分配一遍。也就是说它在整个数据生命周期里是不稳定的。更麻烦的是它跟你想要的顺序没关系。要素创建顺序取决于当初是「先画哪条」跟你现在要的「从北到南」「沿主管线走向」完全是两码事。所以正确做法永远是新建一个业务字段把编号值算出来写进去OBJECTID 只在排序时当一个兜底键用。2.2 排序键怎么定属性、几何、还是混着来排序键决定编号的「阅读顺序」常见三类排序键类型典型字段/几何适用场景注意点属性键图幅号、片区代码、日期、名称已有业务分类顺序由业务定文本排序会按字典序10 2几何键质心 X/Y、四角坐标、包围盒出图好看从上到下从左到右必须在投影坐标系下算混合键先片区、再质心 Y、最后质心 X分幅分片又要求片内有序键越多重复概率越低几何键那条要重点说。地理坐标系WGS84、CGCS2000 经纬度下 X/Y 是度同一纬度带里经度差 0.001 度和纬度差 0.001 度对应的实地距离差好几十个百分点排出来的「从上到下」是歪的。用几何排序前先投影到当地 3 度带或高斯投影坐标单位换成米顺序才对得上眼睛看到的。还有一个隐性要求排序键必须能唯一定序。如果两个要素的排序键完全一样Python 的sort虽然稳定但数据读取顺序本身可能因为压缩、重建索引而变化两次跑出来的结果就不一致。我一般在排序键末尾强制追加OID保证每次结果可复现。2.3 输出字段的类型与格式整型、文本、补零前缀字段类型选错后面全是补丁。对照表如下字段类型取值范围什么时候用备注SHORT短整型-32768 ~ 32767要素数明确不超过 3 万超了直接报溢出LONG长整型约 ±21 亿绝大多数场景的默认选择GDB 推荐TEXT取决于长度要前缀、要补零如A0001shapefile 上限 254DOUBLE浮点基本不用来存编号会出现1.0这种显示补零这件事单独说。业务上经常要PS-0001这种格式很多人先算整数再手动改其实一步到位字段建成 TEXT长度按「前缀 补零位数」算写入时{}{}.format(prefix, str(val).zfill(pad))。pad4就是四位补零超过四位不截断zfill只补不裁。建字段这一步建议封成函数重复跑不会报「字段已存在」# -*- coding: utf-8 -*- import arcpy def ensure_field(fc, name, ftypeLONG, lengthNone, precisionNone): 字段不存在则创建存在则校验类型返回 True 表示新建 exist {f.name.lower(): f for f in arcpy.ListFields(fc)} key name.lower() if key in exist: f exist[key] if f.type.upper() ! ftype.upper(): raise ValueError( 字段 {} 已存在类型为 {}与期望的 {} 不符.format(name, f.type, ftype) ) return False kwargs {in_table: fc, field_name: name, field_type: ftype} if length is not None: kwargs[field_length] length # 仅 TEXT 生效 if precision is not None: kwargs[field_precision] precision # 仅数值型生效 arcpy.management.AddField(**kwargs) return True逻辑说明先ListFields拿现成字段做小写比对避免大小写导致的重复建字段类型不一致直接抛异常而不是硬改因为改字段类型本身要重建字段、丢数据属于该由人来决策的事。length和precision只对特定类型生效按需传传错类型 ArcGIS 会自己忽略。参数说明ftype建议直接用LONG、TEXT、SHORT这种大写常量length对 TEXT 必填比如前缀 3 位加补零 4 位就设 7数值型不要传length传了也没用。注意shapefile 不支持长整型以外的很多类型字段名还有 10 字符上限。要在 shapefile 上跑字段名一律压到 10 字符以内别用中文名。3. 在 ArcGIS Pro 里跑通第一条编号Python 窗口、字段计算器与参数原理清楚之后落到 ArcGIS Pro 里其实就三条路Python 窗口跑脚本、字段计算器填表达式、做成工具箱。日常最顺手的是第一条改一行就能重跑调参不用来回点对话框。3.1 Python 窗口里的最小可运行版本打开 Pro 的「分析 → Python 窗口」先把当前工程的图层路径拿到手。直接写死路径也行但从当前地图取图层更省事换工程不用改代码import arcpy # 从当前工程的活动地图里按名字找图层 aprx arcpy.mp.ArcGISProject(CURRENT) m aprx.activeMap lyr m.listLayers(管线)[0] # 图层名按实际情况改 fc lyr.dataSource # 拿到真实数据源路径 print(fc)拿到fc之后核心编号函数长这样def number_by_order(fc, out_field, order_fieldsNone, start1, step1, prefix, pad0, reverseFalse, whereNone): fc : 要素类路径 out_field : 写入编号的字段名已存在 order_fields: 排序字段列表可为 None start/step : 起始值与步长 prefix/pad : 文本前缀与补零位数pad0 表示不补零 reverse : True 时倒序编号 where : 只对满足条件的要素编号 oid arcpy.Describe(fc).OIDFieldName order list(order_fields or []) # 第一步只读把 OID 和排序键一次性取出来 rows [] with arcpy.da.SearchCursor(fc, [OID] order, where_clausewhere) as cur: for r in cur: # None 单独处理避免排序时 TypeError keys tuple((1, ) if v is None else (0, v) for v in r[1:]) rows.append((r[0], keys)) del cur # 末尾追加 OID 兜底保证排序结果唯一且可复现 rows.sort(keylambda x: x[1] ((0, x[0]),), reversereverse) # 第二步算值编号只跟排序后的位置有关 result {} for i, (o, _) in enumerate(rows): val start i * step if pad 0: result[o] {}{}.format(prefix, str(val).zfill(pad)) else: result[o] val # 第三步写回单独开一个游标避免读写混在一个游标里 n 0 with arcpy.da.UpdateCursor(fc, [OID, out_field], where_clausewhere) as cur: for o, _ in cur: cur.updateRow([o, result[o]]) n 1 del cur return n逐段说清楚。第一步的SearchCursor只读不写把 OID 和排序键全部拉到内存里排序这是关键——如果边读边写游标遍历顺序和写入顺序会互相干扰ArcGIS 官方文档里也明确不建议在同一个游标里做重排序写入。keys那里把None转成(1, )是为了让空值排在所有有效值后面而不是让整个sort抛TypeError。第二步的rows.sort是纯 Python 排序键是「排序字段元组 OID」reverse直接控制升降序。这里没有用 SQL 的ORDER BY原因是ORDER BY在不同数据源上对中文、空值的处理不一致交出去的数据要跨环境复现Python 端自己排更保险。第三步单独开UpdateCursor。where_clause保持一致确保只更新该更新的行不会把未参与编号的要素刷成None。参数上最容易搞错的是pad。如果你要A0001这种格式out_field必须建成 TEXTpad4prefixA如果要纯数字pad0字段建 LONG。混用会直接报类型不匹配。3.2 字段计算器快车道什么时候能用什么时候别碰不是所有编号都值得写脚本。如果只是「全表从 1 到 N顺序无所谓」字段计算器一行就够了# 前提out_field 已存在类型为 LONG arcpy.management.CalculateField( in_table管线, fieldNH, expression!OBJECTID!, # 直接映射内部 ID expression_typePYTHON3 )但这行代码有两个前提一是能接受编号跟内部 ID 绑定删要素后一定断号二是要素数不超过字段类型上限。任何一个不满足就得回到脚本。字段计算器还能用 Code Block 做条件编号比如「按类型分段」arcpy.management.CalculateField( in_table管线, fieldNH, expressionget_num(!TYPE!, !OBJECTID!), expression_typePYTHON3, code_blockoffset {给水: 100000, 排水: 200000, 燃气: 300000} def get_num(t, oid): return offset.get(t, 0) oid )逻辑说明code_block里定义一个字典做段号偏移同一张表里不同专业各占一个号段避免跨专业混编。适合「一眼能从号看出专业」的场景。代价是OBJECTID不保证段内连续删了要素号就断。所以这个写法我只在初版出图时用正式交付前一定换成number_by_order并且把排序键设成业务字段。这种「先快后稳」的切换节奏能省下不少现场等待时间。3.3 参数怎么调起始值、步长、方向、范围四个参数逐个说清楚什么时候改start一般设 1除非要预留号段。比如新老数据合并老数据编到 3200新数据从 3201 起start3201一步到位。step平时都是 1。设成 10 或 100 是给未来插号留空间管线、道路这类会不断新增的要素可以考虑代价是号看上去稀疏。reverse控制方向。按质心 Y 排序时默认是「从下到上」想要「从上到下」就reverseTrue。这个参数建议写在调用处显式给别藏在默认值里否则下一次跑出来的顺序跟上次不一样排查半天。where是范围控制。只给某个片区的要素编号where_clausePIANQU A未命中的要素字段保持原值不动。这一点很重要UpdateCursor带where_clause时只遍历命中的行其他行不会被刷成空。提示第一次跑新数据先where限定到一个片区试跑看一眼属性表顺序对不对再放开全表。全量跑一次几万条要素通常十几秒但跑错了回滚的成本远不止十几秒。4. 分组与空间排序按片区编、按图幅编、沿线路编怎么落地顺序编号只是起点。真实项目里更常见的是「每个片区内部从 1 开始」「每幅图内部从 1 开始」「沿主管线走向从上游编到下游」。这三种需求分别对应分组编号、几何排序编号和沿线编号写法差别不小。4.1 分组编号字典计数与 ORDER BY 的两种写法分组编号的本质是「每个组各有一个独立计数器」。最直观的写法是一个字典键是组值值是当前计数def number_by_group(fc, group_field, out_field, order_fieldsNone, start1, step1, prefix, pad0): 每个 group_field 值内部独立编号编号格式为 组值-序号 order list(order_fields or []) rows [] with arcpy.da.SearchCursor(fc, [OID, group_field] order) as cur: for r in cur: g r[1] if r[1] is not None else __NULL__ keys tuple((1, ) if v is None else (0, v) for v in r[2:]) rows.append((r[0], g, keys)) del cur # 先按组名排再按组内排序键排 rows.sort(keylambda x: ((0, x[1]), x[2], (0, x[0]))) counters {} result {} for o, g, _ in rows: idx counters.get(g, 0) val start idx * step counters[g] idx 1 if pad 0: result[o] {}{}-{}.format(prefix, g, str(val).zfill(pad)) else: result[o] val with arcpy.da.UpdateCursor(fc, [OID, out_field]) as cur: for o, _ in cur: cur.updateRow([o, result[o]]) del cur return len(result)逻辑说明rows.sort的键是三段式——先组名再组内排序键最后 OID 兜底。这样字典计数器只要顺序扫一遍就能保证「组内连续、组间独立」。counters用dict.get(g, 0)而不是defaultdict是为了让逻辑显式一点出问题时好打断点。参数说明group_field建议用短整型或文本别用浮点浮点做分组键会因为精度问题分出莫名其妙的组。pad和prefix同上输出格式是前缀 组值 - 补零序号比如PQ-A-0007。这个写法的边界是数据量。十万条以内Python 排序毫无压力上百万条内存和UpdateCursor的写入时间会明显上去那时候该考虑用arcpy.management.Sort先落一个排好序的临时要素类再顺序扫一遍写号。4.2 按几何位置排序质心、四角、还是包围盒几何排序没有统一答案取决于出图时希望人眼怎么读。三种常见口径按质心 Y 降序、X 升序就是「从上到下、从左到右」的蛇形读法最符合纸质图的阅读习惯def number_by_centroid(fc, out_field, y_descTrue, start1): rows [] with arcpy.da.SearchCursor(fc, [OID, SHAPEXY]) as cur: for oid, (x, y) in cur: # 先按 Y 排再按 X 排Y 降序靠取负实现 rows.append((oid, (-y if y_desc else y, x))) del cur rows.sort(keylambda r: (r[1][0], r[1][1], r[0])) with arcpy.da.UpdateCursor(fc, [OID, out_field]) as cur: for i, (oid, _) in enumerate(rows): pass del cur # 写回略同 number_by_order 的第三步SHAPEXY拿的是几何的质心对线是中点对点是本身比SHAPE快一个数量级十万条以内几乎无感。取负来排序是省掉reverse参数的老写法可读性一般我一般写成key(-y, x)并在注释里标清楚。按四角排序是另一种口径面要素按左上角坐标排适合地籍、地块这类要求「先上后下、先左后右」严格对齐的场景。左上角取extent.XMax, extent.YMax但extent每次都要算包围盒比质心慢。要素数超过二十万时我通常先算一遍质心写进临时字段再拿字段排序比反复访问几何快得多。注意几何排序前务必确认数据框和要素类用的是投影坐标系。经纬度下SHAPEXY返回的是度排序结果跟实地位置对不上这种错在出图时才被发现返工量很大。4.3 沿线路走向编号measureOnLine 的用法与边界沿管线、河道、道路走向编号是管网和水利项目的刚需。核心思路是找一条参考线主管线把每个要素的质心投影到参考线上算出沿线累计距离按距离排序。ArcPy 的Polyline.measureOnLine(point, use_percentage)就是干这个的返回点投影到线上的位置。use_percentageFalse时返回实际长度米True时返回 0~1 的比例def number_along_route(fc, route_fc, out_field, start1, reverseFalse): 按要素质心沿参考线 route_fc 的累计距离编号 # 参考线必须是单条线要素多段线会自动拼成一条 with arcpy.da.SearchCursor(route_fc, [SHAPE]) as cur: route next(iter(cur))[0] del cur rows [] with arcpy.da.SearchCursor(fc, [OID, SHAPE]) as cur: for oid, geom in cur: pt geom.centroid # 要素离参考线太远时 measureOnLine 会投影到最近端点 d route.measureOnLine(pt, False) rows.append((oid, d)) del cur rows.sort(keylambda r: (r[1], r[0]), reversereverse) with arcpy.da.UpdateCursor(fc, [OID, out_field]) as cur: for o, _ in cur: pass del cur # 写回同前第 i 个要素写 start i逻辑说明参考线读成单个Polyline对象而不是逐行处理是因为measureOnLine是几何对象的方法需要完整的线几何。多个要素同时投影到同一条线上才能保证距离值的相对关系一致。参数说明reverseTrue是「从下游往上游编」管网项目里经常要按水流反方向出图这个参数一定要显式给。geom.centroid对面要素是几何中心对线要素是线的中点如果你的要素是分段的管段中点就是这一段的位置合理但如果管段很长、弯折剧烈中点可能偏离主走向这时候改用端点geom.firstPoint更准。边界也很明确measureOnLine是把点垂直投影到线上如果要素离参考线很远比如支线离主管线几百米投影会落到最近的那个端点上多个要素会算出同一个距离排序退化成按 OID。所以这条路线只适合「要素大致沿参考线分布」的场景分叉多的管网建议先按支线分组再在组内沿线排序。5. 排查清单编号重复、顺序乱跳、字段被锁的五种现场脚本跑通只是开始真实现场的问题基本都出在数据状态和运行环境上。下面五条是我这两年改得最多的按出现频率排。第一条编号重复而且重复的位置每次还不一样。现象是跑完属性表里出现两个 17、三个 42重跑一遍重复的位置变了。原因几乎都是排序键不唯一且没有 OID 兜底——sort本身稳定但数据读取顺序可能因为压缩、重建索引变化两次结果就不一致。解决办法是在排序键末尾强制追加OID让每个元素的排序键唯一这样无论读多少次顺序都锁死。这条改完之后重复问题基本绝迹。第二条明明按数字排结果 10 跑到 2 前面去了。现象是按图幅号排序1, 10, 11, 2, 20这么出来。原因是排序字段是 TEXT 类型Python 按字典序排字符1比2小所以10排在2前面。解决有两个方向一是排序前把文本转成intkeylambda r: int(r[1])前提是字段里没有非数字内容二是编号字段本身用整数只在展示时补零。我一般选后者字段类型在建模阶段就定死后面不用反复转换。第三条字段被锁UpdateCursor报「无法获取锁」。现象是脚本跑到写回那一步挂掉提示锁冲突或只读。常见原因有三个属性表窗口开着、编辑会话没保存、要素类被别的程序占用比如同步工具、其他 Pro 工程。逐一排查的顺序是先关属性表再看「编辑」选项卡里有没有未保存的编辑最后看系统里有没有残留的ArcGISPro.exe进程。这三步走完基本都能解。血泪经验是跑批量编号前先保存工程、关掉所有属性表能省掉一大半莫名其妙的失败。第四条脚本中途崩了字段写了一半前功尽弃。现象是跑到第 3000 条报错退出前 3000 条已经写进去了后面的还是旧值数据处于半成品状态。这个问题不是脚本 bug是没有事务保护。解决方式有两种简单的是先arcpy.management.CopyFeatures复制一份备份崩了直接换回来稳妥的是用arcpy.da.Editor开编辑会话把整个编号过程包在startEditing / stopEditing(True)里异常时stopEditing(False)回滚。第二种写法多一点但真正跑大批量时值得。第五条分组字段有空值整组编号从 0 开始或者整组消失。现象是按片区编号PIANQU为空的那批要素要么编号是 0 或者负数要么根本没被编上。原因是None在字典里做键、在排序里做比较都会出问题。解决是在分组前统一把空值映射成一个显式标签比如__NULL__编完再人工确认这批要素该归到哪个片区。这个坑的隐蔽之处在于空值要素往往不多属性表翻到后面才看得见出图时才发现有几条没号。提示这五条里第一、四条的代价最高——一个是结果不可信一个是数据要重来。跑正式数据前拿十条要素的样本全流程走一遍能提前暴露 80% 的问题。6. 把脚本收进 .pyt 工具箱批量跑、带校验、下次直接拖进去脚本在 Python 窗口里跑改参数得回来改代码交给同事还得附一段说明。真正省事的做法是包成.pyt工具箱参数在对话框里填双击就能跑还能扔进模型构建器里串批处理。最小骨架长这样一个工具类对应一个##级功能# -*- coding: utf-8 -*- import arcpy class Toolbox(object): def __init__(self): self.label 编号工具箱 self.alias numbering self.tools [SeqNumberTool] class SeqNumberTool(object): def __init__(self): self.label 顺序编号 self.description 按指定字段排序后连续编号支持前缀与补零 def getParameterInfo(self): p_fc arcpy.Parameter(namein_fc, displayName输入要素类, datatypeDEFeatureClass, parameterTypeRequired, directionInput) p_field arcpy.Parameter(nameout_field, displayName编号字段, datatypeField, parameterTypeRequired, directionInput) p_field.parameterDependencies [in_fc] # 依赖输入表自动列出字段 p_order arcpy.Parameter(nameorder_field, displayName排序字段, datatypeField, parameterTypeOptional, directionInput, multiValueTrue) p_order.parameterDependencies [in_fc] p_start arcpy.Parameter(namestart, displayName起始值, datatypeGPLong, parameterTypeOptional, directionInput) p_start.value 1 return [p_fc, p_field, p_order, p_start] def execute(self, parameters, messages): fc parameters[0].valueAsText out_field parameters[1].valueAsText order parameters[2].values or [] start parameters[3].value or 1 n number_by_order(fc, out_field, [o for o in order], startstart) messages.addMessage(完成编号{} 条.format(n)) # 跑完立刻校验问题当场暴露 check_numbering(fc, out_field)parameterDependencies是关键它让「编号字段」和「排序字段」两个下拉框跟着输入要素类联动用户不用手敲字段名。multiValueTrue允许选多个排序字段顺序就是对话框里的上下顺序。配套的校验函数建议每个项目都跑一遍比肉眼翻属性表靠谱得多def check_numbering(fc, field): 检查空值、重复、断号结果直接打印 vals [r[0] for r in arcpy.da.SearchCursor(fc, [field])] nonnull [v for v in vals if v is not None] dup len(nonnull) - len(set(nonnull)) gaps sum(1 for a, b in zip(sorted(nonnull), sorted(nonnull)[1:]) if b - a ! 1) print(总数 {}空值 {}重复 {}断号 {}.format( len(vals), len(vals) - len(nonnull), dup, gaps)) return dup, gaps批量处理多个要素类时用arcpy.ListFeatureClasses()拿到工作空间下所有要素类循环调number_by_order把out_field统一成一个名字。这一步的前提是所有要素类的字段结构一致不一致的先补字段再跑。我踩过最深的一个坑是校验函数只查了重复没查断号交出去的数据号是连续的但中间缺了 37 号被退回来重做了一遍。从那以后我每次交数据前都强制把check_numbering跑一遍三个数字都归零才敢发出去。这套习惯帮我省掉的返工时间远比写脚本的时间多。希望帮到你。本文还有配套的精品资源点击获取
