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