Plotting basemap on top of other data

# This is a piece of code that I actually used to try to solve the
# Just assumed that I have an int32 array ("data") with satellite data
of dimensions [x,y,3] filled with RGB values ... I have removed the
reading of satellite sample data, which is not relevant here.
# The satellite (geostationary) is located above the equator at zero
degrees east
# in this example I want a plot with 300 dpi resolution.

import numpy
import matplotlib.pyplot as plt

dpi = 300
fig = plt.figure(figsize=[data.shape[0]/dpi,data.shape[1]/dpi])
ax = plt.axes([0,0,1,1], frameon = False)
im = ax.imshow(data, origin='None')

from mpl_toolkits.basemap import Basemap

map = Basemap(projection='ortho', lat_0 = 0, lon_0 = 0,
        resolution = 'l', area_thresh = 1000.)


fig.savefig('figure.png', format='png', dpi=dpi)


-----Original Message-----
From: Patrick Marsh [mailto:patrickmarshwx@…287…]
Sent: zaterdag 10 december 2011 23:41
To: Laat de, Jos (KNMI)
Subject: Re: [Matplotlib-users] Plotting basemap on top of other data


My initial thought was that you need to use the "zorder" keyword
argument and set the zorder to a large value. However, the more I
thought about it, I'm not really sure how you are plotting the
satellite data. Can you provide a code snippet?

Patrick Marsh
Ph.D. Student / Liaison to the HWT
School of Meteorology / University of Oklahoma
Cooperative Institute for Mesoscale Meteorological Studies
National Severe Storms Laboratory