FLAASH大气校正软件使用说明

上机实习容：Flaash大气校正学生王玲学号201420771院系城市与环境学院专业地图学与地理信息系统年级2014级教务处制Flaash大气校正实验报告一、实验目的通过本次实验能够更深一步理解大气校正的原理、方法。

并且熟练掌握Landsat8 OLI 数据的大气校正的流程。

二、实验容1、辐射定标目的：将传感器记录的电压或数字量化值（DN值）转换为绝对辐射亮度值（辐射率）。

原理：L=Gain*DN + Bias步骤：（1）首先，在Envi5.1中打开辐射定标工具，Toolbox/Radiometric Correction/ Radiometric Calibration，并在File Selection对话框中选择数据，如下所示：（2）辐射定标参数设置当选择好辐射定标的数据时，接下来需选择定标参数。

其中，①Calibration Type：辐射定标类型，因Flaash校正要求输入的数据为辐亮度值，因此辐射定标类型选择辐亮度。

当数据的每个波段包含Gain和Offest参数时，Envi会自动从元数据文件中获取这些参数，并按照辐射定标公式进行定标，本实验所使用的Landsat8 OLI 数据的元数据中包含这两个参数。

另外，Envi默认Gain和Offest参数定标单位为W/（m2*sr*μm），因此，计算得到的辐亮度值为W/（m2*sr*μm）。

②Output Interleave：输出数据存储顺序，因Flaash校正要求输入的数据存储类型为BIL或BIP，但因BIL的处理速度快，故在此选择BIL。

③Output Data Type：输出数据类型，辐射定标中可以选择的输出数据类型为三种，分别是：浮点型（Float）、双精度浮点型(Double)和无符号位16整型(Uint)。

本实验中使用的OLI6 原始数据为无符号16位整型，在进行Flaash校正时计算缩放因子是无单位型与浮点型数据之间的缩放关系，因此，该处选择浮点型（Float）。

FLAASH大气校正和黑暗像元法操作指导FLAASH大气校正大气是介于卫星传感器与地球表层之间的一层由多种气体及气溶胶组成的介质层。

在太阳辐射到达地表再到达卫星传感器的过程中，两次经过大气，故大气对太阳辐射的作用影响比较大。

大气校正的目的是消除大气和光照等因素对地物反射的影响，广义上讲是获得地物反射率、辐射率或者地表温度等真实物理模型参数，狭义上是获取地物真实反射率数据。

大气校正可以用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等物质对地物反射的影响，也可以消除大气分子和气溶胶散射的影响。

大多数情况下，大气校正也是反演地物真实反射率的过程。

目前应用广泛的大气辐射传输模型有30多种，常用的辐射传输模型主要有6S、MODTRAN和ATCOR等，各种模型的基本原理都是基本相同，其中MODTRAN模型的精度最高。

MODTRAN模型还可以计算热红外波段。

实验数据实习所用到的数据为TM影像和ETM+影像，其FLAASH大气校正的方法和操作步骤一致，这里以TM影像为例介绍。

对于各个影像所需要的具体参数，将在需要用到的时候说明。

本文采用的实验测试数据为，具体的数据内容列表见下图1：1991年1999年2010年图1该数据包含了7个波段，其中的B6为热红外波段，不在本次Flaash校正范围内，其他剩余波段为可见光波段，需要进行大气校正。

操作步骤1、打开tm原始影像数据ENVI > file > open image file > ‘LT51230321991168BJC00_MTLold.txt’ > 打开文件，如下：可见光波段为选择波段图3辐射定标参数设置对话框3、储存顺序调整Flassh大气校正对于波段存储的要求为：BIL，BIP格式，上述计算得到的存储方式为BSQ，在此进行波段存储顺序的转化，具体操作如下：ENVI > basic tools > convert data (BSQ ,BIL ,BIP)图 4 存放顺序转换4、Flaash校正参数设置大气校正的前期准备工作完毕，现在进行校正参数的设置：ENVI > basic tools > preprocessing > calibration utilities > FLAASH，弹出对话框：的BIL或BIP格式数据，然后会弹出下面对话框（图6），按照下图进行设置。

