<新手教程>—遥感—行星反射率的计算及合并

发一下我以前写的一些envi的简单教程，供新手学习，大神勿喷，很多内容我也不是特别清楚
老师也只是教我怎么操作，并没有很详细讲每一步的意义。
我的目的是抛砖引玉。希望大家把会的东西发出来让大家相互交流学习。可以是遥感方面的，gis方面的或者其他有关软件的应用操作，或者小技巧。也可以从网上转载一些内容但是要注明出处，作者，有转载的字样。回来我单开一贴和大家讨论一下格式。
I

一些简单的说明：
软件：envi4.8
影像：landsat TM5
1. 1.打开landsat5 TM 数据File→Open External File→L andsat→Geo TIFF→文件夹里的图像
如

2, 7 4 2 波段 r g b(只是打开看一下)
如图

3.
单击 basic tools》 band match打开如图所示窗口

4. 在Enter an expression的框里输入(b1-1)*(193.0+1.52)/254-1.52→Add to list→Ok→选定第一波段→Choose建立一个文件夹保存

然后依旧打开 bandmatch 输入!pi*b1*1^2/(1957*cos(!pi*(90-65.3691418)/180))》Add to list→Ok→选定第一波段亮度值。>choose建立一个文件保存。如图

5 重复3。4步骤把 dn值转化为亮度然后转化为反射率所有公式如下
(b1-1)*(193.0+1.52)/254-1.52 这是第一次用band match输入的公式求亮度值
(b2-1)*(365.0+2.84)/254-2.84
(b3-1)*(264.0+1.17)/254-1.17
(b4-1)*(221.0+1.51)/254-1.51
(b5-1)*(30.2+0.37)/254-0.37
(b7-1)*(16.5+0.15)/254-0.15
!pi*b1*1^2/(1957*cos(!pi*(90-65.3691418)/180)) 这是第二次根据亮度求反射率
!pi*b2*1^2/(1829*cos(!pi*(90-65.3691418)/180))
!pi*b3*1^2/(1557*cos(!pi*(90-65.3691418)/180))
!pi*b4*1^2/(1047*cos(!pi*(90-65.3691418)/180))
!pi*b5*1^2/(219.3*cos(!pi*(90-65.3691418)/180))
!pi*b7*1^2/(74.52*cos(!pi*(90-65.3691418)/180))
然后就按照这个顺序做得到最终六个波段反射率I