||
该程序将首先提示用户选择一个DEM文件,读入数据后,同时显示DEM以及与它相关的阴影影像。TOPO_DOIT关键字将调用ENVI地形分析功能来计算获得阴影影像。
pro view_dem
编译器初始化;
compile_opt idl2
envi,/restore_base_save_files
envi_batch_init
选取DEM数据;
dem_file=dialog_pickfile(title='select a DEM')
if(dem_file eq '')then return
打开DEM文件;
envi_open_file,dem_file,r_fid=dem_fid
获取文件信息:
envi_file_query,dem_fid,ns=ns,nl=nl
proj=envi_get_projection(fid=dem_fid,pixel_size=pixel_size)
dims=[-1L,0,ns-1,0,nl-1]
生成阴影影像;
envi_doit,'TOPO_DOIT',azimuth=15.0,bptr=[2],dims=dims,$
elevation=45.0,fid=dem_fid,in_memory=1,pos=[0],$
r_fid=shaded_fid,pixel_size=pixel_size
dem=envi_get_data(fid=dem_fid,dims=dims,pos=[0])
shaded=envi_get_data(fid=shaded_fid,dims=dims,pos=[0])
同时显示DEM和阴影影像;
window,/free,xsize=(2*ns),ysize=nl
tvscl,dem,order=1
tvscl,shaded,ns,0,order=1
envi_batch_exit
end
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-10-20 01:29
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社