CN104050640B - 一种多视角密集点云数据融合方法 - Google Patents

一种多视角密集点云数据融合方法 Download PDF

Info

Publication number
CN104050640B
CN104050640B CN201410266976.3A CN201410266976A CN104050640B CN 104050640 B CN104050640 B CN 104050640B CN 201410266976 A CN201410266976 A CN 201410266976A CN 104050640 B CN104050640 B CN 104050640B
Authority
CN
China
Prior art keywords
point
point data
data
intensive
visual angles
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410266976.3A
Other languages
English (en)
Other versions
CN104050640A (zh
Inventor
史宝全
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201410266976.3A priority Critical patent/CN104050640B/zh
Publication of CN104050640A publication Critical patent/CN104050640A/zh
Application granted granted Critical
Publication of CN104050640B publication Critical patent/CN104050640B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种多视角密集点云数据融合方法,属于光学三维非接触式测量技术领域,该方法的实现步骤包括1)输入多视角密集点云数据;2)拓扑关系构建;3)重叠区域识别;4)重叠区域点数据归属;5)聚类融合;6)输出融合结果。通过所述操作步骤,本发明有效地克服了现有点云融合技术的缺陷,可以将多幅部分区域重叠的多视角密集点云数据一次性融合为一完整的、单层的、光顺的、分布均匀的点云数据。

Description

