Ben, 

I have but no where in there do I see (or at least understand) how to draw a map on top of a current figure with set bounds in data space...

Cheers, 

Brian


On Aug 27, 2010, at 11:24 AM, Benjamin Root wrote:

On Fri, Aug 27, 2010 at 12:11 PM, Brian Larsen <balarsen@lanl.gov> wrote:
Hello all, 

I am plotting satellite orbit files and it looks really nice to plot an Earth in the center with the continents on it to orient people to where the spacecraft is.  Does anyway know how to do this?  All I seem to be able to do is create a whole globe as the figure.

In this simple example:

from pylab import *
r = ndarray(40)
r[:] = 3.3
rad = linspace(0, 2*pi, 40)
figure()
subplot(111, polar=True)
plot(rad, r, lw=3)
draw()


wouldn't it look great to have the earth plotted in the space up to 1.0 (as measured in earth radii)?

Thanks much, 

Brian



Have you checked out the basemap tool?

http://matplotlib.sourceforge.net/basemap/doc/html/

Ben Root






-- 

Brian A. Larsen
Space Science and Applications
Group ISR-1
Los Alamos National Laboratory
PO Box 1663, MS-D466
Los Alamos, NM 87545
USA

(For overnight add:
SM-30, Bikini Atoll Road)

Phone: 505-665-7691
Fax:   505-665-7395
email: balarsen@lanl.gov

Correspondence /
Technical data or Software Publicly Available