[发明专利]一种面向IDL语言的海冰密集度变化趋势求算算法

说明书

本发明涉及一种面向IDL语言的海冰密集度变化趋势求算算法。读取北极海冰密集度数据。以单个数据网格为单位循环行列，求算每个数据网格的海冰密集度在N年的年均值。再以单个数据网格为单位循环行列，对每个网格的在N年的有效年均值进行判断筛选，对于有效年均值数过半数的网格单元，将N年的有效年均值数列与对应的年份数列进行线性拟合，得出海冰密集度的变化率以及此变化率的显著性。将卫星盲区的海冰密集度变化率和显著性数据进行掩膜。写出海冰密集度变化率和显著性数据文件。导入ArcGIS软件进行配准、生成海冰密集度变化趋势图。本发明不仅能反映某区域海冰密集度的变化趋势及空间分布，还能反映海冰密集度变化趋势的显著性。

技术领域

本发明涉及求算海冰密集度变化趋势及利用ArcGIS进行可视化显示，更具体地说是一种面向IDL语言的海冰密集度变化趋势求算算法。

背景技术

IDL程序语言，灵活简洁，功能强大；国际上，IDL语言在海洋科学、遥感地信、地质等领域应用广泛；目前，在我国还没有达到普及水平。卫星具有全方位、多时相特点，打破了现场观测的区域性，时间不同步的缺点，节省了大量的人力物力。卫星数据越来越受到人们的认可，数据所涉及领域也会愈加广泛。目前计算海冰密集度变化趋势的方法主要是先计算某区域整体的年均值，然后对年均值进行线性拟合，求出整个区域的海冰密集度变化趋势值，缺点在于这种结果只能反映整体区域的时间变化性，不能反映海冰密集度变化趋势的区域性。本算法得到的结果既能反映某区域海冰密集度的变化情况，同时也能反映海冰密集度变化趋势空间分布，便于对比不同区域的海冰密集度变化趋势的不同情况等，具有较强的理论价值和应用价值。

发明内容

针对上述现状，本发明的目的提供一种面向IDL语言的海冰密集度变化趋势求算算法。

本发明基于IDL编程语言，求出的海冰密集度变化趋势，不仅能反映某区域海冰密集度的变化情况，同时也能反映海冰密集度变化趋势空间分布，还能得到不同区域海冰密集度变化趋势的显著性。本发明关于海冰密集度变化趋势求算算法是基于IDL程序语言实现，算法简洁，便于操作。

本发明解决其技术问题所采用的技术方案是：一种面向IDL语言的海冰密集度变化趋势求算算法，包括以下步骤：

1)读取北极海冰密集度月均数据；

2)计算当前网格的海冰密集度的当年年均数据；所述计算过程中对北极海冰密集度月均数据的有效值数量进行判断，对当年月均数据的有效值数不过半的网格单元的年均数据设置为无效值；循环行列依次获取每个网格的海冰密集度在N年中的当年年均数据；

3)对当前网格的海冰密集度在N年中的每个年均数据进行筛选，得出有效年均值数列，将有效年均值数不过半数的网格单元设置为无效值；对于当前网格的有效年均值数过半数的网格单元，将这N年中的有效年均值数列与对应的年份数列进行线性拟合，得到海冰密集度的变化率数据及此变化率的显著性数据；依次计算所有网格单元处的海冰密集度的变化率数据及显著性数据；

4)对海冰密集度变化率数据及显著性数据中的卫星盲区数据进行掩膜处理，获取掩膜后的海冰密集度变化率数据及显著性数据；

5)写出掩膜后的海冰密集度变化率数据及显著性数据文件；

6)导入ArcGIS软件进行配准，用不同颜色表示海冰密集度变化率数据、用等值线表示显著性数据，叠加生成海冰密集度变化趋势图。

所述北极海冰密集度月均数据包括12×N个月均数据表：每个月均数据表包含当年该月的按照网格划分的北极海冰密集度月均数据M_i,j，N为统计的年份个数，i为数据表的列数，j为数据表的行数。

所述步骤2)中的计算当前网格的海冰密集度的当年年均数据；所述计算过程中对北极海冰密集度月均数据的有效值数量进行判断，对当年月均数据的有效值数不过半的网格单元的年均数据设置为无效值，包括：