|
I am not sure what you need the lat/lon bands for in the decorrelation stretch. An ASTER 1B product is already registered to a map projection if you are looking to reference the file. If you simply want to run the decorrelation stretch on the TIR bands from an ASTER HDF file, you can use ENVI's ENVI_OPEN_DATA_FILE. The FID contains all 3 FIDs for the 3 datasets in an ASTER file. For example:
file=envi_pickfile()
envi_open_data_file, file, /aster, r_fid=fid
print, fid
envi_file_query, fid[2], nb=nb
print, nb
which returns:
% Compiled module: TEST.
77 73 66
5
This shows I have chosen the FID for the 3rd HDF dataset which has 5 bands - the TIR bands of ASTER. From here you can run decorr_doit on the selected FID.
If you do need to select specific HDF datasets in the file, you can do this using ENVI routines also:
ENVI_OPEN_DATA_FILE, fname, r_fid=fid, /hdf_sd, hdfsd_dataset=0, hdfsd_interleave=0
where hdfsd_dataset is the specific dataset you want to open and are numbered 0 through number of datasets in the file and hdfsd_interleave is the file interleave where 0 = BSQ.
If you want to look at the number of datasets in your ASTER file, you can open the file as generic HDF in ENVI which brings up the HDF file selection dialog.
|