[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: plot, lons, lats overlayed on a map
- To: T Bowers <tbowers(at)nrlssc.navy.mil>
- Subject: Re: plot, lons, lats overlayed on a map
- From: Liam Gumley <Liam.Gumley(at)ssec.wisc.edu>
- Date: Mon, 25 Jan 1999 16:42:08 -0600
- Newsgroups: comp.lang.idl-pvwave
- Organization: Space Science and Engineering Center, University of Wisconsin-Madison
- References: <78iaee$g01$1@filet.nrlssc.navy.mil> <36ACDB45.A741D494@ssec.wisc.edu> <78ioda$hb3$1@filet.nrlssc.navy.mil>
- Xref: news.doit.wisc.edu comp.lang.idl-pvwave:13433
T Bowers wrote:
> Hmmm... it works! But...I've got the problem that the next data station(s)
> may be10 deg. south of this cluster and I need to make sure that when the
> updated set of lats and lons comes through (I just append to the arrays with
> lats = [lats, newLat]
> lon = [lons, newLon]
> ) that the new plot will adjust automatically to encompass ALL the points.
Don't wrestle with the LIMIT keyword - it's too painful. Just pick a
SCALE that is large enough to cover all possibilities. Note that no
matter what your window size, using SCALE creates a map at the same
scale, e.g. try
window, /free, xsize=400, ysize=400
map_set, 35, 125, /ortho, scale=10e6, /cont
window, /free, xsize=800, ysize=800
map_set, 35, 125, /ortho, scale=10e6, /cont
So if you need to, just make a larger image window. And use the mean
value for LAT and LON to center the projection, e.g.
latmean = total( lat ) / float( n_elements( lat ) )
lonmean = total( lon ) / float( n_elements( lon ) )
map_set, latmean, lonmean, /ortho, scale=10e6, /cont
> Thanks, Liam.
No worries.
---
Liam E. Gumley
Space Science and Engineering Center, UW-Madison
1225 W. Dayton St., Madison WI 53706, USA
Phone (608) 265-5358, Fax (608) 262-5974
http://cimss.ssec.wisc.edu/~gumley