ArcGIS 航线 OD 图制作教程

标签：教程 fsub air1 OD list ArcGIS air2 lat 数据

今天我将和大家分享一下通过 ArcGIS 制作航线 OD 图的方法，文章内容包括数据预处理、地图制图及资源分享三部分，所需数据包括机场数据及航班数据。

数据预处理

在获取到数据后需要对其进行整理，提取出各航线起点及终点对应的经纬度坐标。这部分操作可通过编写 Python 代码完成，示例代码如下：

点击查看代码

# 导入所需模块
import pandas as pd
import openpyxl

# 获取机场及航班数据
airports = pd.read_csv("data/airports.csv")
flights = pd.read_csv("data/flights.csv")

# 由于数据量较大，只选取airline列值为"AA"的航班数据
fsub = flights[flights["airline"] == "AA"]
# 重置航班数据表的行索引，使其从零开始顺序排序，以便于后续数据处理
# 代码参考：https://blog.csdn.net/weixin_51098806/article/details/115265573
fsub = fsub.reset_index(drop=True)

air12_id_list = list() # 起始点列表，数据格式："起点 - 终点"
air1_lon_list, air1_lat_list = list(), list() # 起点经纬度
air2_lon_list, air2_lat_list = list(), list() # 终点经纬度

# 对航班数据进行迭代，以提取各航班的起点与终点坐标
for i in range(len(fsub)):
    fsub_i = fsub.loc[i]  # 获取航班表里的第i行数据

    # 获取航班起点及终点的机场数据
    air1 = airports[airports["iata"] == fsub_i["airport1"]]
    air2 = airports[airports["iata"] == fsub_i["airport2"]]

    # 将获取到的机场数据导入对应列表中
    air1_id, air2_id = fsub_i["airport1"], fsub_i["airport2"] 
    air12_id_list.append(f"{air1_id} - {air2_id}")
    air1_lon_list.append(float(air1["long"]))
    air1_lat_list.append(float(air1["lat"]))
    air2_lon_list.append(float(air2["long"]))
    air2_lat_list.append(float(air2["lat"]))

# 将上方的列表合并后导出为Excel格式
data = {"air12_id": air12_id_list, "air1_lon": air1_lon_list, "air1_lat": air1_lat_list,
        "air2_lon": air2_lon_list, "air2_lat": air2_lat_list}
df = pd.DataFrame(data)

# 输出Excel数据
# 输出时须确保输出路径准确无误，如果输出到了一个不存在的文件夹，运行时就会出现报错
# 如果输出数据表本就存在于输出文件夹中，则会直接覆盖源文件，但覆盖过程需要源文件保持关闭状态
df.to_excel("data/data_AA.xlsx")
print("输出成功！")

整理后的部分 Excel 表如下图所示：

注：输出表格的字体、行高及列宽等格式有所调整。

2 地图制图

2.1 导入数据

以世界陆地数据作为底图。

由于先前整理输出的是 xlsx 格式表格，无法直接被 ArcGIS 读取，故需将其转换为 xls 格式后再导入 ArcGIS 中或者通过 ArcGIS 的【Excel 转表】（ArcToolbox：Conversion Tools/Excel/Excel to Table）工具将其转换为 ArcGIS 可正常读取的表格文件。

注：通过【Excel 转表】工具转换时建议将结果表保存到地理数据库中，工具参数设置如下图所示：