底图在纬度经度之间画线答案

【问题标题】：Basemap draw line between latitude longitude底图在纬度经度之间画线
【发布时间】：2017-08-19 08:28:15
【问题描述】：

我正在绘制龙卷风的开始/结束位置。 csv 文件的数据如下：

TouchDownLat TouchDownLong LiftoffLat LiftoffLong
31.53         -97.15       31.74      -96.88
46.45         -100.67      46.67      -100.47
43.1          -83.85       43.17      -83.42

等等……

我所做的是获取每个纬度和经度并将其分成一个 numpy 数组，如下所示：

import matplotlib.pyplot as plt
import csv
import numpy as np
from mpl_toolkits.basemap import Basemap
with open(fname, 'rb') as f:
    w = csv.reader(f, delimiter = ',')
    for i, line in enumerate (w):
      if i == 0 or line[2][0:4] not in str(range(2007,2018)):
          pass
      else:
          lat_td.append(line[27])
          long_td.append(line[28])
          lat_lift.append(line[29])
          long_lift.append(line[30])

touchdown = np.array([[lat_td], [long_td]])
lift = np.array([[lat_lift], [long_lift]])

对于底图，我正在查找所有纬度/经度的最大/最小值，以便制作地图以删除没有龙卷风的州（例如：我不想看到加利福尼亚查看 EF-5 龙卷风位置）

m = Basemap(projection = 'merc', llcrnrlat=float(min(lat_td)) - 2,\
    urcrnrlat=float(max(lat_lift)) + 2, llcrnrlon=float(max(long_td)) - 2,\
    urcrnrlon=float(min(long_lift)) + 2,lat_ts=40,resolution='l')
m.drawcoastlines()
m.fillcontinents(color='white')
m.drawmapboundary(fill_color='white')
m.drawstates(color='black')
m.drawcountries(color='black')
plt.title("#wedgez")

那么，问题来了：我正在寻找使用 numpy 数组中的纬度/经度来绘制龙卷风的轨迹。我该怎么做呢？

【问题讨论】：

不是原始问题的答案；您可以使用genfromtxt 加载 csv 文件。链接到文档here。

标签： python arrays list numpy matplotlib-basemap

【解决方案1】：

如果我理解正确，您正在寻找一种绘制路径的方法，给定路径点的坐标。然后，您可以按如下方式绘制路径：

m = Basemap(projection = 'merc', llcrnrlat=float(min(lat_td)) - 2,\
    urcrnrlat=float(max(lat_lift)) + 2, llcrnrlon=float(max(long_td)) - 2,\
    urcrnrlon=float(min(long_lift)) + 2,lat_ts=40,resolution='l')

lat = [the list of lat coordinates here] 
lon = [the list of lon coordinates here] 

x, y = m(lon, lat)
m.plot(x, y, 'o-', markersize=5, linewidth=1) 

m.drawcoastlines()
m.fillcontinents(color='white')
m.drawmapboundary(fill_color='white')
m.drawstates(color='black')
m.drawcountries(color='black')
plt.title("#wedgez")
plt.show()

【讨论】：