技术领域
[0001] 本
发明属于激光雷达点云数据构建技术领域，具体涉及一种基于机载激光雷达点云数据的高精DEM构建方法。
背景技术
[0002] 目前，数字高程模型(DEM)的获取方法主要包括传统地面测量法、地形图等高线扫描数字化、航空摄影测量法、星载
合成孔径雷达干涉(InSAR)和机载激光雷达技术等。其中，激光雷达是一种主动式遥感测量设备，具有作业周期短、
精度高、成本低等优点，是当前获取高精度DEM的一种新技术。
[0003] 国内主要使用大型机载激光雷达系统获取DEM，但激光雷达系统对
飞行器的要求很高，且价格昂贵，因此应用领域较狭窄。近年来，随着无人机技术的迅速发展，轻巧型激光雷达系统可搭载在无人
旋翼机上，广泛应用于地理信息采集等领域；虽然激光雷达系统可直接测量三维坐标，但点云数据的离散特性不能连续表征地表信息。
发明内容
[0004] 基于上述
现有技术的不足，本发明提出了一种基于机载激光雷达点云数据的高精DEM构建方法。
[0005] 本发明是通过如下技术方案来实现的。
[0006] 一种基于机载激光雷达点云数据的高精DEM构建方法，所述方法为：
[0007] 1)首先进行点云数据预处理，参考现行GBT 24356-2009《测绘成果
质量检查与验收》及相关规范，在满足误差允许的情况下实现点云数据质量检查，避免出现地物扭曲、接边匹配误差影响后续
数据处理精度与后续产品的质量。
[0008] 2)然后采用改进的三
角网
迭代加密滤波法，将每个待判定点赋予一个与距离相关的优先级，保证优先级最高(距三角形
顶点最近)的激光点优先判断，用于地面点云与非地面点云的分离；三角网迭代加密滤波在迭代加密时，采取随机选取待定点的方式，若待定点满足加密要求，则待定点判定为地面点，并更新TIN，再判定下一个待定点；在地形起伏时，低矮植被可能会被误判断为地面点，地面点可能会被误判为非地面点。因此，采用改进三角网迭代加密滤波，将每个待判定点赋予一个与距离相关的优先级，保证优先级最高的激光点优先判断，分离地面点云和非地面点云；
[0009] 3)再使用已知高程点，构建邻近三角形并进行线性内插处理；高程可以表示地表垂直方向上的形态变化，等高线是常用的高程表示模型，表示为一条条互不交叉的闭合线段，通过连续的曲面函数来模拟地表形态；当使用有序数列对地表进行数字化模拟时，此高程模型称为DEM。
[0010] 较佳地，本发明通过不规则三角网(TriangulatedIrregularNetwork，TIN)模型将点云数据构建成DEM；使用Delaunay三角剖分
算法拟合点云数据，首先形成初始三角形，然后寻找最邻近点构成共边图形，直到构建的图形将整个区域
覆盖为止。该方法有三个特点，包括不论从哪三个点开始构建图形网，得到的模拟面都是相同的；每一个图形的外接圆都不含有其他的点；最小角特性，即共边三角形所形成的不含有优角的四边形，两个三角形中最小的两个值不会改变；其具体定义是在平面上有不全共线的两点欧氏距离和两点构成的Voronoi图，比到任意三都近时，此时的Voronoi图即为一个点的
内核，连接所有相邻内核；
[0011] 4)最后采用逻辑分析法与高程误差分析法对DEM精度进行评价，实现基于机载激光雷达点云数据的高精DEM构建。
[0012] 较佳地，所述点云数据大地定向要通过多个
坐标系间的几何关系来确定，包括激光
扫描仪的参考坐标系、瞬时
激光束坐标系、空中平台坐标系、目标所在
水平参考坐标系与垂直坐标系和惯性导航仪的参考坐标系与WGS-84坐标系。
[0013] 较佳地，不规则三角网模型构建完成后，对所述不规则三角网模型进行线性内插处理，线性内插对于点云庞大的数据量来说，可以极大的提高高程插值效率，减少插值函数运行时出现的误差；线性内插处理方法如下：
[0014] 1)在确定目标点高程时，先检索目标点所在的三角形，确定各点的坐标(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)，并得到如下目标点方程：
[0015]
[0016] 2)令两个点与另一顶点间的坐标差值为(x21,y21,z21)、(x31,y31,z31)，则有：
[0017] x21＝x2-x1,x31＝x3-x1,
[0018] y21＝y2-y1,y31＝y3-y1,
[0019] z21＝z2-z1,z31＝z3-z1；
[0020] 进而有目标点的高程h为：
[0021]
[0022] 较佳地，所述逻辑分析法的评价指标包括：
[0023] ①等高线是否存在逻辑错误：即是否存在与等高线特性相悖的情况(例如其的不相交特性，一条高线不能相交、等高线不能打折，等高线必须闭合等，不能有不合理的伪
节点和悬挂点等)；
[0024] ②坡度异常值：对DEM生成的坡度图进行检查，检查坡度异常高的值，对坡度大于50°的地区进行详细检查；
[0025] ③水系与等高线的逻辑关系：指
叠加等高线层与河流层时，河流的
位置在等高线的外凸连线上，且等高线不能两次跨越河流。
[0026] 较佳地，高程误差分析法的评价指标包括：
[0027] ①中误差(即均方根误差，RootMeanSquareError，RMSE)，是指通过统计的方法，客观评价DEM高程的偏移量；假设检验点的高程为Rk(k＝1,2,3……)，与检验点有相同空间坐标的底面点高程为Zk，则DEM的中误差为：
[0028]
[0029] ②平均绝对误差(MeanAbsoluteError，MAE)，是指模拟值与实际高程偏差的绝对值的平均，其公式为：
[0030]
[0031] ③平均误差(MeanRelativeError，MRE)，是指实测高程与DEM高程的误差的平均值，其公式为：
[0032]
[0033] 这三类指标可以衡量DEM产品的整体质量水平，其数值与DEM精度呈负相关，即三类指标的数值越小，DEM的精度越高。
[0034] 机载激光雷达系统是先进的主动式空间遥测系统，其配有全球
定位系统、惯性
导航系统以及激光扫描仪等先进设备，并且在地面还有众多为其提供数据校正的参考站与移动网络，机载激光雷达系统在各个领域中都取得了突破性进展。
[0035] 本发明首次提出采用改进三角网迭代加密滤波法获得地面点云，再基于
随机森林回归法内插DEM提供精
密度，最后通过逻辑分析法与高程误差分析法，对LiDAR点云数据生成的DEM数据进行了精度评价。点云数据作为新一代空间信息数据，极高的精度为各类生产生活提供了准确可靠的副产品，突破了传统摄影测量的
瓶颈。加之社会进步，对精细化空间信息的需求日益增多，也推动LiDAR发展。
附图说明
[0037] 图2为本发明使用Delaunay三角剖分算法拟合点云数据图；
[0038] 图3为机载激光雷达定位的示意图。
具体实施方式
[0040] 如图1所示，一种基于机载激光雷达点云数据的高精DEM构建方法，所述方法为：
[0041] 1)首先进行点云数据预处理，为点云数据质量检查和点云数据大地定向；参考现行GBT 24356-2009《测绘成果质量检查与验收》及相关规范，在满足误差允许的情况下实现点云数据质量检查，避免出现地物扭曲、接边匹配误差影像后续数据处理精度与后续产品的质量。
[0042] 2)然后采用改进的三角网迭代加密滤波法，将每个待判定点赋予一个与距离相关的优先级，保证优先级最高(距三角形顶点最近)的激光点优先判断，用于地面点云与非地面点云的分离；三角网迭代加密滤波在迭代加密时，采取随机选取待定点的方式，若待定点满足加密要求，则待定点判定为地面点，并更新TIN，再判定下一个待定点；在地形起伏时，低矮植被可能会被误判断为地面点，地面点可能会被误判为非地面点。因此，采用改进三角网迭代加密滤波，将每个待判定点赋予一个与距离相关的优先级，保证优先级最高的激光点优先判断，分离地面点云和非地面点云；
[0043] 3)再使用已知高程点，构建邻近三角形并进行线性内插处理；高程可以表示地表垂直方向上的形态变化，等高线是常用的高程表示模型，表示为一条条互不交叉的闭合线段，通过连续的曲面函数来模拟地表形态；当使用有序数列对地表进行数字化模拟时，此高程模型称为DEM。
[0044] 本发明通过不规则三角网(TriangulatedIrregularNetwork，TIN)模型将点云数据构建成DEM；使用Delaunay三角剖分算法拟合点云数据，首先形成初始三角形，然后寻找最邻近点构成共边图形，直到构建的图形将整个区域覆盖为止。该方法有三个特点，包括不论从哪三个点开始构建图形网，得到的模拟面都是相同的；每一个图形的外接圆都不含有其他的点；最小角特性，即共边三角形所形成的不含有优角的四边形，两个三角形中最小的两个值不会改变；其具体定义是在平面上有不全共线的两点欧氏距离和两点构成的Voronoi图，比到任意三都近时，此时的Voronoi图即为一个点的内核，连接所有相邻内核，如图2所示；
[0045] 4)最后采用逻辑分析法与高程误差分析法对DEM精度进行评价，实现基于机载激光雷达点云数据的高精DEM构建。
[0046] 具体地，所述点云数据大地定向要通过多个坐标系间的几何关系来确定，包括激光扫描仪的参考坐标系、瞬时激光束坐标系、空中平台坐标系、目标所在水平参考坐标系与垂直坐标系和惯性导航仪的参考坐标系与WGS-84坐标系。在任意时间内，激光脉冲形成位置为(xSL,ySL,zSL)T的激光点，同时测得此时激光扫描仪与激光点的斜距为ρ，则有[0047] (xSL,ySL,zSL)T＝(0,0,ρ)T
[0048] 激光束坐标系转换到激光扫描仪坐标系时，需要将其坐标绕x轴逆
时针旋转扫描夹角θi，通过转换矩阵Rl可以得知其在
传感器坐标下的位置是(xL,yL,zL)T，公式如下：
[0049] (xL,yL,zL)T＝Rl(xSL,ySL,zSL)T
[0050]
[0051] 激光扫描仪坐标系与惯性导航仪坐标系间存在着不可避免的偏心量激光点在惯性导航仪坐标系内的坐标如下，
[0052]
[0053] 其中，Rm为激光扫描仪与惯性导航仪间x、y、z方向上的安装误差角的转移矩阵乘积。
[0054] 同理，惯性导航仪与GPS天线之间的空间信息也存在差异 与三个方向上的偏差矩阵Rn，所以在当地水平坐标系下的激光点坐标为
[0055]
[0056] 当地的水平坐标系、垂直坐标系与WGS-84坐标系间还存在着偏差，前两者的偏差旋转矩阵为Rg，后两者的偏差旋转矩阵为Rw，结合天线在WGS-84坐标信息 则激光点在WGS-84坐标系下的坐标为
[0057]
[0058] 机载激光雷达定位的向量表示模型如下：
[0059] PWGS-84＝Rw·Rg·Rn(RmRl·r+tl-tg)+APCWGS-84
[0060] 其中，PWGS-84表示数据在地心坐标下的位置信息；r为瞬时激光坐标系下的位置向量；APCWGS-84表示GPS天线在大地坐标系内的位置向量。
[0061] 通过分析机载GPS与CORS获取的数据，可以得到瞬时的差分位置信息(DGPS)；IMU数据与DGPS数据联合解算，参照GPS天线获取的偏心分量，可以得到校正的位置信息；将其与机载LiDAR获取的原始数据相结合，通过计算修正偏差，就可以完成对原始激光数据的大地定向，使激光数据具有了WGS-84大地坐标信息，如图3所示。
[0062] 具体地，本发明不规则三角网模型构建完成后，对所述不规则三角网模型进行线性内插处理，线性内插对于点云庞大的数据量来说，可以极大的提高高程插值效率，减少插值函数运行时出现的误差；线性内插处理方法如下：
[0063] 1)在确定目标点高程时，先检索目标点所在的三角形，确定各点的坐标(x1,y1,z1)、(x2,y2,z2)、(x3,y3,z3)，并得到如下目标点方程：
[0064]
[0065] 2)令两个点与另一顶点间的坐标差值为(x21,y21,z21)、(x31,y31,z31)，则有：
[0066] x21＝x2-x1,x31＝x3-x1,
[0067] y21＝y2-y1,y31＝y3-y1,
[0068] z21＝z2-z1,z31＝z3-z1；
[0069] 进而有目标点的高程h为：
[0070]
[0071] 下面对本发明中DEM评价分析做进一步说明。
[0072] 行数字地形分析的前提是获取精确的数字高程模型，然而由于地形起伏变化与获取过程中的系统误差等因素，数字高程模型的数据质量在空间上呈不均匀分布。因此，需要对整个区域的数字高程模型精度进行评价和分析，较为准确的分析出DEM高程与实际高程之间的误差，并以此来评价数字高程模型的数据质量。已有很多的学者对于空间数据质量的评价体系进行了研究，如空间信息的误差、信息表达的完整性与空间
分辨率等，并且也是评价数字高程模型精度必不可少的指标，为更准确地模拟与分析各种地学过程提供
基础数据。数字高程模型对地表形态与结构的模拟精度非常重要，决定了数字地形分析结果的准确性。为了准确的进行DEM精度评价，不仅需要高程差异的统计，还需要等高线分析、水系分析与地形
可视化分析等辅助分析手段。
[0073] 在
地貌更新编辑完成后，需要进行质量检查，包括图形检查、属性检查和等高线检查，从而发现等高线是否相交，高程点与等高线矛盾，高程点是否落入水面或房屋面，高程点和等高线的属性矛盾等不合理错误，并对应进行
修改，使等高线光滑美观，使高程点达到高精度和高质量，从而满足标准要求，实现数据的高精度和高质。
[0074] 高程误差统计可以对DEM的精度进行估算，但受限于实地
采样的高成本，数量较少的地面实测高程数据不能很好的反应高程误差在空间上的分布情况。为此可视化误差分析正被广泛使用，其可以更形象与直观的定性评价DEM精度与误差的空间分布情况。但可视化误差分析的理论体系还不完善，现阶段主要以DEM反演等高线与地图特征法为主，但这类方法都没有与数字高程模型的
属性信息相结合。在应用领域中，常使用等高线对DEM整体的精度情况进行评价，等高线由DEM数据内插得到，此方法可以定性的评价DEM整体精度。同时，结合实地测量的高程数据，进行点对点的高程误差计算，通过各类高程误差统计模型进行精度评价。
[0075] 进一步地，本发明采用逻辑分析与高程误差分析方法对DEM精度进行评价，具体如下：
[0076] 所述逻辑分析法的评价指标包括：
[0077] ①等高线是否存在逻辑错误：即是否存在与等高线特性相悖的情况(例如其的不相交特性，一条高线不能相交、等高线不能打折，等高线必须闭合等，不能有不合理的伪节点和悬挂点等)；
[0078] ②坡度异常值：对DEM生成的坡度图进行检查，检查坡度异常高的值，对坡度大于50°的地区进行详细检查；
[0079] ③水系与等高线的逻辑关系：指叠加等高线层与河流层时，河流的位置在等高线的外凸连线上，且等高线不能两次跨越河流；
[0080] 所述高程误差分析法的评价指标包括：
[0081] ①中误差(即均方根误差，RootMeanSquareError，RMSE)，是指通过统计的方法，客观评价DEM高程的偏移量；假设检验点的高程为Rk(k＝1,2,3......)，与检验点有相同空间坐标的底面点高程为Zk，则DEM的中误差为：
[0082]
[0083] ②平均绝对误差(MeanAbsoluteError，MAE)，是指模拟值与实际高程偏差的绝对值的平均，其公式为：
[0084]
[0085] ③平均误差(MeanRelativeError，MRE)，是指实测高程与DEM高程的误差的平均值，其公式为：
[0086]
[0087] 这三类指标可以衡量DEM产品的整体质量水平，其数值与DEM精度呈负相关，即三类指标的数值越小，DEM的精度越高。
[0088] 以上所揭露的为本发明的优选实施例，不能以此来限定本发明之权利范围，因此依本发明
申请专利范围所作的等同变化，仍属本发明所涵盖的范围。
