On 8/27/10 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

Sounds like you want the orthographic projection

http://matplotlib.sourceforge.net/basemap/doc/html/users/ortho.html
http://www.scipy.org/Cookbook/Matplotlib/Maps

or the 'near-sided perspective' projection

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


-Jeff
-- 
Jeffrey S. Whitaker         Phone  : (303)497-6313
Meteorologist               FAX    : (303)497-6449
NOAA/OAR/PSD  R/PSD1        Email  : Jeffrey.S.Whitaker@noaa.gov
325 Broadway                Office : Skaggs Research Cntr 1D-113
Boulder, CO, USA 80303-3328 Web    : http://tinyurl.com/5telg