| Re: Overlay images from WMS servers (web mapping servers) on map projections [message #47430 is a reply to message #47426] |
Wed, 08 February 2006 07:41   |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
liamgumley@gmail.com writes:
> I encourage you to pursue the MAP_PROJ_INIT solution, since it is much
> more robust than MAP_SET. I'm not familiar with the projection you
> mention, but here is how I create a UTM projection to match up with a
> GeoTIFF image:
>
> ; Create graphics window
> window, /free, xsize=600, ysize=800
>
> ; Define map projection
> map = map_proj_init('UTM', limit=[10.15, -78.77, 29.87, -95.29],
> zone=16)
>
> ; Configure direct graphics data coordinates to match map projection
> plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], position=[0.0, 0.0, 1.0,
> 1.0], $
> /nodata, /isotropic, xstyle=5, ystyle=5, /noerase
>
> ;- Plot continent outlines
> map_continents, map=map
I'm not sure "robust" is the word I was using the other day,
Liam, when I tried to create a MAP_PROJ_INIT projection
*other* than this one for the GSHHS article I was writing.
I absolutely cannot figure out how this projection space works.
The documentation, as far as I can tell, is hopeless. I presume
you have to look elsewhere to figure out how it works. But
I haven't found the right place yet.
I wanted a map projection of the Great Lakes Region of the US,
positioned in the window according to a POSITION keyword. It
seems straightforward, but a couple of hours of work got me
exactly nowhere. :-(
Can you show me how you would do this robustly? :-)
Cheers,
David
P.S. I think the word I was using had four letters instead
of five. :-)
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
