map_set question [message #6517] |
Tue, 09 July 1996 00:00  |
gfu
Messages: 10 Registered: October 1995
|
Junior Member |
|
|
Hi,
Does anyone know how to decide the default paramter value LIMIT that is
used in MAP_SET so that if I use it that will generate the same result as
the one that I don't use it. For example, how to find 'limit' so that the
following two commands will generate the same results.
MAP_SET, 60, 0, /ORTHO, /GRID, /CONTINENT
and
MAP_SET, 60, 0, /ORTHO, /GRID, /CONTINENT, $
LIMIT=[latmin, lonmin, latmax, lonmax]
Thanks.
Gary
|
|
|
Re: Map_Set Question [message #8252 is a reply to message #6517] |
Fri, 21 February 1997 00:00   |
gennari
Messages: 25 Registered: April 1994
|
Junior Member |
|
|
David Fanning (davidf@dfanning.com) wrote:
: >
: > map_set,0,180,limit=[-20,-280,55,-100],/cont
: >
: Humm, I don't find Australia missing in my view. Have you downloaded
: the map fixes from the RSI anonymous ftp site? Go to gateway.rsinc.com.
: I don't know what subdirectory they are in off-hand, but you should
: be able to find them.
: David
Indeed, the patch fixes a problem dating back to 1995 and three IDL
versions ago ( 4.0.1a, 4.0.1b, 4.0.1c).
ftp://gateway.rsinc.com/pub/idl_4.0.1/updates/map_fix/map_fi x.tar.Z
Shouldn't this have been included in the latest CDROM releases?
The "help yourself to the patch when you stumble on this problem"
approach is not what I had in mind when paying the IDL support
fees This should be included in all CDROM distributions of IDL.
I noticed another problem though using map_continents. Try
filling the continents using this.
map_set,0,180,limit=[-20,-280,55,-100],/cont
map_continents,/fill
This doesn't work. Australia and other Pacific countries get filled
in but then parts of the northeast Pacific Ocean get filled in too.
Neat. Sun Solaris 2.5 with IDL 4.0.1c
Does anyone else get the same output?
And RSI charges us outragous support fees every year....
Scott Gennari
------------------------------------------------------------ ------
University of Hawaii voice (808) 956 5392
Dept. of Information & Computer Sciences fax (808) 956 9399
2565 The Mall, Keller 304A
Honolulu, HI 96822 email gennari@Hawaii.Edu
|
|
|
Re: Map_Set Question [message #8256 is a reply to message #6517] |
Fri, 21 February 1997 00:00   |
gennari
Messages: 25 Registered: April 1994
|
Junior Member |
|
|
David Fanning (davidf@dfanning.com) wrote:
: Siri Flocke <siri@siri.acd.ucar.edu> writes:
: > map_set,0,180,limit=[-20,-280,55,-100],/cont
: >
: > Well, it kind of works, except the outline for Indonesia and Australia are
: > missing. I do get outlines for some small islands in the South Pacific, but
Hi,
map_set,0,180,limit=[-20,-280,55,-100],/cont
produced the same results as Siri's, no south pacific nations. I just
received the latest version on ENVI (2.6) and IDL (4.0.1c) too. I'll
look at the patch that David mentioned but I remember installing a
similar patch about THREE versions ago. This should have been fixied
long ago.
Scott Gennari
------------------------------------------------------------ ------
University of Hawaii voice (808) 956 5392
Dept. of Information & Computer Sciences fax (808) 956 9399
2565 The Mall, Keller 304A
Honolulu, HI 96822 email gennari@Hawaii.Edu
|
|
|
Re: Map_Set Question [message #8265 is a reply to message #6517] |
Fri, 21 February 1997 00:00   |
davidf
Messages: 2866 Registered: September 1996
|
Senior Member |
|
|
Siri Flocke <siri@siri.acd.ucar.edu> writes:
> I am trying to create a map from 20S to 55N and from 100E to 100W (across
the date
> line). After playing around with the limit for Map_Set, I came up with the
> following call:
>
> map_set,0,180,limit=[-20,-280,55,-100],/cont
>
> Well, it kind of works, except the outline for Indonesia and Australia are
> missing. I do get outlines for some small islands in the South Pacific, but
> nothing south of the equator along the left (western) boundary.
Humm, I don't find Australia missing in my view. Have you downloaded
the map fixes from the RSI anonymous ftp site? Go to gateway.rsinc.com.
I don't know what subdirectory they are in off-hand, but you should
be able to find them.
David
-----------------------------------------------------------
David Fanning, Ph.D.
Fanning Software Consulting
2642 Bradbury Court, Fort Collins, CO 80521
Phone: 970-221-0438 Fax: 970-221-4762
E-Mail: davidf@dfanning.com
Coyote's Guide to IDL Programming: http://www.dfanning.com
-----------------------------------------------------------
|
|
|
|
|
Re: Map_Set Question [message #8333 is a reply to message #6517] |
Mon, 24 February 1997 00:00   |
Andy Loughe
Messages: 174 Registered: November 1995
|
Senior Member |
|
|
Scott Gennari wrote:
>
> I noticed another problem though using map_continents. Try
> filling the continents using this.
>
> map_set,0,180,limit=[-20,-280,55,-100],/cont
> map_continents,/fill
>
> This doesn't work. Australia and other Pacific countries get filled
> in but then parts of the northeast Pacific Ocean get filled in too.
> Neat. Sun Solaris 2.5 with IDL 4.0.1c
>
> Does anyone else get the same output?
>
> And RSI charges us outragous support fees every year....
>
> Scott Gennari
Yep, I also have this problem with map_set.
You are preaching to the choir on this one.
IDL often behaves worse than freeware we use here.
Ya' gotta wonder how happy David Stern is about these results,
but what bothers me is that *nothing* ever gets done about it!!!
--
Andrew F. Loughe |
afl@cdc.noaa.gov
University of Colorado, CIRES Box 449 |
http://cdc.noaa.gov/~afl
Boulder, CO 80309-0449 | phn:(303)492-0707
fax:(303)497-7013
------------------------------------------------------------ ---------------
"I do not feel obliged to believe that the same God who has endowed us
with
sense, reason, and intellect has intended us to forego their use."
-Galileo
|
|
|
Re: map_set question [message #17129 is a reply to message #6517] |
Thu, 16 September 1999 00:00  |
Theo Brauers
Messages: 58 Registered: November 1997
|
Member |
|
|
Colin Peterson wrote:
>
> If anyone still has the IDL 4.0 User's Guide (1995) please look at page 19-10 to help in
> answering my question here. The picture for the Oblique
> Cylindrical Equidistant map projection interests me because I would like to see both poles.
> I'm using IDL 5.2 on UNIX so that might make a
> difference, but if I type the command on page 19-11:
>
> map_set, 0, 0, 45, LONDEL = 20, LATDEL = 20, $
> /GRID, /CONTINENT, /CYL, TITLE='...'
>
> to produce this map, it comes up with a simple cylindrical map, but the whole map actually
> rotates 45 degrees, ie the edges of the map are now 45
> degrees from their original horizontal and vertical positions. I would expect the ROT
> value of 45 to create the map on page 19-10 as the book says,
> but it does not. Can anyone explain this discrepancy in the book, and how I can create
> that map on page 19-10? Thanks.
>
> Colin
I dont have IDL 4 on my computer anymore, but I have the manual. I cannot reproduce
Figure 19-4 "oblique.." on WinNT 4.0 / IDL 5.2. Also the "transverse mercator .." one doesnt
work. There is a /TRANSVERSE_MERC keyword, but it doesnt produce Fig 19.4 (upper right).
No idea so far.
Theo
|
|
|
Re: map_set question [message #17268 is a reply to message #6517] |
Thu, 16 September 1999 00:00  |
Craig Markwardt
Messages: 1869 Registered: November 1996
|
Senior Member |
|
|
Colin Peterson <cap21@cornell.edu> writes:
>
> If anyone still has the IDL 4.0 User's Guide (1995) please look at page
> 19-10 to help in answering my question here. The picture for the
> Oblique Cylindrical Equidistant map projection interests me because I
> would like to see both poles. I'm using IDL 5.2 on UNIX so that might
> make a difference, but if I type the command on page 19-11:
>
> map_set, 0, 0, 45, LONDEL = 20, LATDEL = 20, $
> /GRID, /CONTINENT, /CYL, TITLE='...'
>
> to produce this map, it comes up with a simple cylindrical map, but the
> whole map actually rotates 45 degrees, ie the edges of the map are now
> 45 degrees from their original horizontal and vertical positions. I
> would expect the ROT value of 45 to create the map on page 19-10 as the
> book says, but it does not. Can anyone explain this discrepancy in the
> book, and how I can create that map on page 19-10? Thanks.
I have both the book, IDL 4.0 and IDL 5. Under IDL v4, the commands
in book produce the output shown in the book, ie, what appears to be
the correct output. IDL v5 does not.
This looks like a bug, although if the RSI people play lawyer-ball it
might still be interpretted as correct according to the documentation.
Craig
--
------------------------------------------------------------ --------------
Craig B. Markwardt, Ph.D. EMAIL: craigmnet@cow.physics.wisc.edu
Astrophysics, IDL, Finance, Derivatives | Remove "net" for better response
------------------------------------------------------------ --------------
|
|
|
Re: map_set question [message #17272 is a reply to message #6517] |
Thu, 16 September 1999 00:00  |
Harald Frey
Messages: 41 Registered: March 1997
|
Member |
|
|
Colin Peterson wrote:
> If anyone still has the IDL 4.0 User's Guide (1995) please look at
> page 19-10 to help in answering my question here. The picture for the
> Oblique Cylindrical Equidistant map projection interests me because I
> would like to see both poles. I'm using IDL 5.2 on UNIX so that might
> make a difference, but if I type the command on page 19-11:
>
> map_set, 0, 0, 45, LONDEL = 20, LATDEL = 20, $
> /GRID, /CONTINENT, /CYL, TITLE='...'
>
> to produce this map, it comes up with a simple cylindrical map, but
> the whole map actually rotates 45 degrees, ie the edges of the map are
> now 45 degrees from their original horizontal and vertical positions.
> I would expect the ROT value of 45 to create the map on page 19-10 as
> the book says, but it does not. Can anyone explain this discrepancy
> in the book, and how I can create that map on page 19-10? Thanks.
>
The best I could come up with on IDL 5.1 was
map_set, 45, 0., 0, LONDEL = 20, LATDEL = 20,/GRID, $
/CONTINENT, /CYL
but one of the pole ends up at the side of the plot.
Even with the limit keyword I could not shift both poles
into the position of figure 19.4.
=========================================================
Harald U. Frey
Space Sciences Lab phone: 510-643-3323
University of California fax: 510-643-2624
Berkeley, CA 94720-7450 email: hfrey@ssl.berkeley.edu
|
|
|