一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法技术

本发明专利技术公开了一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，属于摄影测量技术领域，包括：多分区最小二乘平面拟合去噪；人工交互提取共有平面；多分区最小二乘平面拟合；对删除离群点后的所有平面点云进行最小中值平方法拟合平面，获取平面模型系数；基于同名平面点云的平面模型系数，利用间接平差法计算激光扫描点云与倾斜实景点云数据之间的转换参数，根据转换参数实现激光扫描点云与倾斜实景扫描点云之间的高精度融合；以平面到平面的距离的均方根误差和迭代次数作为约束条件进行迭代运算，完成点云数据拼接。本发明专利技术使用最小中值平方算法，无需设置过多参数就能够自适应剔除异常点对平面拟合过程的影响，具有较高的工程应用价值。高的工程应用价值。高的工程应用价值。

【技术实现步骤摘要】
一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法


[0001]本专利技术公开了一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，属于摄影测量


技术介绍

[0002]随着我国社会经济的飞速发展，基础设施建设也越来越广泛。在实景三维中国建设的大背景下，对于三维实景模型的需求达到前所未有的高度。桥梁作为交通基础设施的重要组成部分，因其结构复杂，单一的数据获取方式难以实现场景全覆盖，影响后续三维实景模型构建的完整性。因此，开展多源点云数据融合技术方法研究具有重要的意义和应用价值。
[0003]桥梁激光扫描与倾斜实景点云融合的本质就是通过旋转、平移操作将位于不同空间位置处的激光扫描点云和倾斜实景点云转换到同一位置。现有方法多采用ICP算法进行配准，ICP算法实现的基本思路是：对于激光扫描点云中的每一个点，在倾斜实景点云中找到与其欧式距离最近的点，组成匹配点集；根据匹配点集采用SVD分解算法求解旋转、平移矩阵，迭代运算直到所有匹配点对之间的欧式距离平方之和最小。由于道路场景下，桥梁基础设施结构复杂、数据量大直接导致基于点对点的配准方法效率不高；并且不同平台采集的点云数据在扫描精度和密度方面差异较大，不存在严格意义上的同名点对，配准误差较大。

技术实现思路

