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

Home » Public Forums » archive » GeoTiff- grid
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: GeoTiff- grid [message #58257 is a reply to message #58167] Wed, 23 January 2008 05:08 Go to previous messageGo to previous message
robinson.inj is currently offline  robinson.inj
Messages: 32
Registered: August 2007
Member
On Jan 22, 5:31 pm, David Fanning <n...@dfanning.com> wrote:
> Well, I tell you what. Here are the numbers I get when I use a
> LAT=20 and LON=-96.5:
>
> Corner Coordinates:
> Upper Left:  (-167346.019, 3677892.560)
> Lower Left:  (-167346.019, 2982923.485)
> Upper Right: (1222592.132, 3677892.560)
> Lower Right: (1222592.132, 2982923.485)
> Corners in Lon/Lat
> Upper Left:  (-98.555, 51.338)
> Lower Right: (-82.367, 45.408)
>
> Here are the values I get (surprisingly close to
> ENVI's) when I use a LON of 2, instead of 20:
>
> Corner Coordinates:
> Upper Left:  (-167346.019, 3677892.560)
> Lower Left:  (-167346.019, 2982923.485)
> Upper Right: (1222592.132, 3677892.560)
> Lower Right: (1222592.132, 2982923.485)
> Corners in Lon/Lat
> Upper Left:  (-98.037, 33.346)
> Lower Right: (-85.323, 27.832)
>
> Are you sure you did the ENVI numbers correctly?
>
> Here is the way I set up the map projection:
>
>     map = MAP_PROJ_INIT(9, $
>        DATUM=5, $
>        CENTER_LAT=20.0, $
>        CENTER_LON=-96.5)
>     lonlat = MAP_PROJ_INVERSE([xorigin, xend],[yend, yorigin], $
>         MAP_STRUCTURE=map)
>
> Cheers,
>
> David
>
> --
> David Fanning, Ph.D.
> Fanning Software Consulting, Inc.
> Coyote's Guide to IDL Programming (www.dfanning.com)
> Sepore ma de ni thui. ("Perhaps thou speakest truth.")

strange. My map does not have any lat 2 value. converting the tiff
image in a image with header i got:

ENVI
description = {
File Resize Result, x resize factor: 1.000000, y resize factor:
1.000000.
[Thu Jan 17 15:57:16 2008]}
samples = 3000
lines = 1500
bands = 1
header offset = 0
file type = ENVI Standard
data type = 2
interleave = bsq
sensor type = Unknown
byte order = 0
x start = 4380
y start = 1880
map info = {GEOTIFF (Mercator), 1.0000, 1.0000, -167346.0185,
3677892.5566, 4.6331271653e+002, 4.6331271653e+002, WGS-84,
units=Meters}
projection info = {20, 6378137.0, 6356752.3, 20.000000, -96.500000,
0.0, 0.0, WGS-84, GEOTIFF (Mercator), units=Meters}
wavelength units = Unknown
band names = {
Resize (Band 1:MOD09GHK_A2005125_B4.tif)}
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Travel Destination: Holland
Next Topic: t-testing in IDL

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

Current Time: Wed Oct 08 19:59:20 PDT 2025

Total time taken to generate the page: 0.00490 seconds