2017-03-26 3 views
1

Ich zeichne Anfang/Ende Standorte für Tornados. Die CSV-Datei hat die Daten, so dass:Grundkarte zeichnen Linie zwischen Breite Länge

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 

etc ...

Was ich wurde jede geografische Breite und Länge getan genommen und getrennt in eine numpy Array als solche:

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]]) 

Für die Grundkarte, finde ich die Max/Min für alle Lat/Long, so dass es eine Karte erstellt, um die Staaten, die keinen Tornado hatten, auszuschneiden (Beispiel: Ich will Kalifornien nicht sehen, wenn ich es betrachte EF-5 Tornado-Standorte)

So, jetzt kommt die Frage: Ich suche die Spur des Tornados mit dem lat/long in der numply-Array. Wie würde ich das machen?

+0

Keine Antwort auf die ursprüngliche Frage; Sie können CSV-Dateien mit '' 'genfromtxt''' laden. Link zur Dokumentation [hier] (https://docs.scipy.org/doc/numpy/reference/generated/numpy.genfromtxt.html#numpy.genfromtxt). – Niklas

Antwort

3

Wenn ich richtig verstehe, suchen Sie nach einer Möglichkeit, einen Pfad zu zeichnen, gegeben die Koordinaten der Pfadpunkte. Dann können Sie den Pfad wie folgt zeichnen:

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() 
Verwandte Themen