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

Home » Public Forums » archive » Map to data coordinate conversion
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Map to data coordinate conversion [message #83931] Thu, 18 April 2013 07:19
PMan is currently offline  PMan
Messages: 61
Registered: January 2011
Member
Hello all,

I am looking to convert, given a geotff image with map info, the lat/lon data to a pixel location in an image.

What I want to do seems pretty easy, but I am having a hard time coding it in IDL for some reason. I must be missing something.

Using the ENVI() API, this was straight forward:

oEnvi = envi(/headless)
eMask = oEnvi->openRaster(file)
fid = ENVIRasterToFID(eMask)
envi_convert_file_coordinates, fid, xf, yf, lon, lat

with lon, lat being my predefined input and xf, yf being my desired output.

The analogous way to do this in IDL is (I think):

img = image(file, /buffer)
xy = img->convertcoord(lon, lat, /to_data)

however the image() function won't read the geotiff properly (but ENVI() accepts it just fine). Perhaps my geotiff tag is not properly formed or something. Ideas on another way?
[Message index]
 
Read Message
Previous Topic: cgArrow Solid keyword with a cgWindow
Next Topic: issue with CMSAVE

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

Current Time: Wed Oct 08 18:03:45 PDT 2025

Total time taken to generate the page: 0.00278 seconds