FLAASH大气校正参数设置

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

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

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），按照下图进行设置。

FLAASH 大气校正软件使用须知胡顺石 hufrank@（中国科学院遥感应用研究所）1 输入数据要求1．波段范围：卫星图像400—2500nm ，航空图像860—1135nm ；2．数据类型：必须是浮点型、4位有符号整型、2位无符号整型；3．影像存储格式：输入影像必须是BIL 或BIP 格式；4．影像单位：输入影像的单位必须是2/()W cm nm sr μ⋅⋅，如果单位不一致，先转换成所需要的单位。

可以指定为每个波段指定一个缩放因子，这需要一个文本文件，文本文件中包含每个波段对应的缩放因子；也可以为整幅影像的所有波段指定相同的缩放因子，5．水汽获取：如果要获取图像水汽含量，传感器具有1050~1210nm ，770~870nm 或者870~1020nm 范围内的通道，并且这些通道必须具有至少15nm 的光谱分辨率；6．气溶胶获取：如果要获取图像气溶胶含量，传感器必须具有660nm 和2100nm 附近的通道，这些通道主要是用于获取“黑暗像元”，条件为0.662.1 2.1(0.1)&&(0.45)ρρρ≤≈，如果输入图像中还具有800nm 和420nm 附近的通道，可以用于消除阴影和水体，条件为0.880.42 1.0ρρ≤； 7．输入波长信息：对于FLAASH 暂时没有的传感器类型：如果是高光谱数据，需要波长、FWHM 信息，这些信息可以在头文件中，也可以建立一个ASCII 文件进行存储；如果是多光谱数据，由要输入光谱响应函数，这需要在“Multispectral Settings ”中进行设置。

2 多光谱设置1．对于大部分多光谱数据而言，由于其不具备水汽反演通道，并且光谱分辨率没有达到15nm ，水汽反演功能是不能设置的；2．气溶胶反演，如下图所示。

这些参数用于确定黑暗像元，用于气溶胶反演；KT Upper Channel：建议选择2100nm附近的通道；KT Lower Channel：建议选择660nm附近的通道；Maximum Upper Channel Reflectance：建议设置为0.1 ，即：2.10.1ρ≤；Reflectance Ratio：为反射率比值，建议设置为0.45，即：0.662.10.45ρρ≈；Cirrus Channel (optional)：确定云的通道，建议设置为1367—1383nm左右的通道；图1 多光谱设置水汽、气溶胶和云反演通道设置如下图所示：图2 通道设置波长范围3高光谱设置1．Automatic Selection：根据输入高光谱数据的光谱特征信息，自动选择水汽、气溶胶和云反演通道；2．File：选择某一高光谱通道定义文件，然后根据文件中指定的反演通道进行反演。

第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。

本文汇总了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) 分辨率高速度慢精度高，分辨率低，速度快，但是精度差。

本文汇总了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还可以校正垂直成像数据和侧视成像数据。

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”，因此不能处理，我们可以将不能处理的波谱子集裁掉。

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

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

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进行编辑。

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进行编辑。

1.3.2 FLAASH其它参数的设置（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还可以校正垂直成像数据和侧视成像数据。

