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/
Se****e ma de ni thui. ("Perhaps thou speakest truth.")
|
