ArcPy投影转换：定义投影、重投影和空间参考

ArcPy投影转换：定义投影、重投影和空间参考

做 GIS 自动化时，ArcPy投影转换最容易出错的地方不是代码语法，而是把“定义投影”和“重投影”当成同一件事。结果常见表现是图层跑到海上、空间连接全为空、缓冲区距离不对，或者面积计算和 ArcGIS Pro 手工结果不一致。

本文用一个常见项目场景说明：一批道路、地块和 GPS 点来自不同来源，有的缺少坐标系，有的是 WGS 84 经纬度，有的是地方投影坐标。我们要先确认 arcpy空间参考，再决定是否用 DefineProjection 写入正确坐标系，最后用 Project 生成统一坐标系的数据。

问题背景：为什么 ArcPy投影转换 经常跑偏

在真实项目里，同一个工程目录常常混有 Shapefile、File Geodatabase、CAD 转换数据、GPS 点表和第三方平台导出数据。地图上看起来能叠在一起，不代表数据本身的坐标系就正确。ArcGIS Pro 会在地图中动态投影显示图层，但脚本批处理时仍然要依赖数据自己的空间参考。

很多 ArcPy投影转换失败，根源是输入数据的坐标系标签错了。例如坐标值明明是经纬度，却被定义成某个米制投影；或者坐标值是地方投影米制坐标，却被当成 WGS 84。此时直接运行 Project，只会把错误标签当成真实来源继续转换，输出结果一定偏。

正确做法是先判断数据坐标值代表什么，再选择工具。如果数据缺少或写错了坐标系标签，但坐标值本身没有问题，用 DefineProjection；如果数据已有正确空间参考，只是要转换到另一套坐标系，用 Project。

核心原理：arcpy定义投影不是坐标转换

arcpy.management.DefineProjection 的作用是给数据写入或覆盖坐标系信息。它只改空间参考标签，不移动几何坐标。换句话说，输入点坐标是 116.39, 39.90，定义投影后仍然是 116.39, 39.90，只是 ArcGIS 知道这组数字应该按哪个坐标系解释。

arcpy.management.Project 才会创建一个新数据，并把几何坐标从源坐标系换算到目标坐标系。经纬度点重投影到米制投影后，坐标值会明显变化。对面积、长度、缓冲区、叠加分析来说，这一步通常是必须的。

一句话区分：定义投影是“告诉软件这份数据原来是什么坐标系”；重投影是“把数据真正转换成另一个坐标系”。

步骤一：用 arcpy空间参考 先检查输入数据

做 ArcPy投影转换前，第一步不是写目标坐标系，而是读取每个输入数据的空间参考。它通常通过 arcpy.Describe 和 arcpy.SpatialReference 处理。重点看名称、类型、WKID、单位，以及是否为 Unknown。

import arcpy

fc = r"D:\gis_project\work.gdb\road_raw"

desc = arcpy.Describe(fc)
sr = desc.spatialReference

print("数据名称: {}".format(desc.name))
print("空间参考名称: {}".format(sr.name))
print("空间参考类型: {}".format(sr.type))
print("WKID: {}".format(sr.factoryCode))

if sr.name == "Unknown":
    print("坐标系未知：不要直接 Project，先确认真实坐标系")
elif sr.type == "Geographic":
    print("当前是地理坐标系，角度单位: {}".format(sr.angularUnitName))
elif sr.type == "Projected":
    print("当前是投影坐标系，线性单位: {}".format(sr.linearUnitName))

如果输出显示 Unknown，不代表可以随便指定一个目标坐标系。你必须根据坐标值范围、数据来源、元数据、同项目标准图层或原始 .prj 文件判断真实坐标系。比如经纬度坐标通常在经度、纬度范围内，米制投影坐标通常是几十万到几百万的数值。

步骤二：给 Unknown 数据写入正确坐标系

arcpy定义投影适用于“坐标值本身正确，只是数据没有坐标系标签”的情况。典型例子是 GPS 点表转成点图层后显示 Unknown，但你确认经纬度字段来自 WGS 84；或者 CAD 转换数据缺少 .prj，但项目资料明确说明它使用某个地方坐标系。

