# This is a piece of code that I actually used to try to solve the
problem
# 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)
ax.set_axis_off()
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.)
map.drawcoastlines()
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)
Cc: matplotlib-users@lists.sourceforge.net
Subject: Re: [Matplotlib-users] Plotting basemap on top of other data
HI,
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?
PTM
---
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
http://www.patricktmarsh.com