适用于海岛的ICESat-2高程控制点提取方法_戴泽源.pdf

引用格式:戴泽源,张立华,张 林,等.适用于海岛的ICESat-2高程控制点提取方法J.地球信息科学学报,2023,25(8):1559-1569.Dai Z Y,Zhang L H,Zhang L,et al.A method of island elevation control point extraction utilizing ICESat-2 dataJ.Journal of Geo-informa-tion Science,2023,25(8):1559-1569 DOI:10.12082/dqxxkx.2023.220868适用于海岛的ICESat-2高程控制点提取方法戴泽源1,2，张立华1,2*，张林3，刘翔1,2，周寅飞1,2，陈秋41.海军大连舰艇学院军事海洋与测绘系，大连 116018；2.海军大连舰艇学院海洋测绘工程军队重点实验室，大连 116018；3.海军大连舰艇学院导弹与舰炮系，大连 116018；4.91001部队，北京 100071A Method of Island Elevation Control Point Extraction Utilizing ICESat-2 DataDAI Zeyuan1,2,ZHANG Lihua1,2*,ZHANG Lin3,LIU Xiang1,2,ZHOU Yinfei1,2,CHEN Qiu41.Department of Military Oceanography and Hydrography&Cartography,Dalian Naval Academy,Dalian 116018,China;2.KeyLaboratory of Hydrographic Surveying and Mapping of PLA,Dalian Naval Academy,Dalian 116018,China;3.Department ofMissiles and Shipboard Gunnery,Dalian NavalAcademy,Dalian 116018,China;4.91001 Troops,Beijing 100071,ChinaAbstract:To address the challenge of island elevation control information acquisition,a method of islandelevation control point extraction based on ICESat-2 ATL08 data is proposed.Firstly,the sea area points and thegross error points in the data are deleted via the water area mask and the reference data,which reduces thevolume of data waiting to be processed.Then,considering the impact of the atmospheric environment,signal-to-noise ratio,and data integrity on the accuracy of elevation control points,the data are preliminary screened withmultiparameter restrictions.Finally,an adaptive terrain threshold is determined to extract the elevation controlpoints based on the precision index.The accuracy of the extracted elevation control points are validated usinghigh-precision Digital Elevation Model(DEM)data from various islands in Puerto Ricos outlying islands andthe United States Virgin Islands.The experimental results show that,compared to the original data,the MeanAbsolute Error(MAE)and Root Mean Square Error(RMSE)of the elevation control points extracted by theproposed method in the flat,hilly,and mountainous regions are reduced from 2.65 m/7.23 m,3.92 m/7.65 m,and4.93 m/8.29 m to 0.28 m/0.61 m,0.46 m/0.79 m,and 0.63 m/0.89m,respectively.On the premise of ensuringaccuracy,the number of elevation control points extracted by the proposed method is significantly larger thanthat of the existing methods.Our proposed method allows for extraction of a specific number of reliableelevation control points from the island region,which provides a solid foundation for the global nautical chartrenovation and quality control.Key words:island;ICESat-2;ATL08;elevation control point;adaptive threshold;spaceborne LiDAR;photoncounting;nautical chart renovation*Corresponding author:ZHANG Lihua,Email:摘要：针对海岛控制信息获取困难这一现实问题，提出一种基于ICESat-2 ATL08数据的海岛高程控制点提取方法。首先，通Vol.25,No.8Aug.,2023第25卷 第8期2023年8月收稿日期：2022-11-07；修回日期：2022-12-11.基金项目：国家自然科学基金项目（41871369、41871295、41901320、42071439）。Foundation items:National Natural ScienceFoundation of China,No.41871369,41871295,41901320,42071439.作者简介：戴泽源（1994），男，福建长汀人，博士生，主要从事海上作战环境保障研究。E-mail:carl-*通讯作者：张立华（1973），男，湖南桃江人，博士，教授，主要从事海上作战环境保障研究。E-mail:地 球 信 息 科 学 学 报2023年过水域掩膜标记与参考高程数据，去除数据中的海域点与粗差点，降低后续处理的数据量；然后，综合分析大气环境、信噪比、数据完整性等因素对高程值的影响，构建多参数约束初步筛选数据；最后，充分利用数据中的已有参数信息，结合控制点精度指标设计自适应地形阈值对高程控制点进行精细化提取。以波多黎各部分离岛及美属维尔京群岛高精度DEM（Digital Ele-vation Model）数据验证所提取高程控制点的准确性，实验结果表明：相较于原始数据，本文方法提取的高程控制点在平地、丘陵和山地区域的平均绝对误差与中误差分别从2.65 m/7.23 m、3.92 m/7.65 m、4.93 m/8.29 m减小为0.28 m/0.61 m、0.46 m/0.79 m、0.63 m/0.89 m；在保证精度的前提下，本文方法提取的高程控制点数量显著多于现有方法；利用本文方法，可从海岛区域提取一定数量、精度可靠的高程控制点，能为后续全球海图修测与质量控制提供有力数据支撑。关键词：海岛（礁）；ICESat-2；ATL08；高程控制点；自适应阈值；星载激光雷达；光子计数；海图修测1 引言海岛（礁）地形测量是海洋测绘的重要组成部分，高程控制点是地形测量与精度检验的重要基础数据1-2。由于远离大陆、分布零散等现实因素，通过传统外业方式获取高程控制点效率低、成本高、风险大，难以保证数据对现势性的要求，且对于部分位于敏感区域的海岛，实测控制信息的获取几无可能。星载激光雷达（Light Detection and Ranging，LiDAR）凭借其无地域限制、全天时、远距离对地观测等特点，能够快速直接地获取全球范围内地物的三维空间信息3，并能辅助卫星影像区域网平差、提升无地面控制点条件下的影像定位精度4-5，已成为高程控制数据的可靠数据来源之一，可为海图修测与质量控制等海战场环境保障应用提供坚实数据基础6-7。目前已有多颗对地观测卫星搭载了LiDAR测高系统，依其测量体制，可分为全波形LiDAR和光子计数LiDAR共2类。其中，星载全波形LiDAR技术较为成熟、平台众多，已被广泛应用于高程控制点提取8。李国元等9、谭建伟等10、Li等11基于ICESat（Ice,Cloud and land Elevation Satellite）数据，分别通过构建多准则约束、波形分解与拟合、细化理论误差模型等多种方式，提取了部分满足1:50 000控制需求的高程控制点，构建了覆盖全球大陆的高程控制点数据集。但该数据足印直径相对较大（约70 m），数据精度不够稳定，且该卫星已于2009年失效，无法满足海图修测与质量控制对海岛高程控制点现势性的需求。我国于2019年11月发射的高分七号卫星也搭载了全波形LiDAR测高仪，经在轨检校和多参数约束后，所提取的激光测高点平均测量精度可达0.5 m12-13。然而，该数据的足印间距较大（约2.4 km）14，经质量控制后分布较为稀疏，一定程度上限制了其在海岛高程控制点获取中的应用。作为当前唯一搭载光子计数LiDAR的对地观测卫星，ICESat-2（Ice,Cloud and land Elevation Sat-ellite-2）于2018年9月发射，为获取更为精细的地表三维信息提供了新手段。但因测量体制差异，基于星载全波形LiDAR数据的高程控制点提取方法不能直接应用于ICESat-2数据，且ICESat-2发射和接收的均为弱信号，在数据获取过程中不可避免地会受到大气散射、太阳背景及仪器噪声等多种因素的干扰。现有研究表明，采集时间、植被覆盖类型、地形坡度等因素都会对ICESat-2的测高精度造成一定程度的影响15-18。因此，需要对数据进行一定筛选后，才能作为高程控制信息使用。由于具备强大的探测能力和多样化的使命任务，ICESat-2数据按级别可分为Level-1、Level-2、Level-3A、Level-3B四级，其下又细分为ATL00-ATL21（无ATL05）21种数据产品，数据性质与所含属性参数差异较大。其中，Level3A级数据ATL08（陆地和植被高程数据）是在Level2级数据ATL03（全球定位光子数据）的基础上，通过差分回归高斯自适应去噪算法（Differ-ential，Regressive，and Gaussian Adaptive NearestNeighbor，DRAGANN）得到信号光子后，采用迭代滤波方式将信号光子分类为冠层和地表光子，并以100m为统计单元拟合得出中心点的冠层和地表高程值19。相较于记录了沿轨密集分布光子信息的ATL03数据，ATL08作为更高级的数据产品，数据中包含部分测量过程中记录的大气环境参数、高程拟合时的统计特征、粗略的全球参考高程等信息，且数据量相对较小，益于高程控制点的快速提取。依据这些先验知识，王密等20基于参考高程和部分属性参数对基于ICESat-2 ATL08数据的高程控制点提取方法展开了研究。在此基础上，郑迎辉等21进一步扩充了筛选准则，将地面光子数量纳入考量范围，在一定程度上提升了所提取高程控制点的精度。但上述方法均仅基于ICESat-2 ATL08数据内部参数设置固定阈值，不能依据所需的测量精度指标要求自适应地开展高程控制点的提取，在地形复杂区域提取的高程控制点数量往往不能满足需求。同时，由于上述研究大多关注高程控制点在大15608期戴泽源 等：适用于海岛的ICESat-2高程控制点提取方法范围、大比例尺测图中的应用，因此仅在内陆地区对所提取的控制点进行了精度验证，ICESat-2数据是否适用于陆域面积相对较小的海岛高程控制点提取，仍有待进一步探索。为此，本文以ICESat-2ATL08数据为基础，对星载光子计数LiDAR数据高程控制点提取及其在海岛区域的适用性展开研究。首先去除数据中的海域点与粗差点，降低后续处理的数据量；然后构建多参数约束，对数据进行初步筛选；最后依据所需的高程控制精度，设置自适应地形阈值精细化提取高程控制点，并通过实验验证了所提方法的可行性。2 研究方法2.1 粗差点快速去除海岛四面环海，相对大陆面积较小，ATL08数据中绝大部分的高程点位于海部，为了提升后续高程控制点筛选的效率，有必要首先去除数据中的海域点。ATL08 数据中的 segment_landcover 参数提供了基于100 m空间分辨率哥白尼全球土地覆盖（Copernicus Global Land Cover）数据的地表覆盖类型信息。本文利用该参数，结合研究区域公开发行的海图数据构建水域掩膜标记，通过属性分析与叠置分析快速提取数据中的陆域点。同时，仪器误差、极端大气条件等因素都可能造成高程值测量异常，数据中不可避免地存在精度较差的粗差点，这些粗差点可通过同参考高程数据的比较实现快速去除。在统计意义上，中值对离群点的探测更为敏感22，为此，本文提取ATL08数据参数中的地表高程中值（h_te_median）与插值地表高程值（h_te_in-terp），并以式（1）为基础计算同ATL08数据中自带参考高程值（dem_h）的高程差，将高程差大于阈值的高程点视为粗差点予以去除。|Hmed-Href3ref|Hinterp-Href3ref（1）式中：Hmed、Hinterp、Href分别为地表高程中值、插值地表高程值与参考高程值；ref表示参考高程的绝对高程精度。通常，3ref的值可取为25 m22。2.2 多参数约束的高程控制点初步筛选ICESat-2数据的高程测量精度可能受到多种外部因素的复合影响，为了保证所提取高程控制点的有效性、精确性，本文在深入分析ATL08数据算法理论文档22的基础上，结合已有文献对ATL08数据的精度评价结论16-18，对现有研究20-21中所使用的约束指标进行了一定综合与扩充，从大气环境、信噪比、数据完整性3个角度出发，构建多参数约束，以实现对高程控制点的初步筛选。星载LiDAR所发射的激光脉冲在往返大气过程中，复杂的大气环境会使其产生散射，导致能量衰减，影响测高精度，但这为ICESat-2开展大气、云与气溶胶探测提供了理论支撑23，ATL08数据中也包含了观测时大气状态的探测参数。本文首先通过多重散射警告标记（msw_flag）判断当前数据是否受到散射影响：若不存在散射影响，保留该数据；若存在散射影响，进一步分析散射影响来源，当散射 影 响 仅 同 云 层 有 关 时，通 过 云 量 标 记（cloud_flag_atm）判断当前云量情况，并仅保留统计单元内云量小于20%的高程点。现有研究表明，信噪比也是影响ICESat-2测高精度的重要因素24。用于生成ATL08数据的DRA-GANN算法从低信噪比数据中提取的信号光子置信度较差，当原始光子点云信噪比小于1:9时，即便通过人工判读也很难对信号与噪声光子进行准确区分25。为减少信噪比对高程值的影响，本文只保留信噪比参数（snr）大于1:3的高程点。ATL08数据中的高程值是基于统计单元内的信号光子加权插值拟合得出的22。因此，每个统计单元内的信号光子绝对数量与空间分布也对最终所提取的高程控制点精度有很大影响20,26。对光子的绝对数量而言，地表光子数过低可能意味着冠层密集或足印内存在高反射率地物，光子无法到达地表，导致用于计算地形高程的光子数量不足。在ATL08产品生成过程中，当每个统计单元中的信号光子数过少（50）时，将不计算后续的冠层与地表高程值。基于此前提，本文进一步细化约束条件，只保留统计单元中地表光子数（n_te_photons）不少于50的高程点。为了保证高程拟合结果的稳定性，地表光子在理想状态下应在沿轨方向上连续、均匀分布，现有方法大多没有对这一重要因素进行评估。在最新版本（V005）的ATL08数据中，提供了h_te_best_fit_20 m参数来表征每个统计单元内以20 m间隔构建的子统计单元地表光子的最佳拟合高程27，每个100 m统计单元内含有5个20 m子统计单元，这一子统计单元的大小已十分接近ICESat-2原始足印（约17 m）。当子统计单元内的信号光子及地表光子过少（10）时，该值被填充为一个无效值。本文通过对统计单1561地 球 信 息 科 学 学 报2023年元内所有h_te_best_fit_20 m参数进行有效值判定，保留中心点所处及其前后共计3个子统计单元内h_te_best_fit_20 m均为有效值的高程点，间接评估地表光子在沿轨方向上分布的连续性和完整性。2.3 顾及自适应地形阈值的高程控制点精细化提取相较于机/船载LiDAR，星载LiDAR所对应的激光足印较大，已显著区别于传统意义上的“点”，地形对星载 LiDAR 测高精度的影响往往不可忽视。通常认为地形起伏越大、地表越粗糙，ATL08数据中的高程值精度越不稳定。现有方法仅基于ATL08数据中的坡度参数（terrain_slope），采用固定阈值对数据进行筛选，不能依据实际所需精度指标自适应地进行精细化提取，尤其在部分存在山地、丘陵等复杂地形的海岛中可能导致所提取的高程控制点数量较少、空间分布稀疏，而这些高程控制点在辅助影像区域网平差与数据质量控制中恰恰是必不可缺的4。为了自适应地提取高程控制点，进一步分析地表光子在沿轨和高程方向上分布的离散程度，本文首先基于ATL08数据中的坡度参数计算初步筛选所得高程点的坡度，依其坡度情况分别获取平地（2）、丘陵（26）、山地（625）3种地形类型下的高程点，结合ATL08 数据拟合过程中100 m 统计单元内地表光子的高程误差预估值（sigma_atlas_land）、高程标准差（h_te_std）、地表光子高程偏度（h_te_skew）和地表光子数等参数，依据高程控制点所需精度自适应地设定如下式所示的精细化提取准则。landTskewT|hgt-estT（2）est=i=1Nte(iNte-12)H2Nte-1（3）式中：land、skew、hgt和est分别为100 m统计单元内地表光子的高程误差预估值、高程偏度绝对值、高程标准差和基于坡度参数计算的标准差预估值；est的计算形式可如式（3）所示，Nte、H分别为100 m统计单元内地表光子个数和基于坡度参数估算的地形高差；T为所需的高程控制点精度阈值。在海图编绘中，陆部地形要素主要从地形图转绘，1:10 000以上比例尺海图修测与质量控制主要采用1:50 000测图数据为参考28。基于此，本文将T在平地、丘陵和山地地形类型下的值分别取为0.8、1.0和1.2 m29。综上所述，本文方法提取高程控制点的流程可如图1所示，该图中涉及到的变量可参看式（1）式（3）所示，涉及到的ATL08数据参数及其简要释义如表1所示。图1 本文方法海岛高程控制点提取流程Fig.1 Flow chart of the proposed elevation control pointextraction method15628期戴泽源 等：适用于海岛的ICESat-2高程控制点提取方法3 研究区概况与数据来源3.1 研究区概况本文选择的研究区域如图2所示，该区域位于加勒比海北部边缘，包括波多黎各部分离岛及美属维尔京群岛。其中，波多黎各离岛包括别克斯岛、库莱布拉岛及其周边部分海岛（礁），陆域面积合计约163.92 km2，地形以丘陵为主；维尔京群岛包括圣克罗伊岛、圣托马斯岛、圣约翰岛及其周边部分海岛（礁），陆域面积合计约350.58 km2，圣克罗伊岛地势北高南低、面积相对较大，其余两岛则地形起伏较大，多丘陵、山地，面积相对较小。该区域远离敏感海域，存在现势性较强、精度较高的参考DEM数据，可对本文方法进行客观、定量的评价。同时，该区域是典型的群岛区域，而此类区域正是海图修测与质量控制任务中较为关注的重点海岛（礁）区域。不仅如此，该区域内海岛多为火山/大陆岛，海岸陡峭、地形特征复杂、地表覆盖类型丰富，能够进一步验证本文方法在复杂海岛（礁）区域的鲁棒性。3.2 研究数据及来源3.2.1 ICESat-2ATL08数据本文从美国国家航天局Earth-data Search平台（https:/search.earthdata.nasa.gov）收集了研究区域2018年10月至2022年8月的所有ATL08数据作为研究的主要数据，去除海域点后，其陆域高程点分布情况如图3所示。3.2.2 高精度DEM验证数据为验证本文方法所提取的高程控制点的精度，本文从 3DEP Lidar Explorer 平台（https:/apps.na-tionalmap.gov/lidar-explorer/#/）收集了美国地质调查局公开发布的1 m分辨率DEM数据产品作为高精度验证数据，其在研究区域内的数据覆盖范围如图4所示。该数据测量起讫时间为2018年4月至2019年4月，同ATL08数据观测时间相近，现势性较强，且经数据发行方验证，所用高精度DEM数据垂直精度优于0.1 m，可用作精度评价的参考验证数据。3.3 研究数据预处理需要指出的是，本文所采用的研究数据各自具有不同的平面和高程基准：ICESat-2 ATL08数据的表1 本文方法用于高程控制点提取的ATL08数据参数Tab.1 ATL08 data parameter used in the proposedelevation point extraction method参数名segment_landcoverh_te_medianh_te_interpdem_hmsw_flagcloud_flag_atmsnrn_te_photonsh_te_best_fit_20mh_te_skewsigma_atlas_landh_te_stdterrain_slope简要释义统计单元内的地物类型统计单元内地表光子高程中值统计单元内插值地表高程值参考高程值多重散射警告标记云量标记统计单元内的信噪比统计单元内地的地表光子数20 m子统计单元内的最佳拟合高程值统计单元内地表光子的高程偏度统计单元内地表光子的高程误差预估值统计单元内地表光子的高程标准差统计单元内沿轨方向的地面坡度，通过线性拟合得出，数值上等同于沿轨方向上单位长度的高程变化量图2 研究区域海岛（礁）分布示意图Fig.2 Schematic diagram of island distributionin the study area图3 研究区域ATL08陆域高程点分布示意图Fig.3 Schematic diagram ofATL08 land elevation pointsdistribution in the study area1563地 球 信 息 科 学 学 报2023年平面坐标为WGS84，高程为WGS84大地高；高精度DEM验证数据平面坐标为NAD83，高程基准为PRVD02及VIVD09正常高。为保证后续实验结果的准确性，有必要对数据进行基准转换预处理。本文使用了美国国家大气和海洋局发布的基准转换软件VDatum（https:/vdatum.noaa.gov/）将数据转换至相同的参考基准。4 实验与分析利用高精度参考DEM数据对本文研究区域提取前后的高程控制点精度进行验证，将经过海域点去除后的有效陆域点视为原始数据，不同筛选条件下所得的高程点数、高程平均绝对误差（MeanAbsolute Error，MAE）、高 程 中 误 差（Root MeanSquare Error，RMSE）和数据留存率如表2所示；本文方法提取前后的高程残差直方图如图5所示，所提取的高程控制点空间分布如图6所示。综合分析上述精度验证结果，去除海域点后，研究区域内一共含有15 388个陆域高程点，原始高程点MAE和RMSE分别为3.86 m和7.89 m，无法直接用作高程控制。通过本文方法提取后共计保留了1 990个高程控制点，数据留存率为12.93%，MAE和RMSE分别为0.36 m和0.69 m，精度有了较大提升，且所提取的高程控制点残差分布比较集中，92.16%分布在-1 m,1 m区间内，不难看出，通过本文方法提取的高程控制点整体具有较高的高程精度。在高程控制点空间分布上，除了圣约翰岛因地形复杂、植被茂密、数据完整性较差导致高程控制点分布稀疏外，其余四岛经本文方法提取后均有一定高程控制点，能够满足海岛高程控制需求。这表明即便在陆域面积相对较小、地形相对复杂的海岛（礁）地区，本文方法也能够从ICESat-2 ATL08数据中提取一定数量的高程控制点，适用于海岛高程图4 研究区域高精度DEM验证数据示意图Fig.4 Schematic diagram of reference DEM datain the study area表2 高程控制点精度验证结果Tab.2 The accuracy of extracted elevation control point筛选条件原始数据粗差点去除多参数约束自适应地形阈值高程点数/个15 38811 0713 0551 990MAE/m3.861.530.570.36RMSE/m7.892.831.380.69数据留存率/%100.0071.9519.8512.93图5 高程点残差直方图Fig.5 The histogram of residuals for elevation points15648期戴泽源 等：适用于海岛的ICESat-2高程控制点提取方法控制点的提取。进一步分析不同地形类型下高程控制点的提取精度，提取前后高程点的残差直方图如图7图9所示。不失一般性地，在星载LiDAR高程控制点的精度验证过程中，其真实误差可认为是高程控制点同高精度参考数据间的误差与高精度参考数据自身误差的叠加30。考虑到本文所用的高精度参考DEM数据的高程精度约为0.1 m，在平地、丘陵和山地地形下的精度阈值分别为0.8、1.0和1.2 m，依据误差传递原理，本文将高程残差分别在0.806、1.005和1.204范围内的高程点视为有效高程控制点，并统计不同地形类型、不同筛选条件下有效高程点在所提高程控制点中的占比，其结果如表3所示。图6 本文方法提取的高程控制点空间分布示意图Fig.6 The spatial distribution of extractedelevation control points图7 平地地形高程点残差直方图Fig.7 The histogram of residuals for elevation points in the flat region图8 丘陵地形高程点残差直方图Fig.8 The histogram of residuals for elevation points in the hilly region1565地 球 信 息 科 学 学 报2023年对比分析图7图9与表3中的统计结果可知，由于研究区域内的海岛大多为大陆岛或火山岛，多山地丘陵、海岸陡峭、地形较为复杂，经坡度分类后，平地、丘陵和山地类型下的陆域高程点 分别为5 618、4 301和5 138个，数量大致相同。3种地形类型原始高程点的高程残差分布较为分散，且相当一部分高程点残差较大，MAE和RMSE分别为2.65 m/7.23 m、3.92 m/7.65 m、4.93 m/8.29 m，呈现出梯次上升的趋势，精度符合率分别为61.65%、43.87%和32.21%，呈现出梯次下降的趋势，表明地形起伏对ATL08数据的高程精度确实存在显著影响。基于参考高程值去除粗差点后，3种地形类型高程点分别减少18.58%、30.64%和35.36%，精度符合率分别提升12.27%、16.07%和13.05%，这表明在ATL08数据中确实存在一定量的粗差数据，通过同参考高程值的简单比对，减少了参与多参数约束筛选的高程点数量，能提升后续高程点提取的计算效率。经多参数约束后，3种地形类型下的高程点的MAE分别提升了57.14%、56.29%和52.79%，RMSE分别提升了49.76%、47.25%和39.72%，精度符合率分别提升了16.66%、23.94%和29.51%，这表明本文提出的多参数约束能够综合大气环境、信噪比和数据完整性对高程点精度的影响，初步得到精度较为图9 山地地形高程点残差直方图Fig.9 The histogram of residuals for elevation points in the mountainous region表3 不同地形类型下高程控制点精度验证结果Tab.3 The accuracy of extracted elevation control point with different terrain type地形类型平地丘陵山地筛选条件原始数据粗差点去除多参数约束自适应地形阈值原始数据粗差点去除多参数约束自适应地形阈值原始数据粗差点去除多参数约束自适应地形阈值高程点数/个5 6184 5741 8361 2814 3012 9837635325 1383 321444177MAE/m2.650.910.390.283.921.510.660.464.932.331.100.63RMSE/m7.232.111.060.617.652.731.440.798.293.602.170.89数据留存率/%100.0081.4232.6822.80100.0069.3617.7412.37100.0064.648.643.44精度符合率/%61.6573.9290.5893.6043.8759.9483.8889.1032.2145.2674.7786.4415668期戴泽源 等：适用于海岛的ICESat-2高程控制点提取方法可信的数据，且在本文研究区内该约束不受地形环境限制，在本文所述的3种地形类型中均发挥了较大作用。自适应地形阈值数据剔除量较小，说明基于多参数约束初步筛选所得到的高程点已具有了较好的高程精度，但精细化提取后的高程控制点MAE与 RMSE 仍分别提升了 28.20%/42.45%、30.30%/45.14%和 42.72/58.99%，精度符合率分别提升了3.02%、5.22%和11.67%。进一步对比自适应地形阈值提取前后符合精度限差（2倍精度阈值）的高程点占比，如表4所示，精细化提取后限差符合率分别提升了2.46%、2.54%、和6.41%。这表明初步提取的高程点中仍存在部分不合精度阈值的高程粗差点，而本文通过设定自适应地形阈值对高程点进行精细化提取，使不同类型下所提取的高程控制点整体RMSE均能满足精度阈值要求，提升了方法的鲁棒性。5 讨论为进一步比较本文方法同现有方法在高程控制点提取数量与精度上的差异，使用文献20和文献21方法在本文研究区域内提取高程控制点，不同方法的实验结果如表5所示。结果表明，本文方法在研究区域内提取的高程控制点不仅精度较现有方法有了一定提升，还能在保证提取精度的前提上，提取到更多高程控制点。这是因为本文方法能够依据所需测量精度，自适应地在不同地形类型中设定阈值，扩展了高程点的筛选范围，并结合粗差点去除、多参数约束等多种手段，保证了所提取高程点的精度，更适用于面积相对较小、地形相对复杂的海岛高程控制点提取。当然，由于海岛高精度参考数据获取的困难性，本文仅在研究区域内对所提取的高程控制点进行了精度验证与对比分析，更多区域的实验验证还有待后续展开。同时，本文方法在自适应阈值的设计中仅考虑了地形坡度这一基本要素，所提取的高程控制点在ICESat-2沿轨方向和垂轨方向的分布也仍存在一定差异，如何在更多复杂区域（如含有建筑物的人造地表、潮汐影响下的海岸带）结合海图修测与质量控制任务对陆部要素高程控制信息实际需求，扩充自适应阈值的约束要素和选取规则，这些都还有待后续进一步研究。6 结论本文提出了一种适用于海岛的ICESat-2高程控制点提取方法，综合分析现有数据参数，去除数据中的海域点与粗差点，降低后续处理的数据量；然后构建多参数约束，对数据进行初步筛选；最后依据所需的高程控制精度，在不同地形类型中设置自适应阈值精细化提取高程控制点。经过精度验证与对比分析，得结论如下：（1）相较于研究区域内的原始数据，利用本文方法从不同地形类型下所提取的高程控制点整体中误差分别从7.23、7.65、8.29 m减小为0.61、0.79、0.89 m，均符合所需精度阈值要求。（2）与现有方法局限于固定坡度阈值相比，本文所提方法能够实现基于自适应地形阈值的高程控制点提取，扩展了高程控制点的筛选范围，在保证精度的前提下，本文方法在研究区域内提取的高程控制点数量是现有方法的45倍。（3）利用本文方法，可从陆域面积相对较小、地形相对复杂的海岛（礁）地区提取一定数量、精度可靠的高程控制点，能为后续全球海图修测与质量控制提供有力数据支撑。表4 不同地形类型下高程控制点限差符合率Tab.4 Compliance rate of elevation point tolerancebefore and after adaptive terrain threshold filtering地形类型平地丘陵山地筛选条件自适应地形阈值筛选前自适应地形阈值筛选后自适应地形阈值筛选前自适应地形阈值筛选后自适应地形阈值筛选前自适应地形阈值筛选后限差符合率/%95.5998.0593.7096.2490.7797.18表5 不同高程控制点提取方法结果对比分析Tab.5 Comparison of elevation control point extractionresults by different methods提取方法文献20方法文献21方法本文方法（平地）本文方法（丘陵）本文方法（山地）本文方法（总计）高程点数/个3923201 2815321771 990MAE/m0.830.550.280.460.630.36RMSE/m1.771.510.610.790.890.69精度符合率/%74.7486.5693.6089.1086.4491.76限差符合率/%85.4693.1398.0596.2497.1897.491567地 球 信 息 科 学 学 报2023年参考文献(References)：1 申家双,翟国君,黄辰虎,等.海洋测绘学科体系研究(一):总论J.海洋测绘,2021,41(1):1-7.Shen J S,Zhai G J,Huang C H,et al.Research on discipline system of ma-rine surveying and mapping,part I:General introduc