comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » GeoTiff- grid

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: GeoTiff- grid [message #58246 is a reply to message #58167]

Wed, 23 January 2008 07:15

robinson.inj
Messages: 32
Registered: August 2007

Member

On Jan 23, 7:08 am, robinson....@gmail.com wrote:
> 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)}- Hide quoted text -
>
> - Show quoted text -

David, using the GCTP keyword I solved the problem. Thank you very
much for all your assitance.

Report message to a moderator

[Message index]

		GeoTiff- grid By: robinson.inj on Mon, 21 January 2008 13:59
		Re: GeoTiff- grid By: David Fanning on Wed, 23 January 2008 11:00
		Re: GeoTiff- grid By: robinson.inj on Wed, 23 January 2008 23:35
		Re: GeoTiff- grid By: robinson.inj on Wed, 23 January 2008 10:52
		Re: GeoTiff- grid By: David Fanning on Wed, 23 January 2008 07:49
		Re: GeoTiff- grid By: David Fanning on Wed, 23 January 2008 07:29
		Re: GeoTiff- grid By: robinson.inj on Wed, 23 January 2008 07:15
		Re: GeoTiff- grid By: robinson.inj on Wed, 23 January 2008 05:08
		Re: GeoTiff- grid By: David Fanning on Tue, 22 January 2008 15:31
		Re: GeoTiff- grid By: David Fanning on Tue, 22 January 2008 15:17
		Re: GeoTiff- grid By: robinson.inj on Tue, 22 January 2008 15:01
		Re: GeoTiff- grid By: David Fanning on Thu, 24 January 2008 06:48
		Re: GeoTiff- grid By: robinson.inj on Thu, 24 January 2008 06:29
		Re: GeoTiff- grid By: David Fanning on Thu, 24 January 2008 05:40

Previous Topic:	Re: Travel Destination: Holland
Next Topic:	t-testing in IDL

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 07:30:01 PDT 2025

Total time taken to generate the page: 2.15754 seconds