1.程序功能:
打开 ASCII 码格式的光谱文件,读取反射率数据并进行显示。
2.源代码:
2.1 方法1:单一窗口显示光谱数据
#在同一个窗口中用不同线形和颜色显示光谱数据#
pro IDL020101
fn=dialog_pickfile();打开光谱数据
nb=file_lines(fn)
data=fltarr(5,nb)
openr,lun,fn,/get_lun
readf,lun,data
free_lun,lun
wv=transpose(data[0,*]);读取波长
c1=transpose(data[1,*]);读取第1条光谱曲线的反射率
c2=transpose(data[2,*]);读取第2条光谱曲线的反射率
c3=transpose(data[3,*]);读取第3条光谱曲线的反射率
c4=transpose(data[4,*]);读取第4条光谱曲线的反射率
;在同一个窗口中用不同线形和颜色显示光谱数据
p1=plot(wv,c1,xtitle='Wavelength(nm)',ytitle='Reflectance',$
/buffer,xrange=[min(wv),max(wv)],color='red',linestyle=0,$
thick=