ArcPy批量制图怎么做？GIS自动化脚本详解（附：源码）

引言：告别重复劳动，拥抱 GIS 自动化

你是否经历过这样的场景：项目截止日期临近，甲方突然要求对 500 个行政村的土地利用现状图进行批量输出，而且每个村的图名、图例和统计表都需要动态变化？如果你还在使用 ArcGIS 的“导出地图”功能一张张手动操作，或者在“数据驱动页面”（Data Driven Pages）的有限功能中挣扎，那么这篇文章就是为你准备的。

在 GIS 行业，“制图工厂”式的重复劳动是扼杀工程师创造力的最大杀手。手动制图不仅效率低下，而且极易因疲劳产生人为错误。掌握 ArcPy（ArcGIS 的 Python 站点包），是每一位资深 GISer 从“操作员”进阶为“开发者”的必经之路。

本文将深入浅出地解析 ArcPy 批量制图的核心逻辑，不讲晦涩的理论，只教你如何用代码解决实际问题。我们将通过一个真实的案例，带你编写一段能够自动替换文本、调整范围并批量导出 PDF/JPG 的自动化脚本。

为什么选择 ArcPy 而非“数据驱动页面”？

在深入代码之前，我们需要明白为什么 ArcPy 是批量制图的终极解决方案。很多初学者会问：“Data Driven Pages 也能批量导出，为什么还要写代码？”

以下是两者的深度对比，帮你理清思路：

简而言之，数据驱动页面适合简单的“切片式”出图，而 ArcPy 能够应对复杂的、定制化的项目需求。

核心逻辑：ArcPy 批量制图的三部曲

编写一个健壮的批量制图脚本，通常遵循以下三个核心步骤。理解这个逻辑比死记代码更重要。

第一步：模板准备（MXD 设置）

脚本无法凭空创造地图，它需要一个“模具”。在 ArcMap 中配置好一个 .mxd 文档是基础：

• 固定元素：配置好指北针、比例尺等不需要变动的元素。
• 动态元素标识：对于需要代码修改的文本（如标题），需要在其属性中设置 Element Name（元素名称）。例如，将标题的元素名设为 "MapTitle"，这相当于给代码留了一个“接口”。
• 索引图层：准备一个包含所有制图区域（如行政区划）的矢量图层，用于驱动地图范围的变化。

第二步：获取引用与游标遍历

这是脚本的“大脑”。我们需要引用地图文档，并利用搜索游标（SearchCursor）遍历索引图层。

核心代码逻辑解析：

1. 引入 arcpy 模块。

2. 使用 arcpy.mapping.MapDocument 获取 MXD 对象。

3. 使用 arcpy.da.SearchCursor 读取索引图层（例如：某省的县界）。

4. 进入循环：每读取一行记录（一个县），就执行一次制图操作。

第三步：修改元素与导出

在循环内部，我们需要根据当前记录的信息，动态更新地图内容并保存。

• 更新范围：获取当前要素的几何形状（Shape），将数据框（DataFrame）的范围设置为该要素的范围。
• 更新文本：查找名为 "MapTitle" 的文本元素，将其内容修改为当前记录的“县名”。
• 导出文件：使用 arcpy.mapping.ExportToJPEG 或 PDF 函数，将当前状态导出为文件。

实战源码：通用批量制图脚本

以下是一个可以直接复用的 Python 2.7 (ArcMap 环境) 脚本模板。请注意代码中的注释，根据你的实际路径进行修改。

# -*- coding: utf-8 -*-

import arcpy

import os

# 1. 设置工作空间和变量

mxd_path = r"C:ProjectTemplate.mxd" # 模板路径

output_folder = r"C:ProjectOutput" # 输出文件夹

layer_name = "IndexLayer" # 驱动范围的图层名称

name_field = "XZQMC" # 用作文件名的字段（如行政区名称）

# 初始化地图文档对象

mxd = arcpy.mapping.MapDocument(mxd_path)

df = arcpy.mapping.ListDataFrames(mxd)[0] # 获取第一个数据框

