利用MODIS图像反演海岸与海岛的地物光谱反射率

引言

定量遥感初期，人们试图通过植被指数NDVI方法———归一化NDVI和时间序列NDVI最大值合成等方法绕开大气校正问题。目前，大气校正主要采用依据大气参数进行辐射传输计算的方法，如在陆地上使用暗目标方法反演气溶胶光学特性，还有暗体减法等；这些校正方法的目的就是为了减小地面反射率不确定度产生的误差。

由于海洋表面相对陆地表面来说要简单得多，而且在对水体观测时，大气信号所占比重更可达总信号的90%以上，故海洋表面的反射率变化对大气光学特性的影响较小。本文基于这一性质，首先借助AHMAD辐射传输模型，用MODIS图像基于查找表的方法反演出海洋上空的气溶胶的光学特性。若在所选影像为晴空无云条件下，假设一定范围内（中心点15km为半径的圆内）的海岛与海岸上空的大气和水体上空的大气一样，然后借助6S辐射传输模型和反演出的气溶胶模式，计算基于地物光谱反射率的查找表，再由MODIS图像的陆地像元和反演出的气溶胶光学厚度，最后用插值法即可得到地物光谱反射率。

反演方法

2.1 利用MODIS图像反演气溶胶光学特性

用MODIS图像反演气溶胶光学特性，分为陆地与海洋两种反演方法。陆地上空气溶胶的反演是基于暗目标的方法；海洋上空气溶胶的反演是基于查找表的方法。

2.1.1 查找表的计算

有研究人员提出利用MODIS图像，基于查找表的方法来反演海洋上空气溶胶的光学性质。查找表是基于大小粒子模式，用AHMAD辐射传输模型计算大气顶层表观反射率的表。本文采用气溶胶模式，即5种大粒子和4种小粒子模式。海洋上空气溶胶光学特性的反演误差，主要来源于气溶胶粒子的大小与谱分布，因而气溶胶粒子的模型设定有重要的意义。近海的海水辐射变化较大，故海洋表面反射模型采用Cox的海浪模型。把9个粒子模式分别当作9种气溶胶模式，对每一种气溶胶模式，用AHMAD辐射传输模型分别计算在550nm处，6个气溶胶光学厚度（0.0，0.2，0.5，1.0，2.0，3.0），16个观察天顶角（1.5°~88.5°，间隔6°）､16个相对方位角（0°~180°，间隔12°）和9个太阳天顶角（1.5°，12°，24°，36°，48°，54°，60°，66°和72°）在MODIS前7个波段下的表观反射率的值。

2.1.2反演方法和原理

计算生成查找表的方法和利用它们进行反演的基础是按双对数正态模式分布的粒子多次散射，由同样光学厚度下每一个模式粒子产生的辐射加权平均近似计算。当两个模式具有不同的吸收特性时，简化就会遇到困难，但对这里的简化概念是很合适的。这种简化的优点是仅需考虑9种气溶胶模式，而不是所有4×5×11组合（11是表示模式间的相对浓度有11种可能值，即η的所有可能的取值类）。如果卫星观测的总反射率可表示为:

对于海洋上空气溶胶的反演，MODIS的前7个通道能用，但第三通道（470nm）对海洋表面的变化有很大的起伏，而且它会导致在确定还不完全清楚的气溶胶贡献时产生误差，故在反演中没有使用此通道。对于任何η值上的大粒子和小粒子模式的合成，利用式（1）计算550nm通道总的表观反射率，利用550nm通道观测的表观反射率通过在这6个光学厚度下所有大粒子和小粒子模式的组合进行线性内插获得光学厚度，然后在相应的气溶胶模式下用此光学厚度计算其他5个通道的表观反射率。

2.2 地物光谱反射率的反演

2.2.1 6S辐射传输

6S是Vermote等人提出的解辐射传输方程的模型，只要给定几何条件､大气模式､气溶胶模式､光学厚度以及地面反射率情况，利用6S就可以计算大气顶层的表观反射率和大气的其他性质，如大气透过率等。

2.2.2利用6S计算查找表

由6S计算所需的参数，大气模式选择中纬度夏季（根据反演的图像确定）；气溶胶模式用§1.1反演的最佳结果；海拔高度为零；下垫面取朗伯面（反射率分别取0.02到0.18每隔0.02的所有值）。借助6S辐射传输模型，计算在550nm处，6个气溶胶光学厚度（0.0，0.2，0.5，1.0，2.0，3.0）､6个太阳天顶角（18.2°~21.2°，间隔0.5°）､6个相对方位角（0°~30°，间隔5°）和5个观察天顶角（15°，20°，25°，30°，35°）下的在MODIS7个波段处的大气顶层表观反射率，用于反演地物光谱反射率的查找表。查找表在670nm处的结构如表1所示（部分）。

表1计算地物光谱反射率的查找表

2.2.3地物光谱反射率的反演

由MODIS图像，可以读出像元的几何条件即太阳天顶角､观测天顶角､相对观测方位角和表观反射率，再加上§1.1反演出的最佳气溶胶模式和光学厚度，基于表1，利用线性插值法可以求出地物光谱反射率。

2.3 反演的流程

基于查找表反演地物光谱反射率的方法，是建立在海洋表面反射率已知的基础上以及海洋上空的气溶胶模式和海岸上空气溶胶模式相同的条件下。当条件满足时，就可以从MODIS图像反演出地物光谱反射率，其步骤如下。

1）假设9种粒子模式，其中4种小粒子和5种大粒子；

2）用Cox的海浪模型计算海洋表面反射率；

