Flaash大气校正模块常见错误及解决方法
基于TM影像的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）。
Flassh大气校正
[转载]大气校正（转）大气校正是定量遥感中重要的组成部分。
本专题包括以下内容：∙ ●大气校正概述∙∙●ENVI中的大气校正功能1大气校正概述大气校正的目的是消除大气和光照等因素对地物反射的影响，广义上讲获得地物反射率、辐射率或者地表温度等真实物理模型参数；狭义上是获取地物真实反射率数据。
用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等物质对地物反射的影响，消除大气分子和气溶胶散射的影响。
大多数情况下，大气校正同时也是反演地物真实反射率的过程。
很多人会有疑问，什么情况下需要做大气校正，我们购买或者其他途径获取的影像是否做过大气校正。
通俗来讲，如果我们需要定量反演或者获取地球信息、精确识别地物等，需要使用影像上真实反映对太阳光的辐射情况，那么就需要做大气校正。
我们购买的影像，说明文档中会注明是经过辐射校正的，其实这个辐射校正指的是粗的辐射校正，只是做了系统大气校正，就跟系统几何校正的意义是一样的。
常见的绝对大气校正方法有：●基于辐射传输模型∙ ∙♦MORTRAN模型∙ ∙♦LOWTRAN模型∙ ∙♦ATCOR模型∙ ∙♦6S模型等●基于简化辐射传输模型的黑暗像元法●基于统计学模型的反射率反演；相对大气校正常见的是：●基于统计的不变目标法●直方图匹配法等。
既然有怎么多的方法，那么又存在方法选择问题。
这里有一个总结供参考：1、如果是精细定量研究，那么选择基于基于辐射传输模型的大气校正方法。
2、如果是做动态监测，那么可选择相对大气校正或者较简单的方法。
3、如果参数缺少，没办法了只能选择较简单的方法了。
2 ENVI大气校正功能在ENVI中包含了很多大气校正模型，包括基于辐射传输模型的MORTRAN模型、黑暗像元法、基于统计学模型的反射率反演。
基于统计的不变目标法可以利用ENVI一些功能实现。
其中MORTRAN 模型集成在ENVI大气校正扩展模块中。
还有直方图匹配等。
2.1 简化黑暗像元法大气校正黑暗像元法是一种古老、简单的经典大气校正方法。
flaash大气辐射校正
flaash大气辐射校正-李修楠制作一，把DN值转化为亮度1打开landsat5号星图像，7 4 2 波段分别赋值RGB2 单击basic tools >preprocessing>calibration utilities>landsat calibration选择六个反射通道，单击ok3.选择radiance，选择文件保存，等待处理结果处理完得到六个反射通道的亮度值二。
把亮度转化为bip格式1 单击basic tools>convert data(BSQ BIL BIP )然后弹出一个界面选择亮度值，单击ok然后弹出一个界面选择一种格式我选bip，选择文件保存等待处理结果。
三。
FLAASH1 单击spectral>FLAASH2出现以下界面，单击input radiance image选择 bip 格式的那个数据然后ok3再然后出现一个窗口，尺度转换因子，选择单一的转换因子，为10 如图4 。
单击 output reflectance file 选择文件存放位置命名。
5.选择这一景影像的中心经纬度。
如图6，选择传感器类型单击nkown—M> multispectral>landsat TM 57（1）设置获得图像的时间，2006年的8月17日 2时30分（2）大气模型 atmospheric model 选择中纬度夏季mid latitude summer （3）气溶胶选择农村aerosol model 》rural（4）气溶胶反演两个波段k-T aersol retrieval》 2-band（k-T）如图8 单击 Multispetral Settings 出现一个界面单击 kaufman-Tanre aerosol retrieval 选择第一个如图然后okOk完之后单击apply 等待结果结果如图最后，加载校正后的图同样7 4 2 rgb。
进行地理链接，调出光谱曲线比较二者不同。
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大气校正操作流程
利用 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 里面。运行如下：
FLAASH大气校正参数设置
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 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）。
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进行编辑。
EO-1Hyperion高光谱数据FLAASH模块大气纠正及评价
２１ ０ １年 １ ２月
测 绘 与 空 间地 理 信 息
Ｇ ＯＭＡＴ Ｃ ＆ Ｓ ＡＴ Ａ ＮＦ ＭＡＴ Ｏ Ｔ ＨＮＯ ＯＧ Ｅ ＩＳ Ｐ Ｉ Ｌ ｌ ＯＲ Ｉ Ｎ ＥＣ Ｌ Ｙ
Ｖ０ ． ４。 ｏ ６ １３ Ｎ ． Ｄ ｅ ，２ １ ｅ． ０１
Ｅ 一１ ｐ ｒ ｎ高 光 谱 据 Ｆ Ａ Ｈ 模块 Ｏ Ｈｙ ｅｉ ｏ 数 Ｌ ＡＳ 大 气 纠正 及 评 价
Ｏ 引 言
太 阳辐 射在 真空 中传 播 时几 乎 不 受 大气 的 影 响 。然 而， 当它 与 地球 大气 交 互作 用 时 ， 发生 选择 性 散射 和 吸 会 收， 对 于 定 量 研 究 和 分 类 精 度 造 成 很 大 的 影 响。 这 Ｆ Ａ Ｓ 是 目前 精度 最高 的大气 辐射 校正 模 型 ， Ｌ ＡＨ 能有 效 地 消 除大 气 和光 照 等 因素 对地 物 反 射 的 影 响 ， 得 地 物 较 获 为 准确 的反射 率 。本文对 新 疆 地 区 Ｅ Ｏ一１Ｈ ｐｒ ｎ 光 ｙｅ ｏ 高 ｉ 谱 数据 通过 Ｆ Ａ Ｓ 模 块 进 行 大 气 校 正 ， 后 续 的 进 一 Ｌ ＡＨ 为 步研 究 提供 支持 。
周 清 ，祝 民强
（ 东华理工大学 江西省数字国土重点实验室 ， 西 抚 州 ３４０ ） 江 ４００
摘
要： 由于受到大 气的影响 ， 传感 器接收 到 的辐 射信 息不能真 实地反 映地 表反射 光谱信 息 ， 因此 ， 遥感影像 从
中去除大气的影响 ， 即进行 大气校 正， 高光谱遥 感数据 处理 中极 为重要 的环 节。文章介绍 了 Ｅ 是 Ｏ一１ｙｅ ｏ ｈｐｒ ｎ高 ｉ
关键词 ： 高光谱 遥感影像 ； 大气校正 ；Ｌ Ａ Ｈ； ｙｅ ｏ 数据 Ｆ Ａ Ｓ Ｈ ｐｒ ｎ ｉ
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、结果
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进行编辑。
ENVI-Flaash大气校正操作流程
ENVI-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大气校正常见错误及解决方法本文汇总了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）。
flaash 大气校正方法
flaash 大气校正方法相对于常规的校正方法，大气校正在航空、卫星、天文等领域具有广泛的应用。
因为大气对于光的传播和反射都具有很大的影响，如果不考虑大气光学特性，就会产生误差和偏差，从而影响数据处理的结果，降低解析度和精度。
因此，大气校正方法就显得尤为重要。
而FLASH（Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes）算法便是一种可以高效、准确地实现大气校正的方法。
FLASH算法的基本思路FLASH算法是一种基于高光谱数据的大气校正方法，可以高效地实现光学图像数据获取和大气校正。
它是利用大气分子对光的影响，分析高光谱数据中每个波长的传输情况，采用迭代法计算大气反射率和真实地表反射率，从而实现大气校正的目的。
具体而言，FLASH算法需要以下过程：预处理在开始校正之前，需要对光谱数据进行一些预处理操作，包括去除扫描仪噪声、零偏调节、波长校正、光谱辐射校正等。
大气光传输模型FLASH算法采用了光的辐射传输模型来模拟光的透过和反射过程，以计算出大气反射率和地表反射率。
在模型中，光的传输方式可以用以下公式表示：I(l) = I0(e^(-τλ)T0(λ) + (1-e^(-τλ))Tg(λ)ρg(λ)ρs(λ))其中，I(l)是波长为λ时观测值，I0是真实的地物辐亮度值，τλ是指定波长λ处的光学厚度，T0(λ)是大气透过率，Tg(λ)是地表透过率，ρg(λ)是大气反射率，ρs(λ)是地表反射率。
根据传输模型，可以得到以下式子：此外，为了简化计算，该式还可以表示为：其中，S(λ) = I(l) / I0Tg(λ)为观测值与真实值的比例。
迭代计算用上面的模型可以求出大气反射率，但它的准确程度受到不确定因素的影响。
因此，需要进行迭代计算来修正误差。
FLASH采用了快速定量迭代算法（Fast Quantitative Iterative Algorithm，FQIA），以计算整个图像中每个像素的大气校正值。
FLAASH大气校正模块
FLAASH大气校正模块（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中的辐射传输计算方法。
我们可以选取代表研究区的大气模型和气溶胶类型。
FLAASH模型中输入的图像必须是经过辐射定标后的辐射亮度图像，格式为BIL或BIP,数据类型为floating-point，4-byte signed integers，2-byte signed integers，或2-byte unsigned integers，为了能进行大气反演，图像至少包括下面三个范围区间内的15nm分辨率甚至更高的波段，即1050-1210nm，770-870 nm，870-1020nm。
对于已有传感器类型的高光谱遥感图像来说，图像头文件中必须包含波长和波谱带宽（FWHM）。
对于已知的多光谱传感器来说，仅仅需要波长，而未知的多光谱传感器类型，则要求知道波谱响应函数，如论文中ALI数据校正就需要输入波谱响应函数。
环境小卫星多光谱数据FLAASH精确大气校正方法.
环境小卫星多光谱数据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大气校正常见错误及解决方法
本文汇总了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)。
下一页
发帖
打印
123
专题查看
IDL Error：End of input record encountered on file unit:0.
平均海拔高程太大。
注意：填写影像所在区域的平均海拔高程的单位是km：Ground Elevation（Km）。
3.ACC error：avrd:
IDL error：Unable to allocate memory:to make array
Not enough space
ACC_AVRD
为了能处理大数据，ENVI采用分块计算的方式，这个提示是分块（Tile）太大了。
在高级设置里面（Advanced Settings），tile size：它默认是File-> preferences -
>miscellaneous:cache的大小，这个值正常是1~4Mb（0背景很少的情况下）;如果0背景较多，这个值还是需要设置大一些，比如100~200m。
4. ACC error:avrd:No nonblank pixels found
IDL error:
OPENR: Error opening file. Unit: 100,
为了能处理大数据，ENVI采用分块计算的方式，当Tile设置太小，而且有背景值（0），就会出现一个Tile中全部为0的情况，提示这个错误信息。
在高级设置里面（Advanced Settings），tile size：设置稍微大一些，如100~200M等。
5. ACC error:lsmooth2:
IDL error:
ACC_LSMOOTH2:Cannot continue with smoothing calculation
ENVI的FLAASH提供领域纠正功能，但是MODIS、AVHRR等图像分辨率比较低，领域效应区分不出来。
解决方法是在在高级设置里面（Advanced Settings），将领域纠正（Use Adjacency Correction）设置为No。
6. ACC error:modrd5:Nonfinite numbers in coefficient array coef
IDL error:
CDRIVER4V3R2:>>>Wait for MODTRAN4 calculation to finish…
提示传入MODTRAN模型参数有误，常常是由于太阳高度角太小或者太大引起的。
提示这个错误之前会出现以下提示框。
解决方法是确认填写的影像中心位置经纬度信息（西经为负数、南纬为负数）、影像成像时间（格林威治时间）是否正确。
结果错误
1.结果中某一个波段或者多个波段全部为0或者负值
主要产生的原因是输入的辐射亮度数据值偏小。
可能有以下几种情况：
(1)在传感器定标的时候选择的是表观反射率（Reflectance）而不是辐射率数据（Radiance）；
(2)没有做传感器定标，即没有将DN值转换为辐射率数据；
(3)选择了错误的波谱响应函数；
(4)用BandMath做了辐射亮度的单位换算，在FLAASH中导入辐射亮度数据时，Scale Factor选择的不是1.0000。
2.结果中极大值、极小值非常多，也就是0~10000之外的值。
当选择RGB假彩色显示的时候，出现花花绿绿的情况。
辐射定标得到的辐射率数据单位与FLAASH要求的单位不一样。
可能有以下几种情况：
(1)用BandMath做了辐射亮度的单位换算，在FLAASH中导入辐射亮度数据时，Scale Factor选择的不是1.0000。
(2)没有做单位换算。
3.结果中部分像元为负值
这个属于正常现象。
FLAASH是采用MODTRAN辐射传输模型模拟成像中的大气过程，而且很多大气属性都是通过图像来估算，加上大气组成的非均一性，即使MODTRAN4模型精度很高，也不能完全表达大气辐射传输的真实状况。
当影像上有强吸收或者高反射地物时候，就会出现部分像元为负值（如深水、高密度）或者大于10000。
解决方法可以手动修改，如用周围的像元的平均值代替，可以使用ENVI下的The DEM Editing Tool工具。
4.结果图像以RGB显示比原图像视觉效果要差，如模糊。
这个需要了解下遥感软件RGB显示机制。
遥感软件为了让遥感图像显示更加“艳丽”，方便解译，默认会对图像进行拉伸显示，一般是2%的线性拉伸。
也就是我们常常会看到图像值有两个，一个是Scrn值，也就是拉伸之后的显示值，一个是Data值，也就是原始的图像DN值。
另外一个方面，目前我们的RGB加色法显示都是基于8bit显示，也就是0~255。
FLAASH大气校正之后的结果是16bit的整型，而且存在一些极小、极大值，这些对直方图整体形状有一定的影响，影响拉伸效果。
解决方法是利用ENVI下的Interactive stretching工具，选择有效值范围进行拉伸。
还有一种情况是校正图像有很多背景值，比如经过几何校正的整景TM影像。
背景及图像边缘处的像元在大气校正之后变成负值或者0值，由于这样的像元数量多，对图像的整体拉伸影响很大。
解决方法就是将这些像元掩膜掉，如在ROI Tool中，利用Option->Band Threshold to ROI建立一
个ROI进行图像裁剪。
5.结果图像为什么像元值大多是大于1
收藏分享影音更酷！马上访问
没有个人空间？
ESRI 关注
企业级
下一页。