如何在python中的地理数据框架图上叠加箭图

import geopandas as gpd
import matplotlib.pyplot as plt

gdf = gpd.read_file('your_shapefile.shp')  # 替换为你的地理数据框架图文件路径

# 假设你有一些箭头的起点和终点坐标，可以将它们存储为两个列表
start_points = [(x1, y1), (x2, y2), ...]  # 箭头起点坐标
end_points = [(x1, y1), (x2, y2), ...]  # 箭头终点坐标

# 将起点和终点坐标转换为Geopandas的GeoDataFrame
arrows = gpd.GeoDataFrame({'geometry': [LineString([start, end]) for start, end in zip(start_points, end_points)]})

fig, ax = plt.subplots()
gdf.plot(ax=ax)  # 绘制地理数据框架图
arrows.plot(ax=ax, color='red', linewidth=1)  # 绘制箭头数据集

plt.title('Arrow Overlay on Geospatial Data')  # 添加标题
plt.legend(['Geospatial Data', 'Arrows'])  # 添加图例

plt.savefig('arrow_overlay.png')  # 保存图像
plt.show()  # 显示图像