一种基于双边滤波器的多视角密集点云数据融合方法技术

本发明专利技术公开了一种基于双边滤波器的多视角密集点云数据融合方法，属于光学三维非接触式测量技术领域，该方法的实现步骤包括：1)输入多视角密集点云数据；2)拓扑关系构建；3)基于双边滤波器的点数据归属；4)Mean-shift聚类融合；5)输出融合结果。本发明专利技术方法融合多视角点云数据时，引入双边滤波器和Mean-shift聚类，不需要识别重叠区域与非重叠区域，提高了融合多幅多视角点云数据的效率和融合后点云的平滑、光顺程度，有效地克服了现有点云融合技术的缺陷。

【技术实现步骤摘要】
一种基于双边滤波器的多视角密集点云数据融合方法
本专利技术属于光学三维非接触式测量
，涉及一种多视角密集点云数据融合方法，更进一步涉及一种新的基于双边滤波器的多视角密集点云数据融合方法。
技术介绍
光学三维测量技术是集光、机、电和计算机技术于一体的智能化、可视化的高新技术，主要用于对物体空间外形和结构的扫描，以得到物体的三维轮廓，获得物体表面点的三维空间坐标。随着现代检测技术的进步，特别是随着激光技术、计算机技术以及图像处理技术等高新技术的发展，三维测量技术逐步成为人们的研究重点。光学三维测量技术由于具有非接触、快速测量、精度高等优点，使其在航空航天、军工、汽车和装备制造等行业得到了广泛应用。三维测量技术是获取物体表面各点空间坐标的技术，主要包括接触式和非接触式两大类。其中，通过光学三维非接触式测量设备采集的多视角点云数据的处理技术，如融合、降噪、曲面重建等，是当前研究的热点。点云融合是指消除由测量误差和匹配误差等导致的多视角点云重叠区域的噪声、分层和冗余，建立细节特征清晰、表面光顺的单层点云模型。点云融合技术随着三维光学非接触式测量技术的发展而不断发展。按照融合方式可以将现有点云融合方法分为三类:基于隐式曲面重构的点云融合技术，基于显式曲面重构的点云融合技术和基于聚类的点云融合技术。基于隐式曲面重构的点云融合技术一般先通过有向距离场(SignedDistanceField，SDF)或者移动最小二乘(MovingLeastSquares，MLS)等技术构建一光顺的隐式曲面，然后通过三角网格化方法或者B样条曲面构建技术将隐式曲面转换为显式曲面，从而实现多视角密集点云数据的融合。这类点云融合方法在构建隐式曲面的过程中可以有效地消除由测量误差和匹配偏差等导致的分层及噪音，但这类方法在将隐式曲面转换为显式曲面的过程中会占用大量的计算机资源，融合效率低。基于显式曲面重构的点云融合技术在网格化的过程中(通过构建三角网格曲面、泊松曲面或者nurbs曲面等)将多视角点云融合在一起，融合结果为一显式的网格曲面。但这类方法对匹配偏差及噪声比较敏感，融合后重叠区域表面的融合效果一般比较差，会保留融合前的痕迹。此外，这类方法占用计算机资源多、效率低，不适用于大规模点云数据的融合处理。基于聚类的点云融合技术通过聚类的方法将重叠区域的同名点数据聚集在一起并融合，从而实现冗余消除。这类点云融合方法是一种无网格点云处理技术，此类方法占用计算机资源少，但当待融合的点云幅数多于两幅时，需要通过增量式的方法进行融合，即先将其中两幅融合在一起，然后再加进来一幅进行融合，依次类推，直至所有点云融合完毕为止。增量式的融合方式一方面效率低，另外一方面不能保证融合后的点云位于最优的曲面上。
技术实现思路
为了克服现有点云融合技术的缺陷，满足包含重叠区域的多视角密集点云数据的融合处理要求，本专利技术提供了一种基于双边滤波器的多视角密集点云数据融合方法。该方法可一次性将多幅包含匹配偏差、噪声的多视角密集点云数据融合为一幅完整、光顺、分布均匀的单层点云数据。为达到以上目的，本专利技术采取的技术方案是：一种基于双边滤波器的多视角密集点云数据融合方法，包括下述步骤：步骤一，输入多视角密集点云数据输入测量设备采集的待融合的多幅多视角密集点云数据，输入的多幅多视角密集点云数据需同时包含三维坐标信息及法向量信息；步骤二，拓扑关系构建将步骤一所输入的多幅多视角密集点云数据合并在一起，建立一棵k-d树，构建点云数据中点数据间的拓扑关系；步骤三，基于双边滤波器的点数据归属在步骤二所构建的拓扑关系的基础上，查询点云中每一个点数据的局部邻域信息，然后采用双边滤波器更新每一个点数据的三维坐标，从而实现点数据的归属；步骤四，Mean-shift聚类融合对步骤三所归属后的点数据，采用Mean-shift聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据聚集在一起，获得局部模式点数据；并用所获得的局部模式点数据代替所聚集的欧式距离小于多视角密集点云数据的平均点距的点数据，实现冗余消除；步骤五，输出融合结果将步骤四聚类融合后的点数据输出，获得融合点云数据。进一步地，所述基于双边滤波器的多视角密集点云数据融合方法还包括，输入多视角密集点云数据以后，在进行步骤二拓扑关系构建之前，需要计算多视角密集点云数据的平均点距D。所述多视角密集点云数据平均点距D的计算方法如下：2.1)随机抽取一幅待融合的点云数据，为该幅点云数据建立一棵k-d树；2.2)在步骤2.1)所抽取的点云数据中随机抽取若干个点数据，采用k-d树的最邻近点搜索方法查询所抽取的每个点数据的最邻近的一个点数据；2.3)计算步骤2.2)所抽取的每个点数据与其最邻近的一个点数据的欧式距离，并取这些欧式距离的平均值作为多视角密集点云数据的平均点距D。进一步地，所述步骤三中点数据的局部邻域信息是指位于底面半径为r、高度为h的圆柱形区域内，并同时满足欧氏距离判据及法向量判据的点数据集合。所述每一个点数据的局部邻域信息查询方法如下：3.1)由使用者指定圆柱形区域的底面半径r及高度h，并计算该圆柱形区域的外接圆的半径R：3.2)在步骤二所构建的拓扑关系的基础上，采用k-d树的固定半径搜索方法搜索位于以点数据p＝(v,n)为中心，以步骤3.1)计算的半径R为半径的球形区域内的点数据pi＝(vi,ni)；其中，v表示点数据P的三维坐标向量，n表示为点数据P的法向量，vi表示点数据Pi的三维坐标向量，ni表示为点数据Pi的法向量；3.3)对于步骤3.2)所搜索的球形区域内的每个点数据Pi，如果其同时满足所述的欧式距离判据及法向量判据ni·n＞0，则认为该点数据Pi为点数据P的一个邻域点，否则为非邻域点；其中，符号·表示向量内积，符号||表示取绝对值，符号||||表示取三维坐标向量的模；点数据P的邻域点数据Pi的集合构成了点数据P的邻域。进一步地，所述基于双边滤波器的多视角密集点云数据融合方法的步骤三中，采用双边滤波器更新每一个点数据p＝(v,n)的三维坐标的方法如下：其中，v表示待更新的点数据P的三维坐标向量，n表示待更新的点数据P的法向量，v′表示更新后的点数据P的三维坐标向量，hi表示点数据P的第i个邻域点数据Pi至点数据P的距离在点数据P的法向量n方向上的投影：hi＝|(vi-v)·n|，ri表示点数据P的第i个邻域点数据Pi至点数据P的距离在点数据P的切平面上的投影：vi表示点数据P的第i个邻域点数据Pi的三维坐标向量，N为所查询的点数据P的邻域点数据的个数，权值函数r为查询点数据P的邻域时由使用者所指定的圆柱形区域的底面半径，h为查询点数据P的邻域时由使用者本文档来自技高网...

