本扩展工具要求 ENVI 5.7 及以上版本。低版本 ENVI 可以使用如下扩展工具:https://www.cnblogs.com/enviidl/p/16393415.html
自 ENVI 5.7 版本开始,FLAASH 大气校正功能提供了官方 Task 接口,详细信息可查看ENVI帮助内ENVI > Programming > ENVITasks > List of Tasks > FLAASH 章节内容。
由于 FLAASH 工具参数较多,对于初学者不是很友好,所以开发了易用版工具,基于全新 FLAASH Task 实现,隐藏了很多不重要和无需设置的参数,并尽可能多的自动设置了必需参数。
GF1 PMS、GF1 WFV
GF1B/C/D PMS
GF2 PMS
GF4 PMS(如果输入5波段PMS数据,将只输出B/G/R/NIR四波段结果)
GF5 AHSI
GF5/A/B VIMI(仅支持B1~B6波段数据)
GF5/A/B AHSI
GF6 PMS、GF6 WFV
GF7 MUX
ZY3-01 MUX、ZY3-02 MUX、ZY3-03 MUX
ZY1-02C MUX
ZY1-02D(ZY1E) MUX、ZY1-02D(ZY1E) AHSI
ZY1-02E(ZY1F) MUX、ZY1-02E(ZY1F) AHSI
HJ1 CCD、HJ2 CCD
CBERS04 MUX
CBERS04A MSS/MUX/WFI
BJ3N1/2 MS
OHS(ZhuHai-1) CMOSMSS
Landsat-5 TM、Landsat-7 ETM+、Landsat-8 OLI、Landsat 9 OLI-2
Sentinel-2A/B MSI
1 工具安装
方法1:使用 ENVI App Store 进行 ENVI 扩展工具的安装与管理。App Store 下载地址:https://envi.geoscene.cn/appstore
方法2:手动下载安装 https://envi.geoscene.cn/appstore/newlaesuse
2 工具使用
/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]
/辐射校正/大气校正模块/FLAASH 易用版 [新]
在弹出的对话框中,选择输入栅格。可直接选择原始L1级别数据(要求包含定标系数,建议输入此种数据);也可以选择辐射定标后的辐亮度图像。对于存储顺序(Interleave)不再做要求。
Acquisition Time(拍摄时间):自动从元数据中获取,如果没有识别,需手动设置。
Sensor Type(;传感器类型):大多数传感器会自动识别,如果没有识别,需手动设置。
Input Scale(缩放系数):对于包含定标系数的L1级别数据,自动设置为10;对于使用辐射定标工具处理为辐射亮度值图像,自动根据定标的Scale Factor进行设置;其他情况需手动设置。
Aerosol Model(气溶胶模型):需手动设置,根据图像覆盖地物类型选择。
Initial Visibility(能见度):当图像包含SWIR2波段时,将自动进行气溶胶反演,此时本参数将不起作用;其他数据需根据成像天气条件进行设置,默认40km。
Output Raster(输出栅格):设置大气校正输出结果文件名。
3 处理示例
测试环境:
笔记本:ThinkPad P52
操作系统:Windows 11 64-bit
软件版本:ENVI 6.0
处理器:Intel(R) Core(TM) i7-8850H CPU @ 2.60GHz 六核
内存:64.00 GB
硬盘:SSD
3.1 GF1 PMS数据处理
在主菜单中,选择File > Open As > China Satellites > GaoFen > GF-1,打开GF1_PMS_*L1A*_MSS.xml文件。
注:需要安装“中国国产卫星支持工具”。可在ENVI App Store中安装,或手动安装(https://envi.geoscene.cn/china)。
在工具箱中,启动/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]工具,在弹出对话框中选择多光谱数据,如下图所示:
在工具面板,只需设置Aerosol Model(气溶胶模型)、Initial Visibility(能见度)。设置输出路径,点击OK按钮,执行FLAASH大气校正。如下图所示。
处理结果自动加载在视图中,对比大气校正前(左图)后(右图)植被光谱曲线如下图所示。
3.2 GF5B VIMI数据处理
仅支持 B1~B6 波段数据处理。
在主菜单中,选择File > Open As > China Satellites > GaoFen > GF-5,打开GF5B_VIMI_*L1*.xml文件。
注:需要安装“中国国产卫星支持工具”。可在ENVI App Store中安装,或手动安装(https://envi.geoscene.cn/china)。
在工具箱中,启动/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]工具,在弹出对话框中选择多光谱数据,如下图所示:
在工具面板,只需设置Aerosol Model(气溶胶模型)、Initial Visibility(能见度)。设置输出路径,点击OK按钮,执行FLAASH大气校正。如下图所示。
处理结果自动加载在视图中,对比大气校正前(左图)后(右图)植被光谱曲线如下图所示。
3.3 ZY3-03 MUX数据处理
在主菜单中,选择File > Open As > China Satellites > ZiYuan > ZY3-03A,打开ZY303_PMS_*_L1A*-MUX.xml文件。
注:需要安装“中国国产卫星支持工具”。可在ENVI App Store中安装,或手动安装(https://envi.geoscene.cn/china)。
在工具箱中,启动/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]工具,在弹出对话框中选择多光谱数据,如下图所示:
在工具面板,只需设置Aerosol Model(气溶胶模型)、Initial Visibility(能见度)。设置输出路径,点击OK按钮,执行FLAASH大气校正。如下图所示。
处理结果自动加载在视图中,对比大气校正前(左图)后(右图)植被光谱曲线如下图所示。
3.4 ZY1-02E(ZY1F) AHSI数据处理
在主菜单中,选择File > Open As > China Satellites > ZiYuan > ZY1-02E (ZY1F),打开ZY1F_AHSI_*_L1A*.xml文件。
注:需要安装“中国国产卫星支持工具”。可在ENVI App Store中安装,或手动安装(https://envi.geoscene.cn/china)。
在工具箱中,启动/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]工具,在弹出对话框中选择高光谱数据,如下图所示:
在工具面板,只需设置Aerosol Model(气溶胶模型)、Initial Visibility(能见度)。设置输出路径,点击OK按钮,执行FLAASH大气校正。如下图所示。
处理结果自动加载在视图中,对比大气校正前(左图)后(右图)植被光谱曲线如下图所示。
3.5 Landsat-8 OLI数据处理
在主菜单中,选择File > Open,打开*_MTL.txt或xml文件。
在工具箱中,启动/Radiometric Correction/Atmospheric Correction Module/FLAASH Easy-to-Use [New]工具,在弹出对话框中选择多光谱数据,如下图所示:
在工具面板,只需设置Aerosol Model(气溶胶模型)、Initial Visibility(能见度)。设置输出路径,点击OK按钮,执行FLAASH大气校正。如下图所示。
处理结果自动加载在视图中,对比大气校正前(左图)后(右图)植被光谱曲线如下图所示。
4 Task调用示例
file = 'C:\temp\GF1_PMS2_Exxx_L1Axxx-MSS2.xml'
raster = ENVIOpenChinaRaster(file)
Task = ENVITask('NewFLAASHEasyToUse')
Task.INPUT_RASTER = raster
;如果文件名规范,则不需要设置传感器类型
;Task.SENSOR_TYPE = 'GF1 PMS'
;如果包含增益偏移系数,或输入为辐射定标结果,则不需要设置此参数
;Task.INPUT_SCALE = 10.0
;此参数为必选项
Task.AEROSOL_MODEL = 'High-Visibility Rural'
;如果不包含SWIR2波段,需要设置此参数
Task.DEFAULT_VISIBILITY = 40
;输出文件路径
Task.OUTPUT_RASTER_URI = outfile
Task.Execute
OUTPUT_RASTER = Task.OUTPUT_RASTER ;获取输出结果
;如果输入图像不包含拍摄时间,则需手动设置
timeObj = ENVITime(acquisition='2001-01-01T00:30:45.3Z')
Task.ACQUISITION_TIME = timeObj
5 ENVI Modeler中进行批处理
对于L1标准格式或标准命名的辐亮度图像,可以同时进行不同传感器数据的批量处理,输入图像的气溶胶模型需一致。
如下示例为1景GF1 PMS和1景GF2 PMS进行FLAASH批处理。
下载链接:https://pan.baidu.com/s/1jB7XMZB7z5GcXLA3Vcr9gw?pwd=envi
注:需要安装“自定义 ENVITasks 用于 ENVI Modeler”。可在ENVI App Store中安装,或手动安装(https://envi.geoscene.cn/appstore/cusmsfmler)。