3）借助AHMAD辐射传输模型计算用于反演气溶胶光学特性的查找表A；

4）选出晴空的MODIS图，读出海洋部分像元的反射率值，反演出大气气溶胶的模式和光学厚度；

5）由反演出的气溶胶模式和6S辐射传输模型计算用于反演地物光谱反射率的查找表B；

6）从MODIS图中读出海岸和海岛部分像元的反射率值和几何条件，基于查找表B和步骤4）反演的气溶胶光学厚度，利用插值法反演得到地物光谱反射率。

2.4 反演误差的模拟估测

为了验证气溶胶光学厚度的改变对地物光谱反射率的影响，进行如下模拟计算。取太阳天顶角与方位角分别为19°和95°，观测天顶角与方位角分别为28°和90°，气溶胶光学厚度在550nm处的真值为0.28，对查找表B进行插值。

表2模拟了550nm波段气溶胶光学厚度增大10%或减小10%等情况下，对反演地面反射率的影响。从表2中可以看出，随着气溶胶光学厚度误差的加大，对可见光波段的反射率反演影响加大，而对近红外波段的影响相对较小，气溶胶光学厚度变化10%对反演地面各波段光谱反射率产生的误差均在5%以内。如果再考虑插值误差，则总的误差不会超过10%。

表2气溶胶光学厚度的不确定度给地物

通常情况下，陆地上空气溶胶光学厚度总是大于海洋上空的气溶胶光学厚度。本算法中对气溶胶光学厚度的反演误差可能来源于水体､风速和风向的选取，它们带来的误差使气溶胶光学厚度计算值产生偏差；而在一定范围内，认为陆地上空与海洋上空大气条件相同的假设，有可能造成气溶胶光学厚度比真实值小。故气溶胶光学厚度在计算中可能出现偏小，而产生的误差是地面光谱反射率的反演误差的主要来源。

MODIS图像的应用与分析

选用MODIS图像是中国东南海域2003-05-28卫图，所取范围为晴空无云的30像元×30像元（已经归化1km分辨率图，图1（a）圈内的部分）。根据暗目标法，提取图1（b）中的所有水像元，用于反演气溶胶的光学厚度和气溶胶的模式。反演出的气溶胶最佳模式是第三小粒子与第一大粒子按4∶1比例，以双模谱的形式组成；气溶胶光学厚度每行的平均值如图2（a）所示。将每行海水像元计算的气溶胶光学厚度取平均用于该行陆地像元反射率的反演，被海岛隔开的取两个值。

图1MODIS图

图2MODIS图像中每行气溶胶的结果

图2（a）给出了图像中每行气溶胶平均光学厚度（有海岛的分为两个），图2（b）给出气溶胶每行平均光学厚度的相对误差，标准值选自该地区的地面对大气的同步测量。从图2中可以看出该处的大气气溶胶随空间变化很小，所研究范围内的大气均匀状况良好。

图3和图4给出了图1（b）中任意标注的各像元反射率曲线。图3中曲线是海水的光谱曲线，在870nm到近红外波段的反射率近似为零。T1和T2都是把海平面看作朗伯面采用同陆地一样的方法反演出来的反射率，T是用海浪模型计算出来用于反演海洋上空气溶胶光学特性的反射率（风速用6m/s，叶绿素含量用0.5mg/m3，盐份35ppt）。从这三条曲线可以看，出反演的结果与模型计算的结果符合得较好。图4中4个图的曲线是近海岸地物反射率的曲线。这些光谱曲线较为符合植被的光谱特点，在可见光波段反射很小，在670nm~1240nm反射迅速增加，红外波段渐渐减小.

图3图2（a）中所示海洋像元在MODIS7个波段的反射率曲线

图4图1（b）中所示各像元在MODIS7个波段的反射率曲线图

图4（a）中条曲线是大陆沿岸的地物光谱反射曲线，曲线B､F可能含有水像元；图4（b）8条曲线是海岛沿岸地物光谱曲线；曲线G可能是海洋与陆地的复合像元；曲线O是典型的绿色植被反射率曲线；图4（c）､4（d）图中的曲线分别是离海岸较远的内陆与海岛内的地物光谱反射率，从中可以看出内陆的地物较复杂，海岛内的地物较为单一。

4、结论与展望

根据海岸地物近水的特点，本文采用海洋上空反演的气溶胶光学厚度用于地面反射率的反演，对每一个地面像元由卫星接收时的太阳天顶角､方位角，观测天顶角､方位角和反演的大气气溶胶光学厚度，在6S辐射传输模型基于反演的气溶胶模式计算的查找表中，用插值法即可反演出地物光谱反射率。数值模拟显示，在反演地物光谱反射率时气溶胶光学厚度的误差在10%以内，则对地面光谱反射率的反演误差均在10%以内。该方法利用了实时反演的气溶胶模式，反演速度快，反演误差小。但此法仅仅适用于近海岸和小的海岛，对离海洋较远（30km）的陆地反演的结果有待近一步研究。近海岸的水体是Ⅱ类水体，目前对II类水体没有好模型计算其反射率，它对气溶胶反演的精度有影响，从而影响地面反射率的精度。下一步工作的主要任务是改进Ⅱ类水体的反射率模型，并用实测卫星过顶时地物反射率数值与反演结果作对比分析。

欢迎关注公众号：莱森光学，了解更多光谱知识。

莱森光学（深圳）有限公司是一家提供光机电一体化集成解决方案的高科技公司，我们专注于光谱传感和光电应用系统的研发、生产和销售。

审核编辑黄宇