ArcGIS按图斑批量出图：Data Driven Pages与arcpy实战

简介《Arcgis按图斑批量出图》是一份面向GIS从业者与ArcMap用户的实操型文档针对卫片图斑数量多、逐一手动出图耗时低效的问题梳理按图斑自动批量导出JPG图片的完整思路。文档涵盖数据驱动页面配置、布局动态文本设置、arcpy遍历页面导出、按图斑编号命名及结果检查等关键环节并给出可参考的脚本示例与定制方向便于读者理解自动化出图流程。压缩包共1个doc文件约29KB体量轻便适合作为ArcGIS 10.2及以上环境下的操作参考。该文档已有8322人学习适合需要提升大批量图斑出图效率、希望掌握数据驱动页面与arcpy导出方法的初中级GIS用户查阅。1. 一个图斑一张图最磨人的从来不是画图手上有一份遥感监测图斑面层业内口语叫卫片图斑几百上千个要素每个要素要出一张带编号、面积、所在单元的 JPG。手工在 ArcMap 里拖框、改文本框、点导出一天出不了两百张还会串号——第三页的图上写着第七个图斑的属性这种错最难查。ArcGIS 按图斑批量出图这件事桌面端有一套不用装插件的解法Data Driven Pages 管版面索引arcpy.mapping 管循环导出两段接起来就是一条流水线。前提是 ArcGIS 10.2 及以上9.x 没有 Data Driven Pages不用往下看。这套流程适合做图斑核查底图、外业巡查图册、地块台账配图的从业者也适合手里已经有 mxd 模板、只想把它跑成批量的同学。2. 索引图层怎么定Data Driven Pages 的四项配置与取舍Data Driven Pages下称 DDP的模型很朴素挑一个图层当「索引」一条要素记录对应一页页面数量等于要素数量翻一页数据框范围就跳到那条要素身上版面里的动态文本跟着换值。听着简单真正决定出图质量的是三件事——索引图层用什么、名称字段怎么选、比例尺固定还是自适应。这三件事在配置对话框里五分钟就能点完但点错一项后面几百张图全废重跑一遍又是半小时。2.1 索引图层直接拿图斑面层还是单独做一层点最省事的做法是把图斑面层本身设为索引图层要素数就是页数一页对一个图斑。图斑之间不重叠、面积相差不悬殊的时候这么干完全够用也省掉一次数据加工。但只要图斑有相邻、有嵌套、有大面积整片图斑面层做索引就开始出问题。页面范围是以要素的外接矩形为基础算的两个相邻图斑的范围会互相压出图时邻斑大面积入侵版面一个上千公顷的整片图斑会把比例尺撑得极小图上什么都看不清面层里如果有碎面、洞、多部件要素页数和业务上认为的「图斑个数」还会对不上。所以常规做法是单独做一个索引点层每个图斑取一个内部点点层字段里带齐编号、面积、类型、所在单元再拿这个点层当索引。点做索引的好处是页面范围以点为中心、缩放行为稳定版面干净页数和图斑数严格一一对应。生成内部点用「要素转点」Feature to Point勾上 Inside 选项保证点落在面内而不是质心跑到凹多边形外面。这里有个容易被忽略的门槛Feature to Point 属于高级许可级别的工具许可不够时按钮是灰的。许可不够的替代方案是用游标直接读几何真实质心再写点几行代码的事# -*- coding:utf-8 -*- # 许可级别不够时的替代方案读面层真实质心写成索引点层 # 运行环境ArcGIS 10.2 自带 Python 2.7 import arcpy src rD:\GIS_Work\ShpDemo\Patch.shp # 图斑面层 out rD:\GIS_Work\ShpDemo\PatchIndex.shp # 输出索引点层 sr arcpy.Describe(src).spatialReference arcpy.CreateFeatureclass_management( rD:\GIS_Work\ShpDemo, PatchIndex.shp, POINT, spatial_referencesr) arcpy.AddField_management(out, PATCH_ID, TEXT, field_length32) arcpy.AddField_management(out, PATCH_AREA, DOUBLE) # SHAPETRUECENTROID 返回的是几何真实质心不依赖许可级别 with arcpy.da.SearchCursor(src, [PATCH_ID, SHAPETRUECENTROID, SHAPEAREA]) as s: with arcpy.da.InsertCursor(out, [PATCH_ID, SHAPEXY, PATCH_AREA]) as i: for pid, centroid, area in s: i.insertRow([pid, centroid, area])这段脚本的逻辑是先建一个空点层并挂上两个字段再用SHAPETRUECENTROID令牌把每个面的质心读出来配合SHAPEAREA一起写进点层。SHAPEXY是插入游标认的点几何写法直接给元组即可。字段名 PATCH_ID 要和后面 DDP 的名称字段对得上长度给 32 够用了太短会被截断导致重名。注意真实质心对凹多边形可能落在面外做核查底图时如果要求点必须在图斑内部宁可先用「要素折点转点」取顶点再筛或者老老实实上高级许可。2.2 名称字段与排序字段唯一、可读、不含非法字符DDP 配置对话框里第一个下拉是索引图层第二个是 Name Field页面名称字段第三个是 Sort Field排序字段。Name Field 是整条流水线的关键脚本里那张图的文件名就是从它身上取的它决定了输出文件的命名。对名称字段的要求按重要性排第一必须唯一重复了脚本就会拿同一个文件名反复覆盖最后目录里只剩一张图你还以为是脚本没跑第二必须是文本或能被稳定转成文本的短整型浮点型会带小数点出图文件名变成12.0.jpg第三不能含空格和\ / : * ? |这些字符Windows 建不了这种文件第四长度别太长路径总长超过 260 字符时 ArcGIS 报的错往往和真正原因对不上。没有现成合格字段是常态。用字段计算器拼一个Python 解析器下# 字段计算器表达式Python 解析器字段类型选文本 TB str(!OBJECTID!).zfill(5)zfill(5)补前导零好处是按文件名排序时顺序和编号顺序一致不会出现 10 排在 9 前面。这一步不算玄学但确实省掉后面找图的麻烦。排序字段建议按编号排别用 OBJECTID。OBJECTID 在数据导出、合并、修复几何之后会重排今天跑出来的顺序和上周跑出来的不一样对比两批成果时会怀疑人生。字段类型也要留意文本型排序是字典序10 会排在 9 前面数值型才是真按大小排。2.3 显示范围与比例尺固定比例尺和最佳适配的差别DDP 对话框的 Extent 选项卡是第二个决策点它决定每页数据框怎么缩放。常见的三种模式差别很大模式行为适用场景代价Center and maintain current scale以要素为中心保持当前比例尺不变核查底图、需要图与图之间可比小图斑版面留白多大图斑装不下Best fit每页自动缩放到要素刚好填满版面图册、汇报材料每张图比例尺都不同不能横向比较面积Data driven scale按指定字段的值设定每页比例尺已有分级出图规则需要预先算好一份比例尺字段业内做图斑核查我一般建议固定比例尺。比如统一 1:5000 或 1:10000这样翻图册的时候两张图上同一个图斑的大小差异就是真实面积差异看图的人不用每次去找比例尺。代价是小图斑周围留白多看着不饱满但这个代价可以接受。做汇报图册、要给非专业人员看的时候Best fit 更讨喜每张都撑满缺点是比例尺跳来跳去评审时容易被问「这两张能比吗」。Extent 选项卡里还有一个 Margin 百分比控制要素外扩的留白。给 5% 到 10% 比较常见太小了图斑贴边图廓线会把边界压掉太大了图斑缩在中间浪费版面。设完 Margin 点 Apply翻两页看看效果再决定。另外还有一个旋转选项卡按字段旋转数据框图斑走向特殊的时候能救急但一般不碰旋转之后指北针和比例尺的语义会变得别扭。2.4 配置完先翻一遍页再谈脚本配置动作本身是ArcMap 菜单 Customize → Toolbars → Data Driven Pages调出工具栏点第一个按钮进 Setup 对话框选索引图层、名称字段、排序字段切到 Extent 选项卡设范围和比例尺确定。然后把 mxd 存盘。存盘这一步千万别省。脚本读的是磁盘上的 mxd 文件不是 ArcMap 内存里的当前状态。配置改了没存盘跑出来的还是上一版这种低级错误我见过不止一次排查时先怀疑 Python最后发现是没按 CtrlS。存盘后用工具栏上的翻页箭头从第一页翻到最后一页抽几个头、中、尾的页面看范围对不对、图斑有没有跑出图廓、动态文本有没有超框、图例有没有被裁。这一步花五分钟能省掉后面跑完三百张再一张张挑的时间。3. 动态文本与版面让每张图自带图斑身份出图的价值一半在图上一半在图上的字。一张只有图斑轮廓的 JPG交出去对方还得回来问「这是哪个图斑、多大面积、在哪个单元」。版面上把这些信息用动态文本挂上去每页自动换值这才是批量出图相对于手工出图真正省事的地方。3.1 三类动态文本分别对应什么需求ArcMap 的动态文本分三类用途不一样很多人只用了第一类就开始抱怨信息不够。第一类是 DDP 内置标签跟着数据驱动页面走DDP_PAGENAME当前页名称、DDP_CURRENTPAGE当前页码、DDP_PAGECOUNT总页数。做「第 12 页 / 共 386 页」这种页脚用这三个就够。第二类是索引图层的属性字段。插入方式是在 Insert 菜单 → Dynamic Text → Data Driven Page Attribute然后从索引图层里挑字段。选完之后版面上会出现一个标签每翻一页标签内容就换成当前页那条记录对应字段的值。图斑面积、图斑类型、所在单元、核查人都能这么挂上去。第三类是普通动态文本和 DDP 无关但每页都会重算日期、用户名、当前数据框比例尺。比例尺那个尤其实用Best fit 模式下每页比例尺不一样把dyn typedataFrame propertyscale/挂到图廓下方看图的人心里有数。这里有个血泪经验动态文本必须从菜单插入不能手打尖括号。手打出来的DDP_PAGENAME就是一段普通文字翻页时纹丝不动你会以为功能坏了。判断方法很简单选中那个文本框看属性真动态文本在 Properties 里能看到类型定义手打的只有纯字符串。3.2 版面元素哪些会被逐页刷新理解哪些元素刷新、哪些不刷新能避免版面设计踩空。数据框范围每页都会重设这是 DDP 的核心动作。数据框里图层的符号系统默认不变但 Symbology 选项卡提供了一个「按页面名称字段符号化」的开关打开之后当前页对应的图斑可以单独高亮其余图斑用浅色或灰度。做核查图的时候这个功能很值一个图斑一张图但周边环境要保留参考全彩会把视线带偏。图例的行为要单独说。默认图例会把数据框里所有图层都列出来图斑层如果有十几个类别图例能占半张图。在图例属性的 Items 选项卡里勾上「Only show classes that are visible in the current view extent」图例就只列当前范围内的类别版面立刻清爽。静态元素——图廓线、指北针、比例尺条、单位 LOGO 位置、固定的说明文字——每页完全一样改一次管全部不需要在脚本里做任何事。3.3 版面自检清单出批量前的五分钟版面这块的问题跑一张看不出来跑三百张就是灾难。给一份我自己每次都要过的清单检查项怎么看不合格的后果动态文本是否溢出文本框挑名称最长的那个图斑翻到它那页文字被裁图上信息缺失页脚总页数是否正确看最后一页的「共 N 页」页数对不上说明索引层有多余要素比例尺条与数据框比例是否一致固定比例尺下抽查三页比例尺条用的是旧值读图会误判图例是否只列当前范围类别翻到类别最少的图斑图例占版面一半图斑被挤小高亮符号是否生效翻任意一页看当前图斑满屏同色找不到目标图斑图廓外留白是否足够看面积最大的图斑图斑顶到图廓线边界看不清这张表看着琐碎但它对应的是真实翻车场景批量跑了两个小时打开目录一看三成图上的面积文字被裁掉了半个数字只能改版面重跑。五分钟的自检换两小时怎么算都值。4. arcpy.mapping 循环导出脚本骨架、导出参数与命名清洗到这一步mxd 已经配好、版面已经调好剩下的事就一件让程序替你把每一页导出来。ArcGIS 10.x 自带的是 Python 2.7arcpy.mapping模块就在里面不需要额外装任何东西也不需要配环境变量——直接用 ArcGIS 安装目录下那个python.exe跑就行。4.1 最小可跑脚本骨架先给能跑通的最小版本把整条链路走通再谈健壮性。# -*- coding:utf-8 -*- # 运行环境ArcGIS 10.2 自带的 Python 2.7 import arcpy import os MXD rD:\GIS_Work\ShpDemo\IndexMap.mxd # 已配好 DDP 并保存的 mxd OUT_DIR rD:\GIS_Work\Output\JPG # 输出目录必须已存在 mxd arcpy.mapping.MapDocument(MXD) ddp mxd.dataDrivenPages page_count ddp.pageCount print(total pages: %d % page_count) for page_num in range(1, page_count 1): # 上界 1否则漏最后一页 ddp.currentPageID page_num # 切到第 page_num 页 name ddp.pageRow.getValue(ddp.pageNameField.name) # 取当前页名称 out_jpg os.path.join(OUT_DIR, str(name) .jpg) arcpy.mapping.ExportToJPEG(mxd, out_jpg, resolution200) print(%d - %s % (page_num, name)) del mxd # 释放文件锁否则 mxd 一直被占用 print(done)逐行说清楚。arcpy.mapping.MapDocument(MXD)打开存盘的 mxd返回一个地图文档对象注意它读的是磁盘文件ArcMap 界面里没存盘的修改它看不到。mxd.dataDrivenPages拿到 DDP 对象pageCount是总页数等于索引图层的要素数。循环里currentPageID赋值就是翻页这个属性从 1 开始计数赋值超出范围会直接抛异常。ddp.pageRow返回当前页对应的索引图层行对象getValue(字段名)从这行里取值字段名用ddp.pageNameField.name拿这样即使以后改了名称字段脚本也不用动。取到的值可能是 unicode拼路径前用str()转一下更稳。ExportToJPEG(mxd, out_jpg, resolution200)是把当前版面按 200 DPI 导成 JPG。注意它导出的是「当前状态」的 mxd所以必须先设currentPageID再调导出顺序反了就变成三百张一模一样的图。最后del mxd不是可选项。不释放的话 mxd 文件会被 Python 进程锁住下次想改配置打开就提示被占用严重的要重启机器。4.2 ExportToJPEG / ExportToPNG 的参数怎么给导出函数有一串可选参数默认值大部分够用但有两三个必须显式给不然出来的图尺寸或质量不受控。参数作用建议取值map_document地图文档对象用 MapDocument 打开的 mxd不能用图层文件out_jpeg / out_png输出全路径必须带扩展名目录要提前建好函数不会自动创建data_frame指定导出的数据框单数据框可省多数据框版面必须显式指定否则只出第一个df_export_width / df_export_height输出像素宽高留空按版面尺寸走一旦填了就要和分辨率一起算容易把版面文字比例搞乱world_file是否附带世界文件出图给外业看False 即可color_mode色彩模式24-BIT_TRUE_COLOR 常规8-BIT_PALETTE 体积小但渐变色会断层background_color背景色默认白深色底图要单独设resolution输出 DPI屏幕看 150打印 200300分辨率这一项是体积和清晰度的直接权衡。1:5000 的版面150 DPI 出来的 JPG 大约 300KB 到 800KB300 DPI 能到 2MB 以上。三百个图斑跑下来一个是两百多兆一个是六七百兆交付前先想清楚对方是打印还是屏幕看。宽度和高度那两个参数我的建议是能不填就不填。填了之后它按像素定尺寸而版面里的字号、线宽是按页面单位定义的两者不匹配时会出现「图变大了但字没变大」或者反过来看着特别别扭。需要控制清晰度就动 resolution需要控制构图就动版面尺寸两个一起动必然打架。PNG 和 JPEG 的差别也顺手说一句PNG 无损、支持透明背景适合叠在别的底图上JPEG 有压缩、体积小适合发给外业用手机看。业内出核查底图JPEG 更常见。4.3 文件名清洗与单页失败不中断最小脚本能跑通但不能直接上生产。图斑名称里带空格、带斜杠、带括号的情况太常见了一个非法字符就能让整批任务在第 137 页崩掉前面 136 张白跑。先做一层名称清洗# -*- coding:utf-8 -*- import re _BAD re.compile(r[\\/:*?|\s]) # Windows 文件名非法字符 空白 def safe_name(raw, fallback): s _BAD.sub(_, str(raw).strip()) # 非法字符统一替换成下划线 s s.rstrip(.) # Windows 不允许文件名以点结尾 if not s: s fallback # 清洗后为空用页码兜底 return s[:80] # 截断防止路径总长超 260re.sub那行把非法字符全部换成下划线空格也一起换掉——空格本身在 Windows 里合法但传到下游系统、邮件附件、压缩包里经常出问题索性统一处理。rstrip(.)是补 Windows 的一个坑文件名以点结尾时资源管理器不给建报错信息还特别含糊。[:80]的截断长度是按输出目录通常不超过 150 字符估的留出余量。再套一层 try/except让单页失败不影响整批# -*- coding:utf-8 -*- import arcpy, os, csv mxd arcpy.mapping.MapDocument(rD:\GIS_Work\ShpDemo\IndexMap.mxd) ddp mxd.dataDrivenPages OUT_DIR rD:\GIS_Work\Output\JPG log open(rD:\GIS_Work\Output\export_log.csv, wb) log.write(page,name,status\n) for page_num in range(1, ddp.pageCount 1): ddp.currentPageID page_num raw ddp.pageRow.getValue(ddp.pageNameField.name) name safe_name(raw, page_%03d % page_num) out_jpg os.path.join(OUT_DIR, name .jpg) try: arcpy.mapping.ExportToJPEG(mxd, out_jpg, resolution200) log.write(%d,%s,ok\n % (page_num, name)) except Exception as e: # 单页失败写日志继续最后按日志补跑 log.write(%d,%s,fail:%s\n % (page_num, name, str(e).replace(,, ;))) log.close() del mxd日志里 page 列是页码补跑时按这个页码组合成 range 就能只重跑失败的那几页。异常信息里的逗号替换成分号是为了不破坏 CSV 的列结构——这个小细节不处理用 Excel 打开日志就全串列了。5. 避坑与排查五条踩过的记录下面这几条都是实际跑批量时撞过的按「现象 → 原因 → 解决」写撞上了直接对号入座。5.1 最后一页永远不出图现象索引图层有 386 个要素DDP 工具栏也显示 386 页跑完目录里只有 385 张缺的正好是最后一个图斑。原因range(1, pageCount)的上界是开区间。Python 的 range 不包含终点这么写实际只跑到 pageCount - 1。原始脚本里那一行就是这个写法直接抄过来必然少一页。解决写成range(1, pageCount 1)。另外顺手确认一下pageCount取的是 DDP 对象的属性不是索引图层的要素总数——两者在索引图层有筛选、有定义查询的时候会不一致DDP 只认参与分页的要素。5.2 中文路径或中文名称导出报错、乱码现象输出路径放在中文目录下脚本报 IO 错误或者抛一个看不出所以然的异常有时候能跑通但生成的文件名是乱码。原因Python 2.7 里str和unicode混用是经典问题。从getValue拿到的可能是 unicodeos.path.join拼上 str 类型的路径时中间会尝试用 ASCII 解码遇到中文就炸。ArcGIS 的桌面端对中文路径本身是支持的问题出在编码转换这一步。解决文件头加上# -*- coding:utf-8 -*-拼路径前统一转成 str或者干脆全流程用 unicode 字面量uD:\\...输出目录尽量用纯英文这条不是迷信是省时间——真出问题的时候能排除掉一个变量。5.3 文件名带特殊字符整批任务中断现象跑到第 137 页突然抛异常退出报错信息指向文件创建失败前面 136 张都在后面全没有。原因图斑编号里带了/或者:当成文件名就变成了路径分隔符或盘符标记。还有一种情况是名称字段里有回车换行肉眼在属性表里看不出来导出时直接失败。解决走 4.3 节那层safe_name清洗非法字符统一替换。清洗完记得看一眼日志里的 name 列确认替换后的名字仍然能认出来是哪个图斑——全变成下划线也失去了可读性必要时在文件名里保留原编号的数字部分。5.4 导出的图空白或者范围整体偏移现象图能导出来但打开一看是空白版面或者图斑跑到图廓外面去了只剩一个角。原因三种可能。一是导出时currentPageID没设成功导的是默认第一页且第一页的索引要素恰好落在别处二是 DDP 配置改了但 mxd 没存盘脚本读的是旧的版面设置三是多数据框版面下没指定data_frame参数导出的不是你以为的那个框。解决按顺序排查。先在 ArcMap 里手动翻到出问题的那一页确认 DDP 本身没问题再确认 mxd 存盘时间晚于最后一次配置修改最后检查版面里有几个数据框多于一个就在导出函数里显式传data_frame。多数据框的版面还有个隐性坑导出时只有指定数据框会跟着 DDP 翻页其他框保持不动如果那个框里也放了图斑层会显示成上一次的范围。5.5 跑到一半崩了mxd 被锁住打不开现象脚本中途异常退出再去 ArcMap 里打开那个 mxd提示文件被占用或只读。原因Python 进程异常退出时del mxd没执行文件句柄没释放或者脚本还在后台跑着ArcMap 同时打开了同一个文件。解决任务管理器里结束残留的 python 进程养成把导出逻辑包在 try/finally 里的习惯把del mxd放到 finally 里异常也能释放。另外一个实用习惯是跑批量的时候不要用 ArcMap 打开同一个 mxd用副本或者干脆关掉省得互相锁。提示这五条里5.1 和 5.3 出现频率最高写脚本时先把这两处处理掉能省掉大半的返工。6. 从「跑一次」到「可复用」整本 PDF、分批重跑与日志习惯单批跑通之后真正拉开效率差距的是第二件事——把这条流水线变成下次改个路径就能用的东西。这里说三个我常用的手法。第一整本 PDF 一次出图册存档用这个。DDP 对象自带一个exportToPDF方法不需要循环一条语句把全部页面合成一个多页 PDF# -*- coding:utf-8 -*- import arcpy mxd arcpy.mapping.MapDocument(rD:\GIS_Work\ShpDemo\IndexMap.mxd) ddp mxd.dataDrivenPages # ALL 全部页RANGE 配合 page_range_string 只出指定页 ddp.exportToPDF(rD:\GIS_Work\Output\atlas.pdf, page_range_typeALL, resolution200, image_qualityBEST, compress_vectorsTrue, multiple_filesFalse) del mxdpage_range_type有 ALL、CURRENT、RANGE、SELECTED 四个取值配合page_range_string比如1,3,5-20就能只出某几页。multiple_filesTrue会按页拆成多个 PDF出一本还是出一堆看交付要求。PDF 相比 JPG 的好处是矢量要素不糊、文字可选中、体积还小给内部存档很合适。第二失败页分批重跑。上一章那个 CSV 日志就是为这个准备的。把 fail 的行筛出来页码拼成12,45,137-140这样的字符串塞进page_range_string只重跑这几页不用整批重来。三百多页里坏掉十几页的时候这个操作能把返工时间从四十分钟压到两分钟。第三把路径参数抽到脚本头部。下面这些不要散在代码里mxd 路径、输出目录、分辨率、名称字段、清洗规则。抽成一组常量放在文件最上面下次换个项目只改这五六行。我自己的模板里还加了一行前置检查——输出目录不存在就os.makedirs省得每次忘了建目录被函数报一个语焉不详的错。最后说个习惯。从那以后我每次配完 DDP都先憋着不跑全量改成先跑前五页打开目录一张张看完确认版面、动态文本、比例尺、文件命名全对再把range放开到全量。这套动作多花三分钟省下的是跑完三百张再发现文字被裁、全部重来的那两个小时。希望帮到你。本文还有配套的精品资源点击获取