[發明專利]一種面向IDL語言的海冰密集度變化趨勢求算算法

在審

• 申請號: 201810207828.2
• 申請日: 2018-03-14
• 公開（公告）號: CN108509527A
• 公開（公告）日: 2018-09-07
• 發明（設計）人: 王云鶴;畢海波;李運瑋;黃海軍
• 申請（專利權）人: 中國科學院海洋研究所
• 主分類號: G06F17/30
• 分類號: G06F17/30
• 代理公司: 沈陽科苑專利商標代理有限公司 21002
• 代理人: 許宗富

1.一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，包括以下步驟：

1)讀取北極海冰密集度月均數據；

2)計算當前網格的海冰密集度的當年年均數據；所述計算過程中對北極海冰密集度月均數據的有效值數量進行判斷，對當年月均數據的有效值數不過半的網格單元的年均數據設置為無效值；循環行列依次獲取每個網格的海冰密集度在N年中的當年年均數據；

3)對當前網格的海冰密集度在N年中的每個年均數據進行篩選，得出有效年均值數列，將有效年均值數不過半數的網格單元設置為無效值；對于當前網格的有效年均值數過半數的網格單元，將這N年中的有效年均值數列與對應的年份數列進行線性擬合，得到海冰密集度的變化率數據及此變化率的顯著性數據；依次計算所有網格單元處的海冰密集度的變化率數據及顯著性數據；

4)對海冰密集度變化率數據及顯著性數據中的衛星盲區數據進行掩膜處理，獲取掩膜后的海冰密集度變化率數據及顯著性數據；

5)寫出掩膜后的海冰密集度變化率數據及顯著性數據文件；

6)導入ArcGIS軟件進行配準，用不同顏色表示海冰密集度變化率數據、用等值線表示顯著性數據，疊加生成海冰密集度變化趨勢圖。

2.按照權利要求1所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，所述北極海冰密集度月均數據包括12×N個月均數據表：每個月均數據表包含當年該月的按照網格劃分的北極海冰密集度月均數據M_i,j，N為統計的年份個數，i為數據表的列數，j為數據表的行數。

3.按照權利要求1所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，所述步驟2)中的計算當前網格的海冰密集度的當年年均數據；所述計算過程中對北極海冰密集度月均數據的有效值數量進行判斷，對當年月均數據的有效值數不過半的網格單元的年均數據設置為無效值，包括：

對當前網格按年份判斷北極海冰密集度月均數據M_i,j的數值是否有效，對有效值數不過半的年份的網格單元的年均數據設置為無效值；根據有效值過半的年份的網格單元的有效月均數據M_i,j計算年均數據Y_i,j，得到N個按照網格劃分的北極海冰密集度年均數據表。

4.按照權利要求4所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，所述根據有效值過半的年份的網格單元的有效月均數據M_i,j計算年均數據Y_i,j為：年均數據

5.按照權利要求1所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，所述步驟3)中對于當前網格的有效年均值過半數的網格單元，將這N年中的有效年均值數列與對應的年份數列進行線性擬合，得到海冰密集度的變化率數據以及此變化率的顯著性，包括：

對海冰密集度有效年均數據和其對應的年份數列進行線性擬合處理，得到斜率和P值；所述斜率為海冰密集度的變化率數據，所述P值為此網格單元處海冰密集度變化率的顯著性數據。

6.按照權利要求1所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于，所述步驟4)中的掩膜處理為：利用衛星最大盲區制作掩膜，對步驟3)得到的海冰密集度變化趨勢數據及其顯著性數據，濾除盲區內的海冰密集度變化趨勢數據及其顯著性數據。

7.按照權利要求1所述的一種面向IDL語言的海冰密集度變化趨勢求算算法，其特征在于以IDL程序語言編寫實現。