[發明專利]一種基于ICESat-2數據的海冰密集度估算方法及系統

有效

• 申請號: 202110631481.6
• 申請日: 2021-06-07
• 公開（公告）號: CN113360849B
• 公開（公告）日: 2022-03-18
• 發明（設計）人: 劉俊;謝歡;郭亞磊;童小華;李培楠
• 申請（專利權）人: 上海工程技術大學;同濟大學
• 主分類號: G06F17/18
• 分類號: G06F17/18;G06F16/29;G06F16/28;G06F16/2457
• 代理公司: 上海科盛知識產權代理有限公司 31225
• 代理人: 宣慧蘭

1.一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，包括以下步驟：

S1、讀取ICESat-2的ATL10原始數據集，得到海冰數據集，所述海冰數據集中包括多個海冰數據，海冰數據的屬性包括經緯度信息、采集時間、冰水類別參數和段長；

S2、根據海冰數據的采集時間劃分海冰數據集中的海冰數據，得到不同月份的海冰數據集；對于每個月份的海冰數據集，根據海冰數據的經緯度信息，得到各個海冰數據在極地立體投影網格中的格網信息，包括海冰數據所在格網的行列號，根據海冰數據所在格網的行列號劃分海冰數據集中的海冰數據，得到不同格網的海冰數據集；

S3、分別計算極地立體投影網格中各個格網的海冰密集度，得到不同月份各個格網的海冰密集度，計算過程具體為：對于極地立體投影網格中的一個格網，獲取該格網的海冰數據集，計算該格網的海冰密集度：

其中，N表示該格網內ICESat-2的測高軌跡的條數，P_i表示該格網內ICESat-2第i條測高軌跡包含的采樣點的數量，每個采樣點對應海冰數據集中的一個海冰數據，I_track表示沿該格網內ICESat-2測高軌跡的海冰密集度，L_ij表示海冰數據的段長，表示海冰數據的緯度，W_ij表示海冰數據的波形類型，W_ij的取值為0或1，W_ij的大小是根據海冰數據的冰水類別參數的值確定的。

2.根據權利要求1所述的一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，海冰數據的屬性還包括飛行方向參數，所述飛行方向參數用于表征ICESat-2所處的飛行模式，飛行方向參數等于1表示處于前進模式，右側right是強光束，飛行方向參數等于0表示處于后退模式，左側left是強光束，飛行方向參數等于2表示處于過渡模式；步驟S1中還包括基于飛行方向參數的值保留海冰數據集中對應強光束的海冰數據。

3.根據權利要求1所述的一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，冰水類別參數用于區分采樣點是海冰或冰間水道，當冰水類別參數等于0時，表示被云遮蓋的采樣情況，當冰水類別參數等于1時，表示采樣點為沒有冰間水道的冰/雪表面，當冰水類別參數等于2～9時，表示采樣點為冰間水道。

4.根據權利要求3所述的一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，步驟S1中還包括自海冰數據集中剔除冰水類別參數等于0的海冰數據。

5.根據權利要求3所述的一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，步驟S3中，如果一個海冰數據的冰水類別參數等于1，則該海冰數據的波形類型取值為1，如果一個海冰數據的冰水類別參數不等于1，則該海冰數據的波形類型取值為0。

6.根據權利要求1所述的一種基于ICESat-2數據的海冰密集度估算方法，其特征在于，步驟S3中使用的極地立體投影網格是NSIDC的25×25km極地立體網格。