一种多视角密集点云数据融合方法
技术领域
本发明属于光学三维非接触式测量技术领域,涉及一种多视角密集点云数据融合方法,更进一步涉及一种新的结合最小二乘曲面拟合方法及聚类技术的多视角密集点云数据融合方法。
背景技术
光学三维测量技术是集光、机、电和计算机技术于一体的智能化、可视化的高新技术,主要用于对物体空间外形和结构的扫描,以得到物体的三维轮廓,获得物体表面点的三维空间坐标。随着现代检测技术的进步,特别是随着激光技术、计算机技术以及图像处理技术等高新技术的发展,三维测量技术逐步成为人们的研究重点。光学三维测量技术由于具有非接触、快速测量、精度高等优点,使其在航空航天、军工、汽车和装备制造等行业得到了广泛应用。三维测量技术是获取物体表面各点空间坐标的技术,主要包括接触式和非接触式两大类。其中,通过光学三维非接触式测量设备采集的多视角点云数据的处理技术,如融合、降噪、曲面重建等,是当前研究的热点。
点云融合是指消除由测量误差和匹配误差等导致的多视角点云重叠区域的噪声、分层和冗余,建立细节特征清晰、表面光顺的单层点云模型。点云融合技术随着三维光学非接触式测量技术的发展而不断发展。按照融合方式可以将现有点云融合方法分为三类:基于隐式曲面重构的点云融合技术,基于显式曲面重构的点云融合技术和基于聚类的点云融合技术。
基于隐式曲面重构的点云融合技术一般先通过有向距离场(Signed DistanceField,SDF)或者移动最小二乘(Moving Least Squares,MLS)等技术构建一光顺的隐式曲面,然后通过三角网格化方法或者B样条曲面构建技术将隐式曲面转换为显式曲面,从而实现多视角密集点云数据的融合。这类点云融合方法在构建隐式曲面的过程中可以有效地消除由测量误差和匹配偏差等导致的分层及噪音,但这类方法在将隐式曲面转换为显式曲面的过程中会占用大量的计算机资源,融合效率低。
基于显式曲面重构的点云融合技术在网格化的过程中(通过构建三角网格曲面、泊松曲面或者nurbs曲面等)将多视角点云融合在一起,融合结果为一显式的网格曲面。但这类方法对匹配偏差及噪声比较敏感,融合后重叠区域表面的融合效果一般比较差,会保留融合前的痕迹。此外,这类方法占用计算机资源多、效率低,不适用于大规模点云数据的融合处理。
基于聚类的点云融合技术通过聚类的方法将重叠区域的同名点数据聚集在一起并融合,从而实现冗余消除。这类点云融合方法是一种无网格点云处理技术,此类方法占用计算机资源少,但当待融合的点云幅数多于两幅时,需要通过增量式的方法进行融合,即先将其中两幅融合在一起,然后再加进来一幅进行融合,依次类推,直至所有点云融合完毕为止。增量式的融合方式一方面效率低,另外一方面不能保证融合后的点云位于最优的曲面上。
发明内容
为了克服上述现有方法的技术缺陷,满足包含重叠区域的多视角密集点云数据的融合处理要求,本发明提供了一种多视角密集点云数据融合方法。该方法可一次性将多幅包含匹配偏差、噪声的多视角密集点云数据融合为一幅完整、光顺、分布均匀的单层点云数据。
为达到以上目的,本发明采取的技术方案是:
一种多视角密集点云数据融合方法,包括下述步骤:
步骤一,输入多视角密集点云数据
输入测量设备采集的待融合的多幅多视角密集点云数据,输入的多幅多视角密集点云数据需同时包含三维坐标信息及法向量信息。
步骤二,拓扑关系构建
将步骤一所输入的多幅多视角密集点云数据合并在一起,建立一棵k-d树,构建点云数据中点数据间的拓扑关系。
步骤三,重叠区域识别
在步骤二所构建的拓扑关系的基础上,查询每一个点数据的局部邻域信息,如果该点数据的局部邻域包含了其它视角点云中的点数据,则认为该点数据是一个重叠点数据,否则是一个非重叠点数据;重叠点数据的集合构成重叠区域,非重叠点数据的集合构成非重叠区域。
步骤四,重叠区域点数据归属
对重叠区域中的每个点数据,首先依据步骤三所查询的局部邻域信息拟合一最小二乘曲面,然后将该点数据及其微邻域内的点数据归属至所拟合的最小二乘曲面上。
步骤五,聚类融合
对步骤四所归属至最小二乘曲面上的重叠点数据,采用聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据融合在一起。
步骤六,输出融合结果
将步骤三所识别的非重叠区域内的点数据及步骤五所融合后的重叠点数据合并在一起,成为一幅完整的点云数据并输出。
进一步地,所述多视角密集点云数据融合方法还包括,输入多视角密集点云数据以后,在进行步骤二拓扑关系构建之前,需要计算多视角密集点云数据的平均点距。
所述多视角密集点云数据平均点距的计算方法如下:
2.1)随机抽取一幅待融合的点云数据,为该幅点云数据建立一棵k-d树;
2.2)在步骤2.1)所抽取的点云数据中随机抽取若干个点数据,采用k-d树的最邻近点搜索方法查询所抽取的每个点数据的最邻近的一个点数据;
2.3)计算步骤2.2)所抽取的每个点数据与其最邻近的一个点数据的欧式距离,并取这些欧式距离的平均值作为多视角密集点云数据的平均点距。
作为优选,在合并多幅多视角密集点云数据的过程中,为多视角密集点云数据中的每个点数据存储一个用于表明其所属视角的身份信息。
进一步地,所述步骤三中点数据的局部邻域信息是指位于底面半径为r、高度为h的圆柱形区域内,并同时满足欧氏距离判据及法向量判据的点数据集合。
所述每一个点数据的局部邻域信息查询方法如下:
4.1)由使用者指定圆柱形区域的底面半径r及高度h,并计算该圆柱形区域的外接圆的半径R:
R = r 2 + h 2 / 4
4.2)在步骤二所构建的拓扑关系的基础上,采用k-d树的固定半径搜索方法搜索位于以点数据p=(v,n)为中心,以步骤4.1)计算的半径R为半径的球形区域内的点数据pi=(vi,ni);其中,v表示点数据P的三维坐标向量,n表示为点数据P的法向量,vi表示点数据Pi的三维坐标向量,ni表示为点数据Pi的法向量;
4.3)对于步骤4.2)所搜索的球形区域内的每个点数据Pi,如果其同时满足所述的欧式距离判据及法向量判据ni·n>0,则认为该点数据Pi为点数据P的一个邻域点,否则为非邻域点;其中,符号·表示向量内积;点数据P的邻域点数据的集合构成了点数据P的邻域。
进一步地,所述步骤四重叠区域点数据归属中,对重叠区域中的每个点数据,首先依据步骤三所查询的局部邻域信息拟合一加权最小二乘曲面,然后将该点数据及其微邻域内的点数据归属至所拟合的最小二乘曲面上。
进一步地,所述点数据的微邻域是指圆柱形区域底面半径与多视角密集点云数据的平均点距相等时的局部邻域;所述点数据归属是指该点数据沿其法向量方向移动至与所拟合的最小二乘曲面相交的交点处。
进一步地,所述步骤五聚类融合的过程包括如下步骤:
7.1)从某一点数据p出发,采用k-d树的固定半径搜索方法,搜索以该点数据p为中心、半径为r的球形邻域内的点数据。
7.2)对步骤7.1)所搜索的点数据p的球形邻域内的点数据进行加权平均,得到一个新的三维点数据
其中,pi表示点数据p的球形邻域内的第i个点数据,N表示点数据p的球形邻域内的点数据个数。
7.3)判断点数据p′与点数据p是否重合;如果点数据p′与点数据p重合,则该次聚类结束,保留点数据p′,删除步骤7.1)所搜索的点数据p的球形邻域内的所有点数据,转入步骤7.4);如果点数据p′与点数据p不重合,则从点数据p′出发,重复步骤7.1)至7.3)。
7.4)重复步骤7.1)至7.3),直至步骤三所识别的重叠区域内的所有重叠点数据聚类融合完毕为止。
进一步地,所述步骤7.1)中球形邻域半径r的取值范围为多视角密集点云数据的平均点距的δ倍,其中δ为一系数,δ的取值范围为0.5<δ<1。
特别地,所述球形邻域半径r的取值范围为多视角密集点云数据的平均点距的倍。
与现有技术相比,本发明具有以下优点:
(1)本发明方法可一次性融合多幅多视角点云数据,比传统增量式融合方法效率高。
(2)本发明方法融合多视角点云数据时,重叠区域的点数据被归属至具有最小二乘意义的曲面上,因此融合结果更加合理,融合表面更加光顺。
(3)本发明方法融合多视角点云数据时,在点数据归属的过程中,点数据只沿其法向量方向移动,因此可保持原始特征。
(4)本发明方法融合多视角点云数据时,在点数据归属的过程中,每次将一个微小邻域内的点数据同时归属至具有最小二乘意义的曲面上,比传统单点归属的方法效率高。
(5)本发明方法获得的融合点云密度分布较均匀。
以下将结合实施例及附图对本发明做进一步详细阐述。
附图说明
图1是本发明具体操作步骤的流程图。
图2是本发明方法定义的局部邻域示意图。
图3是本发明方法中重叠区域点数据归属示意图。
图4是本发明方法中一次聚类融合示意图。
图5a是待融合的三个视角的部分重叠的人头点云数据。
图5b是图5a中方框圈中部分的放大图。
图5c是融合后的人头点云数据。
图5d是图5c中方框圈中部分的放大图。
具体实施方式
实施例1:
本发明提出一种多视角密集点云数据融合方法,如图1所示。在融合多幅多视角密集点云数据时,第一步,输入多视角密集点云数据。要求输入的多视角密集点云数据同时包含三维坐标信息及法向量信息。也就是说输入的多视角密集点云中的每个点数据p=(v,n),其中v=(vx,vy,vz)表示三维坐标向量,n=(nx,ny,nz)表示法向量。
输入多视角密集点云数据以后,在进行第二步拓扑关系构建之前,需要先计算多视角密集点云数据的平均点距D以备后期使用。多视角密集点云数据平均点距D的计算方法如下:
1)随机抽取一幅待融合的点云数据,采用文献“Multidimensional binarysearch trees used for associative”(JL.Bentley,Commun ACM1975,18:509–17.)所描述的方法为该幅点云数据建立一棵k-d树。
2)在第1)步所抽取的点云中随机抽取若干个点数据(如1000个),然后采用k-d树的最邻近点搜索方法查询所抽取的每个点数据的最邻近的一个点数据。
3)计算第2)步所抽取的每个点数据与其最邻近的一个点数据的欧式距离,并取这些距离的平均值作为多视角密集点云数据的平均点距D。
第二步,拓扑关系构建。将第一步所输入的多幅多视角密集点云数据合并在一起,建立一棵k-d树,构建多视角密集点云中点数据间的拓扑关系。由于将多视角密集点云数据合并在一起会丢失视角信息,为此,在合并的过程中为多视角密集点云中的每个点数据存储一个用于表明其所属视角的身份信息(ID)。k-d树构造方法与第一步中所采用的k-d树构造方法一致。
第三步,重叠区域识别。重叠区域识别流程包括如下步骤:
1)局部邻域查询。在第二步所构建的拓扑关系的基础上,快速查询多视角密集点云中每一个点数据的局部邻域信息。在本发明方法中,任意一点数据p=(v,n)的局部邻域定义为位于底面半径为r、高度为h的圆柱形区域内,并同时满足欧氏距离判据及法向量判据的点数据集合Ω。
所述的欧氏距离判据为 | ( v i - v ) · n | ≤ h / 2 d ( p i ) ≤ r ;
所述的法向量判据为ni·n>0;
所述的点数据集合为Ω={pi|d(pi)≤d(pi+1)},0≤i<N-1;
其中,pi=(vi,ni)为三维空间中与点数据P不重合的点数据,N为点数据集合Ω中包含的点数据个数,符号||表示取绝对值,符号||||表示取三维向量的模,符号·表示向量内积,r为圆柱形区域底面半径,h为圆柱形区域高度,i表示第i个。
如图2所示,图2中点云scan1和点云scan2具有相同方向、相互重叠且点云间隙大小为G。点云scan1中点数据p的局部邻域为包含在以点数据p为中心,以点数据p的法向量n方向为轴线方向,底面半径为r,高度为h的圆柱形区域内的点数据集合Ω。需要强调的是,点数据集合Ω內的点数据需与点数据P同向,即满足所述的法向量判据。底面半径r及高度h的值由使用者指定,要求高度h的值大于多幅点云之间的最大间隙。一般情况下,选择r=3~5D(D为第一步所计算的多视角密集点云数据的平均点距),h=2r即可。
点数据P的局部邻域信息查询步骤如下:
a)由使用者指定底面半径r及高度h这两个参数,并计算该圆柱的外接圆的半径R:
R = r 2 + h 2 / 4 .
b)在第二步所构建的拓扑关系的基础上,采用k-d树的固定半径搜索方法(fixed-radius search,也称为范围搜素方法),搜索位于以点数据P为中心,以第a)计算的半径R为半径的球形区域内的点数据。
c)对于第b)步所搜索的球形区域内的每个点数据Pi,如果其同时满足所述的欧式距离判据及法向量判据ni·n>0,则认为该点数据Pi为点数据P的一个邻域点,否则为非邻域点。其中,符号i表示向量内积。点数据P的邻域点数据的集合构成了点数据P的邻域。
2)重叠点判定。对于第1)步所查询的点数据P的邻域内的任意一个点数据Pi,如果其身份信息与点数据P的身份信息不一致,则点数据P为一重叠点,否则点数据P为一非重叠点。重叠点数据的集合构成了重叠区域,非重叠点数据的集合构成了非重叠区域。
第四步,重叠区域点数据归属。重叠区域点数据归属流程包括以下步骤:
1)对于第三步所识别的重叠区域内的每个点数据P,在其局部邻域内拟合一加权最小二乘多项式曲面(如图3所示)。所拟合的最小二乘多项式曲面的方程如下:
z=a+by+cy2+dy3+ex+fxy+gxy2+hx2+ix2y+jx3
其中,x,y,z为方程的自变量,a,b,c,d,e,f,g,h,i,j为待拟合的多项式曲面的系数。
局部邻域内每个点数据pi=(vi,ni)的权值wi的计算方法如下:
w i = e - d 2 ( p i ) / r 2 .
其中,符号·表示向量内积,r为使用者在第三步所指定的底面半径。
2)将点数据P及其微邻域内的点数据归属至第1)步所拟合的最小二乘曲面上。点数据P的微邻域定义为当底面半径r=1D时的局部邻域(图3),其中,D为第一步所计算的多视角密集点云数据的平均点距。点数据的归属即将该点数据沿其法向量方向移动至与第1)步所拟合的最小二乘曲面相交的交点处。其实质就是求过该点数据并沿其法向量方向的一条直线与第1)步所拟合的三次多项式曲面的交点,并用交点代替该点数据从而实现归属(图3)。
三维空间中过一点p=(v,n),并沿其法向量方向的直线方程可以用参数形式表示如下:
x = v x + n x t y = v y + n y t z = v z + n z t
其中,x,y,z为方程的自变量,v=(vx,vy,vz)为点数据P的三维坐标向量,n=(nx,ny,nz)为点数据P的法向量,t为待求解的未知量。
将直线的参数方程带入第1)步所拟合的多项式曲面方程,可以得到一个关于t的三次多项式方程:
At3+Bt2+Ct+D=0。
其中, A = dn y 3 + jn x 3 + gn x n y 2 + in x 2 n y ;
B = ( h + 3 jv x + iv y ) n x 2 + ( f + 2 iv x + 2 gv y ) n x n y + ( c + 3 dv y + gx ) n y 2 ;
C = ( e + 2 hv x + 3 jv x 2 + fv y + 2 iv x v y + gv y 2 ) n x + ( b + 2 cv y + 3 dv y 2 + fv x + 2 gv x v y + iv x 2 ) n y - n z ;
D = a + ev x + bv y + fv x v y + ( h + iv y ) v x 2 + ( c + gv x ) v y 2 + jv x 3 + dv y 3 - v z .
采用卡尔丹公式法可以求得该方程的所有实根t,带入直线公式便可得到候选的交点。如果候选的交点在点数据的局部邻域内,则认为该点即为所求的交点,否则不是所求的交点。
第五步,聚类融合。对第四步所归属至最小二乘曲面上的重叠点数据,采用聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据聚集在一起并融合,从而实现冗余消除。聚类融合流程包括如下步骤(如图4所示):
1)从某一点p=(v,n)(图4(b))出发,采用k-d树的固定半径搜索方法,搜索以该点为中心,半径r=δD的球形邻域内的点数据。其中,δ为一系数,D为第一步所计算的多视角密集点云数据的平均点距,0.5<δ<1,一般情况下,选
2)对第1)步所搜索的球形邻域内的点数据进行加权平均得到一个新的三维点数据 p ′ = 1 N Σ i = 1 N p i (图4(c))。
其中,pi表示点数据p的球形邻域内的第i个点数据,N表示点数据p的球形邻域内的点数据个数。
3)判断点数据p′与点数据p是否重合。如果重合,则该次聚类结束,保留p′,删除第1)步所搜索的点数据p的球形邻域内的所有点数据,转入第4)步(图4(h));如果不重合,则从p′出发,然后重复步骤1)至3)(图4(d)-(g))。
4)重复步骤1)至3)直至第三步所识别的重叠区域内的所有重叠点数据聚类融合完毕为止。
第六步,输出融合结果。将第三步所识别的非重叠区域内的点数据及第五步所聚类融合的重叠点数据合并在一起成为一幅完整的、单层的、光顺的点云数据并输出。
实施例2:
以下结合模拟实验进行说明,其中本发明方法在VS2010及opengl平台上实现相应的计算流程并在Intel i7-4770CPU3.4GHz、16GB内存的PC机上运行。
三个视角的部分区域重叠的人头点云数据如图5(a)所示,总共包含1023124个三维点数据。为了更清晰的显示重叠区域,图5(b)为图5(a)中方框圈中部分的放大图。图5(c)为融合后点云数据,总共包含574097个点数据,融合过程总共耗时70s。与图5(b)对应,为了更清晰的显示该处重叠区域融合后的效果,图5(d)为图5(c)中方框圈中部分的放大图。通过本例也可以说明,本发明方法可以将多视角密集点云数据快速地融合为一完整的、单层的、光顺的、分布均匀的点云数据。
上面结合实施例对本发明做了进一步的叙述,但本发明并不限于上述实施方式,在本领域的普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。

