Re: [Matplotlib-users] Basemap help, imshow, transform_scaler

[John <was...@gm...>]

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: