高分一号 PMS相机的多光谱数据拥有很好的空间分辨率(8米)和幅宽(60km),包括蓝、绿、红、近红外4个波段。本文介绍利用ENVIFLAASH工具完成高分一号PMS多光谱图像的大气校正,主要分为以下三步:
l数据定标
l波谱响应函数制作
lFLAASH大气校正
注:本文是在ENVI5.0下和ENVIclassic完成,其他版本基本类似。
1 数据定标
利用以下公式可将GF-1卫星各载荷的通道观测值计数值DN转换为卫星载荷入瞳处等效表观辐亮度数据。
式中:Gain为定标斜率;DN为卫星载荷观测值;Bias为定标截距。
定标参数参照表1,来源于中国资源卫星应用中心:http://www.cresda.com/n16/n1115/n1522/n2103/index.html
表1 高分一号星PMS相机的定标系数
卫星载荷 | 波段号 | Gain | Bias |
PMS1 | PAN | 0.1886 | -13.127 |
Band1 | 0.2082 | 4.6186 | |
Band2 | 0.1672 | 4.8768 | |
Band3 | 0.1748 | 4.8924 | |
Band4 | 0.1883 | -9.4771 | |
PMS2 | PAN | 0.1878 | -7.9731 |
Band1 | 0.2072 | 7.5348 | |
Band2 | 0.1776 | 3.9395 | |
Band3 | 0.177 | -1.7445 | |
Band4 | 0.1909 | -7.2053 | |
WFV1 | Band1 | 0.1709 | -0.0039 |
Band2 | 0.1398 | -0.0047 | |
Band3 | 0.1195 | -0.0030 | |
Band4 | 0.1338 | -0.0274 | |
WFV2 | Band1 | 0.1588 | 5.5303 |
Band2 | 0.1515 | -13.642 | |
Band3 | 0.1251 | -15.382 | |
Band4 | 0.1209 | -7.985 | |
WFV3 | Band1 | 0.1556 | 12.28 |
Band2 | 0.1700 | -7.9336 | |
Band3 | 0.1392 | -7.031 | |
Band4 | 0.1354 | -4.3578 | |
WFV4 | Band1 | 0.1819 | 3.6469 |
Band2 | 0.1762 | -13.54 | |
Band3 | 0.1463 | -10.998 | |
Band4 | 0.1522 | -12.142 |
进行数据定标计算时候,我们可以有两种方法选择使用:
一、Bandmath
以下是利用Bandmath工具进行定标,这种方法灵活,不过需要一个一个波段定标。
1)有了绝对定标参数和定标公式,选择Toolbox/Band Ratio/BandMath工具很容易进行传感器定标。
2)由于是单个波段文件的定标,选择Toolbox/Raster Management/LayerStacking将定标后的单波段文件组合成一个多波段文件。
二、Apply Gain and Offset工具
我们可以使用Apply Gain and Offset工具。。
1)启动Toolbox/Radiometric Correction/Apply Gain andOffset。
2)在Gain and Offset Values中填入增益Gain值和偏移Offset。
3)选择输出路径和文件名执行定标。
图1Gain and Offset Values面板
计算好之后,继续执行以下步骤(原数据的储存顺序是BIP):
(1)选择Toolbox/Raster Management/ Edit ENVI Header->Edit Attributes->Wavelength,将每个波段的中心波长输入——b1(503nm),b2(576nm),b3(680nm),b4(810nm)
讨论:由于没有找到官方公布的中心波长信息,这里取与波谱响应中反射率为1的波长。
目前为止,影像数据准备工作完成了。
2 波谱响应函数
波段响应函数,英文名为spectral responsefunction(SRF)或者叫RelativeSpectral Response (RSR) ,与宽波段传感器出现。我们知道每一个波段都有一个波段范围,比如HJ-b1(475nm)波段为520-430纳米,实际上传感器的感光元件在这个波段范围内的每一点所感应的强度都是不一样。在成像中,原则上讲应该根据波段响应函数来进行加权平均,但由于处理起来比较麻烦,而且一般的精度要求不太高,所以大多数图像都是直接取了波段范围内的中点值来运算。波段响应函数是描述一定波长范围内(超出波段范围)的量子效应,当需要精确计算像元响应时候,比如大气校正反演真实地表反射率,就需要使用波谱响应函数。
任何传感器在设计时都会给出严格的波段响应函数,高分一号也不例外。在中国资源卫星应用中心可下载:http://www.cresda.com/n16/n1115/n1522/n2118/186272.html。与资源一号02c和资源三号一样,高分一号也是以.xls格式提供,值得注意的是,第三列和第四例后面有部分无值区,手动填入0。
将数据拷贝到文本文件中。如下图中为波谱响应函数数字表达的一部分,第一列表示波长,后面五列分别表示一个全色和4个多光谱对应波长的波谱响应值。
图2高分一号波谱响应值
在ENVI中,使用波谱曲线来描述波谱响应函数,也就是以波长作为x轴,波谱响应值作为y轴,存储格式为ENVI波谱库文件(.sli)。下面介绍ENVIclassic中的波谱响应函数的制作。
(1)启动ENVI Classic
(2)选择Window->Start New Plot Window,ENVIPlot Window窗口中,选择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)工具箱/RadiometricCorrection/Atmospheric Correction Module/FLAASH AtmosphericCorrection打开FLAASH大气校正模块;
(2)点击InputRadiance Image,前面处理好的数据,在RadianceScale Factors面板中选择Usesingle scale factor for all bands,由于定标的辐射量数据与FLAASH的辐射亮度的单位相差10倍,所以在此Singlescale factor选择默认: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像元大小(pixelSize):8m
l成像区域平均高度可以通过统计DEM数据获取
l成像时间:从数据头文件中读取(.XML),位置为:2013-09-0813:14:07;减去8换算GMT时间
(5)大气模型和气溶胶模型,根据经纬度和影像区域选择。
(6)气溶胶反演方法选择None(缺少短波红外),能见度设置为40km。
图6FLAASH大气校正参数设置
(7)单击MultispectralSetting按钮,在FilterFunction File 导入光谱响应曲线“gf-1pms.sli”,单击OK;
(8)单击AdvancedSettings,在高级设置中
lTile Size 默认的是Cashsize 的大小,手动改为50-100Mb(根据内存大小设定),单击OK;
l从xml文件中得到:97.0867
79.1983
讨论:一般中低分辨率的影像近似天顶角:180和方位角:0,即垂直观测。ENVIFLAASH的天顶角:90-180度,方位角:-180and 180,没有查到高分一号的角度说明,这里选择默认。
图7大气校正高级参数设置
(9)设置好后,在大气校正模块面板中,单击Apply。
(10)大气校正完成后,检查大气校正的结果,分别加载校正前后的图像(选择CIR假彩色方式加载,可以更好的识别植被),查看典型地物的大气校正前后的光谱曲线。
4 总结
从上面可以看到,对于未知的多光谱传感器的FLAASH大气校正,关键是获取波谱响应函数,而波谱响应函数一般都是公开的。类似的方法可用于高分一号宽幅传感器的FLAASH大气校正。