基于像素级的校正，校正由于漫反射引起的连带效应。


1、原理
基于太阳波谱范围内(不包括热辐射)和平面朗伯体(或近似平 面朗伯体)，在传感器处接收的像元光谱辐射亮度公式为：
其中，பைடு நூலகம்

在大气校正参数获取后, 根据该方程逐个波段像元计算地表平均反射率。
2、优点

支持传感器的种类多； 算法精度高；

通过图像像素光谱上的特征来估计大气的属性，
多光谱设置（Multispectral Settings）
校正前
校正后

采用典型地物的光谱曲线特征来检验，如植被、水体。

气溶胶反演（Aerosol Retrieval）

None:选择此项时，初始能见度（Initial Visibility）值将用于气溶胶 反演模型；

2-Band（K-T）:当没有找到合适的黑暗像元时，初始能见度值将用
于气溶胶反演模型；

2-Band Over Water：用于海面上的图像。

尺度转换因子：为了将输入的辐亮度数据变为浮点型，当各波段的辐亮
度单位不一样时，选择第一项；否则选择第二项。

Ground Elevation：影像区域平均海拔，从相应区域的DEM 获取平均值；

大气模型(Atmospheric Model)

气溶胶模型（Aerosol Model）

初始能见度（Initial Visibility Value） 根据天气条件估计。
组员：王玲 王芹 李文苹

模型介绍（原理、优点、条件） 校正过程（步骤、参数说明） 结果评价
1、原理

全称：Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes

第13章高光谱数据FLAASH大气校正本节以AVIRIS高光谱数据为数据源，介绍高光谱数据的FLAASH大气校正过程。

13.1 浏览高光谱数据此AVIRIS高光谱数据为经过传感器定标的辐射亮度数据。

（1）在ENVI主菜单中，选择File→Open Image File，打开JasperRidge98av.img文件。

（2）在波段列表中，选择JasperRidge98av.img，单击右键选择Load True Color，在Display 窗口中显示真彩色合成图像。

（3）在主图像窗口中单击右键，快捷菜单中选择Pixel Locator。

设置Sample:366，Line:179。

此像元为硬质水泥地，吸收特征主要受大气的影响，单击Apply按钮。

（4）在主图像窗口中单击右键，快捷菜单中选择Z Profile，打开Spectral Profile窗口，绘制像素（366,179）的波谱剖面。

（5）在Spectral Profile窗口中，可以看到在760nm，940nm和1135nm处，水汽具有吸收特征，1400nm和1900nm附近基本没有反射能量，二氧化碳在2000nm附近有两个吸收特征。

13.2 AVIRIS数据大气校正（1）在ENVI主菜单中，选择Spectral→FLAASH，打开FLAASH Atmospheric Correction Model Input Parameters对话框。

（2）单击Input Radiance Image按钮，选择JasperRidge98av.img文件。

在Radiance Scale Factors对话框中，选择Read array of scale factors from ASCII file，单击OK按钮。

（3）选择AVIRIS_1998_scale.txt文件，按照默认设置，单击OK按钮。

（4）单击Output Reflectance File按钮，选择输出路径及文件名JasperRidge98av.img。

flaash大气校正Flaash大气校正（IRSP6-08.3.24）大气校正的目的是消除大气和光照等因素对地物反射的影响，获得地物反射率和辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响，消除大气分子和气溶胶散射的影响。

FLAASH可以处理任何高光谱数据、卫星数据和航空数据(860nm/1135nm)，这些数据是由HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)AISA、HARP、DAIS、Probe-1、TRWIS-3、SINDRI、MIVIS、OrbView-4、NEMO 等传感器获得的。

FLAASH还可以校正垂直成像数据和侧视成像数据。

Flaash大气校正使用了MODTRAN 4+辐射传输模型的代码，基于像素级的校正，校正由于漫反射引起的连带效应，包含卷云和不透明云层的分类图，可调整由于人为抑止而导致的波谱平滑。

