Hi
I'm processing external data files and have a problem related to the image coordinates of the tie point within the newly saved files (ENVI-format and geotiff).
This error occurs during batch-processing ALOS AVNIR-2 files, which are available in raw-data format (image centered tie point and rotated) and can be read with the 'read external files routine':
envi_open_data_file, files[i], /alos, r_fid=fid
envi_file_query, fid, ns=ns, nl=nl, nb=nb
My aim is now to save these files as multi-channel geotiffs and ENVI-files, using 'ENVI_OUTPUT_TO_EXTERNAL_FORMAT' and 'envi_doit, 'cf_doit'' by handing over the FIDs.
The new files (geotiffs and ENVI) have the correct map information except for the 'image coordinates of tie point' & 'map coordinate of tie point', which are in the image centre for the original and [1,1] for the newly generated files. Of course it makes sense that he coordinates vary since the images coordinates are shifted. But they are wrong.
One solution would be to somehow adapt the map-information received via the FIDs, by overriding the incorrect map-information with the old one, or, the second approach finding the mistake within my procedure.
I find it in general interesting to know how to change the map-info of existing files (maybe via FIDs), or at least how to change the coordinates of the tie point to [1.5,1.5] (which should be ENVI standard?).
Thank you for any help!
Miguel
PS: The saving/exporting procedures are:
dims = [-1, 0, ns-1, 0, nl-1]
ENVI_OUTPUT_TO_EXTERNAL_FORMAT,fid=fid,dims = dims ,out_name=dir_out+name,/tiff,pos=pos
for j=0,nb-1 do $
envi_doit, 'cf_doit',fid=fid,dims =dims ,out_name=dir_out+name+'test2'+strtrim(j,2),pos=j
|