Jeff,
Thanks for the quick reply, below is my plotting code. Here are the answers
to your question about my arrray:
>>> type(Zdat); type(zdat)
<type 'numpy.ndarray'>
<type 'numpy.ndarray'>
>>> shape(Zdat); shape(zdat)
(180, 360)
(596, 596)
>>> shape(lons); shape(x)
(360,)
(596, 596)
>>> shape(lats); shape(y)
(180,)
(596, 596)
I would like to use contourf, but ultimately I want to plot the log(zdat),
and since my data have so many null or zero values, I can't get that to work
with contourf, whereas imshow seems to handle it.
Thanks!!
PLOTTING FUNCTION HERE:
def fp_plot(datain):
from matplotlib.toolkits.basemap import Basemap, shiftgrid
Zdat=datain;
n,m=N.shape(Zdat)
lons=N.array([i-179.5 for i in range(m)])
lats=N.array([i-89.5 for i in range(n)])
print "The array is bounded: %s : %s lat, %s : %s lon " % (N.min(lats),
N.max(lats),N.min(lons),N.max(lons))
# Set up basmap.
m = Basemap(llcrnrlon=-180.,llcrnrlat=-90.,urcrnrlon=180.,urcrnrlat=90.0
,\
rsphere=(6378137.00,6356752.3142),\
resolution='l',area_thresh=1000.,projection='npstere',\
lat_1=80.,lon_0=0., boundinglat=40.)
# transform to nx x ny regularly spaced native projection grid
nx = int((m.xmax-m.xmin)/20000.)+1; ny = int((m.ymax-m.ymin)/20000.)+1
zdat,x,y = m.transform_scalar(Zdat,lons,lats,nx,ny,returnxy=True)
# create the figure.
fig=figure(figsize=(12,8))
# add an axes, leaving room for colorbar on the right.
ax = fig.add_axes([0.1,0.1,0.7,0.7])
# plot image over map with imshow. Want to use this, but no success yet.
anorm=normalize(.015*N.max(zdat),0.95*N.max(zdat));
# Draw flexpart output
m.imshow(x,y,zdat)
# Draw map components
m.drawcoastlines()
m.drawcountries()
m.drawstates()
meridians=arange(-170.,180.,20.)
m.drawmeridians(meridians,labels=[1,1,0,0])
paralles=arange(-85.,85.,5.)
m.drawparallels(paralles, labels=[0,0,0,1]); show();
Now, to answer your question, Zdat, or datain is a numpy array with the
following characteristics:
-------------------------------------------------------------------------
This SF.net email is sponsored by: Splunk Inc.
Still grepping through log files to find problems? Stop.
Now Search log events and configuration files using AJAX and a browser.
Download your FREE copy of Splunk now >> http://get.splunk.com/
_______________________________________________
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