FLAASH可对Landsat, SPOT, AVHRR, ASTER, MODIS, MERIS, AATSR, IRS等多光谱、高光谱数据、航空影像及自定义格式的高光谱影像进行快速大气校正分析。

能有效消除大气和光照等因素对地物反射的影响，获得地物较为准确的反射率和辐射率、地表温度等真实物理模型参数。

校正过程点击envi——Basic Tools －Preprocessing － Calibration Utilities －FLAASHSpectral －FLAASH.或者点击envi-spectral- FLAASH1、输入数据必须是辐射校正后的数据，对辐射校正数据转成BIL或BIP格式（Basic Tools ——Convert Data）；2、对输入数据进行头文件编辑，主要是对波长wavelenth（即每一波段的波长中心值）和波长宽度fwhm（每一波段的波长范围）的编辑。

不是高光谱数据可以不对fwhm进行编辑。

利用 spectral builder 导入该文件，如下：
FWHM 没有，不用写。 确定后，再次 impor txt 文件，如下：
然后，File 菜单保存为光谱库文件，如下：
此时，可查看光谱响应函数，如下：
2.点击 Mulitispectral Settings 导入刚才的光谱响应函数库，如下：
Flaash 大气校正流程（以 HJ 卫星为例）--于之锋 一．主界面 打开 flaash 主界面， .Input radiance image 必须是辐亮度， 不能是反射率， 需要是 envi standard 格式，并且是 BIL 或者是 BIP，不能是 BSQ 格式。之后还要求输入对应传感器各个波段的 中心波长和半峰全宽 fwhm（Full Wave at Half Maximum”） 。第 j 个波段的中心波长设为 j ， 响应函数设为 f，设光谱响应函数有 n 个间隔，则中心波长计算公式如下：
j
* f
i 1 i
n
i
f
i 1
n
n
i
Fwhm 计算公式为：
Fwhm j
输入的界面为：
f
i 1
i
2
ቤተ መጻሕፍቲ ባይዱ
。
之后选择单一的 scale，scale 设为 10，如下：
1). 2). 3). 4). 5). 6). 7). 8).
Output reflectance file 是输入需要保存反反射率文件的路径。 Out directory for flash file 输入的是保存中间结果的文件夹。 Flight time 是影像成像时间。 Scence center location 可以影像头文件中找到。 Atmospheric Model 可以根据具体地区的温度以及空气湿度找到对应的类型。 Water retrieval 需要有适合水汽反演的波段才行。 Aerosol retrieval 需要？？。 Initial visibility (km) 可以有所测的 AOD 通过软件 modo 进行换算。 大气校正的结果的好 坏可通过调整这个参数，我们所测的 AOD 可以乘以 0.85~0.95 的系数。软件在 E:\大气 校正\flaash 大气校正\modo_install\modo_v3\bin 里面。运行如下：

本文汇总了ENVI FLAASH大气校正模块中常见的错误，并给出解决方法，分为两部分：运行错误和结果错误。

前面是错误提示及说明，后面是错误解释及解决方法。

FLAASH对输入数据类型有以下几个要求：1、波段范围：卫星图像：400－2500nm，航空图像：860nm-1135nm。

如果要执行水汽反演，光谱分辨率<=15nm，且至少包含以下波段范围中的一个：∙∙●1050-1210 nm∙∙●770-870 nm∙∙●870-1020 nm2、像元值类型：经过定标后的辐射亮度（辐射率）数据，单位是：（μW）/（cm2*nm*sr）。

3、数据类型：浮点型（Floating Point）、32位无符号整型（Long Integer）、16位无符号和有符号整型（Integer、Unsigned Int)，但是最终会在导入数据时通过Scale Factor转成浮点型的辐射亮度（μW）/（cm2*nm*sr）。

4、文件类型：ENVI标准栅格格式文件，BIP或者BIL储存结构。

