IDL读取并显示tm影像

学IDL有几天了，今天就调用envi的接口，给他来个tm影像显示。或许明天再加个定标加快速大气校正。

pro open_rs_photo
COMPILE_OPT IDL2
e=envi(/headless);开启envi批处理模式
fn=envi_pickfile(title='select a rs file');选取文件
envi_open_file,fn,r_fid=fid;打开文件
envi_file_query,fid,$
fname=fname,$
dims=dims,$
ns=ns,$
nl=nl,$
nb=nb,$
data_type=dt,$
sensor_type=st;获取图像信息
tm_data=make_Array(ns,nl,nb,type=dt);创建数组储存图像
for i=0,nb-1 do begin
tm_Data[*,*,i]=envi_get_data(fid=fid,dims=dims,pos=i)
endfor;利用循环把图像的三个波段读到数组里面去
dm=dialog_message('Now I will show you the photo!',title='Cyrus tell you');舍友的提醒
window,0,xsize=ns,ysize=nl;创建窗口
tvscl,tm_Data,true=3;将图像拉伸后显示
end

多处参考了夏水_OA的代码，谢谢他。
效果：