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

摘要

本发明公开了一种消除水域遥感数据镜面反射影响的方法，其特征在于，包括以下步骤(A)从图像上提取暗像元，迭代计算大气散射系数、大气透过率、和天空光辐照度；(B)对整幅图像进行大气纠正；(C)消除天空光镜面反射光；(D)消除太阳直射光镜面反射光。该方法有效地消除了光学波段遥感数据中每个水域像元中的水面镜面反射光的影响，为精确的水质遥感、水下地或形地物遥感提供了必要的前提基础。

主权项

一种消除水域遥感数据镜面反射影响的方法，其特征在于，包括以下步骤：(A)从图像上提取山区植被阴影像元或清深水体像元，迭代计算大气散射系数、大气透过率和天空光辐照度；(B)对图像中所有像元进行大气纠正；(C)消除天空光镜面反射光；(D)消除太阳直射光镜面反射光；其中，所述步骤(A)中迭代计算大气散射系数、大气透过率和天空光辐照度的步骤包括：(A1)读取山区植被阴影像元或清深水体像元的表观反射率Rvd，根据山区植被阴影像元或清深水体像元的表观反射率Rvd的计算公式：式中：ω为大气散射系数，T为大气透过率，Rv为地物反射率，P(θ)为大气散射相函数，θ为散射角，首先忽略等式(1)右边第一项的值，计算得到大气散射系数的初值(A2)：根据大气透过率的初值的计算公式：T′＝e‑ω′，计算出大气透过率的初值T′，将大气透过率的初值T′代入等式(1)右边第一项中的T，将大气散射系数的初值ω′代入等式(1)右边第一项中的ω，计算得到大气散射系数(A3)：计算天空光辐照度：将步骤(A2)中计算得到的大气散射系数ω1代入天空光辐照度计算公式计算天空光辐照度Esky1，式中：c为大气纠正前后的大气散射系数比，为测量值；(A4)：将步骤(A2)中计算得到的大气散射系数ω1代入大气透过率的计算公式T＝e‑ω，计算出大气透过率T1；(A5)：将步骤(A2)中计算得到的大气散射系数ω1代入等式(1)右边第一项中的ω、将步骤(A4)中计算得到的大气透过率T1代入等式(1)右边第一项中的T，计算出精度高于大气散射系数ω1的大气散射系数ω2；根据步骤(A3)计算出精度更高的天空光辐照度Esky2；根据步骤(A4)计算出精度高于大气透过率T1的大气透过率T2；依此办法，迭代计算出精度高于大气散射系数ω2的大气散射系数ωk、精度高于大气透过率T2的大气透过率Tk和精度高于天空光辐照度Esky2的天空光辐照度Eskyk；其中，k>2。