用 arcpy 按属性批量导出 shp：Select_analysis 与 SearchCursor 实战配置

1. 从一份海南行政区 shp 说起为什么手动拆分不现实如果你手里有一份带分类字段的 shp比如省级行政区里包含十几个地市、全国区县里包含几百个区县而下游任务又要求每个类别单独成一个 shp 文件那么手动在 ArcMap 或 ArcGIS Pro 里反复「按属性选择 → 数据 → 导出数据」基本是不可行的。字段值一多重复操作几十上百次不仅耗时还容易漏掉某个类别或者命名写错。这类需求在 GIS 数据处理里非常典型按属性批量导出 shp。核心思路其实只有两步第一步用 SearchCursor 把分类字段里的唯一值读出来第二步用 Select_analysis 按每个唯一值做一次条件筛选并输出成独立 shp。arcpy 把这两步串起来就是一个可复用的批量拆分脚本。这篇内容面向已经装了 ArcGIS DesktopArcMap或 ArcGIS Pro、能用 Python 2.7 或 Python 3 跑 arcpy 的读者。我会给出可直接复制的脚本骨架重点讲清楚 Select_analysis 的 where_clause 怎么写、SearchCursor 怎么读字段、字段名怎么替换、导出后怎么校验结果以及几个我实际踩过的坑。整套流程不需要你懂太多 Python照着改字段名和路径就能落地。2. 前置准备环境、数据与 TaoToken 的定位2.1 arcpy 运行环境确认arcpy 不是 pip 能装的普通库它随 ArcGIS 一起安装。ArcMap 对应 Python 2.7ArcGIS Pro 对应 Python 3.x两者语法在字符串和 print 上略有差异。你可以在 ArcMap 的 Python 窗口里直接import arcpy测试或者在命令行运行import arcpy print(arcpy.GetInstallInfo()[Version])能打印出版本号说明环境没问题。如果报ImportError: No module named arcpy通常是用了系统自带的 Python而不是 ArcGIS 自带的解释器。ArcMap 用户请用C:\Python27\ArcGIS10.x\python.exePro 用户用 Pro 安装目录下的arcgispro-py3环境。2.2 数据字段先看清楚以海南行政区 shp 为例属性表里有一个cityName字段存的是各地市名称。你要做的第一件事是确认这个字段的真实名称因为 shp 字段名有时会被截断成 10 个字符以内中文名也可能被转成拼音或英文。在 ArcMap 里右键图层 → 打开属性表看列名或者用 arcpy 直接列字段import arcpy fc rD:/shpData/海南.shp for f in arcpy.ListFields(fc): print(f.name, f.type)把打印出来的字段名记下来后面脚本里的fields和where_clause都要和它完全一致大小写敏感。2.3 TaoToken 在这里能帮什么批量拆分脚本写完后往往还要做字段映射、坐标系转换、批量重命名、甚至把结果接到后续的自动化流程里。这类「写脚本 调模型辅助」的场景我习惯用 TaoToken 来统一管理模型调用。它的官网是 https://taotoken.net/?utm_sourcetaotoken_aicg_blog_endutm_mediumcsdnutm_campaignrewriteutm_content API 入口是 https://taotoken.net/api 。如果你只是想让模型帮你解释一段 arcpy 报错、生成 where_clause 模板用模型对话就够了如果要把脚本辅助能力接进长期的编码或 Agent 工作流可以看 Coding Plan需要自己拿 Key 调接口的去 API Keys 页面创建即可。下面第 3 节先把脚本骨架给全第 4 节再讲怎么验证。3. 可复制配置SearchCursor Select_analysis 脚本骨架3.1 完整脚本ArcMap / Python 2.7 版下面这份脚本可以直接改路径和字段名使用。核心逻辑是先用 SearchCursor 遍历分类字段把唯一值收集到一个集合里用集合去重避免同一个类别导出多次再循环调用 Select_analysis。# -*- coding: utf-8 -*- import arcpy # 工作空间 arcpy.env.workspace rD:/shpData arcpy.env.overwriteOutput True # 允许覆盖同名输出省去手动删除 # 输入要素与分类字段 inFc rD:/shpData/海南.shp splitField cityName # 替换成你实际的分类字段名 outFolder rD:/shpData/out # 输出目录建议单独建一个 # 第一步用 SearchCursor 读取唯一值 uniqueValues set() with arcpy.da.SearchCursor(inFc, [splitField]) as cursor: for row in cursor: if row[0] is not None: uniqueValues.add(row[0]) print(u共读取到 %d 个唯一值 % len(uniqueValues)) # 第二步循环 Select_analysis 导出 for val in uniqueValues: outFc outFolder / str(val) .shp # where_clause 注意字段名用双引号字符串值用单引号 where_clause %s \%s\ % (splitField, val) arcpy.Select_analysis(inFc, outFc, where_clause) print(u成功导出 outFc) print(u全部完成)如果你用的是 ArcGIS ProPython 3把print改成函数形式print(...)u...前缀去掉即可其余逻辑一致。arcpy.da.SearchCursor比老版arcpy.SearchCursor快很多推荐优先用 da 版本。3.2 where_clause 的三种写法对照where_clause 是 Select_analysis 最容易写错的地方字段名和值的引号规则必须分清。下面用表格对照字段类型写法示例说明文本字段cityName 海口字段名双引号值单引号数值字段POP 100000数值不加引号含单引号的文本NAME OBrien值内单引号用两个单引号转义注意在 Python 字符串里写 where_clause 时外层引号要和内层错开。上面脚本用的是%s \%s\外层单引号字段名双引号值单引号这样拼出来就是cityName 海口符合 SQL 表达式要求。3.3 输出命名与中文值处理如果分类字段的值是中文比如「海口」「三亚」直接拼进文件名一般没问题但个别系统对中文路径支持不好。稳妥做法是加一层映射把中文值转成拼音或编号再命名。简单点也可以保留中文只要输出目录路径不含特殊字符即可。另外值里如果带空格或斜杠文件名会非法建议先做一次清洗safeName str(val).replace(/, _).replace(\\, _).replace( , _) outFc outFolder / safeName .shp4. 验证请求与成功结果怎么确认拆分没出错4.1 控制台输出与文件数量核对脚本跑完后控制台会打印每个导出的路径。第一层校验是数量唯一值有几个输出 shp 就应该有几个.shp/.shx/.dbf/.prj 是一组别只数 .shp。你可以在输出目录里用 Python 快速统计import glob files glob.glob(rD:/shpData/out/*.shp) print(u输出 shp 数量%d % len(files))如果数量对不上多半是某个值在 where_clause 里没匹配到或者被 overwriteOutput 覆盖了。4.2 属性与几何双重校验数量对只是第一步还要确认每个 shp 里的要素确实属于对应类别。用 SearchCursor 抽查一个输出文件outFc rD:/shpData/out/海口.shp with arcpy.da.SearchCursor(outFc, [cityName, SHAPEAREA]) as cursor: for row in cursor: print(row[0], row[1])打印出来的 cityName 应该全是「海口」几何面积不为 0。如果出现空要素或面积异常检查输入 shp 的坐标系是否定义正确以及 where_clause 是否误匹配。4.3 用要素计数做批量校验更省事的办法是循环统计每个输出的要素数和原始数据里该值的计数对比for val in uniqueValues: outFc outFolder / str(val) .shp cnt int(arcpy.GetCount_management(outFc).getOutput(0)) print(u%s 要素数%d % (val, cnt))原始数据里每个类别的计数可以用arcpy.Statistics_analysis或直接在属性表里看。两边一致基本可以确认拆分正确。5. 本篇常见错排查Select_analysis 与 SearchCursor 的坑5.1 报错「Parameters are not valid」或 where_clause 无效最常见的原因是字段名写错或引号不匹配。shp 字段名有 10 字符限制cityName可能实际存成了cityName或CITYNAME用第 2.2 节的 ListFields 确认。另一个原因是值里含单引号没转义按 3.2 表格处理。5.2 SearchCursor 读不到值或报「字段不存在」arcpy.da.SearchCursor的第二个参数是字段列表必须和真实字段名完全一致。如果字段名是中文确保脚本文件编码是 UTF-8 且加了# -*- coding: utf-8 -*-。ArcMap 的 Python 2.7 对中文支持较弱必要时把字段名改成英文再跑。5.3 输出 shp 打不开或属性表为空如果输出文件在 ArcMap 里打不开先检查输出路径是否含中文或空格尽量用纯英文路径。属性表为空通常是 where_clause 没匹配到任何要素但 Select_analysis 仍会生成一个空 shp。这时回到 4.3 的计数校验定位是哪个值出了问题。5.4 重复导出与覆盖问题老版脚本里用arcpy.Exists判断再删除其实设了arcpy.env.overwriteOutput True之后就不需要手动删了。但要注意如果两个不同的值清洗后文件名相同比如「A/B」和「A_B」后一个会覆盖前一个。命名清洗时加个序号或哈希可以避免。5.5 性能问题几百个类别跑得慢SearchCursor 读唯一值很快慢的是循环里的 Select_analysis每个类别都要扫一遍全量数据。类别特别多时可以改用arcpy.SplitByAttributes_analysis需要 Advanced 许可或者先按字段建索引。普通几十个类别当前脚本完全够用。6. 把脚本接进你的工作流从单次拆分到长期自动化脚本跑通一次之后真正省时间的是把它变成可复用工具。我的做法是把路径、字段名、输出目录抽成参数用arcpy.GetParameterAsText包一层做成 ArcMap 或 Pro 里的工具箱脚本以后拖进 shp 就能跑。如果你还要在脚本基础上做字段映射、批量投影、结果打包或者让模型帮你生成不同字段类型的 where_clause 模板可以按需分流单纯问模型怎么写表达式用模型对话 https://taotoken.net/api?utm_sourcetaotoken_aicg_blog_endutm_mediumcsdnutm_campaignrewriteutm_content 要把辅助能力接进长期编码流程看 Coding Plan https://taotoken.net/api/coding-plan?utm_sourcetaotoken_aicg_blog_endutm_mediumcsdnutm_campaignrewriteutm_content 需要自己创建 Key 调接口的去 API Keys https://taotoken.net/api/api-keys?utm_sourcetaotoken_aicg_blog_endutm_mediumcsdnutm_campaignrewriteutm_content 接入细节查文档 https://taotoken.net/api/doc?utm_sourcetaotoken_aicg_blog_endutm_mediumcsdnutm_campaignrewriteutm_content 。最后留一个实用技巧批量拆分前先复制一份原始 shp 做备份输出目录单独建别和输入混在一起。跑完用 4.3 的计数校验过一遍确认每个类别的要素数和原始一致再进入下游流程。这套流程我在多个行政区划和土地利用分类数据上用过字段名换一换就能直接迁移。