下面示例把一个 Unknown 点图层定义为 WGS 84。注意，这段脚本不会把坐标转换成米制坐标，只是给数据写入正确空间参考。

import arcpy

unknown_fc = r"D:\gis_project\work.gdb\gps_points_unknown"
correct_sr = arcpy.SpatialReference(4326)

desc = arcpy.Describe(unknown_fc)
current_sr = desc.spatialReference

if current_sr.name == "Unknown":
    arcpy.management.DefineProjection(unknown_fc, correct_sr)
    print("已定义投影: {}".format(correct_sr.name))
else:
    print("当前已有空间参考: {}".format(current_sr.name))
    print("如需覆盖，请先确认原坐标系确实错误，并备份数据")

如果数据已经有正确坐标系，不要为了“转换坐标系”去运行 DefineProjection。这样只会篡改标签，几何坐标不动，后续所有空间分析都会建立在错误解释上。

步骤三：用 arcpy重投影 生成项目标准数据

当输入数据已有正确空间参考后，就可以用 arcpy重投影 生成新数据。项目中更稳的做法是从标准图层读取目标空间参考，而不是手写一个自己不确定的坐标系名称或 WKID。

下面示例把 GPS 点从原始空间参考转换到项目边界图层的空间参考。若源和目标属于不同地理坐标系或基准面，脚本会尝试用 arcpy.ListTransformations 获取可用转换方法，并优先使用返回列表中的第一项。

import arcpy
import os

in_fc = r"D:\gis_project\work.gdb\gps_points_wgs84"
template_fc = r"D:\gis_project\standard.gdb\city_boundary_standard"
out_workspace = r"D:\gis_project\work.gdb"
out_fc = os.path.join(out_workspace, "gps_points_projected")

source_desc = arcpy.Describe(in_fc)
source_sr = source_desc.spatialReference
target_sr = arcpy.Describe(template_fc).spatialReference

if source_sr.name == "Unknown":
    raise ValueError("输入数据缺少坐标系定义，请先确认真实坐标系并运行 DefineProjection")

transformations = arcpy.ListTransformations(source_sr, target_sr, source_desc.extent)

if transformations:
    arcpy.management.Project(
        in_dataset=in_fc,
        out_dataset=out_fc,
        out_coor_system=target_sr,
        transform_method=transformations[0]
    )
else:
    arcpy.management.Project(
        in_dataset=in_fc,
        out_dataset=out_fc,
        out_coor_system=target_sr
    )

print("输出完成: {}".format(out_fc))

这段脚本的关键点有三个：先拒绝 Unknown 输入，目标坐标系来自标准数据，地理转换方法由 ArcPy 根据源、目标和范围列出。实际项目中还应记录使用的转换方法，方便成果复核。

批量 arcpy投影转换 的推荐流程

如果要处理一批要素类，不建议把所有数据直接扔进 Project。批量 arcpy投影转换要先分组：Unknown 数据先人工确认并定义投影，已知坐标系数据再批量重投影，已经是目标坐标系的数据只复制或跳过。

1. 建立标准坐标系。用项目标准边界、标准格网或主管部门提供的模板数据作为目标空间参考。
2. 扫描输入空间参考。用 Describe 输出名称、WKID、类型和单位，保存成检查日志。
3. 处理 Unknown。不要自动猜测。先用元数据、坐标值范围和样本叠加确认真实坐标系，再运行 DefineProjection。
4. 选择地理转换。当源和目标基准面不同，用 ListTransformations 列出可用方法，并记录最终使用项。
5. 运行重投影。对矢量用 Project，对栅格用 ProjectRaster，输出到 File Geodatabase。
6. 做结果质检。检查输出数量、空间范围、坐标单位、面积长度字段和典型样本位置。

批处理时还要避免覆盖原始数据。建议输出命名带上坐标系或用途，例如 roads_prj、parcels_standard，并把原始数据、定义投影后的中间数据和重投影成果分开存放。

常见坑点：这些错误会让重投影结果偏移

