浮石分享 http://blog.sciencenet.cn/u/floatrock

博文

IDL 地形数据阴影影像自动生成

已有 3535 次阅读 2010-11-12 12:40 |个人分类:未分类|系统分类:科研笔记

该程序将首先提示用户选择一个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



http://blog.sciencenet.cn/blog-499915-383011.html


1 刘洋

发表评论 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-9-26 02:50

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部