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

Home » Public Forums » archive » Re: Map_Proj_* Problems Re: Overlay images from WMS servers
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: Map_Proj_* Problems Re: Overlay images from WMS servers [message #47458] Fri, 10 February 2006 11:05
David Fanning is currently offline  David Fanning
Messages: 11724
Registered: August 2001
Senior Member
Liam writes:

> The following works pretty much the way you would expect it to:

Yes, OK. So is my problem really just with the UTM map
projection? Could be, I suppose. I did, at one point,
fool around with other projections, but that may have
been in one of my totally confused periods. These do
seem to work properly.

Thanks, Liam.

Cheers,

David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Re: Map_Proj_* Problems Re: Overlay images from WMS servers [message #47461 is a reply to message #47458] Fri, 10 February 2006 10:52 Go to previous message
liamgumley is currently offline  liamgumley
Messages: 74
Registered: June 2005
Member
The following works pretty much the way you would expect it to:

;- Lambert Azimuthal Equal Area (Great Lakes)
map = map_proj_init(111, center_lat=45, $
center_lon=-85, limit=[40, -95, 50, -75])

;- Default position
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map

;- Default position (isotropic)
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, /isotropic
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map

;- Position specified
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, $
position=[0.1, 0.1, 0.6, 0.6]
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map

;- Position specified (isotropic)
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, $
position=[0.1, 0.1, 0.6, 0.6], /isotropic
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map
  Switch to threaded view of this topic Create a new topic Submit Reply
Previous Topic: Re: more text widget help...
Next Topic: GIF output disappeared?

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

Current Time: Wed Oct 08 15:39:50 PDT 2025

Total time taken to generate the page: 0.00411 seconds