[发明专利]ICESat-2/ATLAS全球高程控制点提取方法及系统

权利要求书

1.一种ICESat-2/ATLAS全球高程控制点提取方法，其特征在于，包括以下步骤：

步骤1，从ICESat-2的ATL08数据产品中提取陆地高程数据，进行预处理，保留属性参数作为高程控制点的筛选条件，以减少数据量，然后基于参考DEM对高程数据进行粗筛选，剔除高程误差较为严重的激光点；

步骤2，结合ICESat-2激光数据的属性信息进行精细化筛选，保留地形平缓且云量低的激光点作为地面高程控制点，导入ICESat-2全球高程控制点数据库。

2.如权利要求1所述ICESat-2/ATLAS全球高程控制点提取方法，其特征在于：所述保留属性参数作为高程控制点的筛选条件，实现方式如下，

从ICESat-2的ATL08数据产品中，提取陆地激光点的参数，包括编号n、经度longitude、纬度latitude、插值陆地高程h_te_interp、均值陆地高程h_te_mean、中值陆地高程h_te_median、最大陆地高程h_te_max、最小陆地高程h_te_min、参考DEM高程dem_h、时间delta_time、坡度terrain_slope、地形标记terrain_flg、城市标记urban_flag、云量标记cloud_flag_atm和夜晚标记night_flag。

3.如权利要求2所述ICESat-2/ATLAS全球高程控制点提取方法，其特征在于：所述基于参考DEM对高程数据进行粗筛选，实现方式如下，

利用ATL08中提取的参数dem_h作为先验参考DEM，对ICESat-2陆地高程数据进行初步筛选，将高程差ΔH大于阈值的高程点作为粗差剔除如下，

ΔH＝|H_ICESat-2-H_{dem_h}|＜3σ

式中，H_ICESat-2表示ICESat-2高程点的插值陆地高程h_te_interp，H_{dem_h}表示参考DEM的高程值dem_h，σ表示参考DEM的绝对高程精度。

4.如权利要求1或2所述ICESat-2/ATLAS全球高程控制点提取方法，其特征在于：所述结合ICESat-2激光数据的属性信息进行精细化筛选，实现方式包括设置基于属性参数的精细化筛选条件如下，

terrain_slopeT_slope

cloud_flag_atmT_cloud

其中，T_slope为坡度阈值，T_cloud为云量阈值；

并支持选择设置night_flag＝1；

将筛选后的陆地高程点载入高程控制点数据库，同时保留相关属性标记，作为高程控制点数据检索时的权重依据。

5.一种ICESat-2/ATLAS全球高程控制点提取系统，其特征在于：用于实现如权利要求1-4任一项所述的一种ICESat-2/ATLAS全球高程控制点提取方法。

6.根据权利要求5所述ICESat-2/ATLAS全球高程控制点提取系统，其特征在于：包括以下模块，

第一模块，用于从ICESat-2的ATL08数据产品中提取陆地高程数据，进行预处理，保留属性参数作为高程控制点的筛选条件，以减少数据量，然后基于参考DEM对高程数据进行粗筛选，剔除高程误差较为严重的激光点；

第二模块，用于结合ICESat-2激光数据的属性信息进行精细化筛选，保留地形平缓且云量低的激光点作为地面高程控制点，导入ICESat-2全球高程控制点数据库。

7.根据权利要求5所述ICESat-2/ATLAS全球高程控制点提取系统，其特征在于：包括处理器和存储器，存储器用于存储程序指令，处理器用于调用存储器中的存储指令执行如权利要求1-4任一项所述的一种ICESat-2/ATLAS全球高程控制点提取方法。

8.根据权利要求5所述ICESat-2/ATLAS全球高程控制点提取系统，其特征在于：包括可读存储介质，所述可读存储介质上存储有计算机程序，所述计算机程序执行时，实现如权利要求1-4任一项所述的一种ICESat-2/ATLAS全球高程控制点提取方法。