代码如下
pro dqsj;打开envi过后读取数据进行计算
cd,'E:\xxsx\idl\fx'
fn='Data_OLI'
;运行前打开envi
envi_open_file,fn,r_fid=fid
envi_file_query,fid,ns=ns,nl=nl,nb=nb,dims=dims
map_info=envi_get_map_info(fid=fid)
r=envi_get_data(fid=fid,dims=dims,pos=3)
nr=envi_get_data(fid=fid,dims=dims,pos=4)
ndvi=(float(nr)-r)/(float(nr)+r);计算ndvi
fns='ndvi'
;写入,记得前面需要打开envi
envi_write_envi_file,ndvi,out_name=fns,map_info=map_info
;利用smooth滤波
smoothndvi5=smooth(ndvi,3)
envi_write_envi_file,smoothndvi5,out_name='smooth_ndvi5',map_info=map_info
end