2.3 对齐实操：定义投影、投影变换与动态对齐

2.3 对齐实操：定义投影、投影变换与动态对齐

本节摘要：图层对齐的标准流程是先诊断、后修复、再验证。修复动作只有两个——给坐标系缺失的图层"定义投影"（贴标签），给坐标系不同但明确的图层做"投影变换"（换算坐标），两者用反会把好数据改坏。本节完整走一遍多源图层对齐案例，界面与 ArcPy 双轨给出。

回到故障现场

章首那条故障：道路图层与行政区图层差了几百公里，前者是 CGCS2000 经纬度，后者是 CGCS2000 高斯 3 度带。现在手里有了前两节的武器，可以把它修到毫米级对齐。先别急着点工具，按诊断、修复、验证三段走——这套节奏适用于一切空间参考问题，比记住任何单个工具重要。

诊断阶段的三个检查动作，顺序不能乱：

1. 查明双方空间参考：用属性页或 Describe 脚本，明确各自坐标系名称与 WKID。Unknown 的先补标签。
2. 看数值范围：小数点后五六位的两位数是经纬度，七八位的整数是米制。确认数值与声明的坐标系相符——声明与实际不符的"伪投影"图层是最隐蔽的坑。
3. 判断偏移性质：整层系统性偏移几十米以内，多半是基准差异（54 对 2000）；偏移几百公里或缩成一小团，是经纬度被当成米制（或反之）硬画。

图：对齐排错决策流

动作A：定义投影——只贴标签，不动数字

坐标系为 Unknown 的图层，几何坐标数字其实是完好的，缺的只是标签。定义投影工具做的事就是把标签贴上：输入图层加目标坐标系，输出还是那份数据，坐标数字一个都没变。

界面路径：ArcToolbox 中的数据管理工具、投影与变换、定义投影。脚本同样干脆：

import arcpy
# 场景：roads 图层数值是 6 位 X 7 位 Y 的米制坐标，源自旧资料，确认是高斯117度带
arcpy.management.DefineProjection(
    "K:/gisdata/incoming/roads_unknown.shp",
    arcpy.SpatialReference(4534)      # CGCS2000 3度带 中央经线117度
)
# 验证：重新读取坐标系名称
sr = arcpy.Describe("K:/gisdata/incoming/roads_unknown.shp").spatialReference
print("已定义为:", sr.name)
# 输出: 已定义为: CGCS2000_3_Degree_GK_Zone_39
​

⚠️ 常见坑：把定义投影当万能修复。如果图层实际是经纬度，你却给它贴上高斯投影的标签，坐标数字被按米解释，图层会缩成原点附近的一小团。定义投影只陈述事实，不创造事实——贴标签前必须人工确认图层真实坐标系（看数值范围、查数据来源、找元数据）。

动作B：投影变换——真正换算坐标数字

两个图层坐标系都明确但不一致时，用投影变换生成一份新坐标系的拷贝。这一步是实打实的数值重算：经纬度换算成米，或基准之间的椭球换算。

import arcpy
# 把经纬度道路图层变换到与行政区一致的高斯117度带
arcpy.management.Project(
    in_dataset="K:/gisdata/incoming/roads_geo.shp",   # CGCS2000 经纬度 4490
    out_dataset="K:/gisdata/aligned/roads_gk117.shp",
    out_coor_system=arcpy.SpatialReference(4534),     # 目标：高斯117度带
    transform_name=""   # 同为 CGCS2000 基准 不需基准变换方法
)
# 快速核对：变换后首个要素的坐标数量级
with arcpy.da.SearchCursor("K:/gisdata/aligned/roads_gk117.shp", ["SHAPE@XY"]) as cur:
    x, y = next(cur)[0]
print(f"首个折点: X={x:.1f} Y={y:.1f}")
# 输出示例: 首个折点: X=449361.8 Y=4417904.3   —— 七位数，符合米制预期
​

跨基准的情形必须多给一个参数。旧 54 基准的地形图要素类转到 CGCS2000：

import arcpy
# 列出本机可用的 54 到 2000 变换方法
methods = [m.name for m in arcpy.ListTransformations(
    arcpy.SpatialReference(4214),    # Beijing1954 地理坐标系
    arcpy.SpatialReference(4490))]   # CGCS2000 地理坐标系
print("可用变换方法:", methods[:3])
# 输出示例: 可用变换方法: ['Beijing_1954_To_CGCS2000_1', 'Beijing_1954_To_CGCS2000_2', ...]
# 执行带基准变换的投影
arcpy.management.Project(
    "K:/gisdata/oldmap/contour_bj54.shp",
    "K:/gisdata/aligned/contour_cgcs2000.shp",
    arcpy.SpatialReference(4490),
    transform_name="Beijing_1954_To_CGCS2000_1"   # 按测区选择合适的方法
)
​

变换方法的选择有地区性：不同方法在不同省区的拟合精度不同，正式项目应采用当地测绘部门推荐的方法，并在元数据中记录所用方法与残差。

动态对齐：不换数据也能叠

其实 ArcGIS 的地图文档本身就有动态投影能力：数据保持原坐标系，显示时软件临时换算到地图的坐标系。这解释了一个新手疑惑——"我没做任何变换，两个图层怎么也叠上了？"那是软件在幕后替你做了动作B的只读版本。

动态投影适合浏览与临时叠加，不适合两种情况：一是分析运算——多数地理处理工具直接取输入数据的坐标数字运算，不做动态换算，混合坐标系的输入会产生错误结果；二是性能敏感场景——每次刷新都重算投影是白白付出的代价。生产纪律一句话：分析前把所有输入统一到同一投影坐标系，一劳永逸。

验证：修复不算完，验证才算完

对齐操作后必须做三件事：目视检查边界咬合（道路应贴着行政边界走，不交叉不悬空）；抽查已知地物（某个显眼的地标应该在正确位置）；量测偏差（选一对同名点计算坐标差，写入项目日志）。第三件事脚本可代劳：

import arcpy
# 同名点对：图上收费站与权威参考坐标
p_layer = arcpy.PointGeometry(arcpy.Point(449382.1, 4417895.5), arcpy.SpatialReference(4534))
p_ref   = arcpy.PointGeometry(arcpy.Point(449385.9, 4417901.2), arcpy.SpatialReference(4534))
print("偏差距离: %.2f 米" % p_layer.angleAndDistanceTo(p_ref, "GEODESIC")[1])
# 输出: 偏差距离: 6.61 米   —— 54转2000 后的典型残差量级，可接受并记录在案
​

本节要点回顾

• 三段节奏：诊断（查参考、看数值、判偏移性质）、修复（定义投影或投影变换）、验证（咬合、地标、偏差量测）
• 两个动作方向相反：定义投影贴标签不改数字，投影变换换算数字产新档
• 跨基准必须给变换方法：54 对 2000 差几十米，方法选择有地区性并须留档
• 动态投影只管看：分析前统一坐标系是硬纪律
• 验证留痕：偏差量测写进项目日志，责任与复现都要靠它

第 2 章到此收束。图层对齐之后，空间参考的纪律将从单张图层上升到整个数据库——第 3 章讲地理数据库如何用要素数据集、属性域与拓扑，让几百张图层长期保持今天对齐好的状态。