使用Cartopy从Excel电子表格绘制纬度/经度

Question

[试图使用Cartopy在地图上绘制测量样点，Cartopy是一个图书馆，与其他图书馆相比，在线信息严重缺乏。我的经/纬度数据位于Excel电子表格的两个单独的列中。

我已经设法建立了底图，我的脚本如下：

import cartopy.feature as cf

map = plt.figure(figsize=(15,15))
ax = plt.axes(projection=ccrs.EuroPP())
ax.coastlines(resolution='10m')

ax.add_feature(cf.LAND)
ax.add_feature(cf.OCEAN)
ax.add_feature(cf.COASTLINE)
ax.add_feature(cf.BORDERS, linestyle=':')
ax.add_feature(cf.LAKES, alpha=0.5)
ax.add_feature(cf.RIVERS)

ax.stock_img()

产生下图

我还从以下电子表格中使用熊猫从Excel列中读取了数据，这些列分别标记为“纬度”（“ L”列）和“经度”（“ M”列）：https://liveplymouthac-my.sharepoint.com/:x:/g/personal/rhodri_irranca-davies_postgrad_plymouth_ac_uk/EYGoEUAO8tBLkdT0r4dwDdsBAijsVChsLCec-DxaYo2Tew?e=rs9cBc

这是使用以下代码完成的：

file = r'C:\Users\Laptop\Documents\Python\Thesis\t2_GeoData.xlsx'
df_lat = pd.read_excel(file, index_col=None, na_values=['NA'], usecols = 'L')
df_lon = pd.read_excel(file, index_col=None, na_values=['NA'], usecols = 'M')

熊猫是正确的选择吗？如果是这样，我现在如何合并创建的两个数据框（“ df_lat”和“ df_lon”），并使用cartopy将其绘制为样条？

提前感谢：）

Answer 1

我不确定您将其绘制为断面时的确切含义，但是您可以使用以下方法在CartoPy中绘制经/纬线：

ax.plot(df_lon, df_lat, linestyle='-', color='orange', transform=ccrs.PlateCarree())