[发明专利]一种基于数字表面模型的卫星激光高度计波形分解的方法

说明书

技术领域

本发明属于遥感对地观测技术领域，特别涉及一种基于数字表面模型的卫星激光高度计波形分解的方法。

背景技术

卫星激光高度计发射的激光脉冲经过大气到达地面，与地物作用后反射回去的能量被接收器接收，即为原始激光回波波形，因卫星激光高度计对激光发射脉冲和地物的后向散射脉冲都进行小间隔采样，因此几乎能记录整个脉冲波形。对接收器记录的原始激光回波波形经过平滑、去噪和大气修正后得到观测波形，观测波形用于波形分解。卫星激光高度计的激光脉冲传播到地面具有一定的发散角，因此激光脉冲在地面照亮的是具有一定面积的光斑。在激光光斑内的地形通常不会是一个简单的平面，激光高度计的回波是对光斑内各点的反射信号按时间先后顺序记录的地物反射率对应的光电压值。光电压值和多个因素相关，传统的波形分解算法基本不利用地面信息，对观测波形直接分解，分解得到的回波波形中只有高度信息，没有与地物的对应关系。本方法利用DSM信息来辅助观测波形分解，首先要利用激光指向数据与激光足印相机数据确定与激光光斑范围对应的DSM数据，之后利用DSM提取相关信息，辅助波形分解，该方法有助于提高波形分解的准确性，分解的波形直接与地物相匹配，从而计算出对应地物的高程。

发明内容

本发明所要解决的问题是：提供一种基于DSM的卫星激光高度计波形分解的方法。

本发明包括以下步骤：

步骤1，对激光高度计的观测波形进行波形分解，得到多个地物回波波形的参数；

步骤2，利用激光光斑内的DSM提取地物斑块，计算各地物斑块的坡度和面积；

步骤3，以各地物斑块到激光高度计的距离、斑块面积和各地物反射率为输入，计算各地物斑块的回波波形参数；

步骤4，将观测波形分解得到的各地物波形参数与计算得到的各地物斑块波形参数进行欧式距离计算，确定回波与地物斑块的对应关系；

步骤5，计算各地物斑块的高程。

所述步骤2还包括以下步骤：

步骤21，利用激光指向数据和足印相机影像数据获取激光光斑内DSM数据；

步骤22，根据垂直方向上的点集拓扑关系以及RANSAC平面提取算法，提取出DSM中多个地物斑块；

步骤23，根据DSM的坐标计算各地物斑块的面积；

步骤24，利用地物斑块的法向量与激光入射方向计算各地物斑块的坡度，坡度在一定阈值内认为其对回波波形无展宽影响；

步骤25，将识别出的斑块按DSM坐标确定的高度进行排序a₁…a_i…a_n，之后首先找到激光光斑内高程最高的地物斑块a₁，以其为基准斑块，以激光高度计设计的高程分辨率为高差判断依据，依次判断后续斑块与基准斑块的垂直距离，直到垂直距离大于激光高度计的高程分辨率，确定该地物斑块a_i，a_i-1之前的地物斑块合并为一个高度a₁，a_i成为其余斑块中高程最高的斑块，递归执行上述过程，得到合并后的各地物斑块并计算各合并斑块的面积。

所述步骤4还包括以下步骤：

步骤41，将观测波形分解的多个地物回波波形参数用向量表示；

步骤42，观测波形分解得到的各地物回波波形参数与计算得到的各地物斑块回波波形参数进行欧式距离计算；

步骤43，欧式距离最小时认为观测波形分解的地物回波与该地物斑块是最匹配的。

所述步骤5还包括以下步骤：

步骤51，根据观测波形分解的地物回波波形对应的时间计算激光高度计到地物斑块的距离；

步骤52，根据步骤51中的距离计算地物斑块的高程。

本发明的优点在于：本发明的方法将DSM信息用于激光高度计波形的分解，得到光斑范围内的回波波形与地物的匹配关系和地物的高程。

附图说明

图1是一种基于数字表面模型的卫星激光高度计波形分解的方法的流程图。

具体实施方式

以下结合附图对本发明的技术方案做进一步详细描述。

如图1，一种基于DSM的卫星激光高度计波形分解的方法，包括以下步骤：