# 2. 获取驱动图层

# 注意：该图层必须要在 MXD 的图层列表中存在

lyr = arcpy.mapping.ListLayers(mxd, layer_name)[0]

# 3. 开始遍历要素

print("开始执行批量制图...")

# 使用 with 语句保证游标正确关闭，字段包括：形状字段 SHAPE@ 和 名称字段

with arcpy.da.SearchCursor(lyr, ["SHAPE@", name_field]) as cursor:

for row in cursor:

feature_geom = row[0]

feature_name = row[1]

# 4. 核心操作：更新地图范围

df.extent = feature_geom.extent

# 可选：将比例尺放大一点，避免图形顶格

df.scale = df.scale * 1.1

# 5. 核心操作：更新标题

# 查找所有文本元素，找到名为 "MainTitle" 的元素

for elm in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):

if elm.name == "MainTitle":

elm.text = feature_name + "土地利用现状图"

# 6. 导出地图

out_name = os.path.join(output_folder, feature_name + ".jpg")

arcpy.mapping.ExportToJPEG(mxd, out_name, resolution=300)

print("已导出: " + out_name)

# 7. 释放内存

del mxd

print("全部完成！")

扩展技巧：高手才知道的避坑指南

在实际生产环境中，仅仅能跑通脚本是不够的，你还需要注意以下两点，这能让你显得更加专业：

1. 内存泄漏与对象释放

ArcMap 的 Python 进程（32位）内存非常有限。如果你需要连续导出上千张图，脚本往往会因为内存溢出而崩溃。技巧： 务必在脚本最后使用 del mxd。对于超大批量的任务，建议每处理 100 张图，就重置一次 mxd 对象，或者利用 Python 的 subprocess 模块调用外部脚本，每张图启动一个新进程（虽然慢，但绝对稳定）。

2. 文本元素的动态布局

当替换的标题长度不一时（例如“北京市” vs “新疆维吾尔自治区伊犁哈萨克自治州”），短标题会居中，长标题可能会偏离中心甚至超出纸张。技巧： 在 MXD 模板中，将文本元素的对齐方式设置为“居中对齐”，并将锚点设置在文本框的中心。这样无论文字长短，它都会向两边延伸，始终保持居中。

FAQ：关于 ArcPy 制图的高频问答

Q1: 这段代码可以在 ArcGIS Pro 中运行吗？

答：不能直接运行。 ArcGIS Pro 使用的是 Python 3.x，且制图模块已由 arcpy.mapping 升级为 arcpy.mp。两者的对象模型有较大差异（例如 Pro 中引入了 Project 对象）。如果你使用的是 Pro，需要将 arcpy.mapping.MapDocument 替换为 arcpy.mp.ArcGISProject，并调整相关方法的调用。

Q2: 导出的图片中中文显示为乱码或方框，怎么办？

答：这通常不是代码的问题，而是字体缺失。 确保你的 MXD 模板中使用的字体在操作系统的字体库中存在。如果是在服务器端进行自动化出图，务必检查服务器是否安装了相应的中文字体（如宋体、黑体）。此外，确保脚本文件本身保存为 UTF-8 编码（并在头部声明），以正确处理中文字符串。

Q3: 如何批量合并导出的 PDF 文件？

答：ArcPy 提供了 PDF 文档管理功能。 你可以使用 arcpy.mapping.PDFDocumentCreate 创建一个新的 PDF，然后利用 appendPages 方法将刚才批量导出的单页 PDF 依次追加进去，最后使用 saveAndClose 保存。这比使用第三方 PDF 工具更高效且易于集成。

总结

ArcPy 批量制图并不是一项高不可攀的技术，它本质上是将我们手动操作的步骤“翻译”给计算机听。通过模板配置、游标遍历、元素更新这三个核心环节，你可以将数天的工作量压缩到几分钟内完成。

不要害怕报错，从最简单的替换标题开始尝试。当你第一次看着文件夹里自动生成的数百张精美地图时，你会发现，GIS 自动化的魅力才刚刚开始。