ArcPy应用详解(含arcpy create feature class实操教程)

很多团队在做 ArcGIS 数据生产与分析时，最头疼的是把“临时应付”的手工步骤沉淀为可复用的自动化流程：新建要素类、统一坐标系、追加字段、写入几何，再到索引与约束。一旦规范了这一链路，数据治理与性能都会显著提升。作为 GIS研习社的 Dr.gis，我用一篇实践导向的指南，带你把思路与代码一次打通。

ArcPy 是什么、能做什么

ArcPy 是 ArcGIS Pro/ArcMap 的 Python 接口，封装了地理处理工具（Geoprocessing）与数据访问（Data Access, arcpy.da）。它的价值在于：

• 自动化与可复现：把点点点的操作转成脚本，确保每次结果一致。
• 规模化：循环处理多个区域、多期数据，减轻重复劳动。
• 可维护：版本化管理脚本，便于团队协作与审计。

类比：如果 ModelBuilder 是“可视化装配线”，ArcPy 则是“可编程装配线”，更灵活、可测试、可集成。

为什么要用脚本创建要素类

用脚本创建要素类的核心诉求是标准化：统一几何类型、坐标系、字段、域与索引，从源头减少后期清洗成本。同时，脚本能把规范变成“执行即检查”的契约，避免口口相传的误差。

选择何种目标存储，取决于数据治理与下游应用需求。生产环境更推荐 .gdb。

核心参数与概念对齐

• out_path/out_name：输出路径与名称，路径可为 .gdb、文件夹（Shapefile）、内存工作空间 in_memory。
• geometry_type：POINT、MULTIPOINT、POLYLINE、POLYGON、MULTIPATCH 等。
• spatial_reference：坐标系，可用 WKID（如 4326）或 .prj 路径。
• has_m / has_z：是否启用 M/Z 值，取 "ENABLED"/"DISABLED"。
• template：用已有要素类或图层克隆结构（字段、域、子类型等）。

提示：生产数据优先使用明确的 spatial_reference。未知坐标系会在空间分析和地图叠置时埋雷。

快速上手：从空白到可用要素类

1. 设置环境与工作空间。
2. 准备目标 .gdb（不存在则创建）。
3. 调用 CreateFeatureclass 创建空要素类。
4. 补充业务字段与索引。
5. 用 arcpy.da.InsertCursor 写入几何与属性。

import arcpy
from pathlib import Path

# 1) 环境与路径
arcpy.env.overwriteOutput = True
root = Path(r"D:gis_projectdata")
gdb = root / "prod.gdb"
arcpy.env.workspace = str(gdb)

# 2) 准备 GDB
if not gdb.exists():
    arcpy.management.CreateFileGDB(str(root), gdb.name)

# 3) 创建要素类（以点为例）
sr = arcpy.SpatialReference(4326)  # WGS 84
fc_name = "pois"
arcpy.management.CreateFeatureclass(
    out_path=str(gdb),
    out_name=fc_name,
    geometry_type="POINT",
    template="",
    has_m="DISABLED",
    has_z="DISABLED",
    spatial_reference=sr
)

# 4) 添加字段
fc = str(gdb / fc_name)
arcpy.management.AddField(fc, "name", "TEXT", field_length=100, field_alias="名称")
arcpy.management.AddField(fc, "category", "TEXT", field_length=30, field_alias="类别")
arcpy.management.AddField(fc, "score", "DOUBLE", field_alias="评分")

# 可选：建立复合索引
arcpy.management.AddIndex(fc, ["category", "score"], "idx_cat_score", "NON_UNIQUE", "ASCENDING")

# 5) 写入数据
rows = [
    {"xy": (116.391, 39.907), "name": "A 咖啡", "category": "cafe", "score": 4.5},
    {"xy": (116.398, 39.909), "name": "B 书店", "category": "book", "score": 4.8},
]
with arcpy.da.InsertCursor(fc, ["SHAPE@XY", "name", "category", "score"]) as cur:
    for r in rows:
        cur.insertRow([r["xy"], r["name"], r["category"], r["score"]])

print("创建与写入完成：", fc)

经验：字段名尽量小写蛇形或统一约定，便于跨库、跨语言互操作（如 name、category）。

进阶：基于模板克隆模式

当已有“标准要素类”时，可把它当模板传给 CreateFeatureclass 的 template 参数，直接克隆字段、域、子类型等结构，大幅减少重复配置。

# 假设 template_fc 为标准要素类
template_fc = str(gdb / "standard_roads")
target_name = "roads_2025q1"

arcpy.management.CreateFeatureclass(
    out_path=str(gdb),
    out_name=target_name,
    geometry_type="POLYLINE",
    template=template_fc,  # 关键：模板
    has_m="ENABLED",
    has_z="DISABLED",
    spatial_reference=""   # 留空则随模板
)

# 现在 target_name 已继承模板的字段/域/子类型等设置

进阶：坐标系与投影注意事项

