ETM蚀变信息提取3-- Flaash大气纠正
来源: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的,所以需要进行数据格式转换。
ENVI功能菜单: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
(std atm-cm)
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合成图(原图):

- 2楼liminlu0314昨天 20:42
- 下次应该讲原理,哈哈
- 1楼xieguanglei前天 22:18
- 沙发