使用ArcGIS制作城市关系强度图

2020.9.28 我又写了篇计算城市引力的文章，点击这里查看。

规划或地理等相关专业的同学可能会遇到城市关系强度图的制作。我们主要介绍制图部分，至于城市间关系强度如何计算，请左转自行查文献或相关资料学习。

ArcMap中有个XY转线（XY to Line）的GP工具，输入两对XY坐标，就可以生成线了。

那我们计算好Excel表中城市A、城市B的坐标放在一张表中，然后使用XY转线工具生成线就可以了。

点击这里下载示例数据。

操作步骤如下：

Step1：计算各城市的XY坐标。

使用GP工具添加XY坐标（Add XY Coordinate），计算点数据的坐标。

打开属性表后，可以看到增加了POINT_X和POINT_Y两个字段。

Step2：将Excel导入到地理数据库中。再分别加入城市A、城市B的XY坐标。

使用GP工具Excel转表（Excel to Table）。

使用GP工具连接字段（Join Field）连接城市A的XY坐标。

打开城市关系强度表的属性表，结果如下：

再添加上城市B的XY坐标。将输入连接字段换成城市B，结果如下：

Step3：进行城市间的连线。

打开GP工具XY转线（XY to Line），输入起终点对应的字段。要保留城市间的强度字段，所以ID（可选中）输入城市强度。

说下测地线类型，ArcGIS提供了以下四种：

GEODESIC（测地线）： 最准确地表示地球表面任意两点之间的最短曲线。
GREAT_CIRCLE（大圆）： 表示地球表面与通过地心的平面的相交线上任意两点之间的路径。下图中的红色线是PQ两点点的大圆距离。

RHUMB_LINE（等角航线）： 表示通过以极点为起点的等方位角所定义的椭球体表面上的任意两点之间的路径。等角航线在墨卡托投影中显示为直线。

NORMAL_SECTION（法截弧）： 表示由椭球体表面与通过椭球体表面上两点并垂直于两点起点处椭球面的平面相交而定义的椭球面上任意两点之间的路径。因此，从 A 点到 B 点与从 B 点到 A 点的法向截面线不同。法截弧是测地线的简版，容易构建，但不是两点间的最短距离。

当研究范围比较小时，使用Geodesic，显示的也是直线。将结果进行分级符号进行渲染，再进行排版后，如下图：

最后再说下类似的案例，依然是研究城市关系强度图。假如两个地市的人合伙开公司，那么就说明这两个地市有关系。

公司名	合伙人来自
A公司	吉林、白山
B公司	四平、吉林、白城
C公司	吉林、白城

组织表格的内容如下，可以看到，B公司和C公司都有来自吉林、白城的合伙人，说明这两个城市合作可能更多一些。

| 公司名 | 合伙人来自 | 合伙人1 | 合伙人2
|------------|---------|
|A公司 | 吉林、白山| 吉林|白山|
|B公司 | 四平、吉林|四平|吉林|
|B公司 | 四平、白城|四平|白城|
|B公司 | 吉林、白城|吉林|白城|
|C公司 | 吉林、白城|吉林|白城|

可以使用GP工具汇总统计（Summary Statistics），对合伙人来自字段统计总和，然后再按照文中的操作步骤进行计算。

这个FRENQUENCY字段可以作为强度字段，进行渲染。

欢迎关注我的公众号~

帮助及参考文档：

XY转线
测地线
大圆
等角航线
法截弧
汇总统计