On Mar 17, 4:06 pm, David Fanning <n...@dfanning.com> wrote:
> David Fanning writes:
>> I suppose I could use the navigation data to establish
>> a limit for a map projection, but what map projection?
>> Satellite? Guess I could try it.
>
> OK, well, I'm in the ballpark with a satellite projection,
> but I'm guessing my satellite is not sitting directly
> over the center of my image. Where is it? How could I
> find out?
>
> ; IDL> help, image
> ; IMAGE BYTE [1200, 862]
>
> ; Corner points.
> ;IDL> print, peru_lat[0,0], peru_lat[0, 861],
> ; peru_lat[1199, 861], peru_lat[1199,0]
> ; -19.4763 12.9533 12.9516 -19.4735
> ;IDL> print, peru_lon[0,0], peru_lon[0, 861],
> ; peru_lon[1199, 861], peru_lon[1199,0]
> ; -99.3065 -98.2939 -51.8151 -50.8084
>
> ; right, top, left, bottom of image for LIMIT.
> ;IDL> print, peru_lat[0,431], peru_lat[600, 861],
> ; peru_lat[1199, 431], peru_lat[600,0]
> ; -3.03467 12.7686 -3.03430 -19.1762
> ;IDL> print, peru_lon[0,431], peru_lon[600, 861],
> ; peru_lon[1199, 431], peru_lon[600,0]
> ; -97.5755 -75.0291 -52.5294 -75.0302
>
> IDL> Window, XSIZE=1200, YSIZE=862
> IDL> Map_Set, -3.03449, -75.0296, /Satellite, SAT_P=[6.5,0.0,0.0], $
> XMargin=0, YMargin=0, $
> LIMIT=[-3.0347, -97.5755, 12.7686, -75.0291, $
> -3.0343, -52.5294, -19.1762, -75.0302]
> IDL> TVIMAGE, image
> IDL> Map_CONTINENTS
> IDL> MAP_GRID
>
> Folks, you *know* being just a little bit off drives me
> to crazy. Please help. I'm completely out of Xanax.
> All I have left is a box of Viagra, left over from
> Coyote's last tete-a-tete. :-(
>
> Cheers,
>
> David
> --
> David Fanning, Ph.D.
> Fanning Software Consulting, Inc.
> Coyote's Guide to IDL Programming:http://www.dfanning.com/
> Sepore ma de ni thui. ("Perhaps thou speakest truth.")
Provide all for corners as 8 numbers in LIMIT will make some
improvements.
From my experience, remapping provides the best results.
Kelly Dean
Fort Collins, CO
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