Claims (8)

1.一种多视角密集点云数据融合方法,其特征在于,包括下述步骤:
步骤一,输入多视角密集点云数据
输入测量设备采集的待融合的多幅多视角密集点云数据,输入的多幅多视角密集点云数据需同时包含三维坐标信息及法向量信息;
步骤二,拓扑关系构建
将步骤一所输入的多幅多视角密集点云数据合并在一起,建立一棵k-d树,构建点云数据中点数据间的拓扑关系;
步骤三,重叠区域识别
在步骤二所构建的拓扑关系的基础上,查询每一个点数据的局部邻域信息,如果该点数据的局部邻域包含了其它视角点云中的点数据,则认为该点数据是一个重叠点数据,否则是一个非重叠点数据;重叠点数据的集合构成重叠区域,非重叠点数据的集合构成非重叠区域;
步骤四,重叠区域点数据归属
对重叠区域中的每个点数据,首先依据步骤三所查询的局部邻域信息拟合一最小二乘曲面,然后将该点数据及其微邻域内的点数据归属至所拟合的最小二乘曲面上;
步骤五,聚类融合
对步骤四所归属至最小二乘曲面上的重叠点数据,采用聚类的方法将欧式距离小于多视角密集点云数据的平均点距的点数据融合在一起;
步骤六,输出融合结果
将步骤三所识别的非重叠区域内的点数据及步骤五所融合后的重叠点数据合并在一起,成为一幅完整的点云数据并输出;
所述步骤三中点数据的局部邻域信息是指位于底面半径为r、高度为h的圆柱形区域内,并同时满足欧氏距离判据及法向量判据的点数据集合;
所述每一个点数据的局部邻域信息查询方法如下:
4.1)由使用者指定圆柱形区域的底面半径r及高度h,并计算该圆柱形区域的外接圆的半径R:
R = r 2 + h 2 / 4
4.2)在步骤二所构建的拓扑关系的基础上,采用k-d树的固定半径搜索方法搜索位于以点数据p=(v,n)为中心,以步骤4.1)计算的半径R为半径的球形区域内的点数据pi=(vi,ni);其中,v表示点数据P的三维坐标向量,n表示为点数据P的法向量,vi表示点数据Pi的三维坐标向量,ni表示为点数据Pi的法向量;
4.3)对于步骤4.2)所搜索的球形区域内的每个点数据Pi,如果其同时满足所述的欧式距离判据及法向量判据ni·n>0,则认为该点数据Pi为点数据P的一个邻域点,否则为非邻域点;其中,符号·表示向量内积;点数据P的邻域点数据的集合构成了点数据P的邻域。
2.根据权利要求1所述的多视角密集点云数据融合方法,其特征在于还包括,输入多视角密集点云数据以后,在进行步骤二拓扑关系构建之前,需要计算多视角密集点云数据的平均点距;
所述多视角密集点云数据平均点距的计算方法如下:
2.1)随机抽取一幅待融合的点云数据,为该幅点云数据建立一棵k-d树;
2.2)在步骤2.1)所抽取的点云数据中随机抽取若干个点数据,采用k-d树的最邻近点搜索方法查询所抽取的每个点数据的最邻近的一个点数据;
2.3)计算步骤2.2)所抽取的每个点数据与其最邻近的一个点数据的欧式距离,并取这些欧式距离的平均值作为多视角密集点云数据的平均点距。
3.根据权利要求1所述的多视角密集点云数据融合方法,其特征在于,在合并多幅多视角密集点云数据的过程中,为多视角密集点云数据中的每个点数据存储一个用于表明其所属视角的身份信息。
4.根据权利要求1所述的多视角密集点云数据融合方法,其特征在于,所述步骤四重叠区域点数据归属中,对重叠区域中的每个点数据,首先依据步骤三所查询的局部邻域信息拟合一加权最小二乘曲面,然后将该点数据及其微邻域内的点数据归属至所拟合的最小二乘曲面上。
5.根据权利要求4所述的多视角密集点云数据融合方法,其特征在于,所述点数据的微邻域是指圆柱形区域底面半径与多视角密集点云数据的平均点距相等时的局部邻域;所述点数据归属是指该点数据沿其法向量方向移动至与所拟合的最小二乘曲面相交的交点处。
6.根据权利要求1所述的多视角密集点云数据融合方法,其特征在于,所述步骤五聚类融合的过程包括如下步骤:
7.1)从某一点数据p出发,采用k-d树的固定半径搜索方法,搜索以该点数据p为中心、半径为r的球形邻域内的点数据;
7.2)对步骤7.1)所搜索的点数据p的球形邻域内的点数据进行加权平均,得到一个新的三维点数据
其中,pi表示点数据p的球形邻域内的第i个点数据,N表示点数据p的球形邻域内的点数据个数;
7.3)判断点数据p′与点数据p是否重合;如果点数据p′与点数据p重合,则该次聚类结束,保留点数据p′,删除步骤7.1)所搜索的点数据p的球形邻域内的所有点数据,转入步骤7.4);如果点数据p′与点数据p不重合,则从点数据p′出发,重复步骤7.1)至7.3);
7.4)重复步骤7.1)至7.3),直至步骤三所识别的重叠区域内的所有重叠点数据聚类融合完毕为止。
7.根据权利要求6所述的多视角密集点云数据融合方法,其特征在于,所述步骤7.1)中球形邻域半径r的取值范围为多视角密集点云数据的平均点距的δ倍,其中δ为一系数,δ的取值范围为0.5<δ<1。
8.根据权利要求7所述的多视角密集点云数据融合方法,其特征在于,所述球形邻域半径r的取值范围为多视角密集点云数据的平均点距的倍。
CN201410266976.3A 2014-06-16 2014-06-16 一种多视角密集点云数据融合方法 Active CN104050640B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410266976.3A CN104050640B (zh) 2014-06-16 2014-06-16 一种多视角密集点云数据融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410266976.3A CN104050640B (zh) 2014-06-16 2014-06-16 一种多视角密集点云数据融合方法

