二类水体遥感数据常见大气校正方法.doc

二类水体遥感数据常见大气校正方法 陈克海 标准方法是在近红外波段的离水辐射率为零的假设下，基于两个近红外波的大气顶层反射率之比估算气溶胶类型和光学厚度，并且外推到可见光部分完成大气校正。该算法的校正精度主要取决于实际海面近红外波段的离水辐射是否可忽略。对于一类水体，由于其气溶胶一般属于非吸收性或者弱吸收性气溶胶，而且其在近红外波段的离水辐射甚微，因而标准校正方法在大部分一类水体处理精度可以达到业务化需求，其校正后得到的离水辐射率误差在5%以内（McClain，1998）。 然而，如果一类水体中存在较高浓度的叶绿素，或者近岸海水中存在较高浓度的泥沙，而高浓度的泥沙和叶绿素即使在近红外波段也具有较高的反射率，从而其在近红外波段的离水辐射率不能忽略。对于二类水体，假设在近红外波段后向散射系数不随波长而变化（Zhang et al，1999），而吸收系数主要取决于水的吸收系数，离水辐射对大气表观反射率的增量近似与成正比，因为，所以，从而标准方法会过高估算，当高过预设阀值，标准方法会报错。的过高估算，会过高估算气溶胶在可见光部分的反射率，波长越短，估算越偏高，所以在蓝光波段（412 nm）校正得到的离水辐射率会出现负值。所以在二类水体中，大气校正标准方法会失效。 另外，大气标准方法的校正精度还取决于实际气溶胶是否在备选气溶胶模型以内。备选气溶胶类型为对流层型气溶胶、海洋型气溶胶和沿岸带气溶胶，都属于非吸收或者弱吸收气溶胶类型。在一类水体，备选气溶胶类型一般能模拟绝大部分气溶胶状态，但对于城市气溶胶类型或者沙尘等吸收性气溶胶，标准方法校正效果较差。近岸水体一般较混浊，而且受到城市气溶胶等吸收性气溶胶的影响，必须对标准方法进行修正或者采用其它思路才能进行大气校正。 到目前为止，科学家提出了多种解决方案，常见的主要有最佳邻近法校正方法、迭代校正方法、光谱匹配校正方法和光谱优化方法、两个短波红外波段组合的校正方法。这些方法在二类水体区域大气校正方面都具有一定的优点和缺点，以下分别介绍。 1 最佳邻近法校正方法 Hu et al（2000） 针对SeaWiFS，提出了一种利用邻近清洁水体上空气溶胶类型对混浊水体进行大气校正的实用方法，能够较好得分离气溶胶和水体对卫星接收信号的影响。该方法基于气溶胶类型中尺度空间范围内（100－1000km）不变的假设（Gordon and Morel，1983），这样，混浊水体的气溶胶类型可以用邻近相对清洁海域的气溶胶类型代替。至于气溶胶浓度可以随像元而发生变化。为计算近红外波段的气溶胶光学厚度，作者进一步假设近红外波段下水体后向散射系数不随波长而变化，则近红外波段离水辐射对大气总反射率的贡献可近似表示为（Carder et al，1999）： 其中为765nm或865nm；为未知常数，相同像元的不随波长而变，但不同像元的值可不同；为水体吸收系数，可近似为水的吸收系数；是大气漫散射透过率，通过最佳邻近法从相邻像元获取。所以，对每个混浊像元的近红外波段，有如下关系 其中为在假设下卫星测量得到的气溶胶单次散射反射率（Gordon and Wang，1994）；为去除离水辐射影响下的气溶胶单次散射反射率。如果该混浊水体气溶胶类型可以从邻近像元获取，则在相同几何观测条件下，气溶胶类型能唯一确定，所以有 根据以上两个式子，可以得到未知常数的估算值 有了气溶胶类型和，就可以计算和相应的光学厚度。 Hu是在集成标准大气校正方法的SeaWiFS软件的基础上整合实现上面方法，首先让标准方法运行一遍，然而根据标准方法计算结果中的标签选定混浊像元。对每一混浊像元，采用最佳邻近法从邻近清洁水体中获取气溶胶类型，并通过上面的计算获取气溶胶光学厚度。最后，利用混浊像元的气溶胶优化参数，接着运行标准方法。Hu方法在墨西哥湾测试并证明它能减小离水辐射率和叶绿素浓度估算误差，并且对于标准方法运行失败的区域，它能较成功的运行。 最佳邻近法利用邻近清洁水体上空气溶胶类别校正混浊水体的大气效应，算法成功运行要求影像中存在清洁像元个数不能太少，否则会影响到校正效果。 2 迭代校正方法 Land and Haigh（1996）对二类水体遥感数据提出一种迭代校正方法，该方法在大气校正的同时能够反演水下物质浓度。该算法需要一定的气溶胶模型和水下光学模型。其中，气溶胶模型为大陆气溶胶、海洋气溶胶和城市气溶胶按照不同比例混合而成。每种气溶胶成分都具特定的单次散射反照率、衰减系数、散射相函数和光