[發明專利]一種結合時空信息的遙感逐日積雪產品去云方法

有效

• 申請號: 201510094932.1
• 申請日: 2015-03-03
• 公開（公告）號: CN104657952B
• 公開（公告）日: 2017-04-12
• 發明（設計）人: 沈煥鋒;付文軒;李星華;張良培
• 申請（專利權）人: 武漢大學
• 主分類號: G06T5/00
• 分類號: G06T5/00
• 代理公司: 武漢科皓知識產權代理事務所(特殊普通合伙)42222
• 代理人: 薛玲

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為積雪單元中非云像元的總數目。