读MOD15A2产品数据的代码

pro read_modis

compile_opt strictarr ;严格编译器要求

envi, /restore_base_save_files

ENVI_BATCH_INIT, LOG_FILE = 'batch.log'

PRINT, 'Starting batch process.'

print, systime()

root_dir = 'C:\Users\Yang\Documents\National Science Fundation\Images\MODIS\LAI\'

;root_dir='c:\users'

fname = 'MOD15A2.A2010169.h28v05.005.2010178151352.hdf'

a = filepath(fname, root_dir=root_dir)

print, a

envi_open_data_file,a,/hdf_sd,hdfsd_dataset=1,r_fid=fid     ;hdfsd_dataset=0、1、2、3 ...

if (fid eq -1) then begin

  envi_batch_exit

  return

endif

envi_file_qry, fid, ns=ns, nl=nl, nb=nb, data_type = data_type, dims = dims

print, 'ns=',ns, ' nl=',nl, ' nb=',nb, ' data_type=',data_type,' dims=', dims

lai = envi_get_data(fid=fid, dims=dims, pos=0)

window, /free,xsize=(ns), ysize=nl

tvscl, lai, order=1

print, 'FINISHED BATCH PROCESS.'

print, systime()

ENVI_BATCH_EXIT

end