ENVI中基于FLAASH模型的大气校正

基于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-大气校正

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

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

1.黑暗像元法描述：黑暗像元法的基本原理是假设待校正遥感图像上存在黑暗像元，地表朗伯面反射和大气性质均一，并忽略大气多次散射辐照作用和邻近像元漫反射作用的前提下，反射率很小（近似0）的黑暗像元由于大气的影响，，使得这些像元的反射率相对增加。

可认为这部分增加的反射率是由于大气影响产生的。

这样，将其他像元减去这些黑暗像元的像元值，就能减少大气（主要是大气散射）对整幅图像的影响，达到大气校正的目的。

步骤：先打开待校正图像，再在主菜单中点击Basic Tools > Preprocessing > General Purpose Utilities > Dark Subtract 选择待校正图像，出现如图1-1，1-2所示对话框，设置参数如图所示。

然后点击OK.图2-2所示为处理之后的图像。

图1-1图1-2 图1-3 2.FLAASH大气校正步骤：先打开待校正的图像。

然后将该文件转换成BIL格式,点击主菜单Basic Tool>convert Data(BSQ,BIL,BIP),选择待校正图像，点击OK。

如图2-1所示。

然后出现如图2-2所示对话框，参数设置如图。

然后是在主菜单点击Basic Tools>Preprocessing>Calibration Utilities>FLAASH.接着就是最重要的一步，对校正对话框的参数设置。

如图2-3所示。

重点参数设置有Lat/Lon （纬度/经度）中输入该待校正图像中心位置的经纬度;获得该数据的方法是：选择主菜单Basic Tools→Preprocessing→Data-Specific Utilities→View HDF Global Attributes,在文件选择对话框中选择要校正的图像文件，打开Global Attributes，将Global Attributes中的内容复制到记事本中，查找scenecenter。

ENVI中基于统计学模型的大气校正方法详解

ENVI中基于统计学模型的大气校正方法详解I i基于统计学模型的反射率反演的方法主要有平场域法( Flat Field , FF)、对数残差法(Log Residuals)、内部平均法(Internal Average Relative Reflectanee，IARR )、经验线性法(Empirical Line)。

集中在Basie Tools->Preprocessing-> Calibration Utilities 菜单下。

1. 平场域定标(Flat Field Calibration)Flat Field定标工具通过选择图像中一块具有高反射率、光谱变化平坦的区域，利用这个区域的平均光谱值来模拟飞行时的大气条件下的太阳光谱。

将每个像元的DN值除以选择区域的平均光谱值得到相对反射率，以此来消除大气的影响。

在使用这个工具前，需要利用ENVI提供的感兴趣区绘制工具(ROI Tool )在被定标图像上选择感兴趣区作为平场域(Flat Field )，感兴趣区可选择沙漠、大块水泥地、沙地等区域。

操作过程如下：(1)在主菜单中，选择Basic Tools->Preprocessing-> Calibration Utilities-〉Flat Field 。

在打开的Calibrati on In put File对话框中，选择输入文件，单击" 0K ”。

(2)在打开的Flat Field Calibration Parameters 面板中(图13.4)，在标有"Select ROI forCalibration ”一栏中，选择感兴趣区(只能选择一个)，作为平场域定标的平均波谱区。

(3)选择输出路径及文件名，单击“0K ”执行定标处理。

