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

Home » Public Forums » archive » map_set question
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
map_set question [message #6517] Tue, 09 July 1996 00:00 Go to next message
gfu is currently offline  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 Go to previous messageGo to next message
gennari is currently offline  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 Go to previous messageGo to next message
gennari is currently offline  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 Go to previous messageGo to next message
davidf is currently offline  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 #8323 is a reply to message #6517] Mon, 24 February 1997 00:00 Go to previous messageGo to next message
Andy Loughe is currently offline  Andy Loughe
Messages: 174
Registered: November 1995
Senior Member
Tim Patterson wrote:
>
> Andy Loughe wrote:
>>
>
>> 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!!!
>>
>
> Perhaps they could email out bug update reports when they
> send out those "Come and Learn IDL With US" flyers every
> few months? Even a section on current bug fixes on their
> web page would be a start.
>
> And while I'm indulging in wishful thinking, perhaps they
> could put up a list of the bug fixes made to Version 5.00 as
> well as the all the new fancy offerings they're pushing.
>
> Tim


You've missed the point.
They don't fix the bugs!

We pay good money for broken software that rarely gets fixed.
RSI is "glitz" happy. "Let's put more functionality into
our product and let the users find work-arounds to old bugs
that are too hard for us to fix. We're too busy to bother."

--
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 #8326 is a reply to message #6517] Mon, 24 February 1997 00:00 Go to previous messageGo to next message
Tim Patterson is currently offline  Tim Patterson
Messages: 65
Registered: October 1995
Member
Andy Loughe wrote:
>

> 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!!!
>

Perhaps they could email out bug update reports when they
send out those "Come and Learn IDL With US" flyers every
few months? Even a section on current bug fixes on their
web page would be a start.

And while I'm indulging in wishful thinking, perhaps they
could put up a list of the bug fixes made to Version 5.00 as
well as the all the new fancy offerings they're pushing.

Tim
Re: Map_Set Question [message #8333 is a reply to message #6517] Mon, 24 February 1997 00:00 Go to previous messageGo to next message
Andy Loughe is currently offline  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 Go to previous message
Theo Brauers is currently offline  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 Go to previous message
Craig Markwardt is currently offline  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 Go to previous message
Harald Frey is currently offline  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
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Save, /Compress (was Re: reading in binary data)
Next Topic: Re: need help on ROI subscript

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

Current Time: Wed Oct 08 19:43:35 PDT 2025

Total time taken to generate the page: 0.00570 seconds