第一作者简介: 张川(1985-),男,硕士研究生,主要从事高光谱遥感数据处理与应用、资源与环境遥感及遥感地质方面的研究。
光谱重建是高光谱数据定量分析的前提,大气校正是光谱重建的关键环节。对于缺少地面同步测量的高光谱数据,大气辐射传输模型是最可行的大气校正方法。应用当前流行的大气辐射传输模型MODTRAN 4.0和6S,分别对试验区EO-1 Hyperion高光谱数据进行大气校正和光谱重建。通过对植被、粘土矿及水体3类地物光谱重建结果的对比分析,并利用统计方法计算光谱重建精度,评价大气辐射传输模型校正的效果和模型的适用性。最后,通过FLAASH模型的试验体现了高光谱数据大气校正集成化处理的优越性。
Spectral reconstruction is a prerequisite for quantitative analysis of hyperspectral data, and atmospheric correction is the key step to spectral reconstruction. The atmospheric radiative transfer model is the most feasible method for hyperspectral data without simultaneous measurement of the ground. In this paper, the prevalent atmospheric radiative transfer models MODTRAN4.0 and 6S were used respectively for atmospheric correction and spectral reconstruction of EO-1’s Hyperion hyperspectral data in the study areas. Then the results of the correction using the atmospheric radiative transfer model and the applicability of the model were evaluated by comparing three types of reconstructed spectra respectively from vegetation, clay mineral and water and by calculating accuracy of spectral reconstruction using the statistical method. The superiority of integrated processing is reflected in the end by the application of FLAASH model.
高光谱遥感是遥感的前沿技术之一。与传统的多光谱遥感相比, 高光谱遥感的显著特点是: 对于每个图像像元, 在光谱维能够形成完整而连续的光谱曲线, 进而能够更为精细地反映像元地物的光谱特征。因此, 运用高光谱遥感技术可以进行更为精细的地物识别和分类[1]。像元光谱特征的准确性是利用高光谱数据进行地物识别的关键和前提条件, 光谱重建的精度决定了高光谱数据的应用效果, 而光谱重建的关键技术是大气校正。国内外研究表明[2], 大气校正方法大致分为基于图像特征的相对校正法、基于地面线性回归模型校正法和基于大气辐射传输模型校正法3类, 而大气辐射传输模型校正法是目前高光谱数据大气校正的常用方法。从20世纪80年代开始, 国外一些学者就开始从事这方面的研究, 已经发展了一系列大气辐射传输模型和相应软件[3, 4], 但是在各模型对高光谱数据大气校正应用效果的评价方面罕有涉及。纵观各类大气辐射传输模型和软件, LOWTRAN~MODTRAN系列和5S~6S系列是最基本的两类模型。本文分别基于MODTRAN 4.0和6S模型对星载高光谱Hyperion数据进行大气校正和光谱重建, 通过对比植被、粘土矿及水体3类地物像元的重建光谱和重建精度, 评价它们的应用效果和适用性。
本研究使用的是Hyperion L1R级辐射定标产品。两个试验区地理位置分别在江西省赣州市中部和东部。试验区Hyperion数据指标如表1所示。
| 表1 试验区Hyperion数据指标 Tab.1 Index of Hyperion in study areas |
大气校正前, 首先去除Hyperion数据集中未经过定标的波段和波长在1.4 μ m和1.9 μ m附近的水气吸收波段, 然后针对USGS对VNIR和SWIR波段所采用过的40和80增益(对DN值进行补偿和改进, 并将数据以整型进行存储), 对VNIR和SWIR数据集分别除以40和80, 将其转化为传感器入瞳辐射值[5, 6]。
坏线和条带是Hyperion图像上经常出现的噪声, 本研究只对明显的坏线和严重的条纹进行平滑处理。图1是第57波段处理前后影像。
大气校正是光谱重建的核心, 基于电磁波在大气中的辐射传输原理建立起来的大气辐射传输模型是理论上精度最高的大气校正方法。在目前的各类模型中, MODTRAN 4.0和6S是最具代表性的, 也是使用最为广泛的两个大气辐射传输模型。
MODTARN(Moderate Resolution Transmission)是美国空军地球物理实验室(AFGL)开发的大气辐射传输模型, 该模型继承了LOWTRAN模型的代码, 并做了适当改进, 目前常用的是MODTRAN 4.0版本[7]; 6S(Second Simulation of the Satellite Signal in the Solar Spectrum)模型是法国大气光学实验室Eric Vemote在5S模式的基础上开发出来的大气辐射传输模型[8]。
从模型的辐射传输方程来看, 在地表条件假设前提相同的情况下, 两者的形式基本一致。式(1)为MODTRAN 4.0在均一朗伯面条件下的辐射传输方程(辐射亮度形式), 即
L(μ v)=L0(μ v)+T(μ v)Fd
式中, L(μ v)为传感器接收的辐射亮度; L0(μ v)为路径辐射, 即程辐射; Fd为太阳下行总辐射; T(μ v)为传感器和目标之间的透过率; ρ 为目标的地表反射率; S为大气半球反照率。
式(2)为同样条件下6S模型的辐射传输方程(反射率形式), 即
ρ T(θ s, θ v, φ s, φ v)=Tg(θ s, θ v)[ρ R+a+T(θ s)T(θ v)
式中, ρ 为地表反射率; ρ T(θ s, θ v, φ s, φ v)为传感器测得的表观反射率; θ s为太阳天顶角; θ v为太阳方位角; φ s为传感器天顶角; φ v为传感器方位角; Tg(θ s, θ v)为大气透过率; T(θ s)为太阳— 目标大气路径透过率; T(θ v)为目标— 传感器大气路径透过率; ρ R+a为分子散射和气溶胶散射所构成的路径辐射反射率; S为大气半球反照率。
2.1.1 共同特点
(1)两个模型包括了透过率计算、多次散射处理和几何路径计算等。
(2)都是以大气和地表条件作为输入参数, 主要输入参数相同。
(3)计算过程中均考虑了水平、垂直、倾斜向上和向下传输等各种复杂的几何关系, 在计算大气倾斜路径及沿着传输路径衰减量时, 都考虑了大气折射和地球的曲率。
(4)大气模型均可以根据理论计算或实测资料, 模型中均包括了具有代表性的大气和气溶胶的模式, 用户可以直接使用标准模式, 也可以自定义大气模式。
(5)在计算地表反射率时, 为了描述地表的均一性和朗伯状态, 都引入了点扩散函数(PSF)与双向反射率分布函数(BRDF)模型。
2.1.2 主要区别
(1)处理谱段范围不同。MODTRAN 4.0的处理谱段范围是0.20 μ m到无穷, 而6S的处理谱段范围是0.25~4.0 μ m。
(2)MODTRAN 4.0可以对云、雨进行处理, 而6S模型没有考虑云、雨的影响。
(3)两者都利用了HITRAN数据库中的基本分子常数, 却采取了各不相同的处置方法把这些常数换算成透过率和辐射度, 因此, 导致两者在理论上具有不同的精确性。
(4)MODTRAN 4.0运行基于Windows平台, 6S模型基于DOS平台, 前者更容易操作。
(5)MODTRAN 4.0的计算过程是正向的, 即由地表反射率→ 传感器入瞳辐射; 6S的计算是反演的过程, 即由传感器入瞳辐射→ 地表反射率。
由于二向反射核驱动模型参数的不确定性, 本文选用朗伯体反射。为了后续的光谱重建评价, 两种模型采用的输入参数类型基本一致。根据现有资料, 以第37波段为例, 两个试验区的参数设置如表2所示。
| 表2 试验区基本参数设置 Tab.2 Parameters of atmospheric correction in study areas |
两种模型的大气校正根据各自的处理顺序逐波段进行计算, 172个Hyperion数据波段处理完毕后, 通过打包处理, 得到试验区的光谱重建数据集。
针对研究区地物类型, 在校正后的Hyperion图像上选取3类地物(植被、粘土、水体)像元反射率曲线进行评价。评价方法是将具有特征光谱的像元反射光谱特征吸收谷波长位置与地面更高光谱分辨率的光谱仪实测同点地物进行特征吸收谷对比, 分析其吸收波长位置的偏移量, 以及偏移量相对于高光谱成像仪光谱分辨率的大小与光谱带宽而言, 是否在一个合理的光谱区间, 如果二者完全吻合, 则图像光谱代表了实际地物的光谱特征[9]。
为了验证大气校正的效果, 利用式(3)计算校正前的表观反射率[10], 即
ρ T=
式中, ρ T为表观反射率; L入为传感器入瞳辐射亮度; D为日地之间距离(天文单位); ESUNλ 为大气层顶的平均太阳光谱辐照度(与波长有关); θ v为太阳的天顶角。
针对植被, 采用USGS标准植被光谱曲线作为光谱重建的评价基准; 针对研究区露头发育的粘土矿, 为了突出地域性, 采集了27个样本, 利用ASD光谱辐射仪(光谱分辨率3~10 nm)进行光谱测试, 用得到的反射率曲线作为粘土矿像元的评价基准。对于试验区光谱重建数据集, 由于同类像元之间也存在光谱差异, 故在光谱重建后的Hyperion图像上寻找面积较大且较均一的植被和粘土矿地区, 提取平均光谱作为光谱重建的评价对象。
按照评价方法进行逐个吸收谷对比评价, 如图2所示。
可以看出, MODTRAN 4.0和6S模型均能将表观反射率转换到地面反射率, 反演的植被、粘土矿反射率与USGS曲线的主要光谱特征基本吻合; 通过图2(b)、(c)和(d)的比较, 粘土矿包含的主要矿物高岭石、蒙脱石在2 165 nm和2 205 nm左右的双吸收特征在两个模型反演结果中均得到较好体现, 两个模型反演的粘土矿反射率在950 nm和1 130 nm附近均出现些异常特征, 6S模型反演的异常特征更为明显。总体上看, MODTRAN 4.0模型的粘土矿反演效果略优于6S模型的; 在植被的反射率反演方面, 两个模型的效果差别不大。
水体的反射率很低, 反射特征基本上都是由水中的其他物质造成的。如图3所示, 两种模型反演的水体反射率曲线整体上较为接近, 基本符合水体的反射特征。
为了反映两种方法的精度, 利用PPI纯净像元指数法分别从两种模型的结果图像上提取植被像元50个、粘土矿像元25个, 根据相同的规则检验光谱曲线的正确性、曲线整体形态的一致性, 每个特征吸收位置的偏移在10 nm以内。结合ENVI平台的匹配滤波方法计算各像元的匹配权重, 逐像元计算得出两种模型重建光谱的相似度分布(表3)。
| 表3 光谱重建像元相似度分布 Tab.3 Index of similarity of spectra after spectral reconstruction |
从表3可以看出, 基于MODTRAN 4.0模型的粘土矿光谱重建精度略优于6S模型的, 而6S模型在植被光谱重建精度上略优。
使用MODTRAN 4.0和6S模型对高光谱数据进行大气校正是一个复杂的过程, 尤其是需要逐波段进行计算和反演, 对于数据量巨大的高光谱数据来说, 是一项庞大的工程。因此, 基于高光谱数据的特点, 发展高效的、能够集成化处理的大气校正软件, 有利于高光谱遥感应用的发展。
FLAASH(Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes)是由美国光谱科技公司(Spectral Sciences Inc)和空气动力研究实验室(Air Force Research Laboratory)共同研制和开发的, 它基于MODTRAN 4.0辐射传输模型的代码, 能对400~2 500 nm波长范围内的遥感影像进行大气校正[11]。FLAASH模型考虑了邻近效应的计算, 基于高光谱数据的特点, 分别采用波段比值法和暗目标法的原理, 引入了水气去除模型和气溶胶去除模型[12]。针对高光谱数据, 设置参数(类似MODTRAN 4.0和6S)之后, 可以进行全谱段的集成运算, 得到大气校正结果。本文利用FLAASH模型对试验区Hyperion数据进行大气校正和光谱重建的结果如图4所示。
可以看出, FLAASH模型的光谱重建效果不亚于MODTRAN 4.0和6S模型的校正效果, 甚至更优, 而且, FLAASH模型只需要运行一次就可得到Hyperion数据的全谱段大气校正结果。
(1)对于缺少地面同步光谱测量的高光谱数据, 大气辐射传输模型是最为可行的大气校正方法。通过对研究区Hyperion数据分别基于MODTRAN 4.0和6S模型的光谱重建, 结合3大类地物像元光谱曲线的反演结果, 说明大气辐射传输模型法能够比较有效地完成高光谱数据的大气校正和光谱重建, 得到地物真实反射率, 为高光谱数据下一步的定量应用打下基础。
(2)MODTRAN 4.0和6S模型都是优秀且成熟的大气校正模型, 从对高光谱Hyperion数据的应用效果来看, 对不同地物类型的处理效果各有千秋, 在实际应用中可适当结合使用, 以达到更好的效果。由于高光谱数据波段很多, 处理效率也是重要的因素之一, FLAASH模型等针对高光谱数据改进的大气校正处理软件能够在保障处理精度的前提下, 大大提高处理效率, 代表了高光谱数据大气校正处理的发展趋势。
The authors have declared that no competing interests exist.
| [1] |
|
| [2] |
|
| [3] |
|
| [4] |
|
| [5] |
|
| [6] |
|
| [7] |
|
| [8] |
|
| [9] |
|
| [10] |
|
| [11] |
|
| [12] |
|
