行星反射率的计算

行星反射率的计算

1.打开ENVI软件→File→Open External File→Landsat→Geo TIFF→文件夹里

的图像

2.选择742波段为RGB值打开图像

3.将六个波段的DN值转化为亮度值：Basic Tools→Band Math→在Enter an

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

4.重复步骤2，将第1至第7（第6除外）波段的DN值全部转化为亮度值，每

次代入相应波段对应的公式（见最后）

5.将六个波段的亮度值转化为反射率：Basic Tools→Band Math→在在Enter

an expression的框里输入!pi*b1*1^2/(1957*cos(!pi*(90-65.3691418)/180))→Add to list→Ok →选定第一波段的亮度值（带公式的那个）→Choose建立一个文件夹保存6.重复步骤4，将第1至第7（第6除外）波段的亮度值全部转化为反射率，

每次代入相应波段对应的公式（见最后）

7.打包：File→Save as→ENVI standard→选择六个波段的反射率→Reorder

重排一下

8.打开保存好的反射率的文件选择742波段为RGB值打开图像 New Display

9.连接俩图像：在图像窗口Tools→Link→Link Displays→Dynamic选择Off

10.在图像上双击左键查看RGB的值

用ENVI的band math功能，将DN值向亮度值转换部分。

(b1-1)*(193.0+1.52)/254-1.52

(b1-1)*(365.0+2.84)/254-2.84

(b1-1)*(264.0+1.17)/254-1.17

(b1-1)*(221.0+1.51)/254-1.51

(b1-1)*(30.2+0.37)/254-0.37

(b1-1)*(16.5+0.15)/254-0.15

亮度值向反射率转换部分

!pi*b1*1^2/(1957*cos(!pi*(90-65.3691418)/180))

!pi*b1*1^2/(1829*cos(!pi*(90-65.3691418)/180))

!pi*b1*1^2/(1557*cos(!pi*(90-65.3691418)/180))

!pi*b1*1^2/(1047*cos(!pi*(90-65.3691418)/180))

!pi*b1*1^2/(219.3*cos(!pi*(90-65.3691418)/180))

!pi*b1*1^2/(74.52*cos(!pi*(90-65.3691418)/180))