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的，所以需要进行数据格式转换。

ＥＮＶＩ功能菜单：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合成图（原图）：