波段比值法统计冰川

pro Part_8
  compile_opt idl2
  e=envi(/headless)
  folder = "D:\Lansat_Data_Composition"
  envi_arr = file_search(folder,'*.dat',count=num)
  for i=0,num-1,1 do begin
    envi_file = envi_arr[i]
    envi_name_length = strlen(envi_file)
    if (envi_name_length eq 32) then begin
      roi_number = strmid(envi_file,27,1)
      roi_number_int = fix(strmid(envi_file,27,1))
    endif else if (envi_name_length eq 33) then begin
      roi_number = strmid(envi_file,27,2)
      roi_number_int = fix(strmid(envi_file,27,2))
    endif
    envi_file_name = "D:\Lansat_Data_Composition\"+roi_number+".dat"
    out_dir = "D:\Lansat_Data_Ratio\"
    out_name = out_dir+roi_number+".dat"
    if (envi_file_name ne "D:\Lansat_Data_Composition\") then begin
      ENVI_OPEN_FILE,envi_file_name,r_fid=fid
      IF (fid[0] EQ -1) THEN CONTINUE
      ENVI_FILE_QUERY,fid,dims=dims
      t_fid=[fid,fid]
      pos=[1,2,3]
      exp='(((float(b2)/b3) ge 1.71)*1)+(((float(b2)/b3) lt 1.71)*0)'
      ENVI_DOIT,'math_doit',fid=t_fid,pos=pos,dims=dims,exp=exp,out_name=out_name,r_fid=r_fid
      print,"Ratio_OK",roi_number
    endif
  endfor
end