[发明专利]一种可见红光通道地表反射率估算方法和装置

说明书

本发明提供了一种可见红光通道地表反射率估算方法和装置，该方法包括：获取待处理遥感图像；基于植被指数确定待处理遥感图像中的亮像元和暗像元；针对亮像元，构建红光波段地表反射率与1.64μm波段表观反射率的二次函数模型，根据二次函数模型计算得到红光波段地表反射率；针对暗像元，构建红光波段地表反射率与1.64μm波段表观反射率的一次函数模型，根据一次函数模型计算得到红光波段地表反射率。本方案提供的可见红光通道地表反射率估算方法同时实现了可见红光通道下暗地表和亮地表的地表反射率的准确估算。

技术领域

本发明涉及遥感技术领域，特别涉及一种可见红光通道地表反射率估算方法和装置。

背景技术

在大气参数反演、地物判别等众多遥感观测研究中，地表反射率是基础参数，其估算精度直接影响着其他遥感参数的测量精度。可见红光通道容易受到大气和气溶胶的影响，因此较难估算该通道天基下视地表反射率。

现有遥感地表反射率估算方法主要包括两种经典算法：深蓝算法和暗目标法。其中深蓝算法假设地表在一定时间段内是稳定的，因此能直接调用基于统计的地表反射率数据库，但引入外部地表反射率数据库会造成不可避免的几何匹配误差。经典的暗目标法具有模型简单、普适性和时效性强、准确度高的优点，但其不适用于干旱、半干旱等亮地表。因此现有的算法均不能达到满足应用需求的效果，阻碍了遥感数据的应用。因此，需提出一种可同时实现可见红光通道下暗地表和亮地表反射率的准确估算模型。

发明内容

本发明提供了一种可见红光通道地表反射率估算方法和装置，该方法能同时实现可见红光通道下暗地表和亮地表的地表反射率的准确估算。

第一方面，本发明实施例提供了一种可见红光通道地表反射率估算方法，包括：

获取待处理遥感图像；

基于植被指数确定所述待处理遥感图像中的亮像元和暗像元；

针对所述亮像元，构建红光波段地表反射率与1.64μm波段表观反射率的二次函数模型，根据所述二次函数模型计算得到红光波段地表反射率；

针对所述暗像元，构建红光波段地表反射率与1.64μm波段表观反射率的一次函数模型，根据所述一次函数模型计算得到红光波段地表反射率。

可选地，所述基于植被指数确定所述待处理遥感图像中的亮像元和暗像元，包括：

获取包括各通道表观反射率的历史遥感图像；

对所述历史遥感图像进行计算分析，得到对应各像元的归一化植被指数；

根据所述归一化植被指数，计算得到在1.64μm波段处的去气溶胶污染植被指数；

针对每个像元，均执行：

判断该像元的去气溶胶污染植被指数是否位于预设去气溶胶污染植被指数范围内；

若是，则确定该像元为亮像元；

若否，则确定该像元为暗像元。

可选地，所述预设去气溶胶污染植被指数范围的取值区间为(0，0.25]。

可选地，所述针对所述亮像元，构建红光波段地表反射率与1.64μm波段表观反射率的二次函数模型，根据所述二次函数模型计算得到红光波段地表反射率，包括：

所述二次函数模型为：

其中，R_Red用于表征所述红光波段地表反射率，R_1.64μm用于表征所述1.64μm波段表观反射率，C用于表征二次拟合系数，D用于表征一次拟合系数，E用于表征常数项拟合系数；其中，C、D和E均与在1.64μm波段处的去气溶胶污染植被指数相关。