将数据拷贝到文本文件中。如下图中为波谱响应函数数字表达的一部分,第一列表示波长,后面五列分别表示一个全色和4个多光谱对应波长的波谱响应值。
图2高分一号波谱响应值
在ENVI中,使用波谱曲线来描述波谱响应函数,也就是以波长作为x轴,波谱响应值作为y轴,存储格式为ENVI波谱库文件(.sli)。下面介绍ENVIclassic中的波谱响应函数的制作。
(1) 启动ENVI Classic
(2) 选择Window->Start New Plot Window,ENVI PlotWindow窗口中,选择File->InputData->ASCII,如图3所示,自动将第一列作为X轴,后面3-6列作为Y轴。单击OK。
(3) 如图所示,生成了4条曲线。选择Edit->DataParameters,更改每一条曲线的名称:b1,b2,b3,b4,便于区分。
(4) 选择File- >Save Plot As->SpectralLibrary,将波谱曲线保存为波谱库文件.sli。

图3 导入ASCII文件

图4 高分一号波谱响应函数的曲线表达
注:如果出现以下情况,原因是没有将多光谱波段1和波段2空白部分填补0。高分一号

图5 错误情况
3 FLAASH大气校正 (1) 工具箱/Radiometric Correction/Atmospheric CorrectionModule/FLAASH Atmospheric Correction打开FLAASH大气校正模块;
(2) 点击Input Radiance Image,前面处理好的数据,在Radiance Scale Factors面板中选择Usesingle scale factor for allbands,由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍,所以在此Single scalefactor选择默认:10,单击OK;
注:定标后的单位是W⋅m−2⋅sr−1⋅μm−1,与FLAASH要求的单位(μW)/(cm2*nm*sr)相差10倍关系,因此在RadianceScale Factors中输入10缩放系数。
(3) 设置输出文件及路径设置;
(4) 传感器基本信息设置:
l 成像中心点经纬度FLAASH自动从影像中获取。
l 传感器高度(Sensor Altitude):645km
注:写这篇文档时没有得到官方正式数据,下面截图是公布之前填写的505,最终公布的数据为645km。
l 像元大小(pixel Size):8m
l 成像区域平均高度可以通过统计DEM数据获取
l 成像时间:从数据头文件中读取(.XML),位置为:2013-09-08 13:14:07;减去8换算GMT时间
(5) 大气模型和气溶胶模型,根据经纬度和影像区域选择。
(6) 气溶胶反演方法选择None(缺少短波红外),能见度设置为40km。
图6FLAASH大气校正参数设置
(7) 单击Multispectral Setting按钮,在Filter Function File导入光谱响应曲线“gf-1pms.sli”,单击OK;
(8) 单击Advanced Settings,在高级设置中
l Tile Size 默认的是Cash size 的大小,手动改为50-100Mb(根据内存大小设定),单击OK;
l 从xml文件中得到: 97.0867
79.1983
讨论:一般中低分辨率的影像近似天顶角:180和方位角:0,即垂直观测。ENVI FLAASH的天顶角:90-180度,方位角:-180and 180,没有查到高分一号的角度说明,这里选择默认。

图7 大气校正高级参数设置
(9) 设置好后,在大气校正模块面板中,单击Apply。
(10)大气校正完成后,检查大气校正的结果,分别加载校正前后的图像(选择CIR假彩色方式加载,可以更好的识别植被),查看典型地物的大气校正前后的光谱曲线。
4 总结从上面可以看到,对于未知的多光谱传感器的FLAASH大气校正,关键是获取波谱响应函数,而波谱响应函数一般都是公开的。高分一号类似的方法可用于高分一号宽幅传感器的FLAASH大气校正。
本文来自电脑杂谈,转载请注明本文网址:
http://www.pc-fly.com/a/shumachanpin/article-32105-2.html
过期商品支接罚50000