ArcGIS如何生成所有面要素间的20000×20000欧氏距离矩阵

方案1：ArcGIS生态内实现 #

• 使用生成邻近表工具（路径：ArcToolbox > 分析工具 > 邻域分析 > 生成邻近表）：该工具支持输出全量要素对的距离结果，使用时将搜索半径设为可覆盖所有要素的足够大值，取消勾选「仅查找最近的要素」选项即可。输出表包含IN_FID（输入要素ID）、NEAR_FID（邻近要素ID）、NEAR_DIST（距离值）三个核心字段，后续可通过透视表工具转换为20000×20000的矩阵格式。
• 工具默认计算两个面要素边界之间的最短平面距离，符合欧氏距离的计算要求，计算前请先将数据转换为平面投影坐标系，避免地理坐标系的球面距离误差。
• 2万要素两两组合总记录数为4亿条，计算前请预留至少10GB磁盘存储空间，也可通过ArcPy批量处理，代码框架参考：

import arcpy
arcpy.env.workspace = r"你的工作空间本地路径"
# 生成全量邻近表
arcpy.GenerateNearTable_analysis(
    in_features="你的面要素类名称",
    near_features="你的面要素类名称",
    out_table="all_near_distance_table",
    search_radius="#",
    location="NO_LOCATION",
    angle="NO_ANGLE",
    closest="ALL", # 核心参数，输出所有匹配的要素对
    closest_count="#",
    method="PLANAR" # 指定为平面欧氏距离计算
)

方案2：开源GIS实现（性能更优） #

• 可通过GeoPandas+Scipy的组合实现，计算速度比ArcGIS工具快3-5倍，内存占用可控。示例代码默认计算面形心的欧氏距离，按需可调整为面边界最短距离：

import geopandas as gpd
import numpy as np
from scipy.spatial.distance import cdist

# 读取面要素并转换为平面投影坐标系
gdf = gpd.read_file("你的面文件路径").to_crs("EPSG:对应平面坐标系的EPSG编码")
# 提取所有面形心的xy坐标
centroid_coords = np.array(list(zip(gdf.geometry.centroid.x, gdf.geometry.centroid.y)))
# 直接生成n×n的欧氏距离矩阵
dist_matrix = cdist(centroid_coords, centroid_coords, metric="euclidean")
# 导出为csv格式，也可存为npy二进制文件大幅减少存储空间
np.savetxt("distance_matrix.csv", dist_matrix, delimiter=",")

• 如果需要计算面边界的最短距离，可搭配rtree空间索引过滤候选配对，再调用geometry.distance()方法计算，避免4亿次逐对计算的性能损耗。

性能优化建议 #

• 若业务不需要无差别全量计算，可设置合理的距离阈值，仅保留小于阈值的距离对，可大幅减少计算量和存储占用。
• 20000×20000的浮点型矩阵存为csv约占用3GB空间，不需要可视化的场景建议直接存为numpy二进制.npy格式，存储空间可压缩到1GB以内。

内容的提问来源于stack exchange，提问作者bram