5、中心波长：数据头文件中（或者单独的一个文本文件）包含中心波长（wavelenth）值，如果是高光谱还必须有波段宽度（FWHM），这两个参数都可以通过编辑头文件信息输入（Edit Header）。

一．高级设置里的选项：1．Aerosol Scale Height大气溶胶高度，用来计算邻近效应的范围，1-2km2．CO2 Mixing Ratio (ppm) 2001年前是370ppm。

2001年以后是390ppm。

3．Use Square Slit Function（是否使用平方函数进行邻近像元亮度的均匀）一般选择no 4．Use Adjacency Correction（进行邻近效应校正）5．Reuse MODTRAN Calculations使用以前的MODTRAN模型计算结果6．Modtran Resolution设置MODTRAN模型的光谱分辨率（推荐值5 cm-1) 分辨率高速度慢精度高，分辨率低，速度快，但是精度差。

1.3.2FLAASH其它参数的设置（1）图像中心点坐标可以从相应的HDF文件中找到，也可以从屏幕上直接读取影像的中心坐标，对反演结果影响不大。

当影像位于西半球时，经度为负值；（2）传感器类型当选择传感器类型时，模块会选择相应的类型的传感器波段响应函数，同时系统一般会自动设置传感器的高度和图像的空间分辨率；（3）海拔高度海拔高度为研究区的平均海拔；（4）数据获取日期和卫星过境时间卫星过境时间为格林尼治时间，可以从相应的HDF文件中找到；（5）大气模型模块提供热带、中纬度夏季、中纬度冬季、极地夏季、极地冬季和美国标准大气模型，研究者根据数据获取时间选择相应的大气模型；（6）水气反演大多数多光谱数据不推荐反演水汽含量；（7）气溶胶模型可供选择的气溶胶模型有无气溶胶、城市气溶胶、乡村气溶胶、海洋气溶和对流层气溶胶模型。

当能见度大于40Km时，气溶胶类型选择对反演没有太多影响，一般情况下利用ASTER 数据不做气胶反演；在高级设置中，①Modtran 分辨率（Modtran resolution）：一般设置成5cm-1；②反射率输出的时尺度系数，默认尺度系数是10000，可以使用默认的尺度系数。

若使用默认的尺度系数，大气校正后得到反射率图像的数值域为：0-10000。

其余参数使用默认值。

大气校正的目的是消除大气和光照等因素对地物反射的影响，获得地物反射率和辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响，消除大气分子和气溶胶散射的影响。

FLAASH 可以处理任何高光谱数据、卫星数据和航空数据(860nm/1135nm)，这些数据是由HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)AISA、HARP、DAIS、Probe-1、TRWIS-3、SINDRI、MIVIS、OrbView-4、NEMO等传感器获得的。

FLAASH还可以校正垂直成像数据和侧视成像数据。

