如何用ArcGIS制作重心迁移轨迹？

Question

在使用ArcGIS制作重心迁移轨迹时，一个常见技术问题是：如何准确计算不同时期地理要素的几何中心并连接形成迁移路径？用户常因数据投影不当导致面积变形，影响重心计算精度；或在处理面要素（如行政区划）时，未能正确应用“要素转点”工具提取几何中心。此外，跨时期重心点的时间序列排序错误，也会导致迁移箭头方向混乱。需注意选择合适的投影坐标系、利用“汇总统计数据”管理时间字段，并通过“XY to Line”工具构建连续迁移轨迹线。

Answer 1

一、问题背景与核心挑战

在地理信息系统（GIS）的空间分析中，重心迁移轨迹是一种用于揭示地理现象时空演变规律的重要可视化手段。广泛应用于人口迁移、经济活动分布、生态变化等领域。其核心在于：对不同时期的地理要素（如行政区划面数据）计算几何中心，并按时间序列连接形成连续迁移路径。

然而，在使用ArcGIS实现该过程时，用户常面临三大技术瓶颈：

1. 投影坐标系选择不当：使用地理坐标系（如WGS84）直接计算面积和质心，导致高纬度区域严重变形，影响重心精度；
2. 面转点操作不规范：未正确使用“要素转点”工具中的“质心”选项，或混淆了“标注点”与“真实几何中心”；
3. 时间序列排序错误：生成的重心点未按时间顺序排列，导致后续构建的迁移线方向混乱，无法反映真实迁移趋势。

二、关键技术流程解析

为确保重心迁移轨迹的科学性与准确性，需遵循以下标准化处理流程：

步骤	操作内容	对应工具	注意事项
1	定义投影坐标系	“投影”工具	选择等积投影（如Albers或Lambert）以保持面积不变形
2	提取各时期面要素质心	“要素转点”工具	勾选“Inside”参数确保点位于多边形内部
3	添加X/Y坐标字段	“添加XY坐标”工具	自动写入质心经纬度（投影后单位为米）
4	按时间字段汇总统计	“汇总统计数据”工具	以时间字段分组，保留每个时期的质心记录
5	排序并生成序列表	字段计算器 + 排序功能	创建“Sequence”字段表示时间先后
6	构建迁移线段	“XY to Line”工具	设置起点与终点为相邻时间的质心坐标
7	符号化与标注	图层渲染设置	使用箭头样式突出迁移方向

三、深度技术实现细节

以下代码片段展示了如何通过Python脚本调用ArcPy自动化完成部分关键步骤：

import arcpy
from arcpy import env

# 设置工作空间
env.workspace = "C:/data/centroid_analysis.gdb"

# 输入面数据（含Time字段）
input_polygon = "AdministrativeRegions"
output_centroids = "Centroids_Raw"

# 步骤1：投影至等积坐标系（例如Asia_North_Albers_Equal_Area_Conic）
projected_fc = "Projected_Regions"
arcpy.Project_management(input_polygon, projected_fc, 
                        arcpy.SpatialReference(102025))

# 步骤2：提取质心
arcpy.FeatureToPoint_management(projected_fc, output_centroids, "CENTROID")

# 步骤3：添加XY坐标
arcpy.AddXY_management(output_centroids)

# 步骤4：按Time字段汇总，获取每期唯一质心
summary_table = "Centroid_Summary"
arcpy.Statistics_analysis(output_centroids, summary_table,
                          [["POINT_X", "MEAN"], ["POINT_Y", "MEAN"]],
                          "Time")
    

四、可视化迁移路径构建

利用“XY to Line”工具将相邻时间点连接成矢量线，是形成动态迁移轨迹的关键。假设已有序列化的时间点数据，包含当前期与下一期的X/Y坐标：

• 起始点字段：MEAN_POINT_X, MEAN_POINT_Y
• 终止点字段：LEAD_MEAN_X, LEAD_MEAN_Y（可通过滞后函数生成）
• 连接字段：SequenceID，用于匹配前后帧

执行命令如下：

arcpy.XYToLine_management(
    in_table="Centroid_Pairs",
    out_feature_class="Migration_Trajectory",
    start_x_field="MEAN_POINT_X",
    start_y_field="MEAN_POINT_Y",
    end_x_field="LEAD_MEAN_X",
    end_y_field="LEAD_MEAN_Y",
    line_type="GEODESIC",
    id_field="SequenceID"
)

五、流程图示：重心迁移轨迹生成逻辑

graph TD
            A[原始面数据] --> B{是否投影?}
            B -- 否 --> C[重投影至等积坐标系]
            B -- 是 --> D[执行要素转点]
            C --> D
            D --> E[添加XY坐标]
            E --> F[按时间字段汇总]
            F --> G[排序生成时间序列]
            G --> H[构造起止点对]
            H --> I[XY to Line生成轨迹]
            I --> J[符号化输出地图]