2017-03-26 65 views
1

我正在绘制龙卷风的开始/结束位置。该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") 

所以,现在问题来了: 我正在寻找使用lat/long在numpy数组中的龙卷风的轨迹。我会怎么做呢?

+0

不是原始问题的答案;你可以使用'genfromtxt'''加载csv文件。链接到文档[这里](https://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html#numpy.genfromtxt)。 – Niklas

回答

3

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

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() 
相关问题