FLAASH使用说明一、为什么要进行大气纠正?（1）太阳辐射通过大气以某种方式入射到物体表面然后再反射回传感器（2）原始影像包含物体表面、大气、以及太阳的信息（3）如果我们想要了解某一物体表面的光谱属性，我们必须将它的反射信息从大气和太阳的信息中分离出来AtmosphericBackscatteringDirect Reflection AdjacencyReflection二、FLAASH简介(1) FLAASH的开发者Spectral Sciences, Inc. (SSI)光谱研究的世界领先者FLAASH的开发者与 AFRL一起进行 MODTRAN模型的研究Air Force Research Labs (AFRL)和SSI 以及Spectral Information Technology Application Center(SITAC) 共同研究开发 FLAASH与 SSI一起进行 MODTRAN模型的研究ITTVIS (RSI)ENVI Integration and FLAASH GUI(2)FLAASH的纠正模型FLAASH 通过光谱特征估计大气的属性采用先进的MODTRAN4+模型光谱反射率方程A, B, S, La 由MODTRAN 的输出结果获得三、如何使用FL Preprocessing / Calibration Utilities FLAASH 或者 FL 有以下六个方面内容 设置 1．多光谱/高光谱传感器数据OT 、A VHRR 、 ASTER 、 MODIS 、 MERIS 、 高光谱HYMAP 、A VIRIS 、HYDICE 、CASI 数据要求带有FWHM 值，这些值可以在头gnd integers ，如果输入的数据是浮点型用户可以选择MODTRAN 模型的光谱分辩率AASH打开FLAASHBasic Tools /Spectral / FLAASH.AASH 模块使用主要1．输入文件准备 2．基本参数设置3．多光谱数据参数4．高光谱数据参数设置 5．高级设置 6．输出文件 7．处理结果 输入文件准备（1）支持多种多光谱： QuickBird 、Ikonos 、 Landsat 、 SP AATSR 、 IRS：HYPERION 、（2）要求数据是经过定标后的辐射率数据（3）数据带有wavelenth 值，如果是高光谱文件里或者单独的ASCII 文件里编写好；（4）数据类型：浮点型,或者2/4byte (un)si 数据,单位：uw/cm 2*nm*sr ，这时缩放因子sacle factor 是1；如果是2或者4byte (un)signd integers 型数据，需要计算缩放因子，缩放因子通过以下换算公式获取，如果缩放因子对每一个波段都不是常量，必须提供一个适当的ASCII 文件。

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-23SCENE_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 TemperatureSub-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 Nov80 SAW SAW SAW M LW M LW S AW70 SAW SAW MLW M LW M LW S AW60 MLW MLW MLW SAS SAS MLW50 MLW MLW SAS SAS SAS SAS40 SAS SAS SAS MLS MLS SAS30 MLS MLS MLS T T MLS20 T T T T T T10 T T T T T T0 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 SASLatitude (°N)Jan March May July Sept Nov-50 SAS SAS SAS MLW M LW SAS-60 MLW MLW MLW M LW M LW M LW-70 MLW MLW MLW M LW M LW M LW-80 MLW MLW MLW M LW M LW M LW处理的数据是10月份的，并且中心纬度是38度多的，所以选择SAS模型（Sub-Arctic Summer）。

利用ENVI对MODIS1B数据进行FLAASH大气纠正1．FLAASH对MODIS数据进行大气纠正的Bug改正：在做大气纠正时需要光谱响应函数，在ENVI4.2和ENVI4.3里有一个BUG：ENVI得到的MODIS的光谱响应函数中波段1和波段2是颠倒的。

可以通过以下方式改正这个BUG：1)在ENVI的波谱库浏览器“Spectral Library Viewer”中打开“modis.sli”文件(/envi42/filt_func directory)。

2)在打开的波谱库波段列表中，首先点击“modis-Band 2”，随即弹出光谱绘图窗口，然后点击“modis-Band 1”加载在绘图窗口中。

(先后顺序不能颠倒)。

3)在绘图窗口中选择“Edit ->Data parameters dialog”，在打开的数据参数编辑对话框中，先点击波段2，将“modis-Band 2”改为“modis-Band 1”，再将“modis-Band 1”改为“modis-Band 2”。

点击“Apply”后再点击“Cancel”关闭对话框。

4)在打开的波谱库波段列表中，继续按顺序点击波段(3－19)，使它们加载在绘图窗口中。

5)把绘图窗口中加载的波谱存储为一个新的波谱库文件：在绘图窗口中选择“File ->Save Plot As->Spectral Library”。

新存储的这个文件就是修改好的MODIS的光谱相应函数。

2．FLAASH输入文件准备1)选择File ->Open External File->EOS->MODIS 打开MODIS1B级500米数据.对于1B级数据，已经经过辐射定标,FLAASH纠正的光谱范围是(400－2500um 可见光与近红外波段数据)，因此只能对MODIS数据的1—19波段进行大气纠正。

在测试过程中，第26波段由于信号强度不够会被标记为“bad band”，因此不能处理，我们可以将不能处理的波谱子集裁掉。