图13.4 Flat Field Calibration Parameters 面板I2. 对数残差(Log Residuals对数残差定标工具将数据除以波段几何均值，后再除以像元几何均值，可以消除光照、大气传输、仪器系统误差、地形影响和星体反照率对数据辐射的影响。

高光谱影像FLAASH大气校正模型参数影响与评价

高光谱影像FLAASH大气校正模型参数影响与评价摘要：图像输出的正确性与否依靠对后期的处理有至关重要的作用，而大气校正模块中参数的正确性是对本文的数据的处理起着重要的作用，在本文中气溶胶模式、能见度、大气模式和水汽相关的参数都有着重要的作用。

本文利用FLAASH模块进行大气校正能很好的提高图像的质量，大气校正模型输入参数的评价是一个重要的工作，也有重要的意义。

关键词:FLAASH大气校正高光谱遥感高光谱影像EO-1是NASA新千年计划（NMP）的第一颗对地观测卫星，是为了接替美国陆地卫星7而研制的新型地球观测卫星，它的发射是对卫星本身和新的遥感器技术进行论证。

E0-1（对地观测卫星）是在2000年11月21日发射成功。

对地观测卫星上搭载了3种传感器，即：2数据预处理数据的前期已经经过过了：回波校正、背景去除、坏线修复等一系列的工作，这些工作做完以后本不应该再有坏像元，但坏的像元依然存在,所以在进行大气校正以前还需要对数据进行一定的处理[1],在这里需要处理的步骤为:（1）未定标及水汽影响严重波段的剔除。

（2）绝对辐射值的转换。

（3）坏线的修复。

（4）垂直条纹的去除。

（5）图像的各发射腔近场非线性效应。

3 FLAASH大气校正模块简述FLAASH(Fast line—of-sight atmospheric analysis of spectral hypercubes是Spectral scienee Inc. (SSI)、U. S. Air force research laboratory(AFRL)禾口Spectral infor—matio n tech no logy applicati on center(SITAC)合作开发的大气订正软件。

MODTRAN4是FLAASH大气校正的核心代码，Hyperspectral和multispectral image的真实地表反射率、Radiation rate和地表的真实温度等可以同大气校正来获得。

高光谱数据FLAASH大气校正

第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模块的大气校正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模型的大气校正

精品课件
辐射定标
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模型对 HJIA-CCD1影像进行校正。
精品课件
基本步骤
• 1、辐射定标 • 2、制作波谱曲线 • 3、FLASSH校正
精品课件
Flaash校正
利用envi中的Flaash模块对影像进行校正，相关参数在头文件中查找，平均海拔高程通过envi 自带的dem中统计获得。
单击Multispectral Setting按钮，在Filter Function File 导入光精谱品课响件应曲线
单击Advanced Settings，在高级设置中， Tile Size 默认的是Cash size 的大小，手动改为50-100Mb（根据内存大小设定），单击OK

ENVI FLAASH大气校正常见错误及解决方法(2013年7月15号更新)

ENVI FLAASH大气校正常见错误及解决方法（2013年7月15号更新）/s/blog_764b1e9d0100pvrk.html本文汇总了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.平均海拔高程太大。

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

ENVIFLAASH大气校正常见错误及解决方法_图文

ENVI FLAASH大气校正常见错误及解决方法（2013年7月15号更新）(2011-03-07 16:55:57)转载▼标签：flaash大气校正分类：ENVI本文汇总了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.平均海拔高程太大。

flassh大气校正ENVI处理方法

大气校正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、结果

大气校正常见错误处理方法及校正后检查

本文汇总了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）。

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大气校正流程

本文汇总了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大气校正

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大气校正

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

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

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

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

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

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

高分二号卫星轨道和姿态控制参数及有效载荷技术指标见下表1、2：表1 高分二号卫星轨道和姿态控制参数参数轨道类型轨道高度倾角降交点地方时侧摆能力指标太阳同步回归轨道631km ° 10:30 AM ±35°，机动35°的时间≦180s 表2 高分二号卫星有效载荷技术指标参数全色1m分辨率全色/4m分辨率多光谱相机—μm —μm 光谱范围多光谱—μm —μm —μm 全色空间分辨率多光谱幅宽重访周期覆盖周期1m 4m 45km 5天69天说明：资料来源中国资源卫星应用中心网站，/n16/n1130/n192766/ 1. 数据打开暂不支持GF2数据.xml 打开方式，但GF2数据为标准TIFF格式，故可直接使用ENVI的Open菜单打开，只是打开后软件不能自动识别元数据信息。

启动；依次File > Open 或直接单击工具栏上的图标，弹出Open 对话框，选择数据文件夹下扩展为.tiff的文件，然后点击Open按钮打开。

2. 辐射定标全色图像没有方法进行大气校正，所以一般在定量遥感中不使用全色图像；这里也仅对多光谱数据进行辐射定标。

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大气校正常见错误及解决方法(2013年7月15号更新)

ENVI FLAASH大气校正常见错误及解决方法（2013年7月15号更新）/s/blog_764b1e9d0100pvrk.html本文汇总了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.平均海拔高程太大。