【技术实现步骤摘要】
本专利技术提供一种基于分区icp算法的地表模型误差分析方法。
技术介绍
1、近年来，以三维激光扫描技术(terrestrial laser scanning，tls)和无人机航空摄影技术(unmanned aerial vehicle，uav)为代表的非接触监测技术迅速发展。相较于传统的监测方式，非接触监测技术具有很大的优势，在效率、安全性以及作业成本方面对比全站仪等传统手段均有很大提升，在工程中得到了广泛的应用。在实际使用中，无人机建模具有很高的完整性，但由于缺乏足够的地面控制点，其形状总是会发生一定的扭曲，为了获得高精度的无人机建模结果，需要比较使用不同方法进行无人机建模的精度，对建模结果进行误差分析。
2、三维激光扫描所得到的点云完整性有限，但是具有很高的精度，因此该误差分析可以通过令无人机建模点云与三维激光扫描点云对比计算来进行。要计算点云所产生的误差是很困难的，因为点云没有除空间性质以外的任何特征，这就导致在分析两个不同点云的差异时，难以找到两个点云之间的对应部分，并正确测量其差距。目前常用的点云误差分析方法有shortest distance(sd)，dem ofdifference(dod)，multiscale model to modelcloud comparison(m3c2)等，但它们都有各自的局限性，不能很好的计算无人机建模扭曲的这种误差。
3、现有技术中，有最短距离法(sd)，选取一个误差计算点，计算该点到参考点云的最短距离，以该距离作为误差计算结果。但是许多情况下，计算所得的
4、还有垂直距离法(dod)，根据dem数据，以计算点处两期dem数据的差值作为该点的误差计算结果。这种方法只在误差产生的方向为垂直方向时比较适用，当误差产生方向包含水平分量时，计算结果会产生偏差。
5、另外还有m3c2算法，该算法首先需要选取核心点云，对于一个核心点，将其与领域指定半径r范围内的点云数据拟合成为一个平面，并沿该平面的法线方向作一个半径为r的圆柱体，分别计算圆柱体内误差计算点云和参考点云数据的平均位置，，并以平均位置之间的距离作为误差计算结果。m3c2算法在点云误差产生方向为沿法线方向时有很好的适用性，但在点云相对于与法线垂直方向出现较大错位时，其计算结果会出现很大的偏差。
技术实现思路
1、为了验证三维地表模型的精度，需要计算模型产生的误差，而这种三维空间中的误差是难以计算的，如何实现三维模型高精度的误差分析。
2、本专利技术提供一种基于分区icp算法的地表模型误差分析方法，该方法使用一种点云匹配算法：iterative closest point(icp)算法，包括以下步骤：
3、s1、对待进行误差分析的点云模型进行分块，令每个点云分块分别与参考点云进行点云匹配算法icp配准，得到一个空间变换矩阵；
4、s2、对点云分块进一步分块，继续进行icp配准，得到进一步的空间变换矩阵。当一次配准过程中点云分块的位移小于设定的阈值时，停止分块，叠加每次得出的空间变换矩阵即可得到点云分块质心处相对于参考点云的误差值，进而完成模型整体的误差分析。
5、具体的，s1使用点云匹配算法icp算法对点云分块与参考点云进行配准，具体包括以下步骤：
6、通过计算点云匹配过程中产生的旋转平移量，从而完成误差的计算。不同点云之间的匹配看作是其中一片点云通过旋转、平移与另一片点云相重合的过程，这种变化用一个旋转矩阵和一个平移向量来表示；通过icp算法，求出待匹配点云p到参考点云q的空间旋转矩阵m与空间平移向量
7、
8、
9、分别在待匹配点云p和参考点云q中，根据一定的约束条件，找到最邻近点(pi，qi)，计算对应点对之间的距离，通过最小化该距离不断迭代旋转矩阵与平移向量，如果对应点距离小于给定阈值设置，则算法结束。
10、每一次迭代过程得到一个空间旋转矩阵mn和一个空间平移向量为了得到原匹配点云p到第n次迭代后的pn的变换，将mn和统一为4×4的齐次矩阵rn和tn：
11、
12、则每一次迭代的空间旋转矩阵和空间平移向量合并为同一个空间变换矩阵mrtn：
13、
14、将每次迭代得到的mrtn相乘即得到所求原匹配点云p到第n次迭代后的pn的空间变换矩阵mrt。
15、
16、通过该空间变换矩阵即得到点云匹配过程中的位移量。
17、s2采用多次分块，叠加每次的计算结果，具体包括以下步骤：
18、使用icp算法使匹配点云与参考点云配准是改变点云整体的位置，配准过程中并不会改变点云的形状，当匹配点云与参考点云之间存在形状差异时，无法做到对点云各个部分进行单独配准，因此最终的配准结果只是一个整体的最优解，无法使得每个局部都精确配准。因此，需要使用分块的方式，使得各个局部分别精确配准，并且分块后仍然可能出现局部无法配准的问题，因此采用多次分块，当一次分块后进行icp配准过程中分块的位移量小于设定的阈值τ时，即认为该分块已经完全配准。此时，叠加每次配准过程中的空间变换矩阵，可得出点云分块从原位置到配准后位置的空间变换矩阵，并由此可以得出点云分块质心处的误差。
19、本专利技术提出了一种基于分区icp算法的误差分析方法，主要可用于以下问题：
20、三维点云模型产生的误差难以计算，传统方法的局限性较大，使用这些方法往往不能得到理想的结果。该技术提出的方法通过对点云分块并进行匹配，能够更好的找出不同点云之间对应点对的关系，精确地计算点云所产生的误差，即使在复杂情况下也有很强的适用性，极大地提高了点云误差分析的精度。
本文档来自技高网...【技术保护点】
1.一种基于分区ICP算法的地表模型误差分析方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于分区ICP算法的地表模型误差分析方法，其特征在于，S1使用点云匹配算法ICP算法对点云分块与参考点云进行配准，具体包括以下步骤：
3.根据权利要求1所述的一种基于分区ICP算法的地表模型误差分析方法，其特征在于，S2采用多次分块，叠加每次的计算结果，具体包括以下步骤：
【技术特征摘要】
1.一种基于分区icp算法的地表模型误差分析方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的一种基于分区icp算法的地表模型误差分析方法，其特征在于，s1使用点云匹配算法icp算法...
【专利技术属性】
技术研发人员：周家文，罗箱陇，蒋楠，李海波，陈俊霖，李从江，胡宇翔，张洁源，杨玉川，刘建，
申请(专利权)人：四川大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。