• 把 DefineProjection 当成 Project。DefineProjection 不会移动坐标，不能用它把 WGS 84 转成米制投影。
• 对错误标签直接 Project。源坐标系写错时，Project 会按错误来源换算，结果通常偏移很远。
• 忽略地理转换。不同基准面之间转换时，缺少合适的 transformation 可能带来位置偏差。尤其是高精度测绘、宗地和工程项目，不要省略检查。
• 在经纬度坐标系下做距离分析。经纬度单位是度，缓冲区、面积和长度统计应优先转到合适的投影坐标系。
• 把地图显示当成数据已转换。ArcGIS Pro 地图能动态显示不同坐标系图层，但原始数据并不会因此变成同一坐标系。
• 矢量和栅格工具混用。矢量用 Project，栅格用 ProjectRaster。栅格还要考虑重采样方法、像元大小和捕捉栅格。
• 只看图形不看坐标值。投影转换后应抽查坐标范围、单位、要素数量和空间范围，而不是只看地图上是否大致重叠。

方法对比：空间参考对象、定义投影和重投影怎么选

空间参考对象不是单独的转换工具，而是脚本判断和传参的基础。你可以用 SpatialReference 创建目标坐标系，用 Describe 读取输入坐标系，再把它们交给 DefineProjection 或 Project。

实践检查清单：运行脚本前后各查一次

• 输入数据是否有 .prj、元数据或项目说明。
• 坐标值范围是否符合声称的坐标系。
• 空间参考检查结果是否为 Unknown。
• DefineProjection 是否只用于修正标签，而不是用于转换坐标。
• Project 的目标坐标系是否来自项目标准数据或可靠 .prj 文件。
• 源和目标基准面不同时，是否记录了 transformation。
• 输出要素数量是否和输入一致，是否出现空几何或异常范围。
• 点位、边界、面积、长度和空间连接结果是否抽样复核。

如果输出仍然偏移，可以把一个样本要素单独导出，分别查看原始坐标、定义后的坐标、重投影后的坐标。小样本能定位问题，再放回批处理脚本中修正，会比盲目改 WKID 更可靠。

FAQ

ArcPy投影转换时，DefineProjection 和 Project 应该先用哪个？

如果输入数据空间参考是 Unknown，但你确认它真实使用的坐标系，先用 DefineProjection 写入正确标签，再用 Project 做坐标转换。如果输入数据已有正确空间参考，直接 Project，不要再 DefineProjection。

arcpy投影转换后图层位置偏了，最先检查什么？

先检查源数据的空间参考是否正确，而不是先改目标坐标系。很多偏移来自源坐标系标签错误、缺少地理转换、经纬度和米制坐标混淆，或者把地图动态投影误认为数据已经转换。

arcpy定义投影会不会改变坐标值？

不会。它只更新数据的坐标系信息，不修改几何坐标。如果你需要把坐标值从经纬度转换成米制投影坐标，应使用 arcpy.management.Project。

arcpy重投影一定要设置 geographic transformation 吗？

不一定。只有源和目标涉及不同地理坐标系或基准面时才需要重点考虑。稳妥做法是用 arcpy.ListTransformations 根据源空间参考、目标空间参考和数据范围列出可用转换，并记录实际使用的方法。

arcpy空间参考可以用 WKID、名称和 prj 文件，哪种更稳？

三种都可以。WKID 简洁但容易被误抄，名称可读但可能存在相似项，.prj 文件适合复用项目标准。工程项目中更推荐从标准要素类或模板数据读取空间参考，减少手写错误。

结论

ArcPy投影转换的稳定流程是：先用 Describe 检查空间参考，再判断是否需要 DefineProjection，最后用 Project 或 ProjectRaster 输出到项目标准坐标系。不要把定义投影当成重投影，也不要对 Unknown 数据直接转换。

只要把空间参考检查、定义投影和重投影的职责分清，投影转换脚本就会更容易复现、排错和交付。对需要长期批处理的 GIS 项目，还应把输入坐标系、目标坐标系、地理转换方法和质检结果写入日志，形成可追踪的数据处理流程。