来源：http://blog.csdn.net/kupe87826/article/details/8130293
参考：http://blog.csdn.net/kupe87826/article/details/8125720
http://www.gisall.com/html/78/6678-5654.html
http://blog.csdn.net/kupe87826/article/details/8125224
软件环境：操作系统win7，ENVI5.0
ETM Flaash大气校正
在ETM蚀变信息提取2文章里，已经得到了辐亮度图像，现在使用ENVI Flaash大气校正功能，进行地表反射率的反演。
1、
数据格式转换
Flaash大气纠正功能需要的辐亮度文件格式为BIL或者BIP，但是目前的文件格式是BSQ的，所以需要进行数据格式转换。
ＥＮＶＩ功能菜单：Basic Tools－》Convert Data(BSQ,BIL,BIP)
选择辐亮度文件
转换参数设置：
再选好输出文件，按OK就行了。
2、
Flaash大气校正
ENVI功能菜单：Spectral->Flaash
下面我们开始进行参数设置。
（1）
输入输出文件设置
Input Radiance Image输入转换格式后的辐亮度文件
Ok后，弹出如下对话框，选择use single scale factor for all bands选项，设置single scale factor为10，具体原因是由于单位不一致。
再选择输出反射率文件的位置，最后设置一下FLAASH输出的一些临时文件的目录。
（2）Flight date栏设置
在*_MTLold.txt找到
ACQUISITION_DATE = 2001-10-23
SCENE_CENTER_SCAN_TIME = 03:50:11.2501468Z
填入即可：
（2）
sensor type栏设置
选择sensor type，MutiSpectral-》landsat TM7
设置好后，sensor Altitude和 pixel size这两项自动就填上了。
还剩下一个Groud Elevation，怎么办呢？google earth啊，哈哈。
Sensor type这栏就填好了，如下图：
（3）
scene center Location栏设置
在earth explore中查看该数据的元文件（具体步骤查看ETM蚀变信息提取1文章）
填上就行了。
（4）
Atmospheric model栏设置
|
Model Atmosphere |
Water Vapor |
Water Vapor (g/cm2) |
Surface Air Temperature |
|
Sub-Arctic Winter (SAW) |
518 |
0.42 |
-16° C (3° F) |
|
Mid-Latitude Winter (MLW) |
1060 |
0.85 |
-1° C (30° F) |
|
U.S. Standard (US) |
1762 |
1.42 |
15° C (59° F) |
|
Sub-Arctic Summer (SAS) |
2589 |
2.08 |
14° C (57° F) |
|
Mid-Latitude Summer (MLS) |
3636 |
2.92 |
21° C (70° F) |
|
Tropical (T) |
5119 |
4.11 |
27° C (80° F) |
|
Latitude (°N) |
Jan |
March |
May |
July |
Sept |
Nov |
|
80 |
SAW |
SAW |
SAW |
MLW |
MLW |
SAW |
|
70 |
SAW |
SAW |
MLW |
MLW |
MLW |
SAW |
|
60 |
MLW |
MLW |
MLW |
SAS |
SAS |
MLW |
|
50 |
MLW |
MLW |
SAS |
SAS |
SAS |
SAS |
|
40 |
SAS |
SAS |
SAS |
MLS |
MLS |
SAS |
|
30 |
MLS |
MLS |
MLS |
T |
T |
MLS |
|
20 |
T |
T |
T |
T |
T |
T |
|
10 |
T |
T |
T |
T |
T |
T |
|
0 |
T |
T |
T |
T |
T |
T |
|
-10 |
T |
T |
T |
T |
T |
T |
|
-20 |
T |
T |
T |
MLS |
MLS |
T |
|
-30 |
MLS |
MLS |
MLS |
MLS |
MLS |
MLS |
|
-40 |
SAS |
SAS |
SAS |
SAS |
SAS |
SAS |
|
-50 |
SAS |
SAS |
SAS |
MLW |
MLW |
SAS |
|
-60 |
MLW |
MLW |
MLW |
MLW |
MLW |
MLW |
|
-70 |
MLW |
MLW |
MLW |
MLW |
MLW |
MLW |
|
-80 |
MLW |
MLW |
MLW |
MLW |
MLW |
MLW |
处理的数据是10月份的，并且中心纬度是38度多的，所以选择SAS模型（Sub-Arctic
Summer）。具体细节请参考我转载的文章Flaash大气纠正
ENVI 。其它的默认。
（5） Aerosol model栏设置
因为数据在中国西部，所以选择rural model
Aerosol retrieaval选择 2-Band（k-T）
Initial visibility设置40km
一般根据天气条件如晴朗、中等雾、厚雾等进行设置，如下图:
具体细节请参考我转载的文章Flaash大气纠正
ENVI，Select an Initial Visibility Value这部分内容，我还没去研究。
（6） MultiSpectral setting
具体为什么还没细的去研究，有些惭愧。
终于设置完了，apply计算吧，哈哈！
对完成后的反射率图像进行波段742合成，如下图：
辐亮度影像的742合成图（原图）：
