[发明专利]一种消除水域遥感数据镜面反射影响的方法

说明书

技术领域

本发明属于水质遥感数据处理方法技术领域，特别涉及一种消除水域遥感数据镜面反射影响的方法。

背景技术

平静的水面可以被近似地看作一面镜子，当太阳光入射时，在与入射方向关于法线对称的反射方向上形成强烈的反射辐射，形成太阳耀光。当水面存在波浪时，水面形成倾角，在海浪的顶端形成闪光点，形成更大范围的太阳耀光区。太阳耀光具有很大的辐射强度，当遥感器直接观测太阳耀光区时，容易出现饱和现象，也即水域遥感数据的镜面反射。

水面镜面反射光基本不含水质和水下地物信息。由于水体反射率低，这部分镜面反射光的信号强度常高于其它水质信息，成为水质遥感或水下信息提取的强的干扰因素。

目前，国内外并无有效的消除水面镜面反射的方法，主要采取“回避”策略，所谓“回避”，就是在传感器探测时尽可能回避存在严重太阳耀光的区域。也即在进行对地观测时，通过设计合理的飞行路线和选择合适的时间，避开水面镜面反射方向，也即强太阳耀斑方向。

但是，由于波浪的作用，水面并非为一平面，镜面反射光在各方向均有分布，无论从哪一方向进行观测，无法完全避开。并且，进入传感器镜头的镜面反射光与水面粗糙度有关，每个像元均不相同，即与水面风浪大小正成相关，在开阔海域尤为严重。因此，这种方法并不能有效消除镜面反射的影响。

发明内容

本发明的目的是提供一种消除水域遥感数据镜面反射影响的方法，该方法有效地消除了光学波段遥感数据中每个水域像元中的水面镜面反射光的影响，为精确的水质遥感、水下地或形地物遥感提供了必要的前提基础。

为了达到上述目的，本发明采用如下技术方案：

一种消除水域遥感数据镜面反射影响的方法，包括以下步骤：

(A)从图像上提取山区植被阴影像元或清深水体像元，迭代计算大气散射系数、大气透过率和天空光辐照度；

(B)对图像中所有像元进行大气纠正；

(C)消除天空光镜面反射光；

(D)消除太阳直射光镜面反射光。

作为一种具体的实施例，步骤(A)中迭代计算大气散射系数、大气透过率和天空光辐照度的步骤包括：

(A1)读取山区植被阴影像元或清深水体像元的表观反射率R_vd，根据山区植被阴影像元或清深水体像元的表观反射率R_vd的计算公式：

式中：ω为大气散射系数，T为大气透过率，R_v为地物反射率，P(θ)为大气散射相函数，θ为散射角，

首先忽略等式(1)右边第一项的值，计算得到大气散射系数的初值

(A2)：根据大气透过率的初值的计算公式：T′＝e^-ω′，计算出大气透过率的初值T′，并将大气透过率的初值T′代入等式(1)右边第一项中的T，将大气散射系数的初值ω′代入等式(1)右边第一项中的ω，计算得到大气散射系数

(A3)：计算天空光辐照度：

将步骤(A2)中计算得到的大气散射系数ω₁代入天空光辐照度计算公式计算天空光辐照度E_sky1，式中：c为大气纠正前后的大气散射系数比，为小于2的测量值；

(A4)：将步骤(A2)中计算得到的大气散射系数ω₁代入大气透过率的计算公式T＝e^-ω，计算出大气透过率T₁；

(A5)：将步骤(A2)中计算得到的大气散射系数ω₁代入等式(1)右边第一项中的ω、将步骤(A4)中计算得到的大气透过率T₁代入等式(1)右边第一项中的T，计算出精度高于大气散射系数ω₁的大气散射系数ω₂；

根据步骤(A3)计算出精度更高的天空光辐照度E_sky2；

根据步骤(A4)计算出精度高于大气透过率T₁的大气透过率T₂；

依此办法，迭代计算出精度高于大气散射系数ω₂的大气散射系数ω_k、精度高于大气透过率T₂的大气透过率T_k和精度高于天空光辐照度E_sky2的天空光辐照度E_skyk；其中，k>2。

进一步地，所述步骤(B)的具体步骤为，对图像中所有像元，均进行以下计算：