一种基于低秩恢复的隧道点云去噪方法技术

本发明专利技术公开了一种基于低秩恢复的隧道点云噪声去除的方法，该方法利用PCA对每个点进行法线计算，构建每个点的局部坐标系；利用双边滤波对法线进行平滑处理，旋转每个点的Z轴使其与法线重合；对每个点构造离散的二维描述子——高度图；对相似的隧道点云块，将其高度图组合为一个矩阵，通过低阶矩阵近似的求解方法得到去噪后的高度图；将去噪后的高度图映射为点的三维坐标得到去噪后的隧道点云，本发明专利技术的方法有效的在去除隧道点云噪声及由于附件结构产生的离群点的基础上并保留了隧道点云模型的局部细小特征，尤其对大噪声具有较强的鲁棒性。

【技术实现步骤摘要】
一种基于低秩恢复的隧道点云去噪方法
本专利技术涉及噪声处理
，具体涉及一种基于低秩恢复的隧道点云去噪方法。
技术介绍
近几年轨道交通的发展随着城市建设的扩展而并行发展。地铁尤其成为了拥有高容纳力的轨道交通系统的重要组成部分。因此有效实现地铁建造过程中的质量检测成为建设一个安全可靠轨道交通系统的必要条件。近几年随着数据测量技术的迅速发展，三维扫描仪采集到的隧道点云越来越密集，并且包含丰富的事物信息。因此现在对于地铁隧道的形变检测基本采用三维激光扫描技术来代替传统的手工检测方法。但即使使用非常高精度的扫描仪，获得的隧道点云总不可避免的受到环境、人为等因素的影响而含有噪声，从而降低了隧道点云质量，影响后续对隧道点云的处理，为了获得良好的扫描点云模型，进行后续的隧道安全质量评估需要实现隧道点云的去噪工作。现有一些点云去噪方法，如RIMLS、MRPCA、EAR、GPF、WLOP等都有可以较好的去除噪声并保留一定的特征，但是，它们都着眼于局部结构，忽略了联合整体点云中的相似结构进行去噪可以增加去噪可信度的这一特点。最近在图像处理领域提出的一些非局部相似点的方法已应用于图像去噪，主要是先识别出图像中目标面片的最相似的几个面片，然后组成一个矩阵，通过低秩矩阵近似的方法恢复出无噪声的矩阵，从而达到去噪的效果。但由于点云的无序性，点云不可以统一的用一个方式来表示点云的结构。另一方面，如何有效地处理低秩模型和从低秩近似中恢复几何信息也是在点云领域存在的难题。目前针对现有的三维扫描仪获取的带有噪声的隧道三维点云数据，如何实现统一方式表达隧道点云结构并恢复隧道几何信息，进而有效地降低噪声对隧道点云质量的影响，尚未提出有效地解决方案。因此，迫切需要一种高效的隧道点云去噪方法。
技术实现思路
(一)解决的技术问题本专利技术的目的在于提供一种基于低秩恢复的隧道点云去噪方法，以解决上述
技术介绍
中提出的实际问题。(二)技术方案为实现上述目的，本专利技术提供如下技术方案：一种基于低秩恢复的隧道点云去噪方法，包括以下步骤：步骤1：对输入隧道点云进行局部投影，生成高度图，并构建相似隧道点云块的高度图相似矩阵Gt；步骤2：利用低阶矩阵近似方法求解每个高度图相似矩阵，得到去噪后矩阵Xt；步骤3：将与一个点相关的每个网格的去噪后的值进行加权平均，得到此点的去噪后的z坐标，将去噪后的高度图映射回局部三维坐标，最后将局部坐标转化为全局坐标，去噪过程结束。进一步的，所述步骤1具体包括以下步骤：步骤1.1：以每个点为球心，半径为r构造一个球，在这个球内的所有点构成了球心的一个局部点云块，利用PCA对每个局部面片估计点云的法线，得到的三个特征向量构成每个点的局部坐标系，其中，法线方向为Z轴；步骤1.2：对每个局部点云块中的点的法线进行法线协方差分析，将得到的三个特征向量作为新的局部坐标系，其中，最大特征值对应的特征向量为Z轴；步骤1.3：对步骤1.1得到的法线进行双边滤波，得到较为连续规整的法线；步骤1.4：将每个局部点云块中的点的坐标转换为在局部坐标系中的坐标，即局部坐标，再构造以每个点为中心，法向为此点局部坐标系的Z轴的切平面，即局部坐标系中的XY平面；步骤1.5：将局部点云块中的点向其对应的切平面投影，将切平面划分为固定大小的网格，每个网格的值由局部点云块中距离网格中心最近的k个点到切平面的加权平均距离决定，因此，高度图确定完毕；步骤1.6：将每个高度图用高维向量表示，通过计算每两个向量间的欧氏距离来判别两个向量间的相似度，找出最相似的K个向量，组成高度图相似矩阵Gt；。进一步的，所述步骤1中利用PCA估计法线的计算方法如下：构造协方差矩阵：其中Pi＝{pj|||pi-pj||2＜＜r}，|Pi|是局部点云块Pi中点的个数；通过分解协方差矩阵，得到三个特征值和特征向量，对应最小特征值的特征向量为点的法线，点pi与其三个特征向量构成点pi的局部坐标系；所述对法线进行协方差分析的计算方法如下：其中，|Ni|为在局部点云块Pi中法线的个数，ni表示点pi的法线，nj表示点pj的法线点pi与其新得到的三个特征向量构成点pi的局部坐标系，最大特征值对应的特征向量作为新局部坐标系的Z轴；所述每个网格的值的计算方法如下：其中，bi是网格的中心，Ti是切平面，height(Ti，p′i)是点p′i到切平面的高度，是高斯权重函数，σd2是空间距离带宽。进一步的，所述步骤2的具体实现步骤如下：步骤2.1：构造优化方程其中，xi，xj为低秩图，ωij为代表两个高度图Hi与Hj之间的相似性的高斯权重，其中λ，β为权衡参数；步骤2.2：利用Ut主方向上的正交性约束，步骤2.1中的优化方程可写为其中，I为单位矩阵，Ut，Vt为去噪后矩阵进行UV分解后得到的分解矩阵；步骤2.3：采用交替最小化方法求解步骤2.2中的优化方程。进一步的，所述步骤3的具体实现步骤如下：步骤3.1：找出与点p′i相关的网格，通过如下公式计算所有相关网格的加权平均值：其中，p′i为去噪后的点云，bi是生成网格的中心，w(bi，p′i)为高斯权重函数，高斯权重函数，是恢复出的中心为bi的网格的高度值，height(T，p′i)是点p′i到切平面的高度，也就是p′i的新z坐标；步骤3.2：将局部坐标转化为全局坐标；步骤3.3：将来自不同局部点云块的同一点的全局坐标取平均，得到最终去除噪声并且保留良好细节特征的隧道点云。(三)有益效果1.本专利技术采用非局部多隧道点云块联合的作用的点云去除噪声的方法，设计了包含每个点的法线信息的具有旋转不变性的矩阵，可以直接作为低阶近似模型的输入；2.提出了一种多个隧道点云块联合作用的方案，能够更好地进行噪声去除，也能够恢复更多的特征；3.在估计每个点法线后，再对其进行双边滤波，使法线更加规整，从而使所建立的局部坐标系对噪声有更强的鲁棒性，保留更多的隧道几何信息。附图说明为了更清楚地说明本专利技术实施例的技术方案，下面将对实施例描述所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1本专利技术流程图；图2是带有噪声的隧道原始点云模型图；图3是带噪声的隧道原始点云模型高度图；图4是隧道点云模型去噪效果对比图；图4(a)是带有噪声的隧道模型；图4(b)本专利技术的方法得到的去除噪声的隧道点云效果图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本发本文档来自技高网...

