卫星激光测高数据辅助提取高程控制点方法与流程

本发明涉及一种高分遥感影像摄影测量、数据匹配和配准、多源遥感数据联合处理等领域，重点涉及激光测高数据与dsm、影像数据的空间匹配和配准方法。
背景技术：
：卫星雷达测高技术研究主要集中在美国和欧洲一些发达国家，而其他国家，尤其是发展中国家该技术的发展则相对滞后，目前世界上已相继实施了10余项卫星测高计划。2003年美国发射的激光测高卫星icesat-1上搭载了地球科学测高传感器glas，其轨道高度590km，计划于2017年发射的icesat-2星上搭载一个先进地形激光测高仪atlas，采用的微脉冲多波束(3×2结构的6束激光)光子计数激光雷达技术克服glas-1快速能量消耗问题，足印大小为10m，测高精度10cm，而计划于2020年发射的list通过阵列探测器以5m分辨率和10cm精度扫描地球，其采用的技术表明未来的激光测高系统将向着高重复、多波束、扫描式、单光子模式发展。目前我国载有测高传感器的卫星有hy-2a和资源三号，计划于2018年发射的高分7号激光测高仪专门为光学影像立体制图而研制。国内外代表性的激光测高卫星见表1。卫星激光测高初期任务以获取海面形状、研究大洋环流和其它海洋学参数为目的。卫星测高数据的应用也逐步从一个个点的变化，逐步扩展到整个面上的变化监测的研究，在海洋学、大地测量学和地球物理学领域也得到了广泛的应用。利用卫星测高所获取的观测量作为边界条件建立海洋动力模型，计算出海洋的深度，从而绘制出海底的地形、地貌。激光测高雷达目标几何定位一般使用严密模型，通过卫星的轨道姿态信息确定的激光波束的空间向量进行构建。卫星测高精度的提高主要依靠高度计传感器的改进和波形重构等数据处理算法的改进。激光测高仪系统的接收信号脉冲可以看作是发射脉冲卷积目标响应函数，地表反射率的变化，会引起回波脉冲能量的强度、峰值发生变化，光斑区地物目标成份、地形差异，会影响其回波波形地表结构的变化，会导致激光脉冲飞行时间的变化，从而导致回波脉冲的波形发生裂变或展宽，波形与地物地貌特征之间，存在密切的联系。表1代表性的激光测高卫星名称国别发射时间高程精度(厘米)足印(千米)jason-1,2,3法/美2001,2008,20164.2,2.5-3.42.2/2.2encisat欧洲20022.51.7icesat-1美国2003150.07cryosat-1,2欧洲2005,20101-31.6hy-2a中国201142zy3-02中国20165730.05理论上，通过光斑中心几何定位处理可以从激光测高数据中提取光斑位置和高程信息，但激光光斑定位平面误差较大、激光光斑直径通常较大，难以从激光测高数据中提取比较复杂地形特定点的高程，同时光斑分布稀疏不均匀，激光测高数据目前还难以直接用于大比例尺地形获取。由于激光测高精度优势，近年来其在光学与sar影像处理中逐渐得到重视。如将icesat/glas数据作为高程控制用于asterdem数据生产中，激光点密集区域直接融合到dem中；利用icesat/glas数据对insar高程进行纠正，获取高精度南极洲dem。针对当前我国高分遥感卫星在完全无控条件下定位精度尤其是高程精度仍难以满足大比例尺制图的难点，文献提出了一种多准则约束的激光高程控制点筛选算法，将筛选后的平地激光足印点数据作为遥感影像高程控制。已有文献对激光高程控制的资三影像定位进行了试验，表明无控条件下可以达到1：5万测图精度，相关学者提出了激光测距数据辅助两线阵影像处理的若干思路，利用激光测距数据与两线阵影像光束法平差仿真试验表明，能有效改善航线模型系统变形。从大光斑激光测高数据中提取不同地形影像上某高分影像像点的高程，目前还鲜有文献涉及。由于影像目标定位、激光遥感目标定位均存在误差，遥感机理的差异使得两者的相对空间关系目前还没有有效的方法实现严格的配准。同时，激光测高数据是光斑区目标的综合反应，遥感影像需要的高程控制点信息包括影像像点坐标及其对应的高程，从激光测高数据中高精度提取像点的高程信息，缺少稳健有效方法。激光测高数据与光学、sar数据间空间联系的缺失，使得激光测高数据在航天摄影测量领域的应用受到严重的制约和局限，激光测高数据辅助精准提取某像点或地物点的高程信息对卫星摄影测量学科发展有着非常重要的意义，也蕴藏了巨大的经济效益。参考文献：[1]李建成,范春波,褚永海,等.2008.icesat卫星确定南极冰盖高程模型研究.武汉大学学报·信息科学版,33(3):226–228.[2]唐新明,李国元,高小明,等.卫星激光测高严密几何模型构建及精度初步验证[j].测绘学报,2016,45(10):1182-1191.[3]杨乐,林明森,张有广,等.中国近岸海域高度计jason-1测量数据的波形重构算法研究[j].海洋学报,2010.32(6):91-100.[4]周辉,李松.激光测高仪接收信号波形模拟器[j].中国激光,2006,33(10):1402-1407.[5]崔云霞,牛燕雄,颜国强,等.自然地物对星载激光测高仪回波特性的影响[j].激光技术,2012,36(4):490-493.[6]wangxy,huanghb,gongp,etal.forestcanopyheightextractioninruggedareaswithicesat/glasdata[j].ieeetransactionsongeoscienceandremotesensing,2014,52(8):4650-4658.[7]李国元,唐新明,张重阳,等.2017.多准则约束的icesat/glas高程控制点筛选.遥感学报,21(1):96–104.[8]鄂栋臣,沈强,徐莹,cheng.2009.基于aster立体数据和icesat/glas测高数据融合高精度提取南极地区地形信息.中国科学d辑:地球科学,39(3):351–359.[9]yamanokuchit,doik,shibuyak.combineduseofinsarandicesat/glasdataforhighaccuracydemgenerationonantarctica[c]//geoscienceandremotesensingsymposium,2007.igarss2007.ieeeinternational.ieeexplore,2007:1229-1231.[10]ligy,tangxm,gaoxm,wanghbandwangy.2016.zy-3blockadjustmentsupportedbyglaslaseraltimetrydata.thephotogrammetricrecord,31(153):88–107.[11]王任享,王建荣.二线阵ccd卫星影像联合激光测距数据光束法平差技术[j].测绘科学技术学报,2014,31(1)1-4.技术实现要素：(1)发明的基本原理首先实现激光光斑与dsm的空间配准，以及光斑区dsm高程系统差提取，进而实现高程控制点提取和其它应用。利用卫星遥感一段时间获取的数据或一定范围内测绘产品误差主要呈现系统性的特点，将同轨一段轨道上的实测激光测高数据当作参考组，并以相同的偏移量对参考组各理论光斑中心周边dsm目标进行搜索和模拟观测，不同偏移量构成不同的模拟观测组，通过整体考察实测组波形与模拟观测组波形的匹配度，寻找dsm上光斑的实际位置。实测组波形作为参考模板时，针对不同的模拟观测搜索目标，有：1)当dsm上模拟观测组目标和实际观测组目标位置上没有误差时，理论上dsm模拟观测组波形与实测组波形的位置和形状一致，图1中的(a)；2)当dsm上模拟观测组目标和实际观测组目标位置一致，仅存在高程系统差时，时间维上整体滑动参考组波形能够使得两种数据达到最优匹配，图1中的(b)；3)当dsm上模拟观测组目标与实际观测组目标位置存在整体偏移时，模拟观测目标地势可能存在升高、下降、地形地势不变等情况，则对应模拟观测波形中心会在时间轴上左移、右移、不变，在地形变化情况下还会产生波形形状的改变，会导致时间维上的整体匹配测度指标显著下降，图1中的(c)。光斑与dsm配准的实现原理，就是寻找与实际观测组最优匹配的模拟观测组波形。即以实测组波形作为模板，以实测组各理论光斑中心为中心，以相同的变化区间和步长改变光斑中心位置，对周边dsm目标进行模拟观测获取系列模拟组波形，并以此为基础通过计算组内波形整体最优匹配条件下获取dsm上光斑的实际位置，如图2。在知道光斑在dsm上位置的情况下，通过模拟观测波形与实测波形在时间维上实现最优匹配时的参考波形在时间轴上滑动的距离，即能获取光斑区dsm的高程系统差和精确高程，进而提取高程控制点。(2)
发明内容本发明提供了一种卫星激光测高数据辅助提取高程控制点方法，通过波形分组匹配实现激光测高遥感目标与dsm平面和高程上的匹配。首先将一段同轨激光测高波形作为整体，将激光测高光斑区与dsm的平面位置配准转化为实测波形与dsm搜索区模拟观测波形的匹配；其次将光斑区dsm的高程系统差计算转化为各测站实现波形与光斑区dsm模拟观测波形间的匹配实现。本发明利用在一定范围内dsm和一段轨道范围内激光测高数据定位平面误差主要表现为系统性的特点，将一段同轨实测激光测高波形作为一组参考模板，以参考模板组中各激光光斑理论位置为中心，以相同的偏移量搜索dsm周边位置，将其作为对dsm数据进行模拟观测的光斑中心，获得多组模拟观测波形；通过实测组与各模拟观测组波形时间轴上进行整体匹配，获取最优匹配模拟观测组，提取最优匹配模拟观测组模拟观测时各光斑中心平面位置，实现光斑与dsm的平面配准；根据实测组各测站波形与最优匹配模拟观测组对应测站波形时间轴上最优匹配时的时间偏移量，获得组内光斑区dsm高程的系统差；利用波形分组匹配获取的光斑区位置和光斑区dsm高程系统差，实现高程控制点提取和测绘应用。本发明主要模拟观测包括以下步骤：步骤1.高分立体影像或insar定向参数精化：对高分立体影像或insar数据进行区域网平差或自由网平差，获取初步精化后的定向参数，用于后继步骤中高分立体影像或insar数据的处理；步骤2.高分立体影像或insardsm生成：高分立体影像或insar生成dsm：利用定向参数，通过测区高分立体影像的密集匹配，或通过insar干涉处理，生成带有地理坐标的dsm；步骤3.激光测高数据实测组划分：将测区激光测高数据按轨道进行分组，在激光数据所在轨道区间过长情况下，将同轨测高数据分为若干组；依次选取每一组实测数据并执行步骤4—步骤10，实现对dsm搜索目标的模拟观测与二次匹配，获取激光测高数据与dsm数据空间目标的平面配准位置、光斑区dsm的高程系统差和高程控制点；步骤4.实测组激光测高波形提取：提取实测组中各个激光测高实测波形，设去除错误和粗差数据后实测组中有n个测高波形数据，对这n个实测波形分别进行激光测距时间因仪器、对流层、潮汐、大气等影响的修正，并进行归一化处理；步骤5.模拟观测参数设计：根据定向参数计算实测组内各激光测高光斑中心的理论位置，作为dsm上光斑配准位置的搜索中心，并设计dsm上光斑中心的搜索区域和步长；以实测组中各测站实际观测参数为准，对各测站统一设计模拟激光测高观测定向参数增量区间和步长，使模拟观测能够对所有设计的dsm搜索区进行观测；具体通过步骤5.1—5.3实现：5.1设计搜索区间和搜索步长：根据光斑与dsm间的平面相对精度和dsm分辨率，设计模拟观测光斑中心的搜索区间和搜索步长；5.2模拟观测变化参数的选择：通过设计模拟激光测高观测定向参数变化区间和步长实现对dsm目标的搜索，以相同的增量同时改变实测组各测站观测姿态、或以相同的增量同时改变实测组各测站轨道位置观测参数、或以相同的增量同时改变各测站轨道和姿态参数，实测组其它观测参数不变，来设置模拟观测参数，实现dsm上模拟观测光斑中心位置相对与实测组理论位置的平面偏移和模拟观测dsm目标的改变；5.3模拟观测参数反算：根据dsm上模拟观测光斑中心的搜索区间和搜索步长，计算各模拟组观测参数相对于实测组观测参数变化的区间和步长，对模拟观测参数进行计算和配置；步骤6.模拟观测获取各模拟观测组波形：利用设计的模拟观测参数对dsm目标进行激光测高模拟观测，设每测站模拟观测数目为m，则可以获取m组模拟观测，每组模拟观测数目为n个；模拟观测时按dsm分辨率大小对地面光斑区进行空间分割，根据分割单元在光斑区的位置计算每一分割光斑单元的能量及其分布，卷积计算分割单元的能量被dsm上对应特定高程的模拟观测目标反射形成的子回波能量及其分布，模拟观测各采样时刻子回波能量的和形成模拟观测波形；对各模拟观测组中所有波形分别进行归一化处理用于后继步骤中的波形匹配；步骤7.实测组分别与各模拟观测组整体匹配：依次选取模拟观测组，将实测组中所有归一化波形整体当作匹配模板，沿时间轴上整体平移实测组各测站归一化波形，计算并记录下实测组分别与m个模拟观测组达到最佳匹配时的m个测度值，完成第一次波形匹配；本步骤中在时间轴上整体平移实测组各测站归一化波形后进行匹配，是为了消除模拟观测组dsm观测目标与实测观测目标高程上存在的系统差；在dsm高程系统差较小或可忽略的情况下，实测组波形与模拟观测组波形的整体匹配可以不在时间轴上整体平移实测组归一化波形实现，而是直接计算两者的匹配测度，用于后继最优匹配组选择；本步骤中实测组与模拟观测组波形匹配以波形对应点间的距离平方和、或模板匹配算法作为匹配测度，计算最佳匹配时的测度值，为了提升一些质量不理想dsm模拟观测波形和实测波形间的匹配效果，匹配时将所有波形脉冲宽度以固定的倍率沿时间轴放大，即以波形高斯均值点为中心，将波形两边分别沿时间轴两方向延展加宽，增加实测波形与模拟观测波形间的重叠度，减小dsm内部非系统差对匹配的影响；步骤8.光斑与dsm平面位置配准：从m个匹配测度值中选取最优匹配测度值，其对应的模拟观测组称为最优模拟观测组，其模拟观测时光斑中心位置及光斑范围，即为光斑区与dsm的平面配准位置，该组模拟观测时光斑中心位置与光斑中心理论位置差异即为配准前后的光斑平面位置偏移量；步骤9.实测组与最优模拟观测组各对应波形分别匹配：分别将实测组中各测站的归一化波形作为参考模板，通过在时间轴上平移参考波形模板，获取实测组各测站归一化波形与最优模拟观测组中对应测站的归一化波形最优匹配时参考波形在时间轴上的偏移量，完成第二次波形匹配；步骤10.各光斑区dsm高程系统差获取：根据时间偏移量和激光传播速度计算最优模拟观测组各光斑区dsm高程的系统差，根据dsm上三维点坐标和dsm高程系统差，精化并提取dsm上特征点作为高程控制点；步骤11.影像高程控制点的提取：根据光斑区dsm上的三维坐标点和生成dsm时高分立体或insar的定向参数，利用影像的几何模型计算该三维点在影像上对应的像点坐标，同时根据dsm上的三维点坐标和该光斑区dsm的高程系统差，计算该三维点修正后的高程值，则影像像点坐标和修正后的高程值组成一对影像高程控制点；步骤12.匹配结果的应用：根据平面和高程匹配结果，直接对dsm产品进行精化，或实现同轨高分传感器与激光测高仪的相对标定，或对遥感影像数据与激光测高数据进一步联合处理；其中，对dsm产品进行精化，根据匹配得到的各激光光斑区dsm高程的系统差，拟合整个测区dsm的高程误差，对dsm高程进行精化处理实现；或根据dsm与激光光斑目标间的高程系统差和提取的高程控制点，精化高分立体影像或insar处理参数；其中，同轨高分传感器与激光测高传感器的相对标定，利用同卫星平台同时期轨道姿态对影像定位和激光测高定位平面位置误差影响相同的特点，根据步骤8中获得的激光光斑中心理论位置和配准位置的相对偏移量以及卫星高度，反算激光测高传感器与高分传感器间的相对姿态修正值。其中，遥感数据与激光测高数据进一步联合处理通过遥感数据与激光测高数据或激光测高数据中提取的高程信息联合平差处理实现。本发明中步骤7—步骤10中，组内二次匹配也可以转化为一次匹配完成，即完成分组和模拟观测后，仍以各模拟观测组和实测组为对象，通过在时间轴上分别滑动组内各实测波形，实现组内各实测波形与模拟观测对应波形的最优匹配，并记录下最优匹配时波形匹配测度值和实测波形在时间维上滑动的距离，计算组内各波形最优匹配测度值均值s1，计算组内各实测波形在时间维上滑动的距离均值及各滑动距离与均值间差值的平方和s2，选取最优匹配测度均值s1和最小s2所在的模拟观测组作为最优模拟观测组，或综合比较s1和s2选取最优模拟观测组，最优模拟观测组光斑位置即为光斑与dsm配准位置，并利用最优模拟观测组中各实测波形与模拟波形最优匹配时在时间轴上滑动的距离即时间差计算对应光斑区高程系统差。本发明提供的激光测高数据与dsm模拟观测数据的波形匹配技术，与生成dsm的源数据没有关系，故不仅能够用于激光测高目标与高分立体影像生成的dsm或insar生成的dsm空间配准，也能实现激光测高目标与其它方法生成的dsm、dem、或点云的空间配准，获取dsm上激光测高光斑区位置和光斑区dsm、dem、或点云高程局部系统差。本发明提供的激光测高数据与dsm模拟观测数据的波形匹配技术，不仅适用于单波束激光测高数据，也能够适用于多波束激光测高数据。附图说明图1是所述不同观测目标的模拟观测组与实际观测组波形、位置关系图；图2是光斑与dsm配准技术实现的示意图；图3是本发明提供的卫星激光测高数据辅助提取高程控制点方法示意图。具体实施方式一种卫星激光测高数据辅助提取高程控制点方法，利用遥感目标定位系统差的特点，将高程控制点的提取分成dsm高程和高程系统差两部分实现，将激光测高目标与dsm的空间配准通过实测波形与搜索区dsm模拟观测波形的分组匹配技术实现，如图3所示，包括以下步骤：步骤1.高分立体影像或insar定向参数精化：对高分立体影像或insar数据进行区域网平差或自由网平差，获取初步精化后的定向参数，用于后继步骤中高分立体影像或insar数据的处理；步骤2.高分立体影像或insar生成dsm：利用定向参数，通过测区高分立体影像的密集匹配，或通过insar干涉处理，生成带有地理坐标的dsm；步骤3.激光测高数据实测组划分：将测区激光测高数据按轨道进行分组，在激光数据所在轨道区间过长情况下，将同轨测高数据分为若干组；依次选取每一组实测数据并执行步骤4—步骤10，实现对dsm搜索目标的模拟观测与二次匹配，获取激光测高数据与dsm数据空间目标的平面配准位置、光斑区dsm的高程系统差和高程控制点；步骤4.实测组激光测高波形提取：提取实测组中各个激光测高实测波形，设去除错误和粗差数据后实测组中有n个测高波形数据，对这n个实测波形分别进行激光测距时间因仪器、对流层、潮汐、大气等影响的修正，并进行归一化处理；步骤5.模拟观测参数设计：根据定向参数计算实测组内各激光测高光斑中心的理论位置，作为dsm上光斑配准位置的搜索中心，并设计dsm搜索区域；以实测组中各测站实际观测参数为准，对各测站统一设计模拟激光测高观测定向参数增量区间和步长，使模拟观测能够对所有设计的dsm搜索区进行观测；具体通过步骤5.1—5.3实现：5.1设计搜索区间和搜索步长：根据光斑与dsm间的平面相对精度和dsm分辨率，设计模拟观测光斑中心的搜索区间和搜索步长；5.2模拟观测参数确定：通过设计模拟激光测高观测定向参数变化区间和步长实现对dsm目标的搜索，可以以相同的增量同时改变实测组各测站观测姿态、或以相同的增量同时改变实测组各测站轨道位置观测参数、或以相同的增量同时改变各测站轨道和姿态参数，实测组其它观测参数不变，来设置模拟观测参数，实现dsm上模拟观测光斑中心位置相对与实测组理论位置的平面偏移、以及模拟观测dsm目标的改变；5.3模拟观测参数反算：根据dsm上模拟观测光斑中心的搜索区间和搜索步长，计算各模拟组观测参数相对于实测组观测参数变化的区间和步长。其中，搜索区间可以是以定向参数初值计算的理论光斑中心为中心的圆形或正方形区域。当搜索区间为正方形，且通过改变姿态来改变搜索区间时，通过步骤(5.4)—(5.6)实现：(5.4)光斑中心搜索范围和步长设计：根据光斑与dsm间的相对定位精度和dsm分辨率，设计模拟观测光斑中心的正方形搜索区间和搜索步长，其边长设为l米，搜索步长设为s米；(5.5)模拟观测姿态变化区间和步长计算：通过改变激光测高传感器侧滚角和俯仰角来实现搜索目标的改变，侧滚角和俯仰角的改变步长为其中hs为卫星相对于地面的高度，针对不同的搜索目标，侧滚角和俯仰角可配置值的数目均设计为m＝m0×2+1，其中int.[]表示对[]内的值取整；(5.6)模拟观测姿态计算：设实测组中第i(i∈[1,n]，n为当前实测组激光测高数)个测高数据的姿态侧滚和俯仰角分别为ωi和以实测组中的每个测站传感器姿态为基准，以不同的姿态侧滚角和俯仰角改变量设置模拟观测参数：整形参数k1依次在区间[-m0,m0]内取值，配置第i个模拟观测姿态侧滚角为ωi+k1δ；并针对每一个设置好的k1值，整数参数k2依次从区间[-m0,m0]内取值，配置第i个模拟观测姿态俯仰角为依次将ωi+k1δ和分别替代实测组第i测站的侧滚角和俯仰角，其余参数不变，作为模拟观测条件，则可获得模拟观测组总数为m＝m×m。通过改变卫星轨道、或通过同时改变卫星轨道和姿态实现对dsm搜索目标的改变，可参考上述的方案实现。采用以理论波束中心为中心的圆形搜索范围等方案，可以在矩形搜索的基础上，对搜索范围超过圆半径的模拟观测直接放弃。步骤6.模拟观测获取各模拟观测组波形：利用设计的模拟观测参数对dsm目标进行激光测高模拟观测，每测站模拟观测数目为m，可以获取m组模拟观测，每组模拟观测数目为n个；模拟观测时按dsm分辨率大小对地面光斑区进行空间分割，根据分割单元在光斑区的位置计算每一分割光斑单元的能量及其分布，卷积计算分割单元的能量被dsm上对应特定高程的模拟观测目标反射形成的子回波能量及其分布，模拟观测各采样时刻子回波能量的和形成模拟观测波形；对各模拟观测组中所有波形分别进行归一化处理用于后继步骤中的波形匹配；步骤7.实测组分别与各模拟观测组整体匹配：依次选取模拟观测组，将实测组中所有归一化波形整体当作匹配模板，沿时间轴上整体平移实测组各测站归一化波形，计算并记录下实测组分别与m个模拟观测组达到最佳匹配时的m个测度值，完成第一次波形匹配；本步骤中在时间轴上整体平移实测组各测站归一化波形后进行匹配，是为了消除模拟观测组目标与实际观测目标高程上存在的系统差。在dsm高程系统差较小或可忽略的情况下，实测组波形与模拟观测组波形的整体匹配可以不在时间轴上整体平移实测组归一化波形实现，而是直接计算两者的匹配测度值，用于后继最优匹配组选择；本步骤中实测组与模拟观测组波形匹配以波形对应点间的距离平方和、或相关模板匹配算法为匹配测度，计算最佳匹配时的测度值，为了提升一些质量不理想dsm模拟观测波形和实测波形间的匹配效果，匹配时将所有波形脉冲宽度以固定的倍率沿时间轴放大，即以波形高斯均值点为中心，将波形两边分别沿时间轴两方向延展加宽，增加实测波形与模拟观测波形间的重叠度，减小dsm内部非系统差对匹配的影响；步骤8.光斑区与dsm平面位置匹配：从m个匹配测度值中选取最优匹配测度值，其对应的模拟观测组称为最优模拟观测组，其模拟观测时光斑中心位置及光斑范围，即为光斑区与dsm的平面配准位置，该组模拟观测时光斑中心位置与光斑中心理论位置差异即为配准前后的光斑平面位置偏移量；步骤9.实测组与最优模拟观测组各对应波形分别匹配：分别将实测组中各测站的归一化波形作为参考模板，通过在时间轴上平移参考波形模板，获取实测组各测站归一化波形与最优模拟观测组中对应测站的归一化波形最优匹配时参考波形在时间轴上的偏移量，完成第二次波形匹配；步骤10.各光斑区dsm高程系统差获取：根据时间偏移量和激光传播速度计算最优模拟观测组各光斑区dsm高程的系统差，根据dsm上三维点坐标和dsm高程系统差，精化并提取dsm上特征点作为高程控制点；步骤11.影像高程控制点的提取：根据光斑区dsm上的三维坐标点和生成dsm时高分立体或insar的定向参数，利用影像的几何模型计算该三维点在影像上对应的像点坐标，同时根据dsm上的三维点坐标和该光斑区dsm的高程系统差，计算该三维点修正后的高程值，则影像像点坐标和修正后的高程值组成一对影像高程控制点；步骤12.根据平面和高程匹配结果，直接对dsm产品进行精化，或实现同轨高分传感器与激光测高仪的相对标定，或对遥感影像数据与激光测高数据进一步联合处理；其中，对dsm产品进行精化，根据匹配得到的各激光光斑区dsm高程的系统差，拟合整个测区dsm的高程误差，对dsm高程进行精化处理实现；或根据dsm与激光光斑目标间的高程系统差和提取的高程控制点，精化高分立体影像或insar处理参数；其中，同轨高分传感器与激光测高传感器的相对标定，利用同卫星平台同时期轨道姿态对影像定位和激光测高定位平面位置误差影响相同的特点，根据步骤8中获得的激光光斑中心理论位置和配准位置的相对偏移量以及卫星高度，反算激光测高传感器与高分传感器间的相对姿态修正值；其中，遥感数据与激光测高数据进一步联合处理通过遥感数据与激光测高数据或激光测高数据中提取的高程信息联合平差处理实现。本发明中步骤7—步骤10中，组内波形二次匹配也可以转化为一次匹配完成，即完成分组和模拟观测后，仍以各模拟观测组和实测组为对象，通过在时间轴上分别滑动组内各实测波形，实现组内各实测波形与模拟观测对应波形的最优匹配，并记录下最优匹配时波形匹配测度值和实测波形在时间维上滑动的距离，计算组内各波形最优匹配测度值均值s1，计算组内各实测波形在时间维上滑动的距离均值及各滑动距离与均值间差值的平方和s2，选取最优匹配测度均值s1和最小s2所在的模拟观测组作为最优模拟观测组，或综合比较s1和s2选取最优模拟观测组，最优模拟观测组光斑位置即为光斑与dsm配准位置，并利用最优模拟观测组中各实测波形与模拟波形最优匹配时在时间轴上滑动的距离即时间差计算对应光斑区高程系统差。本发明提供的激光测高数据与dsm模拟观测数据的波形分组匹配技术，与生成dsm的源数据没有关系，不仅能够用于激光测高目标与高分立体影像生成的dsm或insar生成的dsm空间配准，也能实现激光测高目标与其它方法生成的dsm、dem、或点云间的空间配准，获取激光测高光斑区配准位置和光斑区dsm、dem、或点云高程局部系统差。当前第1页12