《测绘学报》《测绘通报》《测绘工程》《卫星导航（英文）》《测绘科学技术学报》《北京测绘》《Journal of Geodesy and Geoinformation Science》》联合融媒体。

  装备在高分七号卫星上的是我国首个具备全波形记录功能的激光测高仪，大多数都用在获取地面稀疏的高程控制点，提高了同平台立体影像无地面控制点的立体测图精度。高分七号卫星激光测高标准化处理是测绘应用的关键步骤，所生成的激光测高标准产品是后续对外分发和业务化应用的重要前提。本文围绕高分七号卫星的激光数据，研究了激光测高数据处理方法，验证了激光测高标准产品的几何精度。选择几何定标区以及陕西华阴、德国北威州等多个验证区，结合高精度外业测量点和机载LiDAR-DSM数据，对高分七号卫星激光测高标准产品开展精度验证。验证根据结果得出，高分七号SLA03产品定标区两波束激光的平面精度分别为（3.896±1.029）m和（3.286±0.337）m、高程精度分别为（0.018±0.099）m和（-0.017±0.096）m。采用高程控制点质量控制参数ECP_Flag能有效标识出可用于高程控制的激光点，其中陕西华阴验证区两波束激光总体精度分别为（-0.113±2.519）m和（0.191±1.071）m，经质量控制后ECP_Flag标记为1的激光点高程精度为（0.111±0.152）m和（-0.064±0.115）m；德国北威州总体精度为（-0.897±5.485）m和（-0.202±6.207）m，ECP_Flag标记为1的激光点高程精度为（-0.304±0.190）m和（-0.279±0.220）m。目前高分七号卫星激光测高标准产品已在自然资源部国土卫星遥感应用中心实现业务化生产。

  激光测高卫星是对地观测卫星的重要组成部分[1-2]。美国先后于2003年、2018年发射了冰、云和陆地高程卫星(ice、cloud and land elevation satellite，ICESat)及后续卫星ICESat-2，在极地变化监测、湖泊水位测量、高精度地形测量等领域得到普遍应用，成为发展激光测高卫星的成功案例[3-4]，并于2018年12月在国际太空站上成功装备了全球ECO动力学调查(global ecosystem dynamics investigation, GEDI)[5]多波束激光测高仪。虽然我国目前还没有专门的激光测高卫星，但在2016年发射的资源三号02星，成功实现了试验性激光测高载荷对地观测的高程有效测量，并验证了提升立体影像无控高程精度的可行性[6-7]。2019年11月3日成功发射的高分七号(GF-7)卫星上装备的业务化应用的全波形激光测高仪，能为快速获取高程控制点、监测大型湖泊水位变化等提供高精度数据支撑[8]。

  标准化测绘处理是国产卫星激光测高产品走向工程化应用的基础和前提。我国对月观测的“嫦娥”系列卫星上搭载过激光测高系统[1, 9]，但其观测条件、仪器性能、精度指标等与对地观测有很大区别。此外，国外虽然先后建立了ICESat/GLAS、ICESat-2/ATLAS、GEDI等为代表的卫星激光测高产品体系[3-5, 10]，但其载荷与国内存在一定差别，如国产的高分七号卫星激光测高仪同时配备了足印相机而国外没有。因此，研究国产高分七号卫星激光测高标准化处理方法及产品设计，有很重要的现实意义和应用价值。

  围绕激光测高数据处理及产品设计，文献[11]研究了卫星激光测高严密几何模型、对影响激光精度的卫星轨道和姿态、指向角、大气、潮汐、光行差等进行了较深入地分析；文献[12]针对资源三号02星激光测高仪数据处理，并从指向角、观测时间、侧摆等方面做了质量评价分析；文献[13]针对大气折射对激光测距精度影响进行了研究，目前已形成较成熟的方案，但目前大气散射的影响及校正还有待深化[14-15]；文献[16-19]对卫星激光测高误差及精度等做多元化的分析，并精确指出坡度是影响测高精度的一个主要的因素[19]。国产卫星激光测高产品设计相关研究则基本为空白。本文结合高分七号激光测高数据特点，研究了高分七号激光测高产品分级、标准产品数据处理流程、产品结构等内容，然后利用多个验证区外业测量结果验证了标准产品平面和高程精度。所得结论对我国后续激光测高卫星数据处理、精度评定具备极其重大的示范价值。

  高分七号卫星激光测高系统包括2波束激光器, 以3 Hz的工作频率向地面发射1064 nm波长的激光脉冲，在地面形成沿轨间隔约2.4 km、垂轨间隔约12.25 km的离散激光光斑，如图 1所示。对应20 km×20 km幅宽的线 高分七号卫星激光测高数据概况

  高分七号卫星激光测高系统在获取全波形数据的同时，利用足印相机记录激光的发射光斑强度分布及地面落点位置周围的地物影像，通过影像匹配实现激光与线阵遥感影像的几何关联。激光器与足印相机有两种协同工作模式，即同步模式和异步模式。在同步模式下，足印相机对激光光斑和地物同时成像，形成一幅具有激光光斑的足印影像，为避免偏亮的光斑分布影像与实际落点位置地物影像重叠，影响地物的解译与识别，硬件设计时进行了激光波束指向在足印影像上约0.4°的偏移，实现激光光斑和实际落点地物的分离，通过外场几何定标可获得实际落点位置相对于光斑质心的偏移量，如图 2所示。异步模式下，足印相机在激光脉冲发射的前、后曝光形成两幅地物影像，且在激光发射时曝光形成光斑影像，如图 3所示。

  激光足印相机指向记录方式的误差来源于激光导光光路(图 4中红色线部分)，在激光光轴指向与足印相机耦合关系很稳定的前提下，该误差理论上可视为一个常数，对应在足印影像(laser footprint image，LFI)上为一个平移量。在实验室通过严密的测量建立了激光光斑质心与落点位置的转换关系，通过在轨定标能确定激光的落点在足印影像的真实位置，通过监视光斑质心变化能确定实际落点位置是否发生明显的变化。异步模式下可通过激光落点标定位置和发射前后的地物影像内插，获得激光在地物影像上的实际像素位置。

  SLA03作为激光标准测绘产品，目前已在自然资源部国土卫星遥感应用中心实现业务化生产，可以对外分发供相关用户使用[21]。同时随着数据的积累，正在慢慢地构建全球高程控制点数据产品SLA05以及大型湖泊水位产品SLA06。鉴于高分七号激光较低的重频率和极地覆盖能力，目前林业及极地相关这类的产品并未纳入产品体系中。同时考虑到大气在对地观测卫星中的特殊影响，气溶胶、光学厚度、云高等大气专题产品SLA04目前作为备选项进行了保留，后期将结合高分七号同平台的多光谱影像、足印影像及第三方的有关数据，逐步构建大气专题产品，对激光的精细化应用提供支撑。

  高分七号卫星激光测高标准产品SLA03是其他类专题产品的基础，类似ICESat-2卫星的ATLAS系列新产品中的ATL03(全球激光点地理定位产品)。高分七号后续的高程控制点、湖泊水位等专题产品均可由SLA03标准产品经深加工而生成。标准产品SLA03的处理流程如图 6所示，针对高分七号卫星激光测高数据，采用事后处理的精密轨道和姿态数据，并对大气、潮汐等环境影响进行精细改正得到的精确三维坐标[2, 11, 22]、具有精确地理信息的足印影像、标准化的波形特征参数等。处理步骤如下。

  激光测高精度不可避免受大气、地形、地物等因素的影响，质量控制是激光数据处理过程中的一个重要环节，也是非成像类卫星产品的一个难点问题。为尽量满足后续高程控制点应用的自动识别，在高分七号卫星激光测高标准产品处理过程中，综合采用了波形脉宽、波峰数、足印影像基本质量判断、地形坡度等信息[22]，增加了ECP_Flag (elevation control point flag)字段，共标记了8类，除标记为10的代表落水外，剩下的7级用于标识激光点作为高程控制点的可用性。SLA03产品中ECP_Flag标记为1、2、3的激光点表示建议高程控制点使用，其中标记为1的代表地形地物单一、坡度小于2°，高程精度非常高，理论上应优于0.17 m；标记为2的代表坡度小于5°且回波中主峰面积和能量占优，精度较高，理论上优于0.43 m；标记为3的代表坡度小于7.5°，回波中主峰面积和能量基本占优，理论上精度优于0.65 m，但因大光斑内地物影响，精度的可信度可能有一定不足；标记为4~6的不推荐为控制点，主要是地物类型很复杂、波形分解的波峰数不小于3个，且波形有展宽，但可考虑用于林业树高、建筑物高度等特征提取应用；标记为7代表该点信噪非常低、可靠性较差；标记为10的代表根据地理位置和GLC30(global land cover, 地表覆盖分类判断该点落水，是位于水体上的激光点，不能用于高程控制点使用，后续可考虑用于大型湖泊水位测量应用。

  高分七号卫星激光测高标准产品SLA03中包括波形数据、影像数据、激光足印三维坐标以及各种特征参数。主体文件以HDF5格式存储，在HDF5文件内部，每个激光点存储在一个组(Group)中，每个激光点的数据由若干字段组成，根据字段的详细的细节内容分为6个子组(SubGroup)，包括综合信息(Basic_Information)、足印影像信息(LFI_Information)、波形信息(Waveform_Feature)、地形地物信息(TerrainFeature)、地球物理信息(Geophysic)、其他字段(Other)。其中，综合信息中包含激光点的编号、三维坐标、所属激光器等，足印影像信息中包含足印影像分辨率、激光落点在足印影像上的像素位置、足印影像数据体等，波形信息包含经滤波去噪等预处理后的发射和回波波形，经高斯分解后的波形特征参数等，地形地物信息中主要包含基于波形提取的激光光斑内的高程分层值、地表坡度等，地球物理信息中包含大气折射延迟距离改正值以及固体潮、极潮、海潮和负荷潮改正值，其他字段包含基于GLC30获取的激光落点处地物覆盖类别、陆海标识。

  激光测高标准产品主要有数据主体文件、元数据文件、空间分布文件、激光三维坐标文件(表 3)。其中，主体文件为产品主体，存储了激光测高标准产品的所有信息。元数据文件为辅助文件，存储了激光测高标准产品的元信息。空间分布文件为Shape文件，包含了配套的辅助文件，一起存放于LaserRange文件夹中。激光三维坐标文件为文本文件，提供各个激光点的索引、经度、纬度、高程、高程可用性标记等信息。

  高分七号卫星在轨测试期间，自然资源部国土卫星遥感应用中心、中国资源卫星应用中心等多家单位联合，在2020年6月14日、6月15日、7月14日、7月19日先后多次在几何定标场准确捕捉到激光光斑[20]，利用探测器的中心位置可以对定标后经标准处理生产的SLA03产品做平面绝对精度评价，结果见表 4。

  表 4中计算高程值为空的代表该激光点所在区域地面布设了角反射器，波形出现了饱和现象。由表 4能够准确的看出，激光点的平面精度在同一天的同一波束内有很好的一致性，如波束1在7月14日的编号为959536149和959536153的两个激光点平面误差较差为0.12 m，对应激光指向短时间的稳定精度为0.05″；波束2在7月19日的编号为964719118和964719122的两个激光点平面误差较差为0.58 m，对应激光指向短时间的稳定精度约0.24″。试验区波束1的激光点平面绝对误差最大为4.861 m，波束2最大为3.489 m，基于地面探测器统计的两波束激光的平面精度分别为(3.896±1.029)m和(3.286±0.337)m。

  为评价激光点的绝对高程精度，采用定标区，以及陕西华阴地区、德国北威州地区的高分七号多期实际数据来进行评价，高程基准均统一为WGS-84椭球的大地高。陕西华阴地区的高程范围为[299.3 m, 1535.7 m]，平地和山区基本各占一半；德国北威州验证区的高程范围为[73.8 m, 613.2 m]，绝大部分属于城市平坦地区，部分属于丘陵和山区。其中后两个区域的激光点与定标区在时间和空间上均有一定距离，更能反映高分七号激光标准产品的真实精度水平。

  在外场定标区2020年6月9日、6月14日、6月19日、6月24日分别实测了少数的RTK-GPS点，与实际计算的高程值的误差见表 5。波束1和波束2参与统计的个数分别为22个，精度分别为(0.018±0.099)m，(-0.017±0.096)m。即与定标区时间和空间临近的激光点，在坡度小于2°的平坦地区绝对高程精度优于0.10 m。

  选择与内蒙古定标区有一段时间和空间间隔的陕西验证区开展绝对精度验证，数据获取时间分别为2020年4月26日、5月1日的2662轨和2738轨。点位分布如图 7所示，针对激光点的落点位置经纬度，以3~5 m间隔往外扩20 m左右，利用RTK-GPS采集地面点的三维精确坐标，如图 8所示。波束1共45个点、波束2共25个点，经统计波束1和波束2的总体高程精度分别为(-0.113±2.519)m和(0.191±1.071)m，相关结果见表 6。表 6中ECP_Flag字段标记为1的个数均为11个，精度分别为(0.111±0.152)m、(-0.064±0.115)m；标记为2的分别为10个和5个，精度分别为(0.246±0.229)m、(0.122±0.269)m。

  由表 6能够准确的看出，标记为1的激光点绝对高程精度非常高，可完全作为高程控制点使用；标记为2的次之，精度基本在0.3 m以内；标记为3的精度在0.5 m左右、可以考虑使用；剩下的为地形、地物复杂区或信噪比较低的点精度较差，不能作为控制点使用。从统计结果来看，标记为1和2的标记方法完全可信、精度可靠，可当作高精度高程控制点使用。

  进一步地选择境外德国北威州地区共5轨高分七号卫星激光测高数据，如图 9所示。分别为2020年6月21日的第3516轨、7月6日的第3745轨、7月11日的第3820轨、7月21日的第3973轨、8月29日的第4564轨。由于该区域激光点分布范围较大，考虑到激光存在一定间距，因此图中的圆圈要小于上文中的图 7，即圆圈仅代表激光点位置，不代表实际激光光斑大小。图中部分空白段代表该区域因云层较厚激光未达到地面，数据无效。参考高程来自于该区域高精度的机载LiDAR-DSM数据，格网大小为1 m，获取时间为2014-2016年，高程基准为德国采用的DHHN2016，采用EGM2008大地水准面模型转换为WGS-84椭球的大地高。经咨询国外该数据发布方的专家，高程基准转换精度在0.1~0.2 m，综合考虑LiDAR点云本身的精度，该区域参考数据的绝对高程精度约0.25 m。该区域SLA03产品的绝对高程精度评价结果见表 7，两波束激光的总体精度分别为(-0.897±5.485)m和(-0.202±6.207)m。ECP_Flag标记为1的高程精度分别为(-0.304±0.190)m和(-0.279±0.220)m，标记为2的高程精度分别为(-0.110±0.454)m、(0.024±0.501)m，能控制在0.5 m内，即ECP_Flag标记为1和2的可当作高程控制点使用。

  (3) 高分七号卫星两波束激光经质量控制后，在陕西华阴验证区ECP_Flag标记为1的激光点高程精度分别为(0.111±0.152)m和(-0.064±0.115)m，德国北威州验证区ECP_Flag标记为1的激光点高程精度分别为(-0.304±0.190)m和(-0.279±0.220)m。这两个区域标记为2的激光点高程精度波束1分别为(0.246±0.229)m和(-0.110±0.454)m，波束2分别为(0.122±0.269)m和(0.024±0.501)m。说明通过质量控制标记ECP_Flag，可以将精度优于0.5 m能用作高程控制的点进行相对有效识别，标记为1和2的均能用作高程控制点，其中标记为1的精度更高，其中陕西华阴验证区优于0.15 m、德国北威州验证区优于0.3 m。

  在资源三号02星试验性激光测高的基础上，高分七号虽不是专门的激光测高卫星，但其激光测高载荷作为国内首台业务化应用的对地观测激光测高仪，进一步实现了我国卫星激光测高技术的重大进步。与国外同类型的ICESat/GLAS全波形激光测高相比，高分七号激光光斑更小、精度相对略高，但在激光器硬件水平的重频率及应用场景多样性方面还有一定差距。此外，高分七号激光足印影像还存在分辨率偏低、辐射质量不稳定，在大角度侧摆下几何精度下降等问题。在肯定进步的同时，我们也需要正视问题，进一步创新，深入挖掘国产卫星激光测高的精度潜力、应用潜力，同时要加快推进陆海激光卫星的预研攻关[25-26]，争取早日发射国产首颗专业型号的激光测高卫星，在测绘、极地、林业、水利等多个行业实现规模化应用。