[发明专利]一种结合时空信息的遥感逐日积雪产品去云方法

权利要求书

1.一种结合时空信息的遥感逐日积雪产品去云方法，其特征在于，包括以下步骤：

步骤1：将原始的遥感逐日积雪产品重分类为积雪、无雪和云三类，同时根据DEM数据和土地利用数据对遥感逐日积雪产品进行积雪单元划分；

步骤2：利用DEM数据和多时相数据分别计算每个云像元在空间和时间上的积雪概率；

步骤3：利用加权平均方法获得云像元的积雪概率，评判云像元是积雪还是无雪。

2.根据权利要求1所述的结合时空信息的遥感逐日积雪产品去云方法，其特征在于：步骤1中所述的对遥感逐日积雪产品进行积雪单元划分，划分后将云量大于70％的积雪单元进行合并，形成新的积雪单元。

3.根据权利要求1所述的结合时空信息的遥感逐日积雪产品去云方法，其特征在于：步骤2中所述的计算每个云像元在空间和时间上的积雪概率，具体实现过程是：

首选在每个积雪单元的高程中，找到雪线H_S和无雪线H_L；其中，高于雪线H_S的无云像元都为积雪，低于无雪线H_L的无云像元都为无雪；在H_S、H_L之间等间距划分高程，并且高程大于等于H_S和小于等于H_L各为一类；

然后，在每个高程带中，根据南、北和东西三类坡向继续划分高程；

在高程、坡向分类基础上，每个云像元在其所在类中的积雪概率为P_H，具体公式如下：

P_H＝N_s/N_nocloud             (1)；

其中，N_s是该类中积雪像元数目，N_nocloud是该类中非云像元数目；

P_T是当前云像元在时间上前后各5个非云像元中的积雪概率，公式如下：

PT=Σ1/dindexΣi=1101/di,index∈{1,2,3,...,10}---(2);]]>

其中，d_i是非云像元i到当前云像元的时相距离，d_index是非云像元中积雪像元到当前云像元的时相距离。

4.根据权利要求3所述的结合时空信息的遥感逐日积雪产品去云方法，其特征在于：步骤3中的具体实现过程是，根据下列公式得到云像元的总体积雪概率P_s，当P_s大于等于P₀时，将该云像元判断为积雪，反之为无雪：

P_s＝w₁×P_H+w₂×P_T         (3)；

其中，P₀为设定的阈值，w₁和w₂分别是P_H和P_T的权重，且w₁+w₂＝1。

5.根据权利要求4所述的结合时空信息的遥感逐日积雪产品去云方法，其特征在于：步骤3中所述的w₁和w₂，其计算方法是依据积雪单元中的非云像元的信息获取，具体步骤如下：

首先，在每个积雪单元中，利用步骤3中的方法求得非云像元在空间和时间上的积雪概率；

然后，在[0,1]区间内，以0.01为步长设定w₁，以公式(3)求得非云像元的积雪概率，判断其为积雪还是无雪，并将判断结果视为非云像元的预测值；

最后，将非云像元的原始值作为真实值，算出所述的遥感逐日积雪产品去云方法在非云像元上的精度oa，公式如下所示：

oa=Nss+NnsnsNnc---(4);]]>

其中，为积雪单元中非云像元的预测值和真实值都为积雪的数目，为积雪单元中非云像元的预测值和真实值都为无雪的数目，N_nc为积雪单元中非云像元的总数目。