X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 21 Nov 2012 05:19 AM by  anon
batchprocessing external data files; wrong tie points
 1 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
21 Nov 2012 05:19 AM
    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

    Deleted User



    New Member


    Posts:
    New Member


    --
    21 Nov 2012 08:55 AM
    Problem solved. ENVI4.4 ignored the sign of the rotation angle within the map-information. Using ENVI4.6 the problem doesn't occure. Nevertheless it would still interests me, how to influence the map-information having just the FID. Cheers Miguel
    You are not authorized to post a reply.