comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: Finding map SCALE given latitude.
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Switch to threaded view of this topic Create a new topic Submit Reply
Re: Finding map SCALE given latitude. [message #16196] Wed, 07 July 1999 00:00
Liam Gumley is currently offline  Liam Gumley
Messages: 473
Registered: November 1994
Senior Member
Grady Daub wrote:
> I'm making orbit plots. There will be 2 plots, one each of the north and
> south poles.
>
> MAP_SET,/cont,color=0,-90,0,0,pos=[.1,.75,.9,.90],scale=100e 6
>
> would show you what I mean. (Using a window shaped like a sheet of
> paper.)
>
> Is there an easy (or complicated, for that matter) way to determine the
> SCALE= factor to set for the plot to include only up to a certain
> latitude circle?
>
> What I need is a rectangular box (centered on the pole), about 20
> degrees tall and exactly 90 degrees wide (those in latitude). The data
> will be plotted in a similar box, below the map plot, with the data
> corresponding to the orbit position directly above it.
>
> This all could probably be done with LIMIT= , but, then I'd have to
> worry about aspect ratio's. It seems simpler to make a box with
> POSITION= , then scale the earth until it shows the +/- 45 degrees
> latitude circle at the left/right edges. (top/bottom range doesn't
> really matter.)

When you use the SCALE keyword to MAP_SET, the only things defining the
map coverage are the size of your graphics window, and the POSITION
keyword if it is used. Consider the following examples:

window, /free, xsize=300, ysize=300
map_set, 25, -81, scale=10e6, /cont

window, /free, xsize=600, ysize=600
map_set, 25, -81, scale=10e6, /cont

window, /free, xsize=600, ysize=600
map_set, 25, -81, scale=10e6, /cont, position=[0.1, 0.1, 0.5, 0.5]

Notice how in each case, the scale of the map remains exactly the same
(Florida is the same size), but the extent of the mapped area changes.
So you just need to experiment with window size and POSITION, e.g.

window, /free, xsize=714, ysize=924
map_set, -90, 0, scale=50e6, pos=[0.05,0.05,0.95,0.45], /cont
map_grid, /label
map_set, 90, 0, scale=50e6, pos=[0.05,0.55,0.95,0.95], /cont, /noerase
map_grid, /label

Cheers,
Liam.

--
Liam E. Gumley
Space Science and Engineering Center, UW-Madison
http://cimss.ssec.wisc.edu/~gumley
Re: Finding map SCALE given latitude. [message #16200 is a reply to message #16196] Wed, 07 July 1999 00:00 Go to previous message
Martin Schultz is currently offline  Martin Schultz
Messages: 515
Registered: August 1997
Senior Member
Grady Daub wrote:
>
> I'm making orbit plots. There will be 2 plots, one each of the north and
> south poles.
>
> MAP_SET,/cont,color=0,-90,0,0,pos=[.1,.75,.9,.90],scale=100e 6
>
> would show you what I mean. (Using a window shaped like a sheet of
> paper.)
> [...]
> This all could probably be done with LIMIT= , but, then I'd have to
> worry about aspect ratio's.

isn't that what the /isotropic keyword in map_set is for?

Martin

--

|||||||||||||||\\\\\\\\\\\\\-------------------///////////// //|||||||||||||||
Martin Schultz, DEAS, Harvard University, 29 Oxford St., Pierce 109,
Cambridge, MA 02138 phone (617) 496 8318 fax (617) 495 4551
e-mail mgs@io.harvard.edu web http://www-as/people/staff/mgs/

********* ADDRESS CHANGE : AFTER JULY 15, 1999 ***************

Max-Planck-Institut fuer Meteorologie >>> N E W <<<
Bundesallee 55 >>> N E W <<<
20147 Hamburg >>> N E W <<<
Germany >>> N E W <<<
phone (+49 40) 41173 - 0 >>> N E W <<<

email martin.schultz@dkrz.de >>> N E W <<<
(email to mgs@io.harvard.edu will be forwarded)

************************************************************ *******
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Surface Fit request
Next Topic: Bivariate Normal Distribution Function?

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Sun Oct 12 07:02:27 PDT 2025

Total time taken to generate the page: 0.24689 seconds