Publications (2)

Publication Number Publication Date
CN104050640A CN104050640A (zh) 2014-09-17
CN104050640B true CN104050640B (zh) 2017-01-04

Family

ID=51503427

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410266976.3A Active CN104050640B (zh) 2014-06-16 2014-06-16 一种多视角密集点云数据融合方法

Country Status (1)

Country Link
CN (1) CN104050640B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110553624A (zh) * 2018-06-01 2019-12-10 上海辉格科技发展有限公司 一种基于激光扫描仪的隧道模型建立方法
CN109858512B (zh) * 2018-12-10 2021-08-03 北京百度网讯科技有限公司 点云数据的处理方法、装置、设备、车辆及存储介质
CN110415362B (zh) * 2019-07-31 2023-06-06 易思维(杭州)科技有限公司 重叠点云融合方法
CN110765219B (zh) * 2019-08-05 2022-06-28 上海晶赞融宣科技有限公司 地理围栏生成方法、装置、计算机设备和存储介质
CN112580677A (zh) * 2019-09-29 2021-03-30 北京地平线机器人技术研发有限公司 一种点云数据点的分类方法及装置
CN114937124B (zh) * 2022-07-25 2022-10-25 武汉大势智慧科技有限公司 基于倾斜摄影的片状目标物的三维重建方法、装置及设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103729886A (zh) * 2013-12-30 2014-04-16 浙江大学 一种基于顶点重定位的三角网格表面模型概率融合方法
CN103745018A (zh) * 2014-02-11 2014-04-23 天津市星际空间地理信息工程有限公司 一种多平台点云数据融合方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103729886A (zh) * 2013-12-30 2014-04-16 浙江大学 一种基于顶点重定位的三角网格表面模型概率融合方法
CN103745018A (zh) * 2014-02-11 2014-04-23 天津市星际空间地理信息工程有限公司 一种多平台点云数据融合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
移动最小二乘增量式多视点云数据融合算法;曹巨明等;《西安交通大学学报》;20090930;第46-50页 *

