NDVI idl

NDVI的idl算法

pro NDVI
  envi, /restore_base_save_files
  envi_batch_init, log_file='batch.txt'
  envi_select,fid=fid, dims=dims,pos=pos
  num_cols = dims[2]-dims[1]+1
  num_rows = dims[4]-dims[3]+1
  num_bands = n_elements(pos)
  image = fltarr(num_cols,num_rows,num_bands)
  for i=0,num_bands-1 do image(*,*,i) = envi_get_data(fid=fid,dims=dims,pos=pos[i])
  ndvi=fltarr(num_cols,num_rows)
  for i=0,num_cols-1 do begin
    for j=0,num_rows-1 do begin
      ndvi(i,j)=(image(i,j,3)-image(i,j,2))/(image(i,j,3)+image(i,j,2))
    endfor
  endfor
  envi_write_envi_file,ndvi,out_dt=4,out_name='*****'

end