[发明专利]一种实时的数字表面模型生成方法

权利要求书

1.一种实时的数字表面模型生成方法，其特征在于：包括以下步骤：

步骤1：机载相机实时进行地表图像拍摄，并进行实时SLAM处理，得到当前帧的特征点以及构建的地图点云；并对当前帧特征点以及地图点云进行前期预处理：

根据当前帧中特征点的像素坐标，对当前帧进行二维Delaunay三角分割，得到图像平面的二维三角网格；利用地图点云中三维点与当前帧特征点的投影关系生成世界坐标系中以点云为顶点的三维三角网格，并舍弃三维三角网格高程信息，得到世界坐标系中水平面上的二维三角网格；分别利用几何关系对水平面上二维三角网格中的边缘点和非边缘点进行检测过滤，过滤掉边缘点和非边缘点中的噪点；

其中判断边缘点是否为噪点的方法为对每个边缘点进行如下两次判断：

对于某一边缘点，取其所处的每个三角形中的对边组成的折线，得到折线的两个端点，进一步得到两个端点与该边缘点的两条连线，判断这两条连线与所有不包含邻近点的三角形是否相交，如相交，则该边缘点为噪点；所述邻近点指该边缘点所处的每个三角形中的另外两个顶点；

对于某一边缘点，取其所处的每个三角形中的对边组成的折线，得到折线的两个端点，以两个端点连线为直径得到外接圆，取半径为k倍外接圆半径的同心圆，若该边缘点处于同心圆外，则该边缘点为噪点；

判断非边缘点是否为噪点的方法为：对于某一非边缘点，若该非边缘点处于其所处的每个三角形中对边组合成的闭合多边形的外部，则该非边缘点为噪点；

步骤2：生成当前帧的DSM纹理与规则网格：

利用步骤1处理后的当前帧特征点和地图点云，根据图像平面的二维三角网格与水平面上的二维三角网格的映射关系，通过图像平面的二维三角网格将当前帧的图像划分为一系列三角面片，并将划分得到的三角面片投影到水平面上的二维三角网格中，得到当前帧的DSM纹理；

用由若干个矩形子网格组成的矩阵主网格覆盖已添加DSM纹理后的水平面上二维三角网格，并根据水平面上二维三角网格中每个顶点在三维三角网格中的高程信息插值得到矩阵主网格中每个顶点的高程信息；

步骤3：分别进行DSM纹理以及规则网格融合：

步骤3.1：建立当前帧的权值图：对当前帧的水平面上二维三角网格中所有顶点坐标计算算术平均值，得到当前帧的中心坐标，将当前帧的中心处的权值设为255，将当前帧中离中心点最远点的权值设为0，得到权值与距离相关的变化梯度，其他坐标点的权值根据其与中心点的距离以及变化梯度计算得到；

步骤3.2：将步骤2得到的每个矩形子网格作为一个地图瓦片，并根据步骤2确定该地图瓦片在世界坐标系下的坐标，同时在该地图瓦片中添加DSM纹理，以及根据步骤3.1建立的权值图，在地图瓦片中添加该瓦片的权值；

步骤3.3：对于当前帧的每个地图瓦片，根据该地图瓦片在世界坐标系下的坐标，判断瓦片库中是否存在该地图瓦片，若不存在，则将该地图瓦片存入瓦片库中，若存在，则进一步比较当前帧中该地图瓦片与瓦片库中该地图瓦片的权值，以判断该地图瓦片是否为带有拼接线的地图瓦片，如果是带有拼接线的地图瓦片，则对该地图瓦片进行DSM纹理和网格的裂缝拼接，并用裂缝拼接后的该地图瓦片替换瓦片库中该地图瓦片，如果不是带有拼接线的地图瓦片，则比较当前帧中该地图瓦片与瓦片库中该地图瓦片的权值，如果当前帧中该地图瓦片的权值大于瓦片库中该地图瓦片的权值，则用当前帧中该地图瓦片替换瓦片库中的该地图瓦片。

2.根据权利要求1所述一种实时的数字表面模型生成方法，其特征在于：步骤1对水平面上二维三角网格中的边缘点和非边缘点进行检测过滤的过程中，先进行边缘点检测过滤，在进行一次边缘点检测过滤过程中，判断是否存在噪点，如果存在噪点，则在该次边缘点检测过滤过程完成后，更新当前帧特征点以及地图点云，再返回重新对当前帧特征点以及地图点云进行前期预处理；若不存在噪点，则进行非边缘点检测过滤；在一次非边缘点检测过滤过程中，判断是否存在噪点，如果存在噪点，则在该次非边缘点检测过滤过程完成后，更新当前帧特征点以及地图点云，再返回重新对当前帧特征点以及地图点云进行前期预处理。