• WKID vs .prj：WKID 快速、稳定；自定义坐标系用 .prj。
• 水平/垂直坐标系：涉及高程分析时，建议同时指定垂直坐标系。
• 投影转换：源数据坐标系不一致时，统一用 Project 工具。

# 以 .prj 定义坐标系
prj_path = r"D:gis_projectprjCGCS2000_3_Degree_GK_Zone_38.prj"
sr = arcpy.SpatialReference(prj_path)

# 统一投影
arcpy.management.Project(
    in_dataset=str(gdb / "pois"),
    out_dataset=str(gdb / "pois_gk38"),
    out_coor_system=sr
)

进阶：域、子类型与规则

域（Domain）与子类型（Subtype）是数据约束与业务语义的关键，能在源头提高数据质量。

# 创建代码域
arcpy.management.CreateDomain(str(gdb), "RoadClass", "道路等级", "SHORT", "CODED")
for code, desc in [(1, "高速"), (2, "快速"), (3, "主干")]:
    arcpy.management.AddCodedValueToDomain(str(gdb), "RoadClass", code, desc)

# 在要素类上设置子类型字段并定义子类型
fc = str(gdb / "roads_2025q1")
arcpy.management.AddField(fc, "class", "SHORT", field_alias="等级")
arcpy.management.SetSubtypeField(fc, "class")
arcpy.management.AddSubtype(fc, 1, "高速")
arcpy.management.AddSubtype(fc, 2, "快速")

# 绑定域到字段
arcpy.management.AssignDomainToField(fc, "class", "RoadClass")

常见坑与性能优化

• Shapefile 限制：字段名最多 10 字符、日期无时间、字符编码与空值处理受限；生产尽量用 .gdb。
• 坐标系未定义：坐标数值“看着像经纬度”不等于定义正确，务必设置 spatial_reference。
• 路径与锁：避免超长路径与并发写锁；循环写入时注意关闭 Cursor（with 语法）。
• 内存工作空间：中间结果可放 in_memory，最终用 CopyFeatures 落地到 .gdb。
• 批量插入：把属性准备成列表/生成器，减少 Python 与 GP 频繁切换。

# in_memory 加速中间计算，然后落地
temp_fc = r"in_memorytemp_lines"
arcpy.management.CreateFeatureclass("in_memory", "temp_lines", "POLYLINE", spatial_reference=arcpy.SpatialReference(3857))
# ... 中间写入与处理 ...
arcpy.management.CopyFeatures(temp_fc, str(gdb / "lines_final"))

实战：批量生成多区域要素类

示例：按城市批量创建面要素类，统一字段结构，写入示例多边形。

import arcpy
from pathlib import Path

arcpy.env.overwriteOutput = True
root = Path(r"D:gis_projectdata")
gdb = root / "admin.gdb"
if not gdb.exists():
    arcpy.management.CreateFileGDB(str(root), gdb.name)

cities = {
    "beijing": [(116.2,39.7), (116.6,39.7), (116.6,40.0), (116.2,40.0), (116.2,39.7)],
    "shanghai": [(121.2,31.0), (121.8,31.0), (121.8,31.4), (121.2,31.4), (121.2,31.0)]
}
sr = arcpy.SpatialReference(4326)

for name, ring in cities.items():
    fc = str(gdb / f"city_{name}")
    arcpy.management.CreateFeatureclass(str(gdb), f"city_{name}", "POLYGON", spatial_reference=sr)
    arcpy.management.AddField(fc, "city", "TEXT", field_length=30)
    with arcpy.da.InsertCursor(fc, ["SHAPE@", "city"]) as cur:
        arr = arcpy.Array([arcpy.Point(x, y) for x, y in ring])
        poly = arcpy.Polygon(arr, sr)
        cur.insertRow([poly, name])

测试与验证

• 结构验证：Describe/ListFields 检查几何类型、坐标系与字段是否符合预期。
• 规则验证：域与子类型是否生效；插入非法值是否被拦截。
• 性能验证：大样本插入与空间索引对查询/渲染时间的影响。

# 结构检查
fc = str(gdb / "pois")
desc = arcpy.Describe(fc)
print(desc.shapeType, desc.spatialReference.name)
for f in arcpy.ListFields(fc):
    print(f.name, f.type, f.length)

参考文献

• ArcGIS Pro 工具参考：Create Feature Class
• ArcPy Data Access：InsertCursor
• ArcPy 类参考：SpatialReference
• ArcGIS Pro 工具参考：Create File GDB
• ArcGIS Pro 工具参考：Add Field
• ArcGIS Pro 工具参考：Add Index
• ArcGIS Pro 工具参考：Create Domain

结语

本文以“是什么/为什么→怎么做”的脉络，系统梳理了用 ArcPy 创建要素类的关键环节：参数选择、模板克隆、坐标系处理、域与子类型、常见坑与优化。掌握这些要点，你就能把手工流程沉淀为稳健、可复用的生产脚本。

你在团队里最常见的“建库建模难点”是什么？欢迎留言讨论。如果希望获取更多示例脚本与规范模板，关注 GIS研习社(gisyxs.com)。