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

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

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

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

点击这里下载示例数据。  

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

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

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

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

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

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

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

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

GEODESIC（测地线）： 最准确地表示地球表面任意两点之间的最短曲线。

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

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

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

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

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

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

这个FRENQUENCY字段可以作为强度字段，进行渲染。   欢迎关注我的公众号~

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