【技术保护点】
一种基于双边滤波器的多视角密集点云数据融合方法，其特征在于，包括下述步骤：步骤一，输入多视角密集点云数据输入测量设备采集的待融合的多幅多视角密集点云数据，输入的多幅多视角密集点云数据需同时包含三维坐标信息及法向量信息；步骤二，拓扑关系构建将步骤一所输入的多幅多视角密集点云数据合并在一起，建立一棵k‑d树，构建点云数据中点数据间的拓扑关系；步骤三，基于双边滤波器的点数据归属在步骤二所构建的拓扑关系的基础上，查询点云中每一个点数据的局部邻域信息，然后采用双边滤波器更新每一个点数据的三维坐标，从而实现点数据的归属；步骤四，Mean‑shift聚类融合对步骤三所归属后的点数据，采用Mean‑shift聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据聚集在一起，获得局部模式点数据；并用所获得的局部模式点数据代替所聚集的欧式距离小于多视角密集点云数据的平均点距的点数据，实现冗余消除；步骤五，输出融合结果将步骤四聚类融合后的点数据输出，获得融合点云数据。

【技术特征摘要】
1.一种基于双边滤波器的多视角密集点云数据融合方法，其特征在于，包括下述步骤：步骤一，输入多视角密集点云数据输入测量设备采集的待融合的多幅多视角密集点云数据，输入的多幅多视角密集点云数据需同时包含三维坐标信息及法向量信息；步骤二，拓扑关系构建将步骤一所输入的多幅多视角密集点云数据合并在一起，建立一棵k-d树，构建点云数据中点数据间的拓扑关系；步骤三，基于双边滤波器的点数据归属在步骤二所构建的拓扑关系的基础上，查询点云中每一个点数据的局部邻域信息，然后采用双边滤波器更新每一个点数据的三维坐标，从而实现点数据的归属；步骤四，Mean-shift聚类融合对步骤三所归属后的点数据，采用Mean-shift聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据聚集在一起，获得局部模式点数据；并用所获得的局部模式点数据代替所聚集的欧式距离小于多视角密集点云数据的平均点距的点数据，实现冗余消除；步骤五，输出融合结果将步骤四聚类融合后的点数据输出，获得融合点云数据；所述步骤三中点数据的局部邻域信息是指位于底面半径为r、高度为h的圆柱形区域内，并同时满足欧氏距离判据及法向量判据的点数据集合；所述每一个点数据的局部邻域信息查询方法如下：3.1)由使用者指定圆柱形区域的底面半径r及高度h，并计算该圆柱形区域的外接圆的半径R：3.2)在步骤二所构建的拓扑关系的基础上，采用k-d树的固定半径搜索方法搜索位于以点数据p＝(v,n)为中心，以步骤3.1)计算的半径R为半径的球形区域内的点数据pi＝(vi,ni)；其中，v表示点数据P的三维坐标向量，n表示为点数据P的法向量，vi表示点数据Pi的三维坐标向量，ni表示为点数据Pi的法向量；3.3)对于步骤3.2)所搜索的球形区域内的每个点数据Pi，如果其同时满足所述的欧式距离判据及法向量判据ni·n＞0，则认为该点数据Pi为点数据P的一个邻域点，否则为非邻域点；其中，符号·表示向量内积，符号||表示取绝对值，符号||||表示取三维坐标向量的模；点数据P的邻域点数据Pi的集合构成了点数据P的邻域。2.根据权利要求1所述的基于双边滤波器的多视角密集点云数据融合方法，其特征在于还包括，输入多视角密集点云数据以后，在进行步骤二拓扑关系构建之前，需要计算多视角密集点云数据的平均点距D；所述多视角密集点云数据平均点距D的计算方法如下：2.1)随机抽取一幅待融合的点云数据，为该幅点云数据建立一棵k-d树；2.2)在步骤2.1)所抽取的点云数据中随机抽取若干个点数据，采用k-d树的最邻近点搜索方法查询所抽取的每个点数据的最邻近的一个点数据；2.3)计算步骤2.2)所抽取的每个点数据与其最邻近的一个点数据的欧式距离，并取这些欧式距离的平均值作为多视角密集点云数据的平均点距D。3.根据权利要求1所述的基于双边滤波器的多视角密集点云数据融合方法，其特征在于所述步骤三中，采用双边滤波器更新每一个点数据p＝(v,n)的三维坐标的方法如下：其中，v表示待更新的点数据P的三维坐标向量，n表示待更新的点数据P...

【专利技术属性】
技术研发人员：史宝全，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61