使用Python /底圖不在圖上繪製的線

Question

對不起，如果這個問題很簡單，我是使用Python和底圖的新手。無論如何，我試圖在地圖上繪製20條颶風的路徑（圖表）。地圖本身和圖例完美展現，但颶風的路徑卻沒有。此外，我沒有收到任何回溯消息，但我想我有我的問題可能在哪裏的想法。有人能告訴我我哪裏出錯了。

下面是CSV文件的樣本：

Year, Name, Type, Latitude, Longitude 
1957,AUDREY,HU, 21.6, 93.3 
1957,AUDREY,HU,22.0, 93.4 
1957,AUDREY,HU,22.6, 93.5 
1969,AUDREY,HU,28.2,99.6 
1957,AUDREY,HU,26.5,93.8 
1957,AUDREY,HU,27.9,93.8 
1957,AUDREY,HU,29.3,95 
1957,AUDREY,HU,27.9,93.8 
1957,AUDREY,HU,29.3,93.8 
1957,AUDREY,HU,30.7,93.5 
1969,CAMILLE,HU, 21.6,99.3 
1969,CAMILLE,HU,22.0,98.4 
1969,CAMILLE,HU,22.6,90.5 
1969,CAMILLE,HU,23.2,93.6 

這裏是我到目前爲止的代碼：

import numpy as np 
from mpl_toolkits.basemap import Basemap 
import matplotlib.pyplot as plt 
import csv, os, scipy 
import pandas 
from PIL import * 


data = np.loadtxt('louisianastormb.csv',dtype=np.str,delimiter=',',skiprows=1) 
'''print data''' 
fig = plt.figure(figsize=(12,12)) 

ax = fig.add_axes([0.1,0.1,0.8,0.8]) 

m = Basemap(llcrnrlon=-100.,llcrnrlat=0.,urcrnrlon=-20.,urcrnrlat=57., 
      projection='lcc',lat_1=20.,lat_2=40.,lon_0=-60., 
      resolution ='l',area_thresh=1000.) 

m.bluemarble() 
m.drawcoastlines(linewidth=0.5) 
m.drawcountries(linewidth=0.5) 
m.drawstates(linewidth=0.5) 

# Creates parallels and meridians 
m.drawparallels(np.arange(10.,35.,5.),labels=[1,0,0,1]) 
m.drawmeridians(np.arange(-120.,-80.,5.),labels=[1,0,0,1]) 
m.drawmapboundary(fill_color='aqua') 
color_dict = {'AUDREY': 'red', 'ETHEL': 'white', 'BETSY': 'yellow','CAMILLE': 'blue', 'CARMEN': 'green','BABE': 'purple', } 


colnames = ['Year','Name','Type','Latitude','Longitude'] 
data = pandas.read_csv('louisianastormb.csv', names=colnames) 
names = list(data.Name) 
lat = list(data.Latitude) 
long = list(data.Longitude) 
colorName = list(data.Name) 
#print lat 
#print long 
lat.pop(0) 
long.pop(0) 
colorName.pop(0) 
latitude= map(float, lat) 
longitude = map(float, long) 
x, y = m(latitude,longitude) 
#Plots points on map 
for colorName in color_dict.keys(): 
    plt.plot(x,y,linestyle ='-',label=colorName,color=color_dict[colorName], linewidth=5) 
    lg = plt.legend() 
    lg.get_frame().set_facecolor('grey') 
plt.show() 

Answer 1

兩（好吧，我撒謊了，應該有）與您的代碼的問題

我，你輸入的經度應該是否定的是你的底圖定義的邊界內，所以後轉換之前添加此x和y

longitude = [-i for i in longitude] 

二，你的座標轉換線是不對的，你應該換LON和LAT參數列表

x, y = m(longitude, latitude) 

代替

x, y = m(latitude,longitude) 

編輯：

還好，第二個問題是OP張貼在評論，請檢查以下內容的完整代碼，並請注意我做了比較，你的變化

# Last-modified: 21 Oct 2013 05:35:16 PM 

import numpy as np 
from mpl_toolkits.basemap import Basemap 
import matplotlib.pyplot as plt 
import csv, os, scipy 
import pandas 
from PIL import * 


data = np.loadtxt('louisianastormb.csv',dtype=np.str,delimiter=',',skiprows=1) 
'''print data''' 
fig = plt.figure(figsize=(12,12)) 

ax = fig.add_axes([0.1,0.1,0.8,0.8]) 

m = Basemap(llcrnrlon=-100.,llcrnrlat=0.,urcrnrlon=-20.,urcrnrlat=57., 
      projection='lcc',lat_1=20.,lat_2=40.,lon_0=-60., 
      resolution ='l',area_thresh=1000.) 

m.drawcoastlines(linewidth=0.5) 
m.drawcountries(linewidth=0.5) 
m.drawstates(linewidth=0.5) 
# m.bluemarble(ax=ax) 

# Creates parallels and meridians 
m.drawparallels(np.arange(10.,35.,5.),labels=[1,0,0,1]) 
m.drawmeridians(np.arange(-120.,-80.,5.),labels=[1,0,0,1]) 
m.drawmapboundary(fill_color='aqua') 
color_dict = {'AUDREY': 'red', 'ETHEL': 'white', 'BETSY': 'yellow','CAMILLE': 'blue', 'CARMEN': 'green','BABE': 'purple', } 

colnames = ['Year','Name','Type','Latitude','Longitude'] 
data = pandas.read_csv('louisianastormb.csv', names=colnames) 
names = list(data.Name) 
lat = list(data.Latitude) 
long = list(data.Longitude) 
colorNames = list(data.Name) 
#print lat 
#print long 
lat.pop(0) 
long.pop(0) 
colorNames.pop(0) 
latitude= map(float, lat) 
longitude = map(float, long) 
# added by nye17 
longitude = [-i for i in longitude] 
# x, y = m(latitude,longitude) 
x, y = m(longitude,latitude) 
# convert to numpy arrays 
x = np.atleast_1d(x) 
y = np.atleast_1d(y) 
colorNames = np.atleast_1d(colorNames) 
#Plots points on map 
for colorName in color_dict.keys(): 
    plt.plot(x[colorName == colorNames],y[colorName == colorNames],linestyle ='-',label=colorName,color=color_dict[colorName], linewidth=5) 
    lg = plt.legend() 
    lg.get_frame().set_facecolor('grey') 
plt.show() 

Answer 2

我覺得你的困難與其說是在底圖作爲繪圖。您不需要繪製整個x/y數據集，而需要查找與颶風Z相對應的x/y點。然後僅繪製某些顏色的點c。然後找到相應的下一個颶風等點...

下面，雖然沒有使用底圖數據結構應提供一個起點繪製基於某些選擇器向量點的子集。

#given a list of x,y coordinates with a label we'll plot each line individually 

#first construct some points to plot 
x1 = [1,1.1,1.2,1.3, 2.0,2.2,2.3, 4,3.9,3.8,3.7] 
y1 = [5,5.1,5.2,5.3, 6.0,6.2,6.3, 2,2.1,2.2,2.3] 
pointNames = [] 
#generate some labels 
pointNames.extend(['a']*4) 
pointNames.extend(['b']*3) 
pointNames.extend(['c']*4) 

#make things easy by casting to numpy arrays to allow for easier indexing 
x1 = numpy.array(x1) 
y1 = numpy.array(y1) 
pointNames = numpy.array(pointNames) 

for elem in ['a','b','c']: 
    selector = pointNames==elem 
    subsetX = x1[selector] 
    subsetY = y1[selector] 
    #now plot subsetX vs subsetY in color Z 
    plot(subsetX,subsetY,'*-') 

show()