Map to data coordinate conversion [message #83931] |
Thu, 18 April 2013 07:19 |
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?
|
|
|