[0004]本专利技术公开了一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，以解决现有技术中，基于点对点的配准方法效率不高、配准误差较大的问题。
[0005]一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，包括：
[0006]S1.多分区最小二乘平面拟合去噪；
[0007]S1.1.人工交互提取共有平面；
[0008]S1.2.多分区最小二乘平面拟合；
[0009]S2.对删除离群点后的所有平面点云进行最小中值平方法拟合平面，获取平面模型系数；
[0010]S3.基于同名平面点云的平面模型系数，利用间接平差法计算激光扫描点云与倾斜实景点云数据之间的转换参数，根据转换参数实现激光扫描点云与倾斜实景扫描点云之间的高精度融合；
[0011]S4.以平面到平面的距离的均方根误差和迭代次数作为约束条件进行迭代运算，完成点云数据拼接。
[0012]优选地，所述S1包括：在每一个平面点云数据中进行空间多分区的最小二乘平面拟合，对离散点到分区拟合平面的距离进行累计与统计分析，将远离拟合平面的点作为离群点删除；
[0013]通过软件显示出的点云数据，将激光扫描点云和倾斜实景点云数据中同一平面特征目标提取出来。
[0014]优选地，所述S1.2包括：
[0015]对于点云数据中的任意一点P
i
(x
i
,y
i
,z
i
)，以P
i
为中心，r为半径进行球邻域搜索，如果该球形邻域内P
i
的近邻点个数小于3，则将P
i
点标记为离群点删除；如果近邻点个数大于3，则将P
i
其邻域内的点构建协方差矩阵M3×3：式中，m是点P
i
的临近点个数，m>3；是P
i
及其m个临近点的三维质心当i＝0时表示P
i
点本身；
[0016]求解该协方差最小特征值对应的归一化特征向量即为局部平面的单位法向量，进而可得P
i
到局部平面的距离d
i
为：
[0017]遍历计算每个点到局部平面距离的均值μ：
[0018]标准差σ：
[0019]当某一点P
i
到局部平面的距离在(μ
‑
σ,μ+σ)范围内时保留该点，不在该范围则定义为离群点删除。
[0020]优选地，所述S2包括：
[0021]S3.1.随机选取平面内三点，计算模型参数，从平面点云中随机选取3个点，通过这3个点求解平面方程；假设选取的3个点为p1(x1,y1,z1)，P2(x2,y2,z2)，P3(x3,y3,z3)，可得：
[0022][0023]平面的单位法向量n为：
[0024]已知平面的单位法向量和平面上的一点P1，可求得平面方程的系数；
[0025]S3.2.计算各点到S2.1中所得到的平面的距离dis，并统计距离的平方中值Ω，迭代计算k次；
[0026]S3.3.选取k次迭代中Ω的最小值所对应的拟合平面为最佳拟合平面，并根据式(1)、式(2)计算出每个数据对应的权值，通过权值自适应剔除异常点；
[0027]式中，σ为稳健标准差；n为点的个数；r
i
表示第i个点与拟合平面之间的距离；med表示距离的中值；
[0028]式中，w
i
表示第i个点的权重，权重为0则认为是异常点；
[0029]S3.4.在剔除异常点之后的点集中，计算每个点到最佳拟合平面的距离，将距离小于给定阈值的点作为内点；
[0030]S3.5.对内点采用最小二乘法重新计算平面模型的参数。
[0031]优选地，所述S3包括：
[0032]坐标转换模型的获取，基于同名平面点云的平面模型系数，构建间接平差模型，计算激光扫描点云与倾斜实景扫描点云之间的平移向量；
[0033]设同一平面在坐标系O1‑
X1Y1Z1中表示为P(a1,b1,c1,d1)，在坐标系O2‑
X2Y2Z2中表示为Q(a2,b2,c2,d2)，其中，a,b,c,d表示平面的描述参数；由仿射变换知，两平面之间满足如下关系：
[0034]式中：
[0035]其中α为绕X轴的旋转角度；β为绕Y轴的旋转角度；γ为绕Z轴的旋转角度；
[0036]将旋转矩阵线性化，分离九个参数转换为三个参数：
[0037][0038]分别对α,β,γ关于a2,b2,c2求偏导，最终的线性模型如下：V＝Bx
‑
l
[0039]其中：
[0040][0041]其中为a2的估计值，的估计值，的估计值，为b2的估计值，的估计值，的估计值，为c2的估计值，的估计值，的估计值，为d2的估计值，
[0042]按间接平差原理，可计算出变换矩阵中各参数改正值：
[0043]优选地，所述S4包括：
[0044]以平面到平面的距离的均方根误差和迭代次数作为约束条件，对于配准后的任一
同名平面(P
i
,Q
i
)，计算它们之间的欧式距离d
PQ
，则均方根误差rmse为：
[0045]当rmse小于设定阈值时则表示配准成功，同时给定迭代次数约束，避免在均方根误差不满足设定阈值的情况下，出现无限循环迭代；
[0046]用坐标转换参数构建旋转矩阵和平移向量：
[0047][0048]其中R为旋转矩阵，R(1,1)为旋转矩阵中的第一行第一列元素，R(1,2)为旋转矩阵中的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，其特征在于，包括：S1.多分区最小二乘平面拟合去噪；S1.1.人工交互提取共有平面；S1.2.多分区最小二乘平面拟合；S2.对删除离群点后的所有平面点云进行最小中值平方法拟合平面，获取平面模型系数；S3.基于同名平面点云的平面模型系数，利用间接平差法计算激光扫描点云与倾斜实景点云数据之间的转换参数，根据转换参数实现激光扫描点云与倾斜实景扫描点云之间的高精度融合；S4.以平面到平面的距离的均方根误差和迭代次数作为约束条件进行迭代运算，完成点云数据拼接。2.根据权利要求1所述的一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，其特征在于，所述S1包括：在每一个平面点云数据中进行空间多分区的最小二乘平面拟合，对离散点到分区拟合平面的距离进行累计与统计分析，将远离拟合平面的点作为离群点删除；通过软件显示出的点云数据，将激光扫描点云和倾斜实景点云数据中同一平面特征目标提取出来。3.根据权利要求2所述的一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，其特征在于，所述S1.2包括：对于点云数据中的任意一点P
i
(x
i
，y
i
，z
i
)，以P
i
为中心，r为半径进行球邻域搜索，如果该球形邻域内P
i
的近邻点个数小于3，则将P
i
点标记为离群点删除；如果近邻点个数大于3，则将P
i
其邻域内的点构建协方差矩阵M3×3：式中，m是点P
i
的临近点个数，m>3；是P
i
及其m个临近点的三维质心当i＝0时表示P
i
点本身；求解该协方差最小特征值对应的归一化特征向量即为局部平面的单位法向量，进而可得P
i
到局部平面的距离d
i
为：遍历计算每个点到局部平面距离的均值标准差σ：当某一点P
i
到局部平面的距离在(μ
‑
σ，μ+σ)范围内时保留该点，不在该范围则定义为离群点删除。4.根据权利要求3所述的一种桥梁基础设施激光扫描与倾斜实景点云数据拼接方法，其特征在于，所述S2包括：S3.1.随机选取平面内三点，计算模型参数，从平面点云中随机选取3个点，通过这3个点求解平面方程；假设选取的3个点为P1(x1，y1，z1)，P2(x2，y2，z2)，P3(x3，y3，z3)，可得：
平面的单位法向量n为：已知平面的单位法向量和平面上的一点P1，可求得平面方程的系数；S3.2.计算各点到S2.1中所得到的平面的距离dis，并统计距离的平方中值Ω，迭代计算k次；S3.3.选取k次迭代中Ω的最小值所对应的拟合平面为最佳拟合平面，并根据式(1)、式(2)计算出每个数据对应的权值，通过权值自适应剔除异常点；式中，σ为稳健标准差；n为点的个数；r
i
表示第i个点与拟合平面之间的距离；...

【专利技术属性】
技术研发人员：丛波日，辛公锋，李庆营，刘如飞，王飞，陈铮，刘宪明，龙关旭，李想，李子豪，齐麟，刘凯峰，张彤辉，
申请(专利权)人：山东科技大学山东高速集团有限公司创新研究院，
类型：发明
国别省市：