ENVI中基于FLAASH模型的大气校正
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）。
ENVI-Flaash大气校正操作流程
Flash大气校正步骤
1. 打开数据中的MTL文件，可以使数据中的信息全部导入。
2. 辐射定标
此界面选择Radiance，然后点击保存文件“111”
3. 将保存后的文件“111”转换成BIL格式
以下窗口点击BIL，保存文件“222”
4 Flash模块大气校正
Flash校正界面：
其中选择“222”文件，弹出界面如下选择，参数如下填写：
然后选择校正后的保存文件：
下面默认：
下面如下选择：
以下模块，1、3默认，第2个选择该地区高程：
以下选择成像时间和卫星飞行时间，头文件或者下载数据界面可以查找到
以下第3个默认，第一个选择模型，模型选择参考文件“FLAASH大气校正纬度.jpg”
以下默认：
以下选择
界面如下：
选择
界面如下：参数如下：
最后点击APLY即可。
大气校正实验报告
一、实验背景遥感技术作为获取地球表面信息的重要手段，广泛应用于地质、农业、环境、城市规划等领域。
然而，由于大气对太阳辐射的吸收、散射和反射作用，遥感图像中的地物反射率信息受到一定程度的影响。
为了消除大气影响，提高遥感图像的精度和应用价值，大气校正技术应运而生。
本文将针对大气校正实验进行详细报告。
二、实验目的1. 理解大气校正的原理和方法；2. 掌握大气校正实验的操作流程；3. 评估大气校正对遥感图像质量的影响。
三、实验原理大气校正的目的是消除大气对遥感图像的影响，恢复地物真实反射率。
主要原理如下：1. 辐射传输模型：根据遥感成像过程中太阳辐射、大气和地物之间的相互作用，建立辐射传输模型，描述太阳辐射、大气和地物之间的能量传递过程。
2. 大气校正算法：通过分析遥感图像和同步观测的大气参数数据，建立大气校正模型，消除大气影响，恢复地物真实反射率。
3. 大气校正方法：主要包括单窗算法、双窗算法、大气校正模型等。
四、实验数据与工具1. 实验数据：选取Landsat 8卫星的OLI传感器获取的遥感图像作为实验数据。
2. 实验工具：ENVI软件，MODTRAN模型，FLAASH大气校正模型。
五、实验步骤1. 辐射定标：将原始遥感图像进行辐射定标，使其具有物理意义。
2. 大气校正：利用FLAASH大气校正模型对辐射定标后的遥感图像进行大气校正。
3. 结果分析：对比校正前后的遥感图像，分析大气校正对图像质量的影响。
六、实验结果与分析1. 辐射定标对原始遥感图像进行辐射定标，得到具有物理意义的图像数据。
2. 大气校正利用FLAASH大气校正模型对辐射定标后的遥感图像进行大气校正，得到校正后的遥感图像。
3. 结果分析（1）目视效果对比通过目视对比校正前后的遥感图像，可以看出大气校正后的图像清晰度更高，地物信息更丰富。
（2）定量分析通过统计分析校正前后遥感图像的地物反射率，可以发现大气校正后的遥感图像地物反射率更加接近真实值。
环境小卫星高光谱数据FLAASH精确大气校正方法
环境小卫星高多光谱数据FLAASH精确大气校正方法高光谱遥感数据的特点是光谱分辨率高、波段连续性强，能够获得地物在一定范围内连续的、精细的光谱曲线，具有非常好的应用前景。
第一个星载民用成像光谱仪是高光谱成像仪Hyperion，其平均光谱分辨率为10 nm，空间分辨率达到30米。
2008年9 月6日我国HJ-1顺利升空，其中A星搭载了我国自主研制的空间调制型干涉高光谱成像仪（HSI）。
HSI 对地成像幅宽为50 km, 星下点像元地面分辨率为100 m，115个波段，工作谱段：459~ 956nm。
具有30度侧视能力和星上定标功能。
ENVI扩展工具一、HDF5读取补丁从HIS数据中获取元数据信息。
下载地址：/ESRI/viewthread.php?tid=37118&extra=page%3 D1。
二、HJ-1数据预处理补丁为HIS数据添加中心波长信息和波段宽度（FWHM）.下载地址：/ESRI/thread-75575-1-1.html拷贝sav 文件到 ENVI安装目录的save_add 目录下。
步骤：(1)选择ENVI->File->Open External File->HJ-1->HJ-1A /1B Tools工具。
图1(2)选择HIS，单击Input Files选择HIS文件（.xml）(3)选择输出路径，单击Apply执行。
(4)在ENVI主模块中，选择Basic Tools->Convert Data(BSQ,BIL,BIP)，将刚才生成的文件转成BIP储存顺序的文件。
按照FLAASH工具的要求，已经将HSI数据转成BIP储存顺序、带有中心波长信息、波段宽度信息的ENVI格式文件。
说明：如果不用HJ-1数据预处理补丁，可以在波段列表中手动输入中心波长和波段宽度信息。
FLAASH参数设定HSI数据已经经过了定标，单位是100W/(m2*um*sr )，启动FLAASH工具，在输入辐射率数据时候，缩放系数填写：1000（即缩小1000倍），（符合FLAASH对辐射亮度单位的要求（μW）/（cm2*nm*sr））。
flaash大气校正
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进行编辑。
FLAASH大气校正流程
本文汇总了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) 分辨率高速度慢精度高，分辨率低，速度快，但是精度差。
FLAASH模块的大气校正
FLAASH模块的大气校正1.1 FLAASH模块简介FLAASH是由世界一流的光学成像研究所－波谱科学研究所（Spectral Sciences）在美国空气动力实验室支持下开发的大气校正模块。
波谱科学研究所在1989年大气辐射传输模型开发初期就广泛从事MODTRAN的研究工作，已成为大气辐射传输模型开发过程中不可缺少的一员。
FLAASH适用于高光谱遥感数据（如HyMap，AVIRIS，HYIDCE，HYPERION，Probe-1，CASI 和AISA）和多光谱遥感数据（如陆地资源卫星，SPOT，IRS和ASTER）的大气校正。
当遥感数据中包含合适的波段时，用FLAASH还可以反演水气、气溶胶等参数。
ENVI中大气校正模型FLAASH，是高光谱辐射能量影像反射率反演的首选大气校正模型。
FLAASH能够精确补偿大气影响，其适用的波长范围包括可见光至近红外及短波红外，最大波长范围为3μm。
其他的大气校正模型是计算方法基于查找表（Look-up Table）、利用插值方法计算，而FLAASH是直接移植了modtran4中的辐射传输计算方法。
用户可以选取代表研究区的大气模型和气溶胶类型，并且对每景影像，Modtran都有独特的解决方案。
1.2 ASTER数据预处理ASTER L1B数据是记录是DN（Digital Number）值，而基于FLAASH大气校正过程中，需要的是辐射能量值。
因此，需要对ASTER L1B数据辐射定标，即把无量纲的DN值转换成有量纲的分辐辐射亮度值的过程（式1），Radiance=gain*DN+offset （式1）其中，gain是增益,offset是偏差。
经辐射定标后，得到天顶辐射能量值，其量纲为W/(m2.sr.um)。
ASTER数据多以HDF格式储存，利用ENVI软件中Baisc Tools->Preprocessing->Data-Specific Utilities->View HDF Global Attribute功能，读取相应ASTER HDF文件中的增益、偏差、成像时间和中心点坐标信息。
大气校正常见错误处理方法及校正后检查
本文汇总了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.Unable to write to this file.File or directory is invalid or unavailable。
没有设置输出反射率文件名。
解决方法是单击Output Reflectance File按钮，选择反射率数据输出目录及文件名，或者直接手动输入。
2.ACC Error：convert7IDL Error：End of input record encountered on file unit:0.平均海拔高程太大。
注意：填写影像所在区域的平均海拔高程的单位是km：Ground Elevation（Km）。
1.如何使用FLAASH对MODIS数据进行大气纠正准备
利用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大气校正
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进行编辑。
flassh大气校正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大气校正
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。
实验九遥感图像的大气校正2
实验九遥感图像的大气校正（二）目的：基于FLAASH模型的大气校正原理：The ENVI FLAASH ModelThis section is a brief overview of the atmospheric correction method used by FLAASH. FLAASH starts from a standard equation for spectral radiance at a sensor pixel, L, that applies to the solar wavelength range (thermal emission is neglected) and flat, Lambertian materials or their equivalents. The equation is as follows:(1) where:r is the pixel surface reflectancere is an average surface reflectance for the pixel and a surrounding regionS is the spherical albedo of the atmosphereLa is the radiance back scattered by the atmosphereA andB are coefficients that depend on atmospheric and geometric conditions but not on the surface.Each of these variables depends on the spectral channel; the wavelength index has been omitted for simplicity. The first term in Equation (1) corresponds to radiance that is reflected from the surface and travels directly into the sensor, while the second term corresponds to radiance from the surface that is scattered by the atmosphere into the sensor. The distinction between r and re accounts for the adjacency effect (spatial mixing of radiance among nearby pixels) caused by atmospheric scattering. To ignore the adjacency effect correction, set re = r. However, this correction can result in significant reflectance errors at short wavelengths, especially under hazy conditions and when strong contrasts occur among the materials in the scene.The values of A, B, S and La are determined from MODTRAN4 calculations that use the viewing and solar angles and the mean surface elevation of the measurement, and they assume a certain model atmosphere, aerosol type, and visible range. The values of A, B, S and La are strongly dependent on the water vapor column amount, which is generally not well known and may vary across the scene. To account for unknown and variable column water vapor, the MODTRAN4 calculations are looped over a series of different column amounts, then selected wavelength channels of the image are analyzed to retrieve an estimated amount for each pixel. Specifically, radiance averages are gathered for two sets of channels: an absorption set centered at a water band (typically 1130 nm) and a reference set of channels taken from just outside the band. A lookup table for retrieving the water vapor from these radiances is constructed.NoteFor images that do not contain bands in the appropriate wavelength positions to support water retrieval (for example, Landsat or SPOT), the column water vapor amount is determined by the user-selected atmospheric model (see Using Water Retrieval for details).After the water retrieval is performed, Equation (1) is solved for the pixel surface reflectances in all of the sensor channels. The solution method involves computing a spatially averaged radianceimage Le, from which the spatially averaged reflectance re is estimated using the approximate equation:(2) Spatial averaging is performed using a point-spread functionthat describes the relative contributions to the pixel radiance from points on the ground at different distances from the direct line of sight. For accurate results, cloud-containing pixels must be removed prior to averaging. The cloudy pixels are found using a combination of brightness, band ratio, and water vapor tests, as described by Matthew et al. (2000).The FLAASH model includes a method for retrieving an estimated aerosol/haze amount from selected dark land pixels in the scene. The method is based on observations by Kaufman et al. (1997) of a nearly fixed ratio between the reflectances for such pixels at 660 nm and 2100 nm. FLAASH retrieves the aerosol amount by iterating Equations (1) and (2) over a series of visible ranges, for example, 17 km to 200 km. For each visible range, it retrieves the scene-average 660 nm and 2100 nm reflectances for the dark pixels, and it interpolates the best estimate of the visible range by matching the ratio to the average ratio of ~0.45 that was observed by Kaufman et al. (1997). Using this visible range estimate, FLAASH performs a second and final MODTRAN4 calculation loop over water.Input Data RequirementsThe input image for FLAASH must be a radiometrically calibrated radiance image in band-interleaved-by-line (BIL) or band-interleaved-by-pixel (BIP) format. The data type may be floating-point, 4-byte signed integers, 2-byte signed integers, or 2-byte unsigned integers. FLAASH requires input data to be floating-point values in units of mW/cm2 * nm* sr. If the input radiance image is not already in floating-point format, you must also know the scale factor (or factors) used to convert radiance data into these units.要求（1）输入数据应该是辐射亮度（2）辐射亮度的单位应该是[μW/(cm2*sr*nm)]（3）数据的存储格式应该是BIL、BIP而不是我们经常使用的BSQ操作步骤：1、将遥感数据中的DN值转换为辐射亮度L（见实验七自动计算辐射亮度）2、统一单位ENVI’s TM/ETM+ calibration utility outputs data with radiance units of [W/(m2*sr*μm)].However, FLAASH requires radiance in units of [μW/(cm2*sr*nm)]. These two units differ by afactor of 10, so an additional step is required to convert the units. This exercise uses Band Math to divide the radiance units by 10.3、将数据转换为BIL或BIP格式（FLASS输入数据的要求）4、打开FLASSH5、对比大气校正前后的数据。
FLAASH大气校正常见错误及解决方法
FLAASH大气校正常见错误及解决方法本文汇总了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.Unable to write to this file.File or directory is invalid or unavailable。
没有设置输出反射率文件名。
解决方法是单击Output Reflectance File按钮，选择反射率数据输出目录及文件名，或者直接手动输入。
2.ACC Error：convert7IDL Error：End of input record encountered on file unit:0.平均海拔高程太大。
注意：填写影像所在区域的平均海拔高程的单位是km：Ground Elevation（Km）。
ENVI 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。
ENVI常见问题
【转】FLAASH大气校正常见错误及解决方法本文汇总了ENVI FLAASH大气校正模块中常见的错误，并给出解决方法，分为两部分：运行错误和结果错误。
前面是错误提示及说明，后面是错误解释及解决方法。
FLAASH对输入数据类型有以下几个要求：1、波段范围：卫星图像：400－2500nm，航空图像：860nm-1135nm。
如果要执行水汽反演，光谱分辨率<=15nm，且至少包含以下波段范围中的一个：l1050-1210 nml770-870 nml870-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.Unable to write to this file.File or directory is invalid or unavailable。
没有设置输出反射率文件名。
解决方法是单击Output Reflectance File按钮，选择反射率数据输出目录及文件名，或者直接手动输入。
2.ACC Error：convert7IDL Error：End of input record encountered on file unit:0.平均海拔高程太大。
注意：填写影像所在区域的平均海拔高程的单位是km：Ground Elevation（Km）。
FLAASH大气校正常见错误及解决方法
本文汇总了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.Unable to write to this file.File or directory is invalid or unavailable。
没有设置输出反射率文件名。
解决方法是单击Output Reflectance File按钮，选择反射率数据输出目录及文件名，或者直接手动输入。
2.ACC Error：convert7IDL Error：End of input record encountered on file unit:0.平均海拔高程太大。
注意：填写影像所在区域的平均海拔高程的单位是km：Ground Elevation（Km）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
利用envi中的Flaash模块对影像进行校正，相关 参数在头文件中查找，平均海拔高程通过envi 自带的dem中统计获得。
Flaash校正
单击Multispectral Setting按钮，在Filter Function File 导入光谱响应曲线
单击Advanced Settings，在高级设置中，Tile Size 默认的是Cash size 的大小，手动改为 50-100Mb（根据内存大小设定），单击OK
比较校正前后的图像,1为校正前，2为校正后
辐射定标
1、 在envi软件中通过 band math进行计算 定标公式：L=DN/a + L0，式中L为辐亮度， a绝对定标系数增益， L0为偏移量，转换后辐 亮度单位为 W⋅m−2⋅sr−1⋅μm−1
辐射定标
2、将定标后的四个波段合并，并转换为 BIL格式，在头文件中输入每个波段中心 波长 b1(489)b2(568nm)b3(664nm)b4(829nm )，中心波长选择Window>Start New Plot Window，ENVI Plot Window窗口中，选择 File->Input Data>ASCII，生成四条曲 线 2、选择File- >Save Plot As->Spectral Library，将波谱曲线保存为波谱库文件。
基于FLAASH模型的 大气校正
大气校正的目的
大气影响降低了图像的反差比，使 图像可读性降低 ，增加了解译的困难。 为了消除大气的影响，需要对图像进 行大气校正。 本实验是基于flassh模型对HJIACCD1影像进行校正。
基本步骤
• 1、辐射定标 • 2、制作波谱曲线 • 3、FLASSH校正