[发明专利]一种基于ICESat-2数据的海冰密集度估算方法及系统

权利要求书

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极地立体网格。