comp.lang.idl-pvwave archive: archive » map_proj

Home » Public Forums » archive » map_proj_* help

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

Re: map_proj_* help [message #66856 is a reply to message #66805]

Fri, 12 June 2009 08:47

David Fanning
Messages: 11724
Registered: August 2001

Senior Member

David Fanning writes:

> OK, I guess I am giving up on this as a lost cause. I
> have tried enough variations now that I am absolutely
> convinced it is NOT possible to navigate this image
> with the information provided. The numbers don't make
> sense.

I can get closer if I navigate the GeoTiff image you
can download from the same place, but I still can't
line my grid up with the illustration on the web page:

;*********************************************************** *****
file = 'G:\data\pine_2009084_1545_modis_ch02.tif'
image = Read_Tiff(file, Geo=geotag)
image = BytScl(Reverse(image,2))

; Display map in window.
pos = [0.1, 0.1, 0.9, 0.9]
TVImage, image, POSITION=pos, /KEEP_ASPECT, /NOINTERP, /ERASE

; Set up Map Projection. Polar Stereographic, WGS-84 Datum.
map = Map_Proj_Init(106, DATUM=8, CENTER_LAT=-90, CENTER_LON=-70)

s = Size(image, /Dimensions)

xscale = geotag.ModelPixelScaleTag[0]
yscale = geotag.ModelPixelScaleTag[1]
tp = geotag.ModelTiePointTag[3:4]

xOrigin = tp[0]
yOrigin = tp[1] - (yscale * s[1])
xEnd = xOrigin + (xscale * s[0])
yEnd = tp[1]

Plot, [xorigin, xend], [yorigin, yend], $
XStyle=5, YStyle=5, /NoData, POSITION=pos, /NOERASE

; Draw map grids.
Map_Grid, MAP_STRUCTURE=map, LONDEL=5, LATDEL=1, $
LATLAB=-105, LONLAB=-72, /LABEL
END
;*********************************************************** *****

I do notice in the output from listgeo that there is the following
tag in the file:

ProjStraightVertPoleLongGeoKey (Double 1): -100

I wonder if this has something to do with the problem. If so,
there is nothing I can do about it in IDL. :-(

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.")

Report message to a moderator

[Message index]

		map_proj_* help By: mankoff on Sun, 07 June 2009 09:55
		Re: map_proj_* help By: R.G. Stockwell on Fri, 12 June 2009 13:38
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 13:09
		Re: map_proj_* help By: mankoff on Fri, 12 June 2009 13:03
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 12:49
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 12:40
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 12:32
		Re: map_proj_* help By: David Fanning on Sat, 13 June 2009 14:09
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 11:51
		Re: map_proj_* help By: Bruce Bowler on Fri, 12 June 2009 11:48
		Re: map_proj_* help By: Bruce Bowler on Fri, 12 June 2009 11:43
		Re: map_proj_* help By: mankoff on Fri, 12 June 2009 11:28
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 11:16
		Re: map_proj_* help By: David Fanning on Fri, 12 June 2009 08:47
		Re: map_proj_* help By: David Fanning on Mon, 15 June 2009 14:46

Previous Topic:	idl encoder/decoder of json
Next Topic:	Re: file_search: need work around for filenames with square brackets [ ]

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Thu Oct 09 21:52:22 PDT 2025

Total time taken to generate the page: 0.80046 seconds