ENVI辐射校正（FLAASH）自用

一、辐射定标

1. 数据（多光谱数据）

2. 打开envi 5.3 64bit(非经典版)，File—Open—选择元数据“LT51200382005153BJC00_MTL.txt”,使用波段321（真彩色），设置线性拉伸2%。

或者可以这样打开元数据

3.在 Toolbox 中，选择 Radiometric Correction > Radiometric Calibration，对文件

对话框中选择多光谱数据。打开 Radiometric Calibration 面板，选择6波段文件，点击OK设置如下参数：

定标类型（Calibration Type）：辐射率数据 Radiance

储存顺序（Interleave）：BIL

数据类型（Data Type）：Float

辐射率数据单位调整系数（Scale Factor）：0.1。

输出路径不宜过长且不要出现中文，输出格式为.dat（也有书写.img）

PS：经典版/10(输入10) 普通版*0.1（输入0.1）

格式转换：辐射定标后数据格式为BSP，但大气校正需要BIL格式，所以需转换格式。

单击 Apply FLAASH Settings 按钮，自动设置 FLAASH 大气校正工具需要的数据类型

4.打开定标后的数据，选择 Display>Profiles>Spectral 查看波谱曲线，看到

定标后的数值主要集中在 0-10 范围内，单位是 µW/(cm2 * sr * nm)。

左为定标后水体的波谱曲线，右为原始数据水体波谱曲线

二、影像裁剪(本实验不涉及)

若需要进行区域的裁剪，则在辐射定标之后进行，因为裁剪后就没有头文件了，但辐射定标需要头文件，而大气校正不需要。

1. 加载进.shp数据后，在Toolbox中搜索“Subset Data from ROIs”，选择定标后的数据。

2. 然后进入如下界面，选择ROI区域，也就是裁剪区，因为我的适量数据是两个文件组成的，所以有两个，都选上。“Mask pixel outside of ROI?”这个地方选择yes，选yes之后意思就是按照你的不规则边界区域进行裁剪，如果是no就是根据ROI 的那个矩形框进行裁剪，裁剪出来是个矩形。“Mask Background”背景值设置为零，意思是黑色。

3.点击OK，裁剪成功。

三、大气校正

1. 多光谱数据 FLAASH 大气校正

在 Toolbox 中打开 FLAASH 工具：/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction。启动 FLAASH Atmospheric Correction Module Input Parameters 面板，设置参数。

Input Radiance Image：选择辐射定标结果数据或者裁剪数据，在打开的 Radiance Scale

Factors 面板中，设置 Single scale factor：1。（因为普通版前面已经设置过了，所以是1，若是用经典版做的，前面未设置，则此处需要写10）

Output Reflectance File：设置输出路径和文件名；

Output Directory for FLAASH Files：设置其他文件输出目录；（这个参数可能不需要改动）

2.传感器参数设置：

①中心点经纬度 Scene Center Location：如果图像有地理坐标则自动获取；

②选择传感器类型 Sensor Type：Landsat TM5

③Sensor Altitude和Pixel Size对应的TM5自动获取

④设置影像区域的平均地面高程 Ground Elevation：0.033km（详细步骤见附录1）

⑤Flight Date和Flight Time GMT影像成像时间，从元数据文件获取：

【方法1】在 layer manager 中的数据图层中右键选择 View

Metadata，浏览 time 字段获取成像时间，2005 年 6 月 2 号 02:24:48。

【方法2】也可以从元文件"*_MTL.txt"中找到，具体名称：DATE_ACQUIRED 和

SCENE_CENTER_TIME；

⑥大气模型Atmospheric Model根据成像时间和纬度选择，6月和31纬度，选SAS或MLS都可

⑦气溶胶模型 Aerosol Model：Urban（多为农村）

气溶胶反演方法 Aerosol Retrieval：2-band（K-T）

3.多光谱数据参数设置

单击 Multispectral Settings，打开多光谱设置面板；

①K-T 反演选择默认模式：Defaults->Over-Land Retrieval standard（600:2100），

自动选择对应的波段；

②Filter Function File用于设置波普响应函数，当传感器未知时，需手动设置。

（手动设置从安装路径下找文件F:\ENVI\ENVI53\classic\filt_func\tm.sli）

③其他参数选择默认。

4.点击Apply进行大气校正

查看校正前后波谱曲线。

植被：

水体：

影像对比：

附录1

影像区域的平均地面高程 Ground Elevation（km）

【方法1】：File—open—找到envi安装目录下—ENVI—ENVI53—data—打开“GMTED2010.jp2”—工具Statistics—computer statistics—选择“GMTED2010.jp2”文件—stats subset—file—选中区域范围，如选择定标结果—ok—ok—ok—ok—Mean就是区域平均海拔（m），将该值除1000，就是 Ground Elevation值0.033