【技术保护点】
1.一种基于低秩恢复的隧道点云去噪方法，其特征在于，包括以下步骤：/n步骤1：对输入隧道点云进行局部投影，生成高度图，并构建相似隧道点云块的高度图相似矩阵G

【技术特征摘要】
1.一种基于低秩恢复的隧道点云去噪方法，其特征在于，包括以下步骤：
步骤1：对输入隧道点云进行局部投影，生成高度图，并构建相似隧道点云块的高度图相似矩阵Gt；
步骤2：利用低阶矩阵近似方法求解每个高度图相似矩阵，得到去噪后矩阵Xt；
步骤3：将与一个点相关的每个网格的去噪后的值进行加权平均，得到此点的去噪后的z坐标，将去噪后的高度图映射回局部三维坐标，最后将局部坐标转化为全局坐标，去噪过程结束。


2.根据权利要求1所述的一种基于低秩恢复的隧道点云去噪方法，其特征在于，所述步骤1具体包括以下步骤：
步骤1.1：以每个点为球心，半径为r构造一个球，在这个球内的所有点构成了球心的一个局部点云块，利用PCA对每个局部面片估计点云的法线，得到的三个特征向量构成每个点的局部坐标系，其中，法线方向为Z轴；
步骤1.2：对每个局部点云块中的点的法线进行法线协方差分析，将得到的三个特征向量作为新的局部坐标系，其中，最大特征值对应的特征向量为Z轴；
步骤1.3：对步骤1.1得到的法线进行双边滤波，得到较为连续规整的法线；
步骤1.4：将每个局部点云块中的点的坐标转换为在局部坐标系中的坐标，即局部坐标，再构造以每个点为中心，法向为此点局部坐标系的Z轴的切平面，即局部坐标系中的XY平面；
步骤1.5：将局部点云块中的点向其对应的切平面投影，将切平面划分为固定大小的网格，每个网格的值由局部点云块中距离网格中心最近的k个点到切平面的加权平均距离决定，因此，高度图确定完毕；
步骤1.6：将每个高度图用高维向量表示，通过计算每两个向量间的欧氏距离来判别两个向量间的相似度，找出最相似的K个向量，组成高度图相似矩阵Gt。


3.根据权利要求2所述的一种基于低秩恢复的隧道点云去噪方法，其特征在于：所述步骤1中利用PCA估计法线的计算方法如下：
构造协方差矩阵：



其中Pi＝{pj|||pi-pj||2＜＜r}，|Pi|是局部点云块Pi中点的个数；
通过...

【专利技术属性】
技术研发人员：张沅，汪俊，董竟萱，刘树亚，鲁德宁，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32