大气校正ENVI处理方法
吴祎越
1、数据读取
1、需要envi控件 2、解压后打开 ， ，更改数值、调整为CCD、 和融合。
输入和搜寻、输出，并进行辐射定标 3、单击apply，之后等待。
2、修改数据格式
1、在Basic Tools内容下选择Convert Data，出现窗体 ，
选择需要的数据类型（一
般设为BIL），然后选择 输出地方。
2、等待，
结果为所用。
3、构建波谱响应函数
1、选择Window-Start New Plot Window， Window窗口中，选择File-InputData-ASCII ENVI Plot ，选择下载好的
波谱响应函数文件.txt打开。
2、选择需要的项目。
3、出现波谱曲线
，编辑数据，
4、保存，选择存储位置
，
，
可以查找到文件
。
Байду номын сангаас
4、FLAASH大气校正
1、选择Spectral-FLAASH,
2、从
中获取信息填入。
3、
此处加入波谱函数文件。
4、单击apply ，
等待
5、结果

ENVI5.2下高分二号数据FLAASH 大气校正(2015-07-09 15:57:40)转载▼ 标签：envi5.2高分二号gf2flaash 大气校正 分类： ENVI高分二号卫星于2014年8月19日成功发射，搭载有两台高分辨率1米全色、4米多光谱相机，将带来优于1米空间分辨率的光学遥感影像。

高分二号较高分一号来说分辨率提高一倍，同时具有高定位精度和快速姿态机动能力等特点。

从2014年8月21日首次开机成像并下传数据，已逐步被各行业用户使用。

高分二号数据的大气校正与高分一号类似，但由于官方暂时未正式公布卫星的波谱响应函数等参数，ENVI 也未能及时对其进行原生支持。

为方便大家使用，这里以一景GF2-PMS2 L1A 级数据为例，介绍在ENVI5.2下可行的FLAASH 大气校正流程。

注：同样适合在其他ENVI 版本中操作。

高分二号卫星轨道和姿态控制参数及有效载荷技术指标见下表1、2： 表1 高分二号卫星轨道和姿态控制参数启动ENVI5.2；依次File > Open或直接单击工具栏上的图标，弹出Open对话框，选择数据文件夹下扩展为.tiff的文件，然后点击Open按钮打开（本例中为…/GF2_PMS2_E115.7_N42.7_20140928_L1A0000362235-MSS2.tiff）。

说明：GF-2卫星绝对定标系数可从如下地址下载，/n16/n1115/n1522/n2103/193529.html图1使用Apply Gain and Offset工具进行辐射定标3. FLAASH大气校正FLAASH大气校正需要影像的中心波长信息，ENVI暂不能自动识别GF2数据的头文件信息，因此首先需要手动添加中心波长信息。

（1）添加中心波长在Toolbox中，依次Raster Management > Edit ENVI Header，弹出Edit Header Input File对话框，在Select Input File选项卡中选择上一步辐射定标后的结果（本例中为GF2_PMS2_E115.7_N42.7_20140928_L1A0000362235-MSS2_Rad.dat）；弹出Header Info对话框，点击Edit Attributes，选择Wavelengths，弹出Edit Wavelength values对话框（图2），依次填入各波段对应中心波长，点击OK。

环境小卫星多光谱数据FLAASH精确大气校正方法环境与灾害监测预报小卫星星座A、B星（简称环境小卫星，简写HJ-1A /1B）于2009年3月30日开始正式交付使用，HJ-1-A星搭载了CCD相机和超光谱成像仪（HSI），HJ-1-B 星搭载了CCD相机和红外相机（IRS）。

HJ-1A/1B卫星是继我国继气象、海洋、国土资源卫星之后一个全新的民用卫星。

卫星投入使用后，对自然灾害、生态破坏、环境污染进行大范围、全天候、全天时的动态监测，对灾害和环境质量进行快速和科学评估，提高灾害和环境信息的观测、采集、传送和处理能力，为紧急救援、灾后救助及恢复重建和环境保护工作提高科学依据。

HJ-1数据应用于自然灾害、生态环境之前，需要进行几何及光谱方面的预处理。

ENVI 在数据读取、图像配准、精确大气校正等方面提供了非常好的工具。

