1、inputs.pro

calcu,2008144,2008233,'G:tangronglinTangrlmod03'

2、calcu.pro

pro calcu,startday,endday,path;,prefix
a=systime(1)
daystart=long(startday)

dayend=long(endday)

while(daystart le dayend) do begin
arg1=strtrim(string(daystart),2)

suffix='SolarZenith'
;-----------------------------------------------------------
file=file_search(path+'MOD03.A'+arg1+'.'+suffix+'.tif',COUNT=count1);to be modified
;-----------------------------------------------------------
if count1 ne 0 then begin
fill_value=255
example_resize_doit,arg1,path,suffix,fill_value
endif

suffix='View_time'
file=file_search(path+arg1+'.*_'+suffix+'.tif',COUNT=count1);to be modified
;-----------------------------------------------------------
if count1 ne 0 then begin
fill_value=255
example_resize_doit,arg1,path,suffix,fill_value
endif

daystart=daystart+1

endwhile
b=systime(1)
print,b-a
end
stop

3、example_resize_doit.pro

pro example_resize_doit,day,path,suffix,fill_value
  ;
  ; First restore all the base save files.
  ;
  envi, /restore_base_save_files
  ;
  ; Initialize ENVI and send all errors
  ; and warnings to the file batch.txt
  ;
  envi_batch_init, log_file='batch.txt'
  ;
  ; Open the input file
file=file_search(path+'MOD03.A'+day+'.'+suffix+'.tif')

envi_open_file, file, r_fid=fid
  if (fid eq -1) then begin
    envi_batch_exit
    return
  endif
  ;
  ; Set the DIMS and POS to keywords
  ; to processes all spatial and all
  ; spectral data. Output the result
  ; to disk.
  ;
  envi_file_query, fid, ns=ns, nl=nl, nb=nb
  dims = [-1, 0, ns-1, 0, nl-1]
  pos  = lindgen(nb)
  out_name ='MOD03.A'+day+'.'+suffix+'.tif'
  ;
  ; Perform the resize calculation.
  ; Make the output image twice as
  ; large in both X and Y. Use
  ; bilinear interpolation.
  ;
  envi_doit, 'resize_doit', $
    fid=fid, pos=pos, dims=dims, $
    interp=0, rfact=[0.92888890953,0.93055543924], $
    out_name=out_name, r_fid=r_fid; spatial resolution resampled from 1000m to 926m.

 envi_file_mng, id=fid, /remove

close,/all
end