Also Published As

Publication number Publication date
CN104050640A (zh) 2014-09-17

Similar Documents

Publication Publication Date Title
CN104050640B (zh) 一种多视角密集点云数据融合方法
CN111299815B (zh) 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法
US9348877B2 (en) Methods, systems, and data structures for performing searches on three dimensional objects
Mehra et al. Abstraction of man-made shapes
CN104361632B (zh) 一种基于Hermite径向基函数的三角网格补洞方法
CN101719140B (zh) 一种图形检索方法
CN105654483B (zh) 三维点云全自动配准方法
CN106780458B (zh) 一种点云骨架提取方法及装置
CN108171780A (zh) 一种基于激光雷达构建室内真实三维地图的方法
CN107392875A (zh) 一种基于k近邻域划分的点云数据去噪方法
CN101477529B (zh) 一种三维对象的检索方法和装置
Kim et al. Deep-learning-based retrieval of piping component catalogs for plant 3D CAD model reconstruction
CN103886569A (zh) 连续帧多特征点无人机侦察图像并行与匹配精度约束的拼接方法
CN103914837B (zh) 一种适用于多视角点云处理的圆柱形邻域及其搜索方法
CN106874580A (zh) 一种基于点云数据的弯管模型重建方法
Tam et al. Deformable model retrieval based on topological and geometric signatures
CN105844602A (zh) 一种基于体元的机载lidar点云三维滤波方法
CN100585638C (zh) 基于线框的曲面体三维边界表示模型重建方法及其装置
CN113628263A (zh) 一种基于局部曲率及其近邻特征的点云配准方法
CN104504709A (zh) 一种基于特征球的室外场景三维点云数据的分类方法
CN104050639B (zh) 一种基于双边滤波器的多视角密集点云数据融合方法
CN112288859B (zh) 一种基于卷积神经网络的三维人脸建模方法
CN108256218A (zh) 一种基于实测腰线数据的地下连通巷道精细建模方法
Hu et al. Geometric feature enhanced line segment extraction from large-scale point clouds with hierarchical topological optimization
Li et al. Feature-based alignment and comparison between portion and whole of free-form surfaces

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant