2017-03-26 18 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") 

これで質問が出ます: numpy配列のlat/longを使って竜巻の軌跡をプロットしたいと思っています。私はどうしたらいいのですか?

+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()