From: Jeff W. <js...@fa...> - 2010-08-27 17:55:35
|
On 8/27/10 11:39 AM, Brian Larsen wrote: > 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 > Brian: Something like this perhaps? from pylab import * from mpl_toolkits.basemap import Basemap lons = linspace(0, 360, 361) lats = 40.*ones(len(lons)) map = Basemap(projection='ortho',lon_0=270,lat_0=90.) x,y = map(lons,lats) map.plot(x,y,color='y') map.bluemarble(scale=0.5) show() -Jeff > > On Aug 27, 2010, at 11:24 AM, Benjamin Root wrote: > >> On Fri, Aug 27, 2010 at 12:11 PM, Brian Larsen <bal...@la... >> <mailto:bal...@la...>> 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: bal...@la... <mailto:bal...@la...> > > Correspondence / > Technical data or Software Publicly Available > > > > > > ------------------------------------------------------------------------------ > Sell apps to millions through the Intel(R) Atom(Tm) Developer Program > Be part of this innovative community and reach millions of netbook users > worldwide. Take advantage of special opportunities to increase revenue and > speed time-to-market. Join now, and jumpstart your future. > http://p.sf.net/sfu/intel-atom-d2d > > > _______________________________________________ > Matplotlib-users mailing list > Mat...@li... > https://lists.sourceforge.net/lists/listinfo/matplotlib-users -- Jeffrey S. Whitaker Phone : (303)497-6313 Meteorologist FAX : (303)497-6449 NOAA/OAR/PSD R/PSD1 Email : Jef...@no... 325 Broadway Office : Skaggs Research Cntr 1D-113 Boulder, CO, USA 80303-3328 Web : http://tinyurl.com/5telg |