简介ARCGIS编号小工具面向地理信息系统数据处理人员用于为mdb、gdb和shp三类常见地理数据批量添加顺序编号缓解大量要素标识时的手工操作压力可服务于数据管理、地图制作与空间规划等场景。压缩包共74个文件约7.1MB内含32个动态链接库、32个说明文档、4个可执行程序、4个配置文件和2个清单文件动态链接库承载ArcGIS相关功能说明文档辅助理解组件用途可执行程序与配置文件负责工具运行和参数设置。已有3992人学习。工具可先选择数据源再设定编号字段与起始值随后自动写入序列号既支持批量处理也可针对特定图层或要素执行借助完整依赖读者能了解地理数据库读写、几何对象处理和地图显示等环节的调用方式并据此完成编号操作或研究工具结构。1. 图斑编号这件小事两千条数据错位之后我才想明白顺序给要素编号听上去是属性表里拉一列的事真做起来耗掉半天。前段时间接了个图斑整理项目两千多个面要求按行政村分组、每组从 1 开始连续编号、不足四位补零还要按从上到下的空间顺序排。同事先在 Excel 里拉序号再连接回 ArcGIS结果排序字段里躺着一堆空值编号整体错位返工两遍。这类活的共性很清楚顺序、分组、格式化、写回四步里只要有一处靠手工就没有可重复性。这个 ARCGIS 编号小工具把四步压成一次执行输入要素类、排序表达式和分组字段编号直接落进属性表。它适合经常和图斑、宗地、管线、监测点打交道的从业者也适合刚上手 arcpy、想找一个能跑通的最小样例的人。2. 拆开编号工具四个动作和一张参数表2.1 一次编号本质是四个动作串起来很多人对编号的第一反应是“写个循环”但真正决定结果对不对的是循环之前的准备工作。我把整个流程拆成四段确定顺序、划定分组边界、生成序列、写回字段。四段是有严格先后依赖的顺序必须在写回之前一次性确定因为写回动作会改变行状态边写边排只会越排越乱。第一段“确定顺序”要回答的是“谁排第一”。可选依据很多按面积从大到小、按地类编码、按行政村代码、按质心 Y 坐标从上到下。这个依据不能是 FID 或 OBJECTID虽然它们看起来天然有序但编辑、导出、合并、裁剪任意一次操作都会让 OID 重排编号会跟着漂移。血泪经验是临时用 OID 排一次可以但一定要把结果落到实体字段里别指望下次打开还是那个顺序。第二段“划定分组边界”本质是一次分组聚合。分组字段取同一值的连续行计数器归零重新从起始值开始。这里有个容易被忽略的点分组和排序必须同时生效SQL 层面就是ORDER BY 分组字段, 排序键如果只排排序键不排分组字段同一组的数据会被拆散在游标的各个位置计数器会反复归零出来就是 1、2、1、3、1 这种鬼样子。第三段“生成序列”涉及起始值、步长、补零位数、前后缀四个参数。补零是最容易被低估的一个zfill(4)把 1 变成0001看着简单但如果字段类型是整型写进去的0001会被数据库转成 1你打开属性表一看还是没补零——不是脚本错了是字段类型不接受。第四段“写回”用arcpy.da.UpdateCursor逐行更新比CalculateField稳的地方在于它能在同一个游标里带上排序和分组信息不用来回切上下文。顺带说一句CalculateField也能做递增编号写个全局变量在表达式里自增就行但那是玄学写法执行顺序不保证重跑一次结果可能变而且中途报错之后全局变量状态是脏的得重启整个会话。我一般只在一次性、几十条的小数据上这么干。2.2 为什么不直接上字段计算器或者 Excel把几种常见做法摆在一起对比选型理由就很清楚了。做法上手成本能否指定多级排序重复执行主要风险字段计算器Python 解析器低否顺序由 OID 决定每次全表重算编辑后编号漂移报错后状态不可控Excel 拉序列 连接回表低手动排序易错需要重做一遍连接键类型不一致、FID 不稳定、文本被转科学计数模型构建器 迭代器中有限排序要靠前置工具可以参数写死在模型里出错定位困难脚本工具这套思路中任意字段、多级、升降序可以需要处理字段类型、文件锁、空值Excel 那条路最坑的地方不是排序是连接。文本型编号0001在 Excel 里会被吞成 1导出 CSV 再连接回来0001和1对不上你会在属性表里看到一大堆空值然后开始怀疑人生。如果非要用表格当中转务必把编号列先设成文本格式或者干脆在前面拼一个字母前缀。字段计算器的优势是零门槛劣势是不可复现。同一个表达式今天跑出来的编号和下周跑出来的编号可能不一致因为中间有人动过数据、加过字段、改过排序。项目交付前最怕的就是这种“上次是对的这次不对”查起来没有抓手。模型构建器的定位其实是“把流程固定下来”思路是对的但它不擅长表达“按分组重置计数”这种带状态的逻辑。迭代器选中分组、内部再迭代要素能拼出来但每加一个分组就多一层嵌套调试成本直线上升。所以最终的形态还是脚本工具核心逻辑一个函数参数从工具箱传进来既能一键跑也能在 Python 窗口里逐步调试。2.3 参数设计一个排序表达式吃掉所有排序需求工具箱的参数不能无限加动态增减参数在 ArcGIS 里做起来很别扭。我的做法是把所有排序需求压缩成一个字符串表达式格式是字段 方向; 字段 方向分号分隔多级方向写ASC或DESC省略则按升序。比如CUNBM ASC; MJ DESC表示先按行政村代码升序、同组内按面积降序。参数数据类型默认值说明输入要素要素图层 / 要素类无支持 GDB 要素类、Shapefile、SDE 图层编号字段字段无文本型或长整型短整型有 32767 上限排序表达式字符串OID ASC分号分隔多级方向 ASC/DESC分组字段字段可选空为空表示全域连续编号起始值长整型1每组重新从这里开始步长长整型1一般不用改跨图层接续时会用到补零位数长整型00 表示不补零仅对文本字段有效前缀 / 后缀字符串空仅对文本字段有效这个设计牺牲了一点直观性——用户得记住表达式格式但换来的是无上限的排序能力。常见做法是在脚本工具的参数说明里写两行示例再加一层解析容错字段名写成MJ、mj、MJ DESC、MJ desc都能认遇到不存在的字段直接抛一条中文错误别让人对着ERROR 000732猜。还有一个细节值得提前定编号字段的类型检查放在最前面做。文本字段走字符串格式化数值字段直接写整数如果用户给数值字段配了前缀或补零脚本主动报错退出而不是悄悄写个1进去让人以为成功了。这种“早失败”比事后核对省太多时间。3. 把编号跑通从命令行最小闭环到工具箱脚本3.1 先确认 arcpy 能被 import不要一上来就写业务代码。先用一条命令确认解释器找得到 arcpy这一步能把一半的环境问题挡在门外。ArcGIS Pro 自带的是 Python 3 环境ArcMap 10.x 带的是 Python 2.7两者的语法和编码处理不一样脚本要分开写或者做兼容。# ArcGIS Pro直接调用安装目录下自带的 Python 环境路径以实际安装为准 C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe -c import arcpy; print(arcpy.GetInstallInfo()[Version]) # 也可以先激活环境再进交互式窗口调试 propy.bat python -c import arcpy; print(arcpy.env.workspace)# ArcMap 10.x自带 Python 2.7注意 print 是语句不是函数 C:\Python27\ArcGIS10.8\python.exe -c import arcpy; print arcpy.GetInstallInfo()[Version]第一段命令的逻辑很简单绕开系统 PATH 里可能存在的另一个 Python直接指定 ArcGIS 自带的解释器避免出现“装了 arcpy 但 import 报错”的假故障。第二段是给老项目用的Python 2.7 里print后面不能跟括号包多个参数unicode和str的区分也要自己处理中文路径建议统一用u...加.decode(utf-8)。参数上唯一要改的是解释器路径把它换成你机器上的实际安装位置即可。如果这一步就报ImportError先查两件事是不是用了系统 Python是不是 32 位和 64 位混装。ArcMap 后台的 Python 是 32 位Pro 是 64 位同一个脚本在两边跑第三方库基本不能共用。3.2 核心函数排序、分组、补零、写回下面这段是整套工具的主干逻辑全部压在一个函数里方便复制到脚本工具里调用也方便单独调试。写法兼容 Python 2.7 和 3.x没用 f-string。# -*- coding: utf-8 -*- import arcpy def auto_number(fc, num_field, sort_keyOID ASC, group_fieldNone, start1, step1, width0, prefix, suffix, whereNone): 按排序表达式给要素写编号。 fc : 要素类路径或图层名 num_field : 编号字段名文本或长整型 sort_key : 字段 方向; 字段 方向方向省略为 ASC group_field : 分组字段None 表示全域连续编号 start/step : 起始值与步长 width : 补零位数0 表示不补 prefix/suffix: 前后缀仅文本字段有效 where : 可选筛选条件只对满足条件的要素编号 # 1) 拿真实字段名和类型Shapefile 会截断字段名大小写也可能和看到的不一致 desc arcpy.Describe(fc) oid_name desc.OIDFieldName field_map dict((f.name.lower(), f) for f in arcpy.ListFields(fc)) if num_field.lower() not in field_map: raise ValueError(u编号字段不存在{0}.format(num_field)) num_fld field_map[num_field.lower()] is_text (num_fld.type String) if not is_text and (width or prefix or suffix): raise ValueError(u数值型字段不支持补零和前后缀请改用文本字段) # 2) 解析排序表达式拼成 ORDER BY 需要的字段列表 fields [num_fld.name] order_sql [] for seg in (sort_key or ).split(;): seg seg.strip() if not seg: continue parts seg.split() fname parts[0] if fname.lower() in (oid, objectid, fid): fname oid_name # 统一映射到真实 OID 字段名 direction DESC if len(parts) 1 and parts[1][0].upper() D else ASC real field_map.get(fname.lower()) if real is None: raise ValueError(u排序字段不存在{0}.format(fname)) order_sql.append({0} {1}.format(real.name, direction)) fields.append(real.name) if group_field: grp field_map.get(group_field.lower()) if grp is None: raise ValueError(u分组字段不存在{0}.format(group_field)) fields.append(grp.name) sql_clause (None, ORDER BY , .join(order_sql)) if order_sql else (None, None) # 3) 单次游标完成遍历与写回分组变化时计数器归零 seq start last_group object() # 哨兵对象避免和真实的 None / 空串撞车 with arcpy.da.UpdateCursor(fc, fields, where, sql_clausesql_clause) as cur: for row in cur: if group_field: g row[-1] if g ! last_group: seq start last_group g if is_text: body str(seq).zfill(width) if width else str(seq) row[0] u{0}{1}{2}.format(prefix, body, suffix) else: row[0] seq cur.updateRow(row) seq step if __name__ __main__: arcpy.env.workspace rD:\work\demo.gdb auto_number(tb_parcel, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, start1, step1, width4, prefix) print(done)逐段说一下逻辑。开头用Describe拿 OID 真实名字是因为 OID 在不同数据源里可能叫OBJECTID、FID、OBJECTID_1用户写OID的时候必须能映射过去否则排序找不到字段直接报错。字段字典用小写做键是为了兼容 Shapefile 那种大小写不敏感但显示不一致的情况。中间解析排序表达式的循环做了三件事拆字段和方向、把OID别名换成真实字段名、校验字段是否存在。校验这一步很值得写因为arcpy.da.UpdateCursor传了不存在的字段只会给一句RuntimeError: ERROR 999999对新手极不友好自己抛一条中文提示能省半小时。最后一段是核心。sql_clause的第二个元素就是ORDER BY子句游标会按这个顺序吐行这是整套方案能成立的前提。last_group用object()当哨兵是因为分组字段可能真的是空字符串或None用None做初始值会让第一组被误判成“换了组”序号白白重置一次。写回时按字段类型分两支文本走zfill加前后缀数值直接写整数。参数上最常调的是sort_key和group_fieldstep基本只在“接着上次的号往下编”时会变成 1 以外的东西。where参数留了个口子只对满足条件的要素编号比如只编现状图斑、不编历史遗留。3.3 挂进工具箱参数怎么配脚本能跑通之后下一步是把它变成别人也能点的工具。ArcGIS Pro 里用“新建脚本工具”ArcMap 里用“添加脚本”本质一样把 .py 文件挂到工具箱然后按顺序定义参数。参数顺序必须和脚本里GetParameterAsText的索引严格对应错一位就是连环错。序号显示名数据类型方向过滤 / 默认值0输入要素要素图层输入无1编号字段字段输入获得自输入要素过滤文本、长整型2排序表达式字符串输入默认OID ASC3分组字段字段输入获得自输入要素可选4起始值长整型输入默认 15步长长整型输入默认 16补零位数长整型输入默认 07前缀字符串输入可选8后缀字符串输入可选入口代码只做参数转接不放业务逻辑这样在 Python 窗口里调试时还能直接调auto_number# -*- coding: utf-8 -*- import arcpy from number_tool import auto_number # 同目录下的核心模块 if __name__ __main__: p arcpy.GetParameterAsText try: auto_number( fcp(0), num_fieldp(1), sort_keyp(2) or OID ASC, group_fieldp(3) or None, # 工具箱里留空会传空字符串要转成 None startint(p(4) or 1), stepint(p(5) or 1), widthint(p(6) or 0), prefixp(7) or , suffixp(8) or , ) arcpy.AddMessage(u编号完成) except Exception as e: arcpy.AddError(u编号失败{0}.format(e)) raise参数转接里最容易被忽略的是p(3) or None这一句。工具箱里没填的参数传进来是空字符串如果不转成Noneif group_field:判断为真脚本会去找一个名字为空的字段然后报“分组字段不存在”。这个坑我调了不止一次后来干脆在所有可选参数上都加or None或or 默认值。AddError之后要raise否则工具在结果窗口里显示“已完成”但实际什么都没做用户以为成功了问题留到验收才爆出来。工具一旦给别人用报错信息就是唯一的沟通渠道这条不能省。4. 四种编号场景连续、分组、空间、编码型怎么配4.1 全域连续编号参数最少坑最集中最简单的用法就是把排序键设成 OID不分组从 1 编到底。这种场景通常出现在“给要素一个唯一序号”的需求上比如管线点、监测点、抽样点。# 不分组按 OID 保序四位补零前缀加一个字母 auto_number(rD:\work\demo.gdb\pt_monitor, BH, sort_keyOID ASC, width4, prefixJC)这里的逻辑是游标按 OID 顺序吐行编号从 1 递增每行格式化成JC0001这种形式写回BH字段。参数上唯一要留意的是width4如果需要编号过万位数要给到 5否则从 10000 开始zfill不再起作用JC10000和JC0001长度不一致后续按字符串排序会错。排序键换成业务字段也很常见比如按面积降序MJ DESC让大图斑排在前面。但这时必须补一个二级排序键MJ DESC; OID ASC因为面积重复的行之间顺序不稳定每次跑出来的编号可能不一样同一份数据两次执行结果不同验收时没法解释。4.2 分组重编号先洗干净分组字段再谈编号按行政村、按图层、按权属单位分组重新编号是最常见的需求也是最容易翻车的。翻车点不在编号逻辑在分组字段本身有的值带前后空格有的是全角空格有的None和空串混用结果就是同一组被拆成三组编号出来是 1、2、1、3。先做一遍清洗再编号# 清洗分组字段去空格含全角、None 转空串避免同组被拆散 fc rD:\work\demo.gdb\tb_parcel with arcpy.da.UpdateCursor(fc, [CUNBM]) as cur: for row in cur: v row[0] or u row[0] v.replace(u\u3000, u).strip() cur.updateRow(row) # 清洗完再编号先按分组字段排组内按面积降序 auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, width4)第一步把CUNBM里的全角空格和首尾空白干掉None统一成空串这一步必须单独做不能塞进编号函数里因为清洗是破坏性操作出问题要能单独回滚。第二步的排序表达式把分组字段放在第一位这是分组编号的硬性要求ORDER BY里分组字段必须在前编号才能保证同组数据连续出现计数器归零的时机才对得上。参数上的对应关系是group_field决定什么时候归零sort_key的第一段必须和group_field一致方向也一致否则游标吐出来的行序和分组边界对不上会出现“组 A 编到 5突然跳到组 B又跳回组 A 从 1 开始”的情况。这个约束我在函数文档里写了但还是有人踩后来干脆加了一条检查如果sort_key的第一段不是group_field直接抛错。4.3 空间顺序编号从上到下、从左到右规划类项目经常要求“按从上到下的顺序编号”本质是按质心 Y 坐标降序、X 坐标升序。SHAPEXY能直接取到要素的质心坐标但要注意坐标系地理坐标系下经纬度排序会歪因为经度一度的实际距离随纬度变化同一行要素的 Y 值差异会被放大。正确做法是先投影到投影坐标系再取坐标。# 1) 加两个双精度字段存质心坐标已存在则跳过 fc rD:\work\demo.gdb\tb_parcel exist set(f.name.lower() for f in arcpy.ListFields(fc)) add [[n, DOUBLE] for n in (CX, CY) if n.lower() not in exist] if add: arcpy.management.AddFields(fc, add) # 2) 写入质心坐标SHAPEXY 返回 (x, y) with arcpy.da.UpdateCursor(fc, [SHAPEXY, CX, CY]) as cur: for row in cur: row[1], row[2] row[0] cur.updateRow(row) # 3) 先按 Y 降序从上到下同一行内按 X 升序从左到右 auto_number(fc, BH, CY DESC; CX ASC, group_fieldCUNBM, width4)第一步用ListFields判断字段是否已存在避免重复执行时报“字段已存在”。第二步把SHAPEXY元组拆成两个数值字段这一步是必须的因为ORDER BY不支持直接对几何排序只能对属性字段排。第三步的排序表达式里CY DESC控制从上到下CX ASC控制同一水平带内从左到右。参数上有个经验值如果要素在空间上呈明显的条带状分布比如沿道路两侧的地块纯按质心排序会出现“上一行还没编完就跳到下一行”的锯齿这时可以在排序表达式里加一个分级字段比如ROW_NO ASC; CX ASC先用人工划定的行号分组再在行内按 X 排。工具本身不生成行号行号通常靠“按位置选择”或人工辅助图层得到。4.4 编码型编号前缀、年份、流水号的拼接政务和规划类项目里编号往往不是纯数字而是XZ-2024-0001这种带业务含义的编码。这类需求的关键不在拼接在字段长度和字符集。# 文本字段长度要提前算够前缀 3 分隔符 1 年份 4 分隔符 1 流水 4 13 # 实际建字段时留余量直接给 32 或 50 prefix uXZ-2024- auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, start1, step1, width4, prefixprefix)逻辑上就是把前缀交给prefix参数流水号由width补零拼出来是XZ-2024-0001。参数上要注意两点一是字段长度GDB 的文本字段有长度上限长度不够会被静默截断XZ-2024-0001变成XZ-2024-00打开属性表才发现二是字符集前缀里有中文的话Shapefile 的编码问题会让人抓狂建议一律用 GDB且前缀只用字母数字和连字符。如果编码需要按年份分组重置做法是把年份字段也放进分组字段的组合里但单字段分组满足不了。常见做法是先把年份和行政村拼成一个辅助字段GRP YEAR || CUNBM用这个字段当group_field编号完再删掉辅助字段。多一步操作但比改脚本逻辑稳妥。5. 避坑与排查五个现场记录和它们背后的原因5.1 字段与数据层面的两个坑现象一脚本报“字段不存在”但属性表里明明有这个字段。原因基本是两个Shapefile 的字段名上限是 10 个字符你建的是BIANHAO实际存的是BIANHAO刚好 7 位没事但如果叫PARCL_CODE落到 shp 里会变成PARCL_COD另一个是大小写ListFields返回的名字可能和你输入的完全一致也可能因为数据源不同有差异。解决方式是所有字段查找都走小写字典并且在报错信息里把可用字段列表打出来一眼就能看出真名是什么。现象二编号跑到一半报溢出或者大数字变成了负数。原因是编号字段用了短整型SmallInteger上限 32767超过就溢出。这个坑在测试数据上永远不会出现因为测试数据只有几十条一上生产就炸。解决办法是建字段时就选长整型Integer上限约 21 亿或者干脆用文本字段。工具里加一条前置检查字段类型是 SmallInteger 且要素数可能超过 3 万时直接提示换字段。5.2 排序与分组层面的两个坑现象三分组编号出现 1、2、1、3 这种乱跳。原因在 4.2 节提过分组字段里有隐藏的空白字符或全半角差异同一组被当成多组。还有一种是None和空串混用游标里None ! 每次遇到都会触发归零。解决办法是先跑一遍清洗再用“汇总统计”按分组字段计数看看组的数量和业务上预期的是否一致对不上就别急着编号。现象四同一份数据跑两次编号结果不一样。原因是排序键有重复值重复行之间的相对顺序由数据库决定不保证稳定。这个问题最阴因为第一次跑完看着是对的第二次跑完也“看着是对的”只有把两次结果做对比才会发现差异。解决办法是排序表达式最后永远补一个OID ASC兜底让排序在全序上唯一。我现在的习惯是任何排序表达式不写 OID 都不许提交。5.3 运行环境层面的一个坑现象五报“无法获取锁”或者“要素类正在被使用”。原因有三种编辑会话没停、属性表窗口开着、文件被别的程序占着比如同步盘、杀毒软件扫描。ArcGIS 对 GDB 的锁比较严格一个打开的属性表就足以让UpdateCursor拿不到写锁。解决办法是按顺序排查先停编辑并保存再关掉所有属性表窗口和目录面板还不行就重启 Pro最后才怀疑是不是同步盘在后台抓文件。生产数据建议放在本地盘或网络盘的非同步目录里跑别放在会自动上传的文件夹下。顺带提一个相关的Shapefile 在编号过程中会生成.lock文件异常中断后锁文件可能残留导致后续操作一直失败。这时不要删.lock硬来正确做法是关掉所有 ArcGIS 进程再重新打开让系统自己释放。6. 编号之后的验收重复值、断号与增量维护的检查清单编号跑完不等于交付完成验收这一步能挡住大部分返工。我一般做三件事。第一件是查重复和空值。用Counter快速扫一遍重复的直接列出来from collections import Counter fc rD:\work\demo.gdb\tb_parcel with arcpy.da.SearchCursor(fc, [BH, CUNBM]) as cur: rows [(r[0], r[1]) for r in cur] dup [k for k, v in Counter(x[0] for x in rows).items() if v 1] empty [i for i, x in enumerate(rows) if not x[0]] print(u重复编号 {0} 个空编号 {1} 个.format(len(dup), len(empty)))第二件是查断号。按分组统计最小值和最大值和组内要素数对比max - min 1 count说明连续不等就是中间跳号了。这个检查用“汇总统计”工具做也行用几行脚本做也行关键是每组都查别只抽查一组。第三件是接续编号。新增要素要接着已有的最大号往下编不能从头开始。做法是先扫一遍现有编号取最大值把它加一作为startdef next_start(fc, num_field): 取现有编号的最大流水值 1兼容带前缀的编码型编号 mx 0 with arcpy.da.SearchCursor(fc, [num_field]) as cur: for (v,) in cur: digits u.join(ch for ch in str(v or ) if ch.isdigit()) if digits: mx max(mx, int(digits[-4:])) # 取末四位流水 return mx 1 auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, startnext_start(fc, BH), width4, prefixuXZ-2024-)这段的取数逻辑是“把编号里的数字挑出来取末四位当流水号”对XZ-2024-0123这种格式能正确还原出 123。参数上唯一要改的是末四位这个长度如果流水位是 5 位就改成digits[-5:]和width保持一致。这个函数不完美——如果前缀里本身带数字比如年份会被一起挑出来所以取的是末位而不是全部这是一种取舍。从那以后我每次跑完编号都强制走一遍“重复 断号 抽样核对”三连哪怕只有两百条数据。这三步加起来不到两分钟但能挡住交付后被打回来的尴尬。希望帮到你。本文还有配套的精品资源点击获取
简介ARCGIS编号小工具面向地理信息系统数据处理人员用于为mdb、gdb和shp三类常见地理数据批量添加顺序编号缓解大量要素标识时的手工操作压力可服务于数据管理、地图制作与空间规划等场景。压缩包共74个文件约7.1MB内含32个动态链接库、32个说明文档、4个可执行程序、4个配置文件和2个清单文件动态链接库承载ArcGIS相关功能说明文档辅助理解组件用途可执行程序与配置文件负责工具运行和参数设置。已有3992人学习。工具可先选择数据源再设定编号字段与起始值随后自动写入序列号既支持批量处理也可针对特定图层或要素执行借助完整依赖读者能了解地理数据库读写、几何对象处理和地图显示等环节的调用方式并据此完成编号操作或研究工具结构。1. 图斑编号这件小事两千条数据错位之后我才想明白顺序给要素编号听上去是属性表里拉一列的事真做起来耗掉半天。前段时间接了个图斑整理项目两千多个面要求按行政村分组、每组从 1 开始连续编号、不足四位补零还要按从上到下的空间顺序排。同事先在 Excel 里拉序号再连接回 ArcGIS结果排序字段里躺着一堆空值编号整体错位返工两遍。这类活的共性很清楚顺序、分组、格式化、写回四步里只要有一处靠手工就没有可重复性。这个 ARCGIS 编号小工具把四步压成一次执行输入要素类、排序表达式和分组字段编号直接落进属性表。它适合经常和图斑、宗地、管线、监测点打交道的从业者也适合刚上手 arcpy、想找一个能跑通的最小样例的人。2. 拆开编号工具四个动作和一张参数表2.1 一次编号本质是四个动作串起来很多人对编号的第一反应是“写个循环”但真正决定结果对不对的是循环之前的准备工作。我把整个流程拆成四段确定顺序、划定分组边界、生成序列、写回字段。四段是有严格先后依赖的顺序必须在写回之前一次性确定因为写回动作会改变行状态边写边排只会越排越乱。第一段“确定顺序”要回答的是“谁排第一”。可选依据很多按面积从大到小、按地类编码、按行政村代码、按质心 Y 坐标从上到下。这个依据不能是 FID 或 OBJECTID虽然它们看起来天然有序但编辑、导出、合并、裁剪任意一次操作都会让 OID 重排编号会跟着漂移。血泪经验是临时用 OID 排一次可以但一定要把结果落到实体字段里别指望下次打开还是那个顺序。第二段“划定分组边界”本质是一次分组聚合。分组字段取同一值的连续行计数器归零重新从起始值开始。这里有个容易被忽略的点分组和排序必须同时生效SQL 层面就是ORDER BY 分组字段, 排序键如果只排排序键不排分组字段同一组的数据会被拆散在游标的各个位置计数器会反复归零出来就是 1、2、1、3、1 这种鬼样子。第三段“生成序列”涉及起始值、步长、补零位数、前后缀四个参数。补零是最容易被低估的一个zfill(4)把 1 变成0001看着简单但如果字段类型是整型写进去的0001会被数据库转成 1你打开属性表一看还是没补零——不是脚本错了是字段类型不接受。第四段“写回”用arcpy.da.UpdateCursor逐行更新比CalculateField稳的地方在于它能在同一个游标里带上排序和分组信息不用来回切上下文。顺带说一句CalculateField也能做递增编号写个全局变量在表达式里自增就行但那是玄学写法执行顺序不保证重跑一次结果可能变而且中途报错之后全局变量状态是脏的得重启整个会话。我一般只在一次性、几十条的小数据上这么干。2.2 为什么不直接上字段计算器或者 Excel把几种常见做法摆在一起对比选型理由就很清楚了。做法上手成本能否指定多级排序重复执行主要风险字段计算器Python 解析器低否顺序由 OID 决定每次全表重算编辑后编号漂移报错后状态不可控Excel 拉序列 连接回表低手动排序易错需要重做一遍连接键类型不一致、FID 不稳定、文本被转科学计数模型构建器 迭代器中有限排序要靠前置工具可以参数写死在模型里出错定位困难脚本工具这套思路中任意字段、多级、升降序可以需要处理字段类型、文件锁、空值Excel 那条路最坑的地方不是排序是连接。文本型编号0001在 Excel 里会被吞成 1导出 CSV 再连接回来0001和1对不上你会在属性表里看到一大堆空值然后开始怀疑人生。如果非要用表格当中转务必把编号列先设成文本格式或者干脆在前面拼一个字母前缀。字段计算器的优势是零门槛劣势是不可复现。同一个表达式今天跑出来的编号和下周跑出来的编号可能不一致因为中间有人动过数据、加过字段、改过排序。项目交付前最怕的就是这种“上次是对的这次不对”查起来没有抓手。模型构建器的定位其实是“把流程固定下来”思路是对的但它不擅长表达“按分组重置计数”这种带状态的逻辑。迭代器选中分组、内部再迭代要素能拼出来但每加一个分组就多一层嵌套调试成本直线上升。所以最终的形态还是脚本工具核心逻辑一个函数参数从工具箱传进来既能一键跑也能在 Python 窗口里逐步调试。2.3 参数设计一个排序表达式吃掉所有排序需求工具箱的参数不能无限加动态增减参数在 ArcGIS 里做起来很别扭。我的做法是把所有排序需求压缩成一个字符串表达式格式是字段 方向; 字段 方向分号分隔多级方向写ASC或DESC省略则按升序。比如CUNBM ASC; MJ DESC表示先按行政村代码升序、同组内按面积降序。参数数据类型默认值说明输入要素要素图层 / 要素类无支持 GDB 要素类、Shapefile、SDE 图层编号字段字段无文本型或长整型短整型有 32767 上限排序表达式字符串OID ASC分号分隔多级方向 ASC/DESC分组字段字段可选空为空表示全域连续编号起始值长整型1每组重新从这里开始步长长整型1一般不用改跨图层接续时会用到补零位数长整型00 表示不补零仅对文本字段有效前缀 / 后缀字符串空仅对文本字段有效这个设计牺牲了一点直观性——用户得记住表达式格式但换来的是无上限的排序能力。常见做法是在脚本工具的参数说明里写两行示例再加一层解析容错字段名写成MJ、mj、MJ DESC、MJ desc都能认遇到不存在的字段直接抛一条中文错误别让人对着ERROR 000732猜。还有一个细节值得提前定编号字段的类型检查放在最前面做。文本字段走字符串格式化数值字段直接写整数如果用户给数值字段配了前缀或补零脚本主动报错退出而不是悄悄写个1进去让人以为成功了。这种“早失败”比事后核对省太多时间。3. 把编号跑通从命令行最小闭环到工具箱脚本3.1 先确认 arcpy 能被 import不要一上来就写业务代码。先用一条命令确认解释器找得到 arcpy这一步能把一半的环境问题挡在门外。ArcGIS Pro 自带的是 Python 3 环境ArcMap 10.x 带的是 Python 2.7两者的语法和编码处理不一样脚本要分开写或者做兼容。# ArcGIS Pro直接调用安装目录下自带的 Python 环境路径以实际安装为准 C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe -c import arcpy; print(arcpy.GetInstallInfo()[Version]) # 也可以先激活环境再进交互式窗口调试 propy.bat python -c import arcpy; print(arcpy.env.workspace)# ArcMap 10.x自带 Python 2.7注意 print 是语句不是函数 C:\Python27\ArcGIS10.8\python.exe -c import arcpy; print arcpy.GetInstallInfo()[Version]第一段命令的逻辑很简单绕开系统 PATH 里可能存在的另一个 Python直接指定 ArcGIS 自带的解释器避免出现“装了 arcpy 但 import 报错”的假故障。第二段是给老项目用的Python 2.7 里print后面不能跟括号包多个参数unicode和str的区分也要自己处理中文路径建议统一用u...加.decode(utf-8)。参数上唯一要改的是解释器路径把它换成你机器上的实际安装位置即可。如果这一步就报ImportError先查两件事是不是用了系统 Python是不是 32 位和 64 位混装。ArcMap 后台的 Python 是 32 位Pro 是 64 位同一个脚本在两边跑第三方库基本不能共用。3.2 核心函数排序、分组、补零、写回下面这段是整套工具的主干逻辑全部压在一个函数里方便复制到脚本工具里调用也方便单独调试。写法兼容 Python 2.7 和 3.x没用 f-string。# -*- coding: utf-8 -*- import arcpy def auto_number(fc, num_field, sort_keyOID ASC, group_fieldNone, start1, step1, width0, prefix, suffix, whereNone): 按排序表达式给要素写编号。 fc : 要素类路径或图层名 num_field : 编号字段名文本或长整型 sort_key : 字段 方向; 字段 方向方向省略为 ASC group_field : 分组字段None 表示全域连续编号 start/step : 起始值与步长 width : 补零位数0 表示不补 prefix/suffix: 前后缀仅文本字段有效 where : 可选筛选条件只对满足条件的要素编号 # 1) 拿真实字段名和类型Shapefile 会截断字段名大小写也可能和看到的不一致 desc arcpy.Describe(fc) oid_name desc.OIDFieldName field_map dict((f.name.lower(), f) for f in arcpy.ListFields(fc)) if num_field.lower() not in field_map: raise ValueError(u编号字段不存在{0}.format(num_field)) num_fld field_map[num_field.lower()] is_text (num_fld.type String) if not is_text and (width or prefix or suffix): raise ValueError(u数值型字段不支持补零和前后缀请改用文本字段) # 2) 解析排序表达式拼成 ORDER BY 需要的字段列表 fields [num_fld.name] order_sql [] for seg in (sort_key or ).split(;): seg seg.strip() if not seg: continue parts seg.split() fname parts[0] if fname.lower() in (oid, objectid, fid): fname oid_name # 统一映射到真实 OID 字段名 direction DESC if len(parts) 1 and parts[1][0].upper() D else ASC real field_map.get(fname.lower()) if real is None: raise ValueError(u排序字段不存在{0}.format(fname)) order_sql.append({0} {1}.format(real.name, direction)) fields.append(real.name) if group_field: grp field_map.get(group_field.lower()) if grp is None: raise ValueError(u分组字段不存在{0}.format(group_field)) fields.append(grp.name) sql_clause (None, ORDER BY , .join(order_sql)) if order_sql else (None, None) # 3) 单次游标完成遍历与写回分组变化时计数器归零 seq start last_group object() # 哨兵对象避免和真实的 None / 空串撞车 with arcpy.da.UpdateCursor(fc, fields, where, sql_clausesql_clause) as cur: for row in cur: if group_field: g row[-1] if g ! last_group: seq start last_group g if is_text: body str(seq).zfill(width) if width else str(seq) row[0] u{0}{1}{2}.format(prefix, body, suffix) else: row[0] seq cur.updateRow(row) seq step if __name__ __main__: arcpy.env.workspace rD:\work\demo.gdb auto_number(tb_parcel, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, start1, step1, width4, prefix) print(done)逐段说一下逻辑。开头用Describe拿 OID 真实名字是因为 OID 在不同数据源里可能叫OBJECTID、FID、OBJECTID_1用户写OID的时候必须能映射过去否则排序找不到字段直接报错。字段字典用小写做键是为了兼容 Shapefile 那种大小写不敏感但显示不一致的情况。中间解析排序表达式的循环做了三件事拆字段和方向、把OID别名换成真实字段名、校验字段是否存在。校验这一步很值得写因为arcpy.da.UpdateCursor传了不存在的字段只会给一句RuntimeError: ERROR 999999对新手极不友好自己抛一条中文提示能省半小时。最后一段是核心。sql_clause的第二个元素就是ORDER BY子句游标会按这个顺序吐行这是整套方案能成立的前提。last_group用object()当哨兵是因为分组字段可能真的是空字符串或None用None做初始值会让第一组被误判成“换了组”序号白白重置一次。写回时按字段类型分两支文本走zfill加前后缀数值直接写整数。参数上最常调的是sort_key和group_fieldstep基本只在“接着上次的号往下编”时会变成 1 以外的东西。where参数留了个口子只对满足条件的要素编号比如只编现状图斑、不编历史遗留。3.3 挂进工具箱参数怎么配脚本能跑通之后下一步是把它变成别人也能点的工具。ArcGIS Pro 里用“新建脚本工具”ArcMap 里用“添加脚本”本质一样把 .py 文件挂到工具箱然后按顺序定义参数。参数顺序必须和脚本里GetParameterAsText的索引严格对应错一位就是连环错。序号显示名数据类型方向过滤 / 默认值0输入要素要素图层输入无1编号字段字段输入获得自输入要素过滤文本、长整型2排序表达式字符串输入默认OID ASC3分组字段字段输入获得自输入要素可选4起始值长整型输入默认 15步长长整型输入默认 16补零位数长整型输入默认 07前缀字符串输入可选8后缀字符串输入可选入口代码只做参数转接不放业务逻辑这样在 Python 窗口里调试时还能直接调auto_number# -*- coding: utf-8 -*- import arcpy from number_tool import auto_number # 同目录下的核心模块 if __name__ __main__: p arcpy.GetParameterAsText try: auto_number( fcp(0), num_fieldp(1), sort_keyp(2) or OID ASC, group_fieldp(3) or None, # 工具箱里留空会传空字符串要转成 None startint(p(4) or 1), stepint(p(5) or 1), widthint(p(6) or 0), prefixp(7) or , suffixp(8) or , ) arcpy.AddMessage(u编号完成) except Exception as e: arcpy.AddError(u编号失败{0}.format(e)) raise参数转接里最容易被忽略的是p(3) or None这一句。工具箱里没填的参数传进来是空字符串如果不转成Noneif group_field:判断为真脚本会去找一个名字为空的字段然后报“分组字段不存在”。这个坑我调了不止一次后来干脆在所有可选参数上都加or None或or 默认值。AddError之后要raise否则工具在结果窗口里显示“已完成”但实际什么都没做用户以为成功了问题留到验收才爆出来。工具一旦给别人用报错信息就是唯一的沟通渠道这条不能省。4. 四种编号场景连续、分组、空间、编码型怎么配4.1 全域连续编号参数最少坑最集中最简单的用法就是把排序键设成 OID不分组从 1 编到底。这种场景通常出现在“给要素一个唯一序号”的需求上比如管线点、监测点、抽样点。# 不分组按 OID 保序四位补零前缀加一个字母 auto_number(rD:\work\demo.gdb\pt_monitor, BH, sort_keyOID ASC, width4, prefixJC)这里的逻辑是游标按 OID 顺序吐行编号从 1 递增每行格式化成JC0001这种形式写回BH字段。参数上唯一要留意的是width4如果需要编号过万位数要给到 5否则从 10000 开始zfill不再起作用JC10000和JC0001长度不一致后续按字符串排序会错。排序键换成业务字段也很常见比如按面积降序MJ DESC让大图斑排在前面。但这时必须补一个二级排序键MJ DESC; OID ASC因为面积重复的行之间顺序不稳定每次跑出来的编号可能不一样同一份数据两次执行结果不同验收时没法解释。4.2 分组重编号先洗干净分组字段再谈编号按行政村、按图层、按权属单位分组重新编号是最常见的需求也是最容易翻车的。翻车点不在编号逻辑在分组字段本身有的值带前后空格有的是全角空格有的None和空串混用结果就是同一组被拆成三组编号出来是 1、2、1、3。先做一遍清洗再编号# 清洗分组字段去空格含全角、None 转空串避免同组被拆散 fc rD:\work\demo.gdb\tb_parcel with arcpy.da.UpdateCursor(fc, [CUNBM]) as cur: for row in cur: v row[0] or u row[0] v.replace(u\u3000, u).strip() cur.updateRow(row) # 清洗完再编号先按分组字段排组内按面积降序 auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, width4)第一步把CUNBM里的全角空格和首尾空白干掉None统一成空串这一步必须单独做不能塞进编号函数里因为清洗是破坏性操作出问题要能单独回滚。第二步的排序表达式把分组字段放在第一位这是分组编号的硬性要求ORDER BY里分组字段必须在前编号才能保证同组数据连续出现计数器归零的时机才对得上。参数上的对应关系是group_field决定什么时候归零sort_key的第一段必须和group_field一致方向也一致否则游标吐出来的行序和分组边界对不上会出现“组 A 编到 5突然跳到组 B又跳回组 A 从 1 开始”的情况。这个约束我在函数文档里写了但还是有人踩后来干脆加了一条检查如果sort_key的第一段不是group_field直接抛错。4.3 空间顺序编号从上到下、从左到右规划类项目经常要求“按从上到下的顺序编号”本质是按质心 Y 坐标降序、X 坐标升序。SHAPEXY能直接取到要素的质心坐标但要注意坐标系地理坐标系下经纬度排序会歪因为经度一度的实际距离随纬度变化同一行要素的 Y 值差异会被放大。正确做法是先投影到投影坐标系再取坐标。# 1) 加两个双精度字段存质心坐标已存在则跳过 fc rD:\work\demo.gdb\tb_parcel exist set(f.name.lower() for f in arcpy.ListFields(fc)) add [[n, DOUBLE] for n in (CX, CY) if n.lower() not in exist] if add: arcpy.management.AddFields(fc, add) # 2) 写入质心坐标SHAPEXY 返回 (x, y) with arcpy.da.UpdateCursor(fc, [SHAPEXY, CX, CY]) as cur: for row in cur: row[1], row[2] row[0] cur.updateRow(row) # 3) 先按 Y 降序从上到下同一行内按 X 升序从左到右 auto_number(fc, BH, CY DESC; CX ASC, group_fieldCUNBM, width4)第一步用ListFields判断字段是否已存在避免重复执行时报“字段已存在”。第二步把SHAPEXY元组拆成两个数值字段这一步是必须的因为ORDER BY不支持直接对几何排序只能对属性字段排。第三步的排序表达式里CY DESC控制从上到下CX ASC控制同一水平带内从左到右。参数上有个经验值如果要素在空间上呈明显的条带状分布比如沿道路两侧的地块纯按质心排序会出现“上一行还没编完就跳到下一行”的锯齿这时可以在排序表达式里加一个分级字段比如ROW_NO ASC; CX ASC先用人工划定的行号分组再在行内按 X 排。工具本身不生成行号行号通常靠“按位置选择”或人工辅助图层得到。4.4 编码型编号前缀、年份、流水号的拼接政务和规划类项目里编号往往不是纯数字而是XZ-2024-0001这种带业务含义的编码。这类需求的关键不在拼接在字段长度和字符集。# 文本字段长度要提前算够前缀 3 分隔符 1 年份 4 分隔符 1 流水 4 13 # 实际建字段时留余量直接给 32 或 50 prefix uXZ-2024- auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, start1, step1, width4, prefixprefix)逻辑上就是把前缀交给prefix参数流水号由width补零拼出来是XZ-2024-0001。参数上要注意两点一是字段长度GDB 的文本字段有长度上限长度不够会被静默截断XZ-2024-0001变成XZ-2024-00打开属性表才发现二是字符集前缀里有中文的话Shapefile 的编码问题会让人抓狂建议一律用 GDB且前缀只用字母数字和连字符。如果编码需要按年份分组重置做法是把年份字段也放进分组字段的组合里但单字段分组满足不了。常见做法是先把年份和行政村拼成一个辅助字段GRP YEAR || CUNBM用这个字段当group_field编号完再删掉辅助字段。多一步操作但比改脚本逻辑稳妥。5. 避坑与排查五个现场记录和它们背后的原因5.1 字段与数据层面的两个坑现象一脚本报“字段不存在”但属性表里明明有这个字段。原因基本是两个Shapefile 的字段名上限是 10 个字符你建的是BIANHAO实际存的是BIANHAO刚好 7 位没事但如果叫PARCL_CODE落到 shp 里会变成PARCL_COD另一个是大小写ListFields返回的名字可能和你输入的完全一致也可能因为数据源不同有差异。解决方式是所有字段查找都走小写字典并且在报错信息里把可用字段列表打出来一眼就能看出真名是什么。现象二编号跑到一半报溢出或者大数字变成了负数。原因是编号字段用了短整型SmallInteger上限 32767超过就溢出。这个坑在测试数据上永远不会出现因为测试数据只有几十条一上生产就炸。解决办法是建字段时就选长整型Integer上限约 21 亿或者干脆用文本字段。工具里加一条前置检查字段类型是 SmallInteger 且要素数可能超过 3 万时直接提示换字段。5.2 排序与分组层面的两个坑现象三分组编号出现 1、2、1、3 这种乱跳。原因在 4.2 节提过分组字段里有隐藏的空白字符或全半角差异同一组被当成多组。还有一种是None和空串混用游标里None ! 每次遇到都会触发归零。解决办法是先跑一遍清洗再用“汇总统计”按分组字段计数看看组的数量和业务上预期的是否一致对不上就别急着编号。现象四同一份数据跑两次编号结果不一样。原因是排序键有重复值重复行之间的相对顺序由数据库决定不保证稳定。这个问题最阴因为第一次跑完看着是对的第二次跑完也“看着是对的”只有把两次结果做对比才会发现差异。解决办法是排序表达式最后永远补一个OID ASC兜底让排序在全序上唯一。我现在的习惯是任何排序表达式不写 OID 都不许提交。5.3 运行环境层面的一个坑现象五报“无法获取锁”或者“要素类正在被使用”。原因有三种编辑会话没停、属性表窗口开着、文件被别的程序占着比如同步盘、杀毒软件扫描。ArcGIS 对 GDB 的锁比较严格一个打开的属性表就足以让UpdateCursor拿不到写锁。解决办法是按顺序排查先停编辑并保存再关掉所有属性表窗口和目录面板还不行就重启 Pro最后才怀疑是不是同步盘在后台抓文件。生产数据建议放在本地盘或网络盘的非同步目录里跑别放在会自动上传的文件夹下。顺带提一个相关的Shapefile 在编号过程中会生成.lock文件异常中断后锁文件可能残留导致后续操作一直失败。这时不要删.lock硬来正确做法是关掉所有 ArcGIS 进程再重新打开让系统自己释放。6. 编号之后的验收重复值、断号与增量维护的检查清单编号跑完不等于交付完成验收这一步能挡住大部分返工。我一般做三件事。第一件是查重复和空值。用Counter快速扫一遍重复的直接列出来from collections import Counter fc rD:\work\demo.gdb\tb_parcel with arcpy.da.SearchCursor(fc, [BH, CUNBM]) as cur: rows [(r[0], r[1]) for r in cur] dup [k for k, v in Counter(x[0] for x in rows).items() if v 1] empty [i for i, x in enumerate(rows) if not x[0]] print(u重复编号 {0} 个空编号 {1} 个.format(len(dup), len(empty)))第二件是查断号。按分组统计最小值和最大值和组内要素数对比max - min 1 count说明连续不等就是中间跳号了。这个检查用“汇总统计”工具做也行用几行脚本做也行关键是每组都查别只抽查一组。第三件是接续编号。新增要素要接着已有的最大号往下编不能从头开始。做法是先扫一遍现有编号取最大值把它加一作为startdef next_start(fc, num_field): 取现有编号的最大流水值 1兼容带前缀的编码型编号 mx 0 with arcpy.da.SearchCursor(fc, [num_field]) as cur: for (v,) in cur: digits u.join(ch for ch in str(v or ) if ch.isdigit()) if digits: mx max(mx, int(digits[-4:])) # 取末四位流水 return mx 1 auto_number(fc, BH, CUNBM ASC; MJ DESC, group_fieldCUNBM, startnext_start(fc, BH), width4, prefixuXZ-2024-)这段的取数逻辑是“把编号里的数字挑出来取末四位当流水号”对XZ-2024-0123这种格式能正确还原出 123。参数上唯一要改的是末四位这个长度如果流水位是 5 位就改成digits[-5:]和width保持一致。这个函数不完美——如果前缀里本身带数字比如年份会被一起挑出来所以取的是末位而不是全部这是一种取舍。从那以后我每次跑完编号都强制走一遍“重复 断号 抽样核对”三连哪怕只有两百条数据。这三步加起来不到两分钟但能挡住交付后被打回来的尴尬。希望帮到你。本文还有配套的精品资源点击获取
