ArcPy入门1:Python环境与ArcPy概览

要注意的是，并不是整个窗口都能输入代码，窗口上层是运行结果显示界面，下面小文标注绿色框的位置才是我们需要的代码输入框（有点像AutoCAD的命令窗口）。这个划分跟老版本的Python窗口不大一样，老版本转Pro的童鞋注意一下。

代码输入后，按两次回车键可执行代码。

import arcpyprint(arcpy.GetInstallInfo())

import sysprint(sys.version)

ArcPy采用了模块化的设计，核心模块就是arcpy，按功能又细分为了多个子模块。

模块总览：

arcpy├── arcpy.mp        → 地图文档与制图自动化├── arcpy.da        → 数据访问（游标、字段操作）├── arcpy.sa        → 空间分析（栅格计算、表面分析）├── arcpy.na        → 网络分析（路径、服务区、位置分配）├── arcpy.management → 数据管理工具集├── arcpy.analysis  → 叠加分析、邻域分析、提取分析├── arcpy.ia        → 影像分析（ArcGIS Image Analyst 扩展）├── arcpy.editing   → 高级编辑操作└── arcpy.geoprocessing → 地理处理服务与工具箱管理

这个是负责地图文档（.aprx）的自动化操作，可以打开、保存、导出地图文档，能管理图层、布局、图例、比例尺和打印布局自动化等。我们常说的批量出图，一般就是用这个模块配合地图系列来实现的，可以批量把图件导出为pdf、png、jpeg等格式。

比如批量生成不同区域的专题地图代码如下：

import arcpyaprx = arcpy.mp.ArcGISProject("CURRENT")layout = aprx.listLayouts()[0]layout.exportToPDF(r"E:\\output\\map.pdf")

这个模块提供了高性能的数据读写接口，可以读取查询数据（SearchCursor）、插入新记录（InsertCursor）、更新记录和删除现有记录（UpdateCursor），还能够在Numpy数组与要素类之间相互转换（NumPyArrayToFeatureClass）。

相比于传统的ArcPy游标，arcpy.da游标的速度提高了5-10倍！

示例代码：

import arcpywith arcpy.da.SearchCursor("cities", ["NAME", "POPULATION"]) as cursor:for row in cursor:        print(f"{row[0]}: {row[1]}")

这个模块需要启用Spatial Analyst扩展许可才能使用，可以使用Spatial Analyst的所有栅格数据分析功能。如表面分析、栅格计算、水文分析、插值分析和距离分析等能力。

示例代码：

import arcpy.sa as saslope = sa.Slope("dem.tif", "DEGREE")slope.save("slope.tif")

需要启用Network Analyst 扩展许可，可以执行最短路径、服务器分析、最近设施点、位置分配、起源-目的地成本矩阵等网络分析。

示例代码：

import arcpy.na as na# 创建网络分析图层并进行路径求解

此模块提供地理数据库与图层的基础管理操作，可以复制要素到新位置（CopyFeatures）、坐标系的投影变换（Project）、创建空要素类（CreateFeatureclass）、批量追加数据（Append）和数据融合与聚合（Merge/Dissolve）等。数据管理是自动化GIS数据预处理的第一站，覆盖了90%以上的日常数据整理需求。

import arcpy# 复制并投影到 CGCS2000arcpy.management.CopyFeatures("cities", "cities_copy")arcpy.management.Project("cities_copy", "cities_cgcs2000", 'CGCS2000_3_Degree_GK_Zone_38')

提供矢量数据的核心空间运算，如我们常用的缓冲区分析（Buffer）、裁剪提取（Clip）、叠加分析（Intersect/Union）、按属性选择（SelectLayerByAttribute）和邻近度分析（Near/GenerateNearTable）等等。以上所有的分析工具输出的结果均为新的要素，原始数据不会改动，支持链式GP工作流。

import arcpy# 对城市做 5km 缓冲区，再与行政区求交集arcpy.analysis.Buffer("cities", "city_buffer", "5 Kilometers")arcpy.analysis.Intersect(["city_buffer", "districts"], "city_district_overlap")

提供遥感影像与深度学习分析能力，需要启用Image Analyst 扩展许可，可以进行像素级影像分类（ClassifyPixelsUsingDeepLearning）、目标检测（DetectObjectsUsingDeepLearning）、影像分割（Segment）、波段提取（ExtractBand）、植被指数计算（NDVI）和水体指数计算（NDWI）等。支持GPU加速，深度学习推理速度比传统方法提升5-10倍。

示例代码：

import arcpy# 使用 Deep Learning 检测建筑物arcpy.ia.DetectObjectsUsingDeepLearning(    in_raster="orthophoto.tif",    out_detected_objects="buildings.shp",    in_model_definition="building_detection.dlpk")

提供了地理数据库的事务级编辑的能力，可以开启编辑会话（Editor.startEditing()）、开始编辑操作（Editor.startOperation()）、批量更新及删除记录（UpdateCursor）、调用几何修复工具（Snap/ExtendLine）等操作，支持回滚操作。

import arcpy# 在 GDB 中安全更新属性edit = arcpy.da.Editor("C:/data.gdb")edit.startEditing(False, False)edit.startOperation()with arcpy.da.UpdateCursor("cities", ["NAME", "POPULATION"]) as cursor:for row in cursor:        row[1] = int(row[1] * 1.05)        cursor.updateRow(row)edit.stopOperation()edit.stopEditing(True)

提供工具发现、环境控制与执行监控，如枚举可用工具（ListTools）、获取数据属性和元数据（Describe）、检查数据是否存在（Exists）、环境设置（env）和日志与消息管理（AddMessage）等操作。通过脚本控制GP环境，可以实现跨工程、跨工具的自动化编排。

示例代码：

import arcpy# 设置 GP 环境arcpy.env.workspace = "C:/data.gdb"arcpy.env.overwriteOutput = True# 工具发现与数据检查tools = arcpy.ListTools("*Buffer*")print(f"找到 {len(tools)} 个 Buffer 相关工具")if arcpy.Exists("cities"):    arcpy.AddMessage("数据就绪，开始处理...")

# -*- coding: utf-8 -*-"""第一个ArcPy脚本：创建要素缓冲区作者: NeverLandv日期: 2026-08-19"""import arcpy# === 步骤1: 导入模块并设置环境 ===# 设置工作空间（存放输入/输出数据的目录）arcpy.env.workspace = r"E:\\GIS_Data\\project"# 设置输出覆盖（开发调试阶段可开启）arcpy.env.overwriteOutput = True# === 步骤2: 验证数据存在 ===input_fc = "cities.shp"ifnot arcpy.Exists(input_fc):    raise FileNotFoundError(f"输入数据不存在: {input_fc}")# === 步骤3: 调用地理处理工具 ===# 使用 arcpy.analysis 子模块中的 Buffer 工具output_fc = "cities_buffer.shp"arcpy.analysis.Buffer(    in_features=input_fc,    out_feature_class=output_fc,    buffer_distance_or_field="1000 Meters",    line_side="FULL",    line_end_type="ROUND",    dissolve_option="NONE")# === 步骤4: 输出结果信息 ===count = arcpy.management.GetCount(output_fc)print(f"缓冲区创建成功！")print(f"输出路径: {arcpy.env.workspace}\\{output_fc}")print(f"生成要素数量: {count}")print(f"ArcPy版本: {arcpy.GetInstallInfo()['Version']}")

下期预告：

GIS+AI给中华白海豚当"保镖"：厦门怎么做到的