FLAASH大气校正常见错误及解决方法
- 格式:doc
- 大小:301.00 KB
- 文档页数:8
[转载]大气校正(转)大气校正是定量遥感中重要的组成部分。
本专题包括以下内容:∙ ●大气校正概述∙∙●ENVI中的大气校正功能1大气校正概述大气校正的目的是消除大气和光照等因素对地物反射的影响,广义上讲获得地物反射率、辐射率或者地表温度等真实物理模型参数;狭义上是获取地物真实反射率数据。
用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等物质对地物反射的影响,消除大气分子和气溶胶散射的影响。
大多数情况下,大气校正同时也是反演地物真实反射率的过程。
很多人会有疑问,什么情况下需要做大气校正,我们购买或者其他途径获取的影像是否做过大气校正。
通俗来讲,如果我们需要定量反演或者获取地球信息、精确识别地物等,需要使用影像上真实反映对太阳光的辐射情况,那么就需要做大气校正。
我们购买的影像,说明文档中会注明是经过辐射校正的,其实这个辐射校正指的是粗的辐射校正,只是做了系统大气校正,就跟系统几何校正的意义是一样的。
常见的绝对大气校正方法有:●基于辐射传输模型∙ ∙♦MORTRAN模型∙ ∙♦LOWTRAN模型∙ ∙♦ATCOR模型∙ ∙♦6S模型等●基于简化辐射传输模型的黑暗像元法●基于统计学模型的反射率反演;相对大气校正常见的是:●基于统计的不变目标法●直方图匹配法等。
既然有怎么多的方法,那么又存在方法选择问题。
这里有一个总结供参考:1、如果是精细定量研究,那么选择基于基于辐射传输模型的大气校正方法。
2、如果是做动态监测,那么可选择相对大气校正或者较简单的方法。
3、如果参数缺少,没办法了只能选择较简单的方法了。
2 ENVI大气校正功能在ENVI中包含了很多大气校正模型,包括基于辐射传输模型的MORTRAN模型、黑暗像元法、基于统计学模型的反射率反演。
基于统计的不变目标法可以利用ENVI一些功能实现。
其中MORTRAN 模型集成在ENVI大气校正扩展模块中。
还有直方图匹配等。
2.1 简化黑暗像元法大气校正黑暗像元法是一种古老、简单的经典大气校正方法。
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。
进行地理链接,调出光谱曲线比较二者不同。
第13章高光谱数据FLAASH大气校正本节以AVIRIS高光谱数据为数据源,介绍高光谱数据的FLAASH大气校正过程。
13.1 浏览高光谱数据此AVIRIS高光谱数据为经过传感器定标的辐射亮度数据。
(1)在ENVI主菜单中,选择File→Open Image File,打开JasperRidge98av.img文件。
(2)在波段列表中,选择JasperRidge98av.img,单击右键选择Load True Color,在Display 窗口中显示真彩色合成图像。
(3)在主图像窗口中单击右键,快捷菜单中选择Pixel Locator。
设置Sample:366,Line:179。
此像元为硬质水泥地,吸收特征主要受大气的影响,单击Apply按钮。
(4)在主图像窗口中单击右键,快捷菜单中选择Z Profile,打开Spectral Profile窗口,绘制像素(366,179)的波谱剖面。
(5)在Spectral Profile窗口中,可以看到在760nm,940nm和1135nm处,水汽具有吸收特征,1400nm和1900nm附近基本没有反射能量,二氧化碳在2000nm附近有两个吸收特征。
13.2 AVIRIS数据大气校正(1)在ENVI主菜单中,选择Spectral→FLAASH,打开FLAASH Atmospheric Correction Model Input Parameters对话框。
(2)单击Input Radiance Image按钮,选择JasperRidge98av.img文件。
在Radiance Scale Factors对话框中,选择Read array of scale factors from ASCII file,单击OK按钮。
(3)选择AVIRIS_1998_scale.txt文件,按照默认设置,单击OK按钮。
(4)单击Output Reflectance File按钮,选择输出路径及文件名JasperRidge98av.img。
flaash大气校正Flaash大气校正(IRSP6-08.3.24)大气校正的目的是消除大气和光照等因素对地物反射的影响,获得地物反射率和辐射率、地表温度等真实物理模型参数,用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧对地物反射的影响,消除大气分子和气溶胶散射的影响。
FLAASH可以处理任何高光谱数据、卫星数据和航空数据(860nm/1135nm),这些数据是由HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)AISA、HARP、DAIS、Probe-1、TRWIS-3、SINDRI、MIVIS、OrbView-4、NEMO 等传感器获得的。
FLAASH还可以校正垂直成像数据和侧视成像数据。
Flaash大气校正使用了MODTRAN 4+辐射传输模型的代码,基于像素级的校正,校正由于漫反射引起的连带效应,包含卷云和不透明云层的分类图,可调整由于人为抑止而导致的波谱平滑。
FLAASH可对Landsat, SPOT, AVHRR, ASTER, MODIS, MERIS, AATSR, IRS等多光谱、高光谱数据、航空影像及自定义格式的高光谱影像进行快速大气校正分析。
能有效消除大气和光照等因素对地物反射的影响,获得地物较为准确的反射率和辐射率、地表温度等真实物理模型参数。
校正过程点击envi——Basic Tools -Preprocessing - Calibration Utilities -FLAASHSpectral -FLAASH.或者点击envi-spectral- FLAASH1、输入数据必须是辐射校正后的数据,对辐射校正数据转成BIL或BIP格式(Basic Tools ——Convert Data);2、对输入数据进行头文件编辑,主要是对波长wavelenth(即每一波段的波长中心值)和波长宽度fwhm(每一波段的波长范围)的编辑。
不是高光谱数据可以不对fwhm进行编辑。
ETM Flaash大气校正在ETM蚀变信息提取2文章里,已经得到了辐亮度图像,现在使用ENVI Flaash大气校正功能,进行地表反射率的反演。
1、数据格式转换Flaash大气纠正功能需要的辐亮度文件格式为BIL或者BIP,但是目前的文件格式是BSQ的,所以需要进行数据格式转换。
ENVI功能菜单: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大气校正下高分二号数据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按钮打开。
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进行编辑。
ENVI5.2下高分二号数据FLAASH 大气校正(2015-07-09 15:57:40)转载▼ 标签:envi5.2高分二号gf2flaash 大气校正 分类: ENVI高分二号卫星于2014年8月19日成功发射,搭载有两台高分辨率1米全色、4米多光谱相机,将带来优于1米空间分辨率的光学遥感影像。
高分二号较高分一号来说分辨率提高一倍,同时具有高定位精度和快速姿态机动能力等特点。
从2014年8月21日首次开机成像并下传数据,已逐步被各行业用户使用。
高分二号数据的大气校正与高分一号类似,但由于官方暂时未正式公布卫星的波谱响应函数等参数,ENVI 也未能及时对其进行原生支持。
为方便大家使用,这里以一景GF2-PMS2 L1A 级数据为例,介绍在ENVI5.2下可行的FLAASH 大气校正流程。
注:同样适合在其他ENVI 版本中操作。
高分二号卫星轨道和姿态控制参数及有效载荷技术指标见下表1、2: 表1 高分二号卫星轨道和姿态控制参数启动ENVI5.2;依次File > Open或直接单击工具栏上的图标,弹出Open对话框,选择数据文件夹下扩展为.tiff的文件,然后点击Open按钮打开(本例中为…/GF2_PMS2_E115.7_N42.7_20140928_L1A0000362235-MSS2.tiff)。
说明:GF-2卫星绝对定标系数可从如下地址下载,/n16/n1115/n1522/n2103/193529.html图1使用Apply Gain and Offset工具进行辐射定标3. FLAASH大气校正FLAASH大气校正需要影像的中心波长信息,ENVI暂不能自动识别GF2数据的头文件信息,因此首先需要手动添加中心波长信息。
(1)添加中心波长在Toolbox中,依次Raster Management > Edit ENVI Header,弹出Edit Header Input File对话框,在Select Input File选项卡中选择上一步辐射定标后的结果(本例中为GF2_PMS2_E115.7_N42.7_20140928_L1A0000362235-MSS2_Rad.dat);弹出Header Info对话框,点击Edit Attributes,选择Wavelengths,弹出Edit Wavelength values对话框(图2),依次填入各波段对应中心波长,点击OK。
flaash大气校正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.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 大气校正方法相对于常规的校正方法,大气校正在航空、卫星、天文等领域具有广泛的应用。
因为大气对于光的传播和反射都具有很大的影响,如果不考虑大气光学特性,就会产生误差和偏差,从而影响数据处理的结果,降低解析度和精度。
因此,大气校正方法就显得尤为重要。
而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精确大气校正方法环境与灾害监测预报小卫星星座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 nm
2、像元值类型:经过定标后的辐射亮度(辐射率)数据,单位是:(μ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:convert7
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<?xml:namespace prefix = o ns =
"urn:schemas-microsoft-com:office:office" />
IDL error:
CDRIVER4V3R2:>>>Wait for MODTRAN4 calculation to finish…
<?xml:namespace prefix = v ns = "urn:schemas-microsoft-com:vml" /> 提示传入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
ENVI FLAASH考虑到数据储存和后续处理,将大气校正得到的反射率结果乘以10000变成16bit整型。
如果想让反射率结果在0~1范围,可用BandMath,表达式为b1/10000.0。