CN104123746B - 一种三维扫描点云中实时法向量的计算方法 - Google Patents
一种三维扫描点云中实时法向量的计算方法 Download PDFInfo
- Publication number
- CN104123746B CN104123746B CN201410327054.9A CN201410327054A CN104123746B CN 104123746 B CN104123746 B CN 104123746B CN 201410327054 A CN201410327054 A CN 201410327054A CN 104123746 B CN104123746 B CN 104123746B
- Authority
- CN
- China
- Prior art keywords
- point
- normal vector
- neighborhood
- data
- points
- 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.)
- Expired - Fee Related
Links
- 239000013598 vector Substances 0.000 title claims abstract description 142
- 238000000205 computational method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims description 17
- 238000011156 evaluation Methods 0.000 claims description 14
- 238000000513 principal component analysis Methods 0.000 claims description 10
- 235000014161 Caesalpinia gilliesii Nutrition 0.000 abstract 1
- 244000003240 Caesalpinia gilliesii Species 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- RGCLLPNLLBQHPF-HJWRWDBZSA-N phosphamidon Chemical compound CCN(CC)C(=O)C(\Cl)=C(/C)OP(=O)(OC)OC RGCLLPNLLBQHPF-HJWRWDBZSA-N 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种三维扫描点云实时法向量的计算方法,其步骤为:(1)采用kinect摄像机进行实物扫描,读取点云数据,采用KD树搜索点云中点的邻域点;(2)根据主元分析法(PCA)对上述查找的邻域点拟合平面,以这拟合平面的法向量作为点云数据的各个点的法向量;(3)通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值;(4)设置点的法向量评估置信度,利用上述步骤(3)所述的数据点的各个邻域点的法向量加权均值进行评估;(5)设定点的法向量置信度的阀值a,判断各个点的法向量置信度,修正该点的法向量。该方法降低计算点云数据的法向量估计的开销时间,还能对点的法向量具有修正作用,避免对点云法向量重新定向计算,减小了计算的复杂度。
Description
技术领域
本发明涉及一种三维扫描点云中实时法向量计算方法,属于计算机三维建模技术领域。
背景技术
在逆向工程中,主要通过三维扫描仪获得物体表面的点云数据,将这些点云数据输入文件进行存储,称为三维点云模型,通常,三维点云模型不能够直接应用在3D造型。应将三维点云模型转化为曲面模型,即实现三维点云模型的表面重建。
由点云模型重建出完整的曲面模型,其中,点云法向量的计算是目前表面重建方法中的重要步骤。现有的点云法向量计算方法可分为三类:
.基于局部表面拟合的方法:该方法首先假设点云的采样表面处处光滑,对点的局部邻域用平面进行拟合,然后对于点云中的每个点p采用主元分析法(PrincipalComponent Analysis,PCA)求解出拟合平面的法向量作为点p的法向量。但是,主元分析法PCA对外点噪声较为敏感。
.基于Delaunay/Voronoi方法:该方法仅适用于无噪声的点云模型的曲面重建。
.基于鲁棒统计的方法:该方法首先处理点云模型中的噪声、外点和尖锐特征,将每个点的邻域进行分类,为每个模型设置合适的最小邻域和距离阈值,但是该方法对法向量估计的计算时间开销很大,无法直接应用到大规模的点云模型。
为提高对点云法向量计算的精度、速度、稳定性, Alexandre Boulch和RenaudMarlet在2012届欧洲几何处理研讨会(Eurographics Symposium on GeometryProcessing 2012)上,针对具有尖锐特征的模型的法向量计算提出了基于随机变换(RHT)鲁棒性法向量计算方法, 该方法对尖锐特征敏感,对噪声和外点具有稳定性,并获得良好的结果。2012年,S.Holzer、R.B.Rusu和M.Dixon等在智能机器人与***国际会议(International Conference on Intelligent Robots and systems)上,提出了基于有序点云的实时法向量计算方法,该方法没有处理边界的一些点的法向量计算的错误问题。例如,经过现有技术的文献检索发现,中国专利号为:200810239327.9,名称为:点云数据的主曲率和主方向估计方法,该方法针对现实世界中由激光扫描得到的离散点云数据,提供一个对主曲率和主方向的估计方法。从该方法中的步骤2中所述的“对于点云数据的每一个点,利用点云数据的kd树查找15个或30个邻近点, 根据最小二乘方法把这些近邻点拟合出一个平面,以这个平面的法向量作为点云法向量的初始估计值,然后通过加权平均算法修正点云数据的各个点的法向量估计,可以看出,该方法计算点云数据的各点的法向量估计需要进行两次计算点云的法向量,其计算时间开销大。
针对点云数据的各点的法向量的重定向,经检索,目前是由Hoppe等在该文提出的方法将法向量的重定向问题转化为一个求解图的最大代价问题。还检索到,在2009年,RaduBogdan Rusu在其博士论文Semantic 3D Object Maps for Everyday Manipulation inHuman Living Environments 中提出,将点云质心作为视点,重新设定点云的法向量的方法,但是,该方法在尖锐特征复杂模型表面法向量测量时,可能产生局部法向量错误。
发明内容
本发明的目的是克服已有技术存在的不足,提出一种三维扫描点云中实时法向量的计算方法,该方法降低计算点云数据的法向量估计的开销时间,还能修正点的法向量,避免对点云法向量重新定向计算,减小了计算的复杂度。
本发明的目的是通过下述技术方案实现的。
一种三维扫描点云中实时法向量计算的方法,其具体步骤为:
(1).采用kinect摄像机进行实物扫描,读取点云数据,对点云数据按KD树划分,得到每个点云数据中每个数据点的k个邻域点;
(2).对于点云数据的每点,利用点云数据的KD树查找i个邻域点,其中的取值为区间[5~20]中的正整数,其中,,根据主元分析法(PCA)对上述查找的邻域点拟合出一个平面,以这个拟合平面的法向量作为点云数据的各个点的法向量,以kinect摄像机位置为视点,法向量朝向视点;
(3).通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值;
(4).设置数据点的法向量评估置信度,利用上述步骤(3)所述的数据点的各个邻域点的法向量加权均值对数据点的法向量进行评估;
(5).设定评估点的法向量置信度的阀值a, 判断数据点的法向量置信度是否大于阈值,即是否满足,>a,如果是大于阈值a,则认为该点的法向量可信,否则认为该点的法向量不可信,将该邻域点法向量的加权均值替换为该点的法向量。
上述步骤(3)所述的通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,其具体如下:
(3-1).设点云数据中的第i个邻域点的法向量,记为,其表达式为:
= (1)
其中,为第i个邻域点的法向量x分量,为第i个邻域点的法向量y分量,为第i个邻域点的法向量z分量,i是指第i个邻域点;
(3-2).计算点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,记为,其表达式为:
(2)
其中,为i个邻域点的法向量加权均值的x分量,为i个邻域点的法向量加权均值的y分量, 为i个邻域点的法向量加权均值的z分量;
所述的i个邻域点的法向量加权均值的x分量,其表达式为:
, (3)
其中,为第i个邻域点的权重,为第i个邻域点的法向量x分量,n为在半径r内邻域点的个数,i的范围为[1,n];
所述的i个邻域点的法向量加权均值的y分量,其表达式为:
(4)
其中,为第i个邻域点的权重,为第i个邻域点的法向量y分量,n为在半径r内邻域点的个数,i的范围为[1,n];
所述的i个邻域点的法向量加权均值的z分量,其表达式为:
(5)
其中,为第i个邻域点的权重,为第i个邻域点的法向量z分量,n为在半径r内邻域点的个数,i的范围为[1,n];
(3-3).计算第i个邻域点的权重,其计算表达式为:
(6)
其中,为第i个邻域点的权重, 为第i个邻域点与数据点之间的距离,为i个邻域点与数据点之间的最大距离。
使用KD树搜寻每一个数据点在半径r内的数据点,r的大小由人为根据点的疏密确定,取r为3mm~10mm中的任意一个值,由公式(4)得到每点的邻域点法向量的加权均值。
上述步骤(4)所述的设置数据点的法向量评估置信度,利用上述步骤(3)所述的数据点的各个邻域点的法向量加权均值对数据点的法向量进行评估,其评估值为置信度,记为,其评估公式(7)如下:
(7)
本发明的一种三维扫描点云中实时法向量的计算方法与现有技术相比较具有的优点如下:该方法降低计算点云数据的法向量估计的开销时间,能对噪声具有一定的抵抗性;还能对点云的法向量计算具有修正作用,避免对点云法向量重新计算,减小了计算的复杂度。
附图说明
图1为本发明的一种三维扫描点云中实时法向量的计算方法的流程图;
图2为kinect摄像机进行实物扫描的物体模型图;
图3为kinect扫描实物得到的一帧点云数据图;
图4为实时计算扫描点云的法向量的图;
图5为采用本发明的方法计算后的点云模型的法向量的效果图;
图6为采用本发明的方法法向量计算后的泊松重建的曲面模型图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
本发明的一种三维扫描点云中实时法向量的计算方法,如图1所示,其步骤为:
(1).采用kinect摄像机进行实物扫描,如图2所示,读取点云数据,如图3所示。,对点云数据按KD树划分,得到每个点云数据的k个邻域点;
(2).对于点云数据的每点,利用点云数据的KD树查找i个邻域点,其中的取值为区间[5~20]中的正整数,其中,,根据主元分析法(PCA)对上述查找的邻域点拟合出一个平面,以这个拟合平面的法向量作为点云数据的各个点的法向量,以kinect摄像机位置为视点,法向量朝向视点;
(3).通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,其具体如下:
(3-1).设点云数据中的第i个邻域点的法向量,记为,其表达式为:
= (1)
其中,为第i个邻域点的法向量x分量,为第i个邻域点的法向量y分量,为第i个邻域点的法向量z分量,i是指第i个邻域点
(3-2).计算点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,记为,其表达式为:
(2)
其中,为i个邻域点的法向量加权均值的x分量,为i个邻域点的法向量加权均值的y分量, 为i个邻域点的法向量加权均值的z分量;
所述的i个邻域点的法向量加权均值的x分量,其表达式为:
, (3).
其中,为第i个邻域点的权重,为第i个邻域点的法向量x分量,n为在半径r内邻域点的个数,i的范围为[1,n];
所述的i个邻域点的法向量加权均值的y分量 ,其表达式为:
(4)
其中,为第i个邻域点的权重, 为第i个邻域点的法向量y分量,n为在半径r内邻域点的个数,i的范围为[1,n];
所述的i个邻域点的法向量加权均值的z分量 ,其表达式为:
(5).
其中,为第i个邻域点的权重,为第i个邻域点的法向量z分量,n为在半径r内邻域点的个数,i的范围为[1,n];
(3-3).计算第i个邻域点的权重,其计算表达式为:
(6)
其中,为第i个近邻点的权重, 为第i个邻域点与数据点之间的距离,为i个邻域点与数据点之间的最大距离;
使用KD树搜寻每一个数据点在半径r内的数据点,r的大小由人为根据点的疏密确定,取r为3mm~10mm中的任意一个值,由公式(4)得到每点的邻域点法向量的加权均值;
(4).设置数据点的法向量评估置信度,利用上述步骤(3)所述的数据点的各个邻域点的法向量加权均值对数据点的法向量进行评估,其评估值为置信度,记为,其评估公式(7)如下:
(7)
(5).设定评估点的法向量置信度的阀值a, 判断点的法向量置信度是否大于阈值,即是否满足>a,如果是大于阈值a,则认为该点的法向量可信,否则认为该点的法向量不可信,将该邻域点法向量的加权均值替换为该点的法向量,如图4所示。
为了验证本发明的一种三维扫描点云中实时法向量的计算方法的效果,点云经过法向量计算后,配准得到点云模型,如图5所示,将该点云模型泊松重建的曲面模型,如图6所示,图6所显示的模型中未出现不闭合的现象。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (2)
1.一种三维扫描点云中实时法向量计算的方法,其特征在于,其具体步骤为:
(1).采用kinect摄像机进行实物扫描,读取点云数据,对点云数据按KD树划分,得到每个点云数据中每个数据点的k个邻域点;
(2).对于点云数据的每点,利用点云数据的KD树查找i个邻域点,其中K的取值为区间[5~20]中的正整数,其中,i=1,2,…,根据主元分析法(PCA)对上述查找的邻域点拟合出一个平面,以这个拟合平面的法向量作为点云数据的各个点的法向量,以kinect摄像机位置为视点,法向量朝向视点;
(3).通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值;
(4).设置数据点的法向量评估置信度,利用上述步骤(3)每个数据点在半径为r内各个邻域点的法向量加权均值对数据点的法向量进行评估;
(5).设定评估点的法向量置信度的阀值a,判断数据点的法向量置信度是否大于阈值,即是否满足,cosθ>a,如果是大于阈值a,则认为该点的法向量可信,否则认为该点的法向量不可信,将该邻域点法向量的加权均值替换为该点的法向量;
上述步骤(3)所述的通过加权平均算法计算出点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,其具体如下:
(3-1).设点云数据中的第i个邻域点的法向量,记为ni,其表达式为:
ni=(Nix,Niy,Niz) (1)
其中,Nix为第i个邻域点的法向量x分量,Niy为第i个邻域点的法向量y分量,Niz为第i个邻域点的法向量z分量,i是指第i个邻域点;
(3-2).计算点云数据的每个数据点在半径为r内各个邻域点的法向量加权均值,记为n',其表达式为:
n'=(n′ix,n′iy,n′iz) (2)
其中,n′ix为i个邻域点的法向量加权均值的x分量,n′iy为i个邻域点的法向 量加权均值的y分量,n′iz为i个邻域点的法向量加权均值的z分量;
所述的i个邻域点的法向量加权均值的x分量n′ix,其表达式为:
其中,αi为第i个邻域点的权重,Nix为第i个邻域点的法向量x分量,m为在半径r内邻域点的个数,i的范围为[1,m];
所述的i个邻域点的法向量加权均值的y分量n′iy,其表达式为:
其中,αi为第i个邻域点的权重,Niy为第i个邻域点的法向量y分量,m为在半径r内邻域点的个数,i的范围为[1,m];
所述的i个邻域点的法向量加权均值的z分量n′iz,其表达式为:
其中,αi为第i个邻域点的权重,Niz为第i个邻域点的法向量z分量,m为在半径r内邻域点的个数,i的范围为[1,m];
(3-3).计算第i个邻域点的权重,其计算表达式为:
其中,αi为第i个邻域点的权重,li为第i个邻域点与数据点之间的距离,lmax为i个 邻域点与数据点之间的最大距离。
2.根据权利要求1所述的一种三维扫描点云实时法向量计算方法,其特征在于,上述步骤(4)所述的设置数据点的法向量评估置信度,利用上述步骤(3)所述的数据点的各个邻域点的法向量加权均值n'对数据点的法向量n进行评估,其评估值为置信度,记为cosθ,其评估公式(7)如下:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410327054.9A CN104123746B (zh) | 2014-07-10 | 2014-07-10 | 一种三维扫描点云中实时法向量的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410327054.9A CN104123746B (zh) | 2014-07-10 | 2014-07-10 | 一种三维扫描点云中实时法向量的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104123746A CN104123746A (zh) | 2014-10-29 |
CN104123746B true CN104123746B (zh) | 2017-07-25 |
Family
ID=51769143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410327054.9A Expired - Fee Related CN104123746B (zh) | 2014-07-10 | 2014-07-10 | 一种三维扫描点云中实时法向量的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104123746B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106327558B (zh) * | 2015-07-08 | 2019-11-19 | 深圳市腾讯计算机***有限公司 | 点云立面提取方法和装置 |
CN105184855B (zh) * | 2015-08-25 | 2018-09-11 | 广州市城市规划勘测设计研究院 | 基于三维点云的特征面构建方法及装置 |
KR101888963B1 (ko) * | 2017-03-06 | 2018-08-17 | (주)오앤드리메디컬로봇 | 레이저 치료를 위한 영역 분할 방법, 그를 이용한 레이저 치료 방법 및 장치 |
CN107463871A (zh) * | 2017-06-19 | 2017-12-12 | 南京航空航天大学 | 一种基于角特征加权的点云匹配方法 |
CN107292956B (zh) * | 2017-07-12 | 2020-09-22 | 杭州电子科技大学 | 一种基于曼哈顿假设的场景重建方法 |
CN108520550B (zh) * | 2018-04-10 | 2021-09-24 | 南京信息工程大学 | 基于噪声分类与mls的点云数据泊松曲面重建方法 |
CN109002418B (zh) * | 2018-06-20 | 2020-08-11 | 厦门大学 | 基于体素生长和地面激光点云的树木胸径自动计算方法 |
CN110070611B (zh) * | 2019-04-22 | 2020-12-01 | 清华大学 | 一种基于深度图像融合的人脸三维重建方法和装置 |
CN110782531A (zh) * | 2019-09-16 | 2020-02-11 | 华为技术有限公司 | 处理三维点云数据的方法和计算设备 |
CN112102178A (zh) * | 2020-07-29 | 2020-12-18 | 深圳市菲森科技有限公司 | 一种点云保特征去噪方法、装置、电子设备及存储介质 |
CN114930402A (zh) * | 2020-11-12 | 2022-08-19 | 深圳元戎启行科技有限公司 | 点云法向量计算方法、装置、计算机设备和存储介质 |
CN113129402B (zh) * | 2021-04-19 | 2024-01-30 | 中国航发沈阳发动机研究所 | 一种截面数据云图绘制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101751695A (zh) * | 2008-12-10 | 2010-06-23 | 中国科学院自动化研究所 | 点云数据的主曲率和主方向估计方法 |
CN102890828A (zh) * | 2012-06-15 | 2013-01-23 | 北京理工大学 | 基于法向夹角的点云数据精简方法 |
CN103065354A (zh) * | 2012-12-24 | 2013-04-24 | 中国科学院深圳先进技术研究院 | 点云优化方法及其装置 |
-
2014
- 2014-07-10 CN CN201410327054.9A patent/CN104123746B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101751695A (zh) * | 2008-12-10 | 2010-06-23 | 中国科学院自动化研究所 | 点云数据的主曲率和主方向估计方法 |
CN102890828A (zh) * | 2012-06-15 | 2013-01-23 | 北京理工大学 | 基于法向夹角的点云数据精简方法 |
CN103065354A (zh) * | 2012-12-24 | 2013-04-24 | 中国科学院深圳先进技术研究院 | 点云优化方法及其装置 |
Non-Patent Citations (1)
Title |
---|
On the normal vector estimation for point cloud data from smooth surfaces;DaoShan OuYang et al.;《Computer-Aided Design》;20051231;第1071-1079页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104123746A (zh) | 2014-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104123746B (zh) | 一种三维扫描点云中实时法向量的计算方法 | |
JP5772821B2 (ja) | 顔特徴点位置補正装置、顔特徴点位置補正方法および顔特徴点位置補正プログラム | |
US20190065885A1 (en) | Object detection method and system | |
JP5297403B2 (ja) | 位置姿勢計測装置、位置姿勢計測方法、プログラムおよび記憶媒体 | |
JP6816058B2 (ja) | パラメータ最適化装置、パラメータ最適化方法、プログラム | |
US9317735B2 (en) | Information processing apparatus, information processing method, and program to calculate position and posture of an object having a three-dimensional shape | |
JP6349418B2 (ja) | 高精度単眼移動によるオブジェクト位置特定 | |
CN104517289B (zh) | 一种基于混合摄像机的室内场景定位方法 | |
US20170161911A1 (en) | System and method for improved distance estimation of detected objects | |
CN111144213B (zh) | 一种对象检测方法和相关设备 | |
TW201635198A (zh) | 用於人臉邊緣特徵點定位的方法及裝置 | |
JP2013045433A (ja) | 学習装置、学習装置の制御方法、検出装置、検出装置の制御方法、およびプログラム | |
CN105678241B (zh) | 一种级联的二维图像人脸姿态估计方法 | |
WO2021051526A1 (zh) | 多视图3d人体姿态估计方法及相关装置 | |
WO2019196476A1 (zh) | 基于激光传感器生成地图 | |
CN111998862B (zh) | 一种基于bnn的稠密双目slam方法 | |
US20140168204A1 (en) | Model based video projection | |
JP2013205175A (ja) | 3次元対象面認識装置および方法、ならびにプログラム | |
JP6922348B2 (ja) | 情報処理装置、方法、及びプログラム | |
KR102126722B1 (ko) | 3차원 피사체 인식의 안티-스푸핑 방법 | |
JP6086491B2 (ja) | 画像処理装置およびそのデータベース構築装置 | |
US10440350B2 (en) | Constructing a user's face model using particle filters | |
CN104182971B (zh) | 一种高精度图像矩定位方法 | |
CN113793379A (zh) | 相机姿态求解方法及***、设备和计算机可读存储介质 | |
Rhee et al. | Two-pass icp with color constraint for noisy rgb-d point cloud registration |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170725 Termination date: 20200710 |
|
CF01 | Termination of patent right due to non-payment of annual fee |