一种结合立体像对的激光高度计高程控制点生成方法
【技术领域】
[0001]本发明属于遥感对地观测技术领域，具体涉及一种结合立体像对的激光高度计高程控制点生成方法。
【背景技术】
[0002]为了提高立体测图卫星的精度，亟需研究和解决双线阵点云生产过程中所面临的一系列关键技术问题，其中一项重要内容就是使用国产激光高度计提供精确的广义高程控制点，作为立体测图卫星的高程控制点，从而提高双线阵图像的测图精度。
[0003]要实现国产激光高度计在不同地形之下的每一次都能够得到正确的波形分解结果，并最终准确得到每一个像点/地物所对应的高程，对足印相机与双线阵图像的匹配技术、双线阵图像点云生产技术、波形分解技术、波形分解结果与地面点对应技术等均提出了较高的要求。由于卫星姿态测量的误差，激光出射点的测量误差，直接得到的激光指向存在较大误差，需要通过足印相机图像与双线阵图像进行匹配得到激光的精确指向，以及其光斑内对应的地物。并且由于国产卫星激光高度计只有一个全波形数据，无法知道光斑内每一个点的高程，需要在辅助数据的支持下确定光斑内每一个像点的高程，从而得到广义高程控制点。
[0004]由全波形的激光回波数据生成控制点的技术至关重要，若直接使用波形分解方法对激光回波进行分解，只能知道每一种高程在光斑内部所占的比重，无法知道每一个回波具体针对的是哪一个面片/像点。甚至可能有多个面片的回波耦合到同一个分解后的波形信息当中，无法进行区分。
[0005]传统获取高程控制点的方法为水准测量的方法，虽然水准测量方法能够比较高精度的获取高程控制点，但是人工测量的方法不能满足大面积测图的需要，另外人工测图需要花费大量人力、物力和财力，同时实地测量的方法不能在全球范围内获取高程控制点。由于激光雷达能够比较准确的获取地物的高程信息(星载激光雷达高程定位误差〈0.3m)，因此利用星载激光雷达技术能够在全球范围内进行广义高程控制点的获取，不仅可以节省大量人力物力，同时解决了全球范围内高程控制点提取的问题。
[0006]相比于小光斑激光雷达(光斑直径0.2m-2m)，大光斑激光雷达能够更大范围的获取地表反射数据，更加适用于大范围的观测，特别是在全球广义高程控制点提取过程中，大光斑激光雷达更有优势。同时由于大光斑激光雷达光斑覆盖地面面积较大，在地物比较复杂的区域波形也会比较复杂，造成波形分解难度增加。由于采用激光雷达同轴足印相机与立体相机影像配准方法获取激光光斑地面位置，因此平面定位精度由立体相机的定位精度决定。
【发明内容】
[0007]本发明所要解决的问题是:提供一种结合立体相机影像的激光高度计广义高程控制点生成方法。
[0008]结合立体像对的激光高度计高程控制点生成方法，包括以下步骤:
[0009](I)将立体相机影像与大光斑激光雷达同轴的足印相机影像进行配准，得到激光的指向信息;具体步骤如下:
[0010](1-1)对同轴相机影像进行重采样，使其分辨率与立体相机影像一致；
[0011](1-2)采用图像匹配算法进行影像之间的匹配，获取得到影像之间的对应关系；
[0012](1-3)通过最小二乘方法剔除匹配误差比较大的点获取较准确的影像匹配点；
[0013](1-4)通过影像匹配点获取影像匹配关系；
[0014](1-5)结合轨道的位置，根据匹配到的激光指向位置，计算得到激光的指向信息；
[0015](2)根据双线阵立体相机获取密集高程点云数据，具体步骤如下:
[0016](2-1)根据立体相机影像与大光斑激光雷达同轴的足印相机影像进行配准的结果，以匹配位置为中心，设置一个光斑大小的缓冲区；
[0017](2-2)对立体相机影像的前视和后视图像进行匹配，对包含于缓冲区内的任意一个像点，对前视和后视的图像进行匹配；
[0018](2-3)根据匹配得到的像对在各自扫描行上的瞬时成像姿态和相机参数，计算像点的粗略高程；
[0019](2-4)重复步骤(2-1 )-(2-3),得到激光光斑范围内每一个像点的高程；
[0020](3)通过对在足印内的点云数据按高程进行分类，得到光斑内不同高程地物类别数目；方法步骤如下:
[0021](3-1)根据各个像点的高程大小，将其按照高程的大小分为五等；
[0022](3-2)根据每一类所有像点的高程计算类中心的高程；
[0023](3-3)根据任意两类类中心高程的差异判断是否合并，若高程差小于0.1m，则对两类点云进行合并，若高程大于等于0.1m，则不合并；
[0024](3-4)根据每一类3-点云的高程的方差分布决定是否，若方差大于0.1m，则将该类拆分为两个新的类，若方差小于等于0.1m，则不进行拆分；
[0025](3-5)反复重复步骤(3-2)-(3-4)步，直到类别中心不发生变化为止；
[0026](4)通过类别数目指导进行大光斑激光雷达波形分解，通过波形分解得到每一高程类别地物准确的高程信息;方法步骤如下:
[0027](4-1)由于结合地类数目的高斯波形分解方法对噪声不敏感，因此只采用高斯低通滤波器对波形进行平滑，除去激光雷达波形中的毛刺现象；
[0028](4-2)通过形态学方法获取激光雷达波形的拐点，波峰和波谷的位置；
[0029](4-3)通过获取的波峰，波谷和拐点的位置，从激光雷达波形数据中获取组成激光雷达波形的高斯组分；
[0030](4-4)通过高斯波形的标准差和峰值大小对分解出来的高斯波形根据其重要进行排序，标准差越大，峰值越大的高斯波形越重要；
[0031](4-5)选取与点云按高程分类数目获取的类别数目相同的前几个重要的高斯波形作为激光雷达波形分解的波形；
[0032](5)通过点云计算每个高程类别信息高程的标准差，取标准差最小的一类作为参考提取广义高程控制点;方法步骤如下:
[0033](5-1)采用LM方法调整每一个高斯波形分量的均值和方差，使得选取的高斯波形能够最佳拟合激光雷达获取的波形，由此得到与点云数据得到的不同高程类别地物对应的高斯分量；
[0034](5-2)根据不同高程类别地物的高斯分量，计算出地物的高程；
[0035](5-3)将高程与点云所在区域的影像进行配准，得到影像上各个面片的高程信息，即为广义高程控制点。
[0036]本发明具有较高的控制点生成精度，能够生成多种形态的高程控制点信息，可满足国产激光高度计高精度广义高程控制点生产这一业务需求。
【附图说明】
[0037]图1是结合立体像对的激光高度计高程控制点生成方法流程图。
【具体实施方式】
[0038]本发明提供的技术方案采用如下步骤:
[0039]SI足印相机几何校