ENVI FLAASH大气校正今天忙了一下午，终于将Hyperion数据的FLAASH大气校正搞定了，下面是FLAASH User's GUIDE 的一个例子的步骤，共享一下Start FLAASH and Set the Input and Output Parameters1. From the main ENVI menu select Spectral →FLAASH or Basic Tools →Calibration Utilities →FLAASH.2. Click on the Input Radiance Image button and select theJasperRidge98av.img file from ENVI’s tandard Input File dialog and click OK.3. In the Radiance Scale Factors dialog, click the Read array of scale factors(one per band) from ASCII file button, then click OK.4. When the file selection dialog appears, navigate to theenvidata/flaash/hyperspectral/input_files/ directory on Tutoral Data CD #3 and select the A VIRIS_1998_scale.txt file.5. In the Input ASCII File dialog, accept all of the default values and click OK.Recall that the input radiance image has been scaled into two-byte signed integers. In order for FLAASH to compute the atmospheric correction, these data must be converted into floating-point radiance values in units of. The 1998 A VIRIS scale factors (which are valid for all A VIRIS data collected between 1995 and 2003) are 500 for the first 160 bands and 1000 for the remainder.in the FLAASH Atmospheric Correction Model Input Parameters dialog, the default path and output reflectance file name for the FLAASH-correctedreflectance result are displayed in the Output Reflectance File text box.6. Change the output path to one on a writable drive by typing in the text box, or by clicking on the Output Reflectance File button to navigate to the desired output directory and define the output file name.7. In the Output Directory for FLAASH Files text box, enter the full path of the directory where you want to have all other output FLAASH files written. You may also click on the Output Directory for FLAASH Files button to the left of the text box to navigate to the desired directory.8. In the Rootname for FLAASH Files text box, enter the name you want to use as a prefix for the FLAASH Output Files.ENVI will automatically add an underscore character to the rootname that youenter.Restore a Template and Review the Model Parameters1. In the bottom right-hand corner of the FLAASH Atmospheric Model InputParameters dialog, click the Restore button.The FLAASH run parameters for the sample Jasper Ridge A VIRIS image havebeen provided in a template file on the Tutorial Data CD #3.2. Navigate to the envidata/flaash/hyperspectral/input_files/directory on Tutorial Data CD #3 in the file selection dialog, select the JasperRidge98av_template.txt file, and click Open.The FLAASH Atmospheric Correction Model Input Parameters dialog fills in automatically and looks like Figure 4-2 (except for the file input and outputfields, which will be specific to your machine).NoteIf a message appears, warning you that the paths don’t exist, click OK to dismiss it.3. Review the scene collection details and model parameters for the Jasper Ridge scene.4. Click the Advanced Settings button at the bottom of the dialog window.The parameters in the Advanced Settings dialog allow you to adjust additional controls for the FLAASH model. Note that the default setting for Automatically Save Template File is Yes and Output Diagnostic Files is No.While you may find it excessive to save a template file for each FLAASH run,this file is often the only way to determine the model parameters that were used to atmospherically correct an image after the run is complete, and access to it can be quite important. The ability to output diagnostic files is offered solely as an aid for RSI Technical Support engineers to help diagnose problems. For more details about Advanced Settings see “FLAASH Advanced Settings”on page 39.5. Click Cancel in the Advanced Settings dialog to return to the input parameters dialog.Start the Processing and View the Corrected Image1. In the FLAASH Atmospheric Model Input Parameters dialog, click Apply to begin the FLAASH processing. A standard ENVI status report dialog will be displayed. •You may cancel the processing at any point, but be aware that there are some FLAASH processing steps that can’t be interrupted, so the response to the Cancel button may not be immediate.2. When FLAASH completes, the output reflectance image, as well as the column water vapor image and the cloud classification map, will be entered into the Available Bands List. You should also find the journal file and the template file in the FLAASH output directory.3. Using standard ENVI procedures, display the reflectance result into a new imagedisplay.4. In the Main Image window, right-click and select Z Profile from the shortcut menu.5. Move the cursor in the Main Image window and observe the reflectance spectra in the Spectral Profile window.Verify the Model ResultsThe results you produce with the Jasper Ridge files should be identical to the data found in the envidata/flaash/hyperspectral/flaash_results/ directory on Tutorial Data CD #3.Compare Images by Linking1. From the ENVI main menu, select File →Open Image File and open the perRidge98av_flaash_refl.img file from the envidata/flaash/hyperspectral/flaash_results/ directory on Tutorial Data CD #3.The image bands are listed in the Available Bands List.2. In the Available Bands List, highlight one of the image bands, then right-click and select Load True Color to <new> from the shortcut menu (or display an RGB combination of bands 28, 19, 10, spectively).3. Display the same RGB combination from your FLAASH reflectance result into another image window.4. Right-click in one of the image display windows and choose Link Displays from the shortcut menu. In the Link Displays dialog, set the Dynamic Overlay to Off and click OK.5. Double-click in one of the Main Image display windows to display the Cursor Location/Value window.6. Roam around one of the images using your mouse and note the data values int he Cursor Location/Value window. You should see that the data values are identical for corresponding bands in both images.Compute a Difference Image Using Band MathFor a more quantitative verification of the reflectance results, compute a difference image using Band Math.1. From the ENVI main menu, select Basic Tools →Band Math. The Band Math dialog appears.2. In the Enter an Expression text box, enter the following expression: float(b1) –b23. Click OK.4. In the Variables to Band Pairings dialog, click on B1 to highlight it and click the Map Variable to Input File button.5. When the input file selection dialog appears, click on the JasperRidge98av_flaash_refl.img file and click OK.6. Use the same procedure to assign B2 to your FLAASH reflectance image.7. Enter or choose an output file name for the result and click OK.NoteThe file size for this difference image will be twice as large as the FLAASH reflectance image file, so be sure you have sufficient disk space for this Band Math result.8. Every value in the difference image should be zero. To ensure that the results are identical, select Basic Tools →Statistics →Compute Statistics from the ENVI main menu to calculate the basic statistics for the difference image.Note the Max and Min columns in the statistics report window.NoteDue to differences in computer machine precision, your FLAASH reflectanceimage result may differ from those in the verification directory by approximately1-5 DNs, or 0.0001 to 0.0005 reflectance units。

环境小卫星多光谱数据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数据可以直接输出定标结果（辐射亮度）。

应用FLAASH模块对TM数据进行大气校正目前大气校正的方法主要有基于辐射传输方程法和经验线性法，其中6S，MODTRAN、LOWTRAN、ATCOR等软件主要是基于辐射传输方程理论；经验线性法主要是基于野外实测地物波谱，建立经验线性方程完成大气校正。

ATCOR和FLAASH都是基于MODTRAN 发展起来的，FLAASH和ATCOR2主要针对平原区，A TCOR3针对有DEM的山区。

FLAASH适用于高光谱遥感数据和多光谱遥感数据的大气校正。

当遥感数据中包含合适的波段时，用FLAASH还可以反演水气、气溶胶等参数。

FLAASH能够精确补偿大气影响，其适用的波长范围包括可见光至近红外及短波红外，最大波长范围为3μm。

FLAASH 直接移植了modtran4中的辐射传输计算方法。

用户可以选取代表研究区的大气模型和气溶胶类型，并且对每景影像，Modtran都有独特的解决方案。

一、数据预处理1 打开数据：File/Open External File/Landsat/Fast或者GeoTIFF2 编辑头文件：数据层右键快捷菜单Edit/Header可编辑头文件。

添加Wavelength（输入每个波段的中心波长，单位为微米）、FWHM（输入与波长同单位的数据，可近似用每个波段的带宽代替）、Gains和Offsets值（使用USGS发布的数值），MapInfo信息（一般在图像打开时，会进行自动识别。

这一项一定要定义正确，否则FLAASH模块不能正常运行）、Sensor Type。

这些信息编辑好以后，以后生成的数据都具有相同的头文件信息。

表1 Landsat5中心波长和半幅全宽数值Landsat 7 中心波长和半幅全宽数值表2 Landsat5的增益和偏置值3辐射定标：Baisc Tools/Preprocessing/Calibration Utilities/Landsat TM定标公式为Radiance=gain*DN+offset。

上机实习容：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）。