IndexError’ too many indices’ in basemap contour
Hi!
I try to plot some interpolated data on a map and get an error saying there are too many indices. When I use contour in matplotlib without basemap I don’t get the error. Also the map without a contour plot on it works. Maybe some of you know what I do wrong?
Here is my code:
···
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.mlab import griddata
from mpl_toolkits.basemap import Basemap
yi = np.linspace(min(Lat),max(Lat),100)
xi = np.linspace(min(Lon),max(Lon),50)
lon,lat,var = np.array(Lon), np.array(Lat), np.array(Var)
zi = griddata(lon,lat,var,xi,yi)
map = Basemap(projection=‘cyl’, llcrnrlat=67, urcrnrlat=73,\
llcrnrlon=4, urcrnrlon=20, resolution='h')
map.drawcoastlines()
map.drawmeridians(np.arange(0,360,1))
map.drawparallels(np.arange(-90,90,1))
xlon, ylat = map(xi,yi)
cs = map.contour(xlon,ylat,zi)
plt.show()
and this is the error message:
cs = map.contour(xlon,ylat,zi)
File “/usr/local/lib/python2.6/dist-packages/mpl_toolkits/basemap/init.py”, line 2820, in contour
xx = x[x.shape[0]/2,:]
IndexError: too many indices
Any help appreciated!
Sünnje