1、数据读取目前，网上免费获取的HJ-1A/1B卫星CCD和HSI影像的分发的格式主要有两种：CCD为Geotiff，每一个波段为一个Geotiff文件，并提供一个元数据说明（.XML）; HSI为HDF5格式，也提供一个元数据说明（.XML）。

CCD相机的Geotiff格式直接可以在ENVI软件下打开，利用ENVI->Basic Tools->Layer Stacking工具将各个波段组合成一个文件输出。

HIS高光谱数据是以HDF5格式提供，安装ENVI的HDF5读取补丁后直接在ENVI中打开。

更为方便的方法是直接使用HJ-1数据读取补丁，下载网址为：/ESRI/thread-83044-1-3.html。

直接双击运行.sav 或拷贝sav 文件到ENVI安装目录的save_add 目录下，1）启动ENVI->File->Open External File->HJ-1->HJ-1A/1B Tools工具（下图）。

直接读取CCD、HIS、IRS数据，输出结果为一个多波段的ENVI标准栅格文件，并带有中心波长等信息，其中CCD数据可以直接输出定标结果（辐射亮度）。

ENVI-Flaash大气校正操作流程
Flash大气校正步骤
1. 打开数据中的MTL文件，可以使数据中的信息全部导入。

2. 辐射定标
此界面选择Radiance，然后点击保存文件“111”
3. 将保存后的文件“111”转换成BIL格式
以下窗口点击BIL，保存文件“222”
4 Flash模块大气校正
Flash校正界面：
其中
选择“222”文件，弹出界面如下选择，参数如下填写：
然后选择校正后的保存文件：
下面默认：
下面如下选择：
以下模块，1、3默认，第2个选择该地区高程：
以下选择成像时间和卫星飞行时间，头文件或者下载数据界面可以查找到
以下第3个默认，第一个选择模型，模型选择参考文件“FLAASH 大气校正纬度.jpg”
以下默认：
以下选择
界面如下：
选择
界面如下：参数如下：
最后点击APLY即可。

F L A A S H大气校正参数设置1.3.2FLAASH其它参数的设置（1）图像中心点坐标可以从相应的HDF文件中找到，也可以从屏幕上直接读取影像的中心坐标，对反演结果影响不大。

当影像位于西半球时，经度为负值；（2）传感器类型当选择传感器类型时，模块会选择相应的类型的传感器波段响应函数，同时系统一般会自动设置传感器的高度和图像的空间分辨率；（3）海拔高度海拔高度为研究区的平均海拔；（4）数据获取日期和卫星过境时间卫星过境时间为格林尼治时间，可以从相应的HDF文件中找到；（5）大气模型模块提供热带、中纬度夏季、中纬度冬季、极地夏季、极地冬季和美国标准大气模型，研究者根据数据获取时间选择相应的大气模型；（6）水气反演大多数多光谱数据不推荐反演水汽含量；（7）气溶胶模型可供选择的气溶胶模型有无气溶胶、城市气溶胶、乡村气溶胶、海洋气溶和对流层气溶胶模型。

当能见度大于40Km时，气溶胶类型选择对反演没有太多影响，一般情况下利用ASTER数据不做气胶反演；在高级设置中，①Modtran 分辨率（Modtran resolution）：一般设置成5cm-1；②反射率输出的时尺度系数，默认尺度系数是10000，可以使用默认的尺度系数。

若使用默认的尺度系数，大气校正后得到反射率图像的数值域为：0-10000。

其余参数使用默认值。

大气校正的目的是消除大气和光照等因素对地物反射的影响，获得地物反射率和辐射率、地表温度等真实物理模型参数，用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响，消除大气分子和气溶胶散射的影响。

FLAASH 可以处理任何高光谱数据、卫星数据和航空数据(860nm/1135nm)，这些数据是由HyMAP、AVIRIS、CASI、 HYDICE、HYPERION(EO-1)AISA、HARP、DAIS、Probe-1、TRWIS-3、SINDRI、MIVIS、 OrbView-4、NEMO等传感器获得的。