CN106651752A - 三维点云数据配准方法及拼接方法 - Google Patents
三维点云数据配准方法及拼接方法 Download PDFInfo
- Publication number
- CN106651752A CN106651752A CN201610855932.3A CN201610855932A CN106651752A CN 106651752 A CN106651752 A CN 106651752A CN 201610855932 A CN201610855932 A CN 201610855932A CN 106651752 A CN106651752 A CN 106651752A
- Authority
- CN
- China
- Prior art keywords
- cloud data
- centre
- point cloud
- cloud
- target
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000000007 visual effect Effects 0.000 claims abstract description 57
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 42
- 230000009466 transformation Effects 0.000 claims description 32
- 230000011218 segmentation Effects 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 6
- 238000003384 imaging method Methods 0.000 abstract description 2
- 230000004075 alteration Effects 0.000 abstract 1
- 241001269238 Data Species 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000013450 outlier detection Methods 0.000 description 4
- 230000007704 transition Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000004888 barrier function Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明涉及一种三维点云数据配准方法及拼接方法。包括:获取三维扫描装置对设有包括多个标靶球的标定靶目标空间分别以第一视角、第二视角扫描对应得到的第一点云数据、第二点云数据,三维扫描装置分别以第一视角、第二视角扫描时标定靶均放置于第一视角和第二视角之间的公共扫描区域的同一个位置;从第一点云数据、第二点云数据中分别找出第一球心点云、第二球心点云,第一球心点云、第二球心点云均包含标定靶内所有标靶球的球心坐标;根据第一球心点云与第二球心点云之间的转换关系将第一点云数据与第二点云数据配准到同一坐标系下。该方法利用标定靶作为两个视角扫描数据匹配的桥梁,标靶球面的成像畸变程度小于平面,提高了配准的精确度。
Description
技术领域
本发明涉及土木工程技术领域,特别是涉及一种三维点云数据配准方法及拼接方法。
背景技术
房屋测量是建筑工程、家装设计工程重要的组成部分,是一项专业性很强的工作,如果在建筑面积计算上出现问题,将对预算造价的准确性产生直接影响,从而带来大量人力、物力的损失。三维激光扫描仪是一种非接触式海量高精度数据的获取手段,可以获得房屋内部的数字模型,然而房屋内部各房间被墙体隔断,每个房间是相对独立的空间,使得三维激光扫描仪得到的各房间的三维点云数据之间的配准比较困难。
传统的三维点云配准方法通常都是通过被测量对象表面粘贴标记点的方法完成,然而此类方法存在的缺点为:平面标记在二维空间存在很大畸变,导致识别难度较大,定位精度不高,从而降低配准的精确度。
发明内容
基于此,有必要针对如何改善传统三维点云配准方法配准精确度较低的问题,提供一种三维点云数据配准方法及拼接方法。
一种三维点云数据配准方法,包括:
获取三维扫描装置对目标空间分别以第一视角、第二视角扫描对应得到的第一点云数据、第二点云数据;所述目标空间设有包括多个标靶球的标定靶;所述三维扫描装置分别以第一视角、第二视角扫描时,所述标定靶均放置于所述第一视角和第二视角之间的公共扫描区域的同一个位置;
从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云;所述第一球心点云、所述第二球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到同一坐标系下。
在其中一个实施例中,还包括:
获取所述目标空间在所述第一视角、所述第二视角下对应的第一二维图像数据、第二二维图像数据;
同时,从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云的步骤包括:
从所述第一二维图像数据、所述第二二维图像数据中根据设定的标靶球判断准则分别检测出所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标;其中,从所述第一二维图像数据中检测的所有所述疑似标靶球、从所述第二二维图像数据中检测的所有所述疑似标靶球均至少包括所述标定靶内的所有标靶球;
根据所有所述第一二维圆心坐标从所述第一点云数据找出包含所述标定靶内所有所述标靶球球心坐标的第一球心点云;
根据所有所述第二二维圆心坐标从所述第二点云数据中找出包含所述标定靶内所有所述标靶球球心坐标的第二球心点云。
在其中一个实施例中,所述第一二维图像数据包括第一灰度图像和第一深度图像;所述第二二维图像数据包括第二灰度图像和第二深度图像;
同时,从所述第一二维图像数据、所述第二二维图像数据中根据设定的标靶球判断准则分别检测出所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标的步骤包括:
分别对所述第一灰度图像、所述第二灰度图像进行阈值分割,对应得到第一灰度二值化图像、第二灰度二值化图像;
分别计算所述第一灰度二值化图像、所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标;
判断所述第一灰度二值化图像内疑似标靶球的数量小于所述标定靶内标靶球的总数量时,利用局部阈值分割的方法对所述第一深度图像进行分割得到第一深度二值化图像,并计算所述第一深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第一二维圆心坐标;
判断所述第二灰度二值化图像内疑似标靶球的数量小于所述标靶球的总数量时,利用局部阈值分割的方法对所述第二深度图像进行分割得到第二深度二值化图像,并计算所述第二深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第二二维圆心坐标。
在其中一个实施例中,分别计算所述第一灰度二值化图像、所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标的步骤包括:
分别对所述第一灰度二值化图像、所述第二灰度二值化图像进行滤波;
分别计算滤波后的所述第一灰度二值化图像、滤波后的所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标。
在其中一个实施例中,根据所有所述第一二维圆心坐标从所述第一点云数据找出包含所述标定靶内所有所述标靶球球心坐标的第一球心点云的步骤包括:
根据各所述第一二维圆心坐标从所述第一点云数据找出相应的第一三维球面点云,并拟合出各所述第一三维球面点云的球心;
对所有所述第一三维球面点云的球心进行离群点处理,得到包含所述标定靶内所有所述标靶球球心坐标的第一球心点云。
在其中一个实施例中,根据所有所述第二二维圆心坐标从所述第二点云数据中找出包含所述标定靶内所有所述标靶球球心坐标的第二球心点云的步骤包括:
根据各所述第二二维圆心坐标从所述第二点云数据找出相应的第二三维球面点云,并拟合出各所述第二三维球面点云的球心;
对所有所述第二三维球面点云的球心进行离群点处理,得到包含所述标定靶内所有所述标靶球球心坐标的第二球心点云。
在其中一个实施例中,根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到同一坐标系下的步骤包括:
利用配准算法计算将所述第二球心点云配准到所述第一球心点云所在坐标系下的第一变换矩阵;
根据所述第一变换矩阵对所述第二点云数据进行转换;
利用配准算法计算转换后的所述第二点云数据配准到所述第一点云数据所在坐标下的第二变换矩阵;
根据所述第一变换矩阵及所述第二变换矩阵将所述第二点云数据配准到所述第一点云数据所在坐标系下。
在其中一个实施例中,根据所述第一变换矩阵对所述第二点云数据进行转换的步骤之前还包括:
分别对所述第一点云数据、所述第二点云数据进行抽稀;
同时,根据所述第一变换矩阵对所述第二点云数据进行转换的步骤为:
根据所述第一变换矩阵对抽稀后的所述第二点云数据进行转换,得到转换后的第二抽稀点云数据;
同时,利用配准算法计算转换后的所述第二点云数据配准到所述第一点云数据所在坐标下的第二变换矩阵的步骤为:
利用配准算法计算所述转换后的第二抽稀点云数据配准到抽稀后的所述第一点云数据所在坐标系下的第二变换矩阵。
在其中一个实施例中,所述标定靶还包括底座及多根连接杆;各所述连接杆的一端安装于所述底座上,且各所述连接杆的另一端分别与相应的所述标靶球连接;其中,所述连接杆的颜色为黑色,所述标靶球的颜色为白色。
一种三维点云数据拼接方法,包括:
获取三维扫描装置分别对依次相邻的第一目标空间、第二目标空间、第三目标空间、第四目标空间扫描得到的第一点云数据、第二点云数据、第三点云数据、第四点云数据;其中,所述第一点云数据与所述第二点云数据分别由所述三维扫描装置以不同视角扫描得出,且均包括位于同一位置的标定靶的点云数据,其中,所述标定靶包括多个标靶球;所述第二点云数据和所述第三点云数据为相同的点云;所述第三点云数据和所述第四点云数据分别由所述三维扫描装置以不同视角扫描得出,且均包括位于同一位置的所述标定靶的点云数据;
从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云;所述第一球心点云、所述第二球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到第一坐标系下;
从所述第三点云数据、所述第四点云数据中分别找出第三球心点云、第四球心点云;所述第三球心点云、所述第四球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第三球心点云与第四球心点云之间的转换关系,将所述第四点云数据配准到所述第三点云数据所在的坐标系下;
利用所述第三点云数据与所述第二点云数据之间的关系,将配准后的所述第四点云数据配准到所述第一坐标系下;
将配准至所述第一坐标下的第一点云数据、配准至所述第一坐标下的第二点云数据及配准至所述第一坐标下的第四点云数据拼接在一起。
上述三维点云数据配准方法及拼接方法具有的有益效果为:其中,三维扫描装置分别以第一视角、第二视角扫描时,标定靶均放置于第一视角和第二视角的公共扫描区域的同一个位置,因此第一点云数据和第二点云数据分别包括位于各自坐标系下且实际空间位置相同的两组标定靶的点云数据。基于上述关系,之后从第一点云数据、第二点云数据中分别找出第一球心点云、第二球心点云,即可根据第一球心点云与第二球心点云之间的转换关系将第一点云数据与第二点云数据配准到同一坐标系下。因此该三维点云数据配准方法及拼接方法利用标定靶作为两个视角扫描数据进行匹配的桥梁,由于标靶球面的成像畸变程度远小于平面,因此提高了配准的精确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他实施例的附图。
图1为一实施例提供的三维点云数据配准方法的执行流程图;
图2为图1所示实施例的三维点云数据配准方法涉及的标定靶的立体示意图;
图3(a)、图3(b)、图3(c)、图3(d)、图3(e)、图3(f)分别为图2所示实施例的标定靶的仰视图、俯视图、右视图、左视图、正视图、后视图;
图4为图1所示实施例的三维点云数据配准方法的其中一种具体实施方式流程图;
图5为图4所示实施例的三维点云数据配准方法中步骤S310的其中一种具体实施方式流程图;
图6为图5所示实施例的三维点云数据配准方法中步骤S312的其中一种具体实施方式流程图;
图7为图4所示实施例的三维点云数据配准方法中步骤S320的其中一种具体实施方式流程图;
图8为图4所示实施例的三维点云数据配准方法中步骤S330的其中一种具体实施方式流程图;
图9为图1所示实施例的三维点云数据配准方法中步骤S400的其中一种具体实施方式流程图;
图10为另一实施例提供的三维点云数据拼接方法的执行流程图。
具体实施方式
为了便于理解本发明,下面将参照相关附图对本发明进行更全面的描述。附图中给出了本发明的较佳实施例。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
除非另有定义,本文所使用的所有的技术和科学术语与属于发明的技术领域的技术人员通常理解的含义相同。本文中在发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
一实施例提供了一种三维点云数据配准方法,用于将多视角的点云数据配准至同一坐标系下。该三维点云数据配准方法可以由电脑或其他具有数据处理能力的智能设备来执行,如图1所示,包括以下步骤。
步骤S100.获取三维扫描装置对目标空间分别以第一视角、第二视角扫描对应得到的第一点云数据、第二点云数据。其中,该目标空间设有包括多个标靶球的标定靶。三维扫描装置分别以第一视角、第二视角扫描时,该标定靶均放置于第一视角和第二视角之间的公共扫描区域的同一个位置。
其中,三维扫描装置例如为三维激光扫描仪。目标空间例如为房屋。第一视角和第二视角为两个不同的视角,故第一点云数据和第二点云数据分别处于各自的坐标系下。在具体扫描过程中,可以先将三维扫描装置置于目标空间内与第一视角对应的位置,同时将标定靶置于第一视角和第二视角之间的公共扫描区域,这时三维扫描装置扫描得到的数据即为第一点云数据;之后,保持标定靶的位置不变,将三维扫描装置移动至与第二视角对应的另一个位置,这时三维扫描装置扫描得到的数据即为第二点云数据。在上述两次扫描过程中,由于标定靶均放置于第一视角和第二视角的公共扫描区域的同一个位置,因此三维扫描装置在该两次扫描过程中均能扫描到位于同一位置的标定靶,故第一点云数据、第二点云数据分别包括位于各自坐标系下且实际位置相同的两组标定靶的点云数据。
标定靶的结构如图2所示,本实施例中,标定靶包括4个标靶球,即第一标靶球111、第二标靶球112、第三标靶球113、第四标靶球114。各标靶球的直径可以为16cm或其他数值。此外,标定靶还包括底座130及多根连接杆。其中,各连接杆的一端安装于底座130上,且各连接杆的另一端分别与相应的标靶球连接。本实施例中共包括4跟连接杆,即第一连接杆121、第二连接杆122、第三连接杆123、第四连接杆124,且分别对应连接第一标靶球111、第二标靶球112、第三标靶球113、第四标靶球114。
本实施例中,标靶球的颜色与连接杆及底座130的颜色不同,以提高标靶球的辨识度,更便于从第一点云数据和第二点云数据中提取出各标靶球的点云数据。例如,各标靶球的颜色均为白色,连接杆及底座130的颜色均为黑色。
上述4根连接杆的其中一种空间位置关系为,请参考图3(a)至图3(f):第一连接杆121、第二连接杆122、第三连接杆123、第四连接杆124的长度分别为400mm、250mm、140mm、200mm。进一步的,在图3(a)中,设定该仰视图对应XOZ平面,则第一连接杆121沿X轴负方向顺时针旋转的角度a为85°,第二连接杆122沿X轴正方向逆时针旋转的角度b为84°,第三连接杆123沿X轴正方向逆时针旋转的角度c为30°,第四连接杆124沿X轴负方向逆时针旋转的角度d为60°。在图3(d)中,设定该左视图对应ZOY平面,则第二连接杆122沿Z轴正方向顺时针旋转的角度f为30°,其余的连接杆均与Z轴平行。在图3(e)中,设定该正视图对应XOY平面,则第二连接杆122沿X轴正方向顺时针旋转的角度e为80°,其余的连接杆均与X轴平行。因此,标定靶的上述设置方式能够减少各标靶球之间的相互遮挡部位,提高了从不同角度都能看到4个标靶球的可靠性,从而进一步增强了配准的准确度。可以理解的是,标定靶的结构设置方式不限于上述情况,只要能够减少各标靶球之间的相互遮挡部位以满足配准的精度要求即可。
步骤S300.从上述第一点云数据、第二点云数据中分别找出第一球心点云、第二球心点云。其中,第一球心点云、第二球心点云均包含标定靶内所有标靶球的球心坐标。
该步骤换言之,从第一点云数据中找出第一球心点云,且第一球心点云包括标定靶内所有标靶球在第一视角对应坐标系下的球心坐标。从第二点云数据中找出第二球心点云,且第二球心点云包括标定靶内所有标靶球在第二视角对应坐标系下的球心坐标。同时,第一球心点云中任一个标靶球的球心坐标与第二球心点云中相同标靶球的球心坐标之间的区别仅在于对应的坐标系不同,而两者的实际空间位置相同。
步骤S400.根据第一球心点云与第二球心点云之间的转换关系将第一点云数据与第二点云数据配准到同一坐标系下。
由于第一球心点云与第二球心点云为实际空间位置相同的同一特征(即标定靶)分别在两个坐标系下的数据,因此将第一球心点云与第二球心点云进行匹配,可得到相应的旋转和平移矩阵,进而即可将第一点云数据与第二点云数据配准到同一坐标系下。
因此,上述三维点云数据配准方法利用标定靶作为两个视角扫描数据进行匹配的桥梁,由于标靶球面的成像畸变程度远小于平面,因此提高了配准的精确度。另外,在进行多次测量时,只需移动标定靶的位置即可,操作简单。
进一步的,在上述三维点云数据配准方法的其中一种具体实施方式中,还包括以下步骤,请参考图4。
步骤S200.获取上述目标空间在上述第一视角、上述第二视角下对应的第一二维图像数据、第二二维图像数据。
其中,第一二维图像数据、第二二维图像数据均为二维图像数据,例如灰度图像。第一二维图像数据、第二二维图像数据可以直接由该智能设备从三维扫描装置的扫描数据中提取出来,或者由其他装置另外对目标空间进行扫描得出,并将得出的第一二维图像数据、第二二维图像数据发送至该智能设备。
基于步骤S200,上述步骤S300具体包括以下步骤,请继续参考图4。
步骤S310.从上述第一二维图像数据、上述第二二维图像数据中根据设定的标靶球判断准则分别检测出所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标。其中,从第一二维图像数据中检测的所有疑似标靶球至少包括标定靶内的所有标靶球。从第二二维图像数据中检测的所有疑似标靶球至少包括标定靶内的所有标靶球。
该步骤换言之,从第一二维图像数据中根据设定的标靶球判断准则检测出所有疑似标靶球的第一二维圆心坐标,且第一二维圆心坐标为二维坐标数据;从第二二维图像数据中根据设定的标靶球判断准则检测出所有疑似标靶球的第二二维圆心坐标,且第二二维圆心坐标同样为二维坐标数据。其中,从第一二维图像数据中检测出的所有疑似标靶球除了包括标定靶内所有真实的标靶球之外,还可能包括其他非球类对象。同理,从第二二维图像数据中检测出的所有疑似标靶球除了包括标定靶内所有真实的标靶球之外,还可能包括其他非球类对象。
步骤S320.根据所有第一二维圆心坐标从第一点云数据找出包含标定靶内所有标靶球球心坐标的第一球心点云。
由于步骤S310中检测出的疑似标靶球中可能包括其他非球类对象,因此该步骤中,可以对从第一二维图像数据中检测出的所有疑似标靶球,在三维空间内通过拟合、离群点处理等措施,从而最终找出准确的标靶球球心坐标,完成三维空间内标靶球的准确定位。
步骤S330.根据所有第二二维圆心坐标从第二点云数据中找出包含标定靶内所有标靶球球心坐标的第二球心点云。
由于步骤S310中检测出的疑似标靶球中可能包括其他非球类对象,因此该步骤中,可以对从第二二维图像数据中检测出的所有疑似标靶球,在三维空间内通过拟合、离群点处理等措施,从而最终找出准确的标靶球球心坐标,完成三维空间内标靶球的准确定位。
由于二维数据与三维数据相比,运算量较少,因此步骤S300的上述实施方式根据二维数据来检测标靶球,减少了运算量,从而加快了运算速度。
具体的,上述第一二维图像数据包括第一灰度图像和第一深度图像。上述第二二维图像数据包括第二灰度图像和第二深度图像。其中,第一灰度图像和第二灰度图像均为灰度图像,灰度图像反映了物体的部分纹理和空间关系。第一深度图像和第二深度图像均为深度图像,深度图像反映了物体的空间关系和物体到三维扫描装置的距离。在此基础上,上述步骤S310的其中一种具体实施方式包括以下步骤,请参考图5。
步骤S311.分别对第一灰度图像、第二灰度图像进行阈值分割,对应得到第一灰度二值化图像、第二灰度二值化图像。
该步骤换言之,对第一灰度图像进行阈值分割得到第一灰度二值化图像;对第二灰度图像进行阈值分割得到第二灰度二值化图像。其中,阈值分割的基本原理是:通过设定不同的特征阈值,把图像象素点分为若干类。第一灰度二值化图像和第二灰度二值化图像均为二值化图像,二值化图像是指图像上的像素点的灰度值为0或255,呈现出明显的只有黑和白的视觉效果。本实施例中,由于标靶球在灰度图像中具有较高的亮度,因此通过阈值分割的方式易于将标靶球与其他特征分割开来。
步骤S312.分别计算上述第一灰度二值化图像、第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标。
该步骤换言之,计算上述第一灰度二值化图像内各连通区域的圆形度,以检测出各疑似标靶球的第一二维圆心坐标;计算上述第二灰度二值化图像内各连通区域的圆形度,以检测出各疑似标靶球的第二二维圆心坐标。
其中,连通区域是指没有“洞”的区域。圆形度的计算公式为面积乘上4π除以周长的平方。上述根据圆形度而采取的检测方式可以为:设定与标靶球对应的合适阈值,当某一连通区域的圆形度大于该阈值时,则认为是疑似标靶球,进而提取出相应的圆心坐标,即第一二维圆心坐标。
具体的,步骤S312还包括以下步骤,请参考图6。
步骤S312a.分别对第一灰度二值化图像、第二灰度二值化图像进行滤波。
其中,可以利用形态学滤波或其他滤波方式对第一灰度二值化图像、第二灰度二值化图像进行滤波,以去除噪声和干扰点。
步骤S312b.分别计算滤波后的第一灰度二值化图像、滤波后的第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标。
该步骤换言之,计算滤波后的第一灰度二值化图像内各连通区域的圆形度,以检测出各疑似标靶球的第一二维圆心坐标。计算滤波后的第二灰度二值化图像内各连通区域的圆形度,以检测出各疑似标靶球的第二二维圆心坐标。
接下来请继续参考图5,步骤S310还包括以下具体步骤:
步骤S313.判断上述第一灰度二值化图像内疑似标靶球的数量是否小于标定靶内标靶球的总数量,若是,执行步骤S314;否则,执行步骤S320。
如果疑似标靶球的数量小于标定靶内标靶球的总数量,代表第一灰度二值化图像中的干扰因素较多,导致丢失了一些标靶球特征。因此,本实施例针对这种情况,进一步通过步骤S314来提高检测的精确性。
步骤S314.利用局部阈值分割的方法对上述第一深度图像进行分割得到第一深度二值化图像。
其中,局部阈值分割是指将图像分割成若干子图像进行处理。由于深度图像反映了物体到三维扫描装置的距离,且标靶球所在的局部空间干扰较少,那么如果去除标靶球所在空间位置外的前景和背景,则能去除大部分干扰,因此可以利用步骤S313检测到的疑似标靶球圆心到三维扫描装置的距离信息,设置一个阈值distance1,然后将第一深度图像中距离大于distance1的像素设置为0,从而仅保留了包含疑似标靶球的局部子图像。之后对该局部子图像进行阈值分割,即可得出第一深度二值化图像。
步骤S315.计算上述第一深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第一二维圆心坐标。
其中,其余的疑似标靶球是指第一灰度二值化图像由于干扰因素而丢失的标靶球特征。
步骤S316.判断上述第二灰度二值化图像内疑似标靶球的数量是否小于标定靶内标靶球的总数量,若是,执行步骤S317;否则,执行步骤S330。
如果疑似标靶球的数量小于标定靶内标靶球的总数量,代表第二灰度二值化图像中的干扰因素较多,导致丢失了一些标靶球特征。因此,本实施例针对这种情况,进一步通过步骤S317来提高检测的精确性。
步骤S317.利用局部阈值分割的方法对上述第二深度图像进行分割得到第二深度二值化图像。
其中,局部阈值分割是指将图像分割成若干子图像进行处理。由于深度图像反映了物体到激光束的距离,且标靶球所在的局部空间干扰较少,那么如果去除标靶球所在空间位置外的前景和背景,则能去除大部分干扰,因此可以利用步骤S316检测到的疑似标靶球圆心到三维扫描装置的距离信息,设置一个阈值distance2,然后将第二深度图像中距离大于distance2的像素设置为0,从而仅保留了包含疑似标靶球的局部子图像。之后对该局部子图像进行阈值分割,即可得出第二深度二值化图像。
步骤S318.计算上述第二深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第二二维圆心坐标。
其中,其余的疑似标靶球是指第二灰度二值化图像由于干扰因素而丢失的标靶球特征。
因此,在上述步骤S310的具体实施方式中,通过灰度图像和深度图像来共同检测疑似标靶球,从而能够去除干扰点,提高抗干扰能力,增强了该三维点云数据配准方法的鲁棒性。可以理解的是,步骤S310的具体实现方式不限于上述情况,只要能够检测所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标即可。
具体的,上述步骤S320的其中一种具体实施方式包括以下步骤,请参考图7。
步骤S321.根据各第一二维圆心坐标从第一点云数据找出相应的第一三维球面点云,并拟合出各第一三维球面点云的球心。
其中,第一三维球面点云是对应某一疑似标靶球的球面点云。不同的第一二维圆心坐标可以找出不同的第一三维球面点云。具体的,可以先将各第一二维圆心坐标映射到三维空间,再利用标靶球的半径约束原理从第一点云数据中找出相应的球面点云。具体方式例如:对于某一第一二维圆心坐标,以其映射到三维空间的点为中心,在第一点云数据中寻找位于(r-1)cm到(r+1)cm范围内的点云(其中,r为标靶球的半径),找出的点云即为该第一二维圆心坐标对应的第一三维球面点云。具体的,可以利用采样一致性算法拟合出各第一三维球面点云的球心,从而可以得到相应的球心坐标。
步骤S322.对所有第一三维球面点云的球心进行离群点处理,得到包含标定靶内所有标靶球球心坐标的第一球心点云。
由于步骤S310中得出的所有第一二维圆心坐标中可能有错误值,这些错误的第一二维圆心坐标通过步骤S321拟合出的第一三维球面点云的球心则不是真正的标靶球球心,即为离群点,因此通过步骤S322即可去掉这些离群点。离群点处理的具体实现方式例如为:先采用基于统计的离群点检测、基于距离的离群点检测或其他离群点检测方法检测出各离群点,然后去除这些离群点,从而得出包含标定靶内所有标靶球球心坐标的第一球心点云。
可以理解的是,步骤S320的具体实现方式不限于上述情况,只要能够根据所有所述第一二维圆心坐标从第一点云数据找出包含标定靶内所有标靶球球心坐标的第一球心点云即可。
具体的,上述步骤S330的其中一种具体实施方式包括以下步骤,请参考图8。
步骤S331.根据各第二二维圆心坐标从第二点云数据找出相应的第二三维球面点云,并拟合出各第二三维球面点云的球心。
其中,第二三维球面点云是对应某一疑似标靶球的球面点云。不同的第二二维圆心坐标可以找出不同的第二三维球面点云。具体的,可以先将各第二二维圆心坐标映射到三维空间,再利用标靶球的半径约束原理从第二点云数据中找出相应的球面点云。具体方式例如:对于某一第二二维圆心坐标,以其映射到三维空间的点为中心,在第二点云数据中寻找位于(r-1)cm到(r+1)cm范围内的点云(其中,r为标靶球的半径),找出的点云即为该第二二维圆心坐标对应的第二三维球面点云。具体的,可以利用采样一致性算法拟合出各第二三维球面点云的球心,从而可以得到相应的球心坐标。
步骤S332.对所有第二三维球面点云的球心进行离群点处理,得到包含标定靶内所有标靶球球心坐标的第二球心点云。
由于步骤S310中得出的所有第二二维圆心坐标中可能有错误值,这些错误的第二二维圆心坐标通过步骤S331拟合出的第二三维球面点云的球心则不是真正的标靶球球心,即为离群点,因此通过步骤S332即可去掉这些离群点。离群点处理的具体实现方式例如为:先采用基于统计的离群点检测、基于距离的离群点检测或其他离群点检测方法检测出各离群点,然后去除这些离群点,从而得出包含标定靶内所有标靶球球心坐标的第二球心点云。
可以理解的是,步骤S330的具体实现方式不限于上述情况,只要能够根据所有第二二维圆心坐标从第二点云数据找出包含标定靶内所有标靶球球心坐标的第二球心点云即可。
具体的,步骤S400的其中一种具体实施方式包括以下步骤,请参考图9。
步骤S410.利用配准算法计算将第二球心点云配准到第一球心点云所在坐标系下的第一变换矩阵。
该步骤相当于进行点云初始配准,用于缩小第一点云数据与第二点云数据之间的旋转和平移错位,使得后续的精确配准不至于趋向错误的方向。其中,配准算法例如为ICP(Iterative Closest Point,迭代就近点)配准算法。第一变换矩阵是指相应的旋转和平移矩阵。
步骤S420.根据上述第一变换矩阵对第二点云数据进行转换。
其中,转换后的第二点云数据为第二点云数据与第一变换矩阵的乘积。
步骤S430.利用配准算法计算转换后的第二点云数据配准到第一点云数据所在坐标下的第二变换矩阵。
该步骤相当于进行精确配准,即在上述初始配准的基础上使第二点云数据配准到第一点云数据所在坐标系的精度更高。其中,配准算法例如为ICP配准算法。
步骤S440.根据上述第一变换矩阵及第二变换矩阵将第二点云数据配准到第一点云数据所在坐标系下。
因此,第二点云数据配准到第一点云数据所在坐标系的配准转换矩阵为第一转换矩阵与第二转换矩阵的乘积,故最终配准后的第二点云数据为第二点云数据与该配准转换矩阵的乘积。
进一步的,上述步骤S400的具体实现方式中,在执行步骤S420之前,还可以先分别对第一点云数据、第二点云数据进行抽稀。抽稀是指在保证矢量曲线形状不变的情况下,最大限度地减少数据点个数,以去除重复的点。对第一点云数据、第二点云数据进行抽稀处理后,减少了运算量,从而进一步提高了运算速度。可以理解的是,若第一点云数据、第二点云数据中没有重复的点或重复点的个数小于设定阈值,也可省略上述抽稀的过程。
在经过上述抽稀处理后,步骤S420即为:根据上述第一变换矩阵对抽稀后的第二点云数据进行转换,得到转换后的第二抽稀点云数据。设置第一点云数据为pCloud1、第二点云数据为pCloud2、抽稀后的第一点云数据为pCloud1_1,抽稀后的第二点云数据为pCloud2_1,那么转换后的第二抽稀点云数据pCloud2_2=pCloud2_1*T1,其中,T1为第一变换矩阵。
之后,步骤S430即为:利用配准算法计算上述转换后的第二抽稀点云数据pCloud2_2配准到抽稀后的第一点云数据pCloud1_1所在坐标系下的第二变换矩阵T2。因此,pCloud2_1配准到pCloud1_1的最终转换矩阵为:T1*T2,相当于第二点云数据pCloud2配准到第一点云数据pCloud1的最终转换矩阵为:T1*T2。
可以理解的是,步骤S400的具体实现方式不限于上述情况,只要能够根据第一球心点云与第二球心点云之间的转换关系将第一点云数据与第二点云数据配准到同一坐标系下即可。例如,也可以将第一点云数据配准到第二点云数据所在的坐标系下。
另外,基于上述的三维点云数据配准方法,另一实施例提供了一种三维点云数据拼接方法,包括以下步骤,请参考图10。
步骤S510.获取三维扫描装置分别对依次相邻的第一目标空间、第二目标空间、第三目标空间、第四目标空间扫描得到的第一点云数据、第二点云数据、第三点云数据、第四点云数据。其中,第一点云数据与第二点云数据分别由三维扫描装置以不同视角扫描得出,且均包括位于同一位置的标定靶的点云数据,标定靶包括多个标靶球。第二点云数据和第三点云数据为相同的点云。第三点云数据和第四点云数据分别由三维扫描装置以不同视角扫描得出,且均包括位于同一位置的标定靶的点云数据。
该步骤中,第一目标空间、第二目标空间、第三目标空间、第四目标空间分别可以为依次相邻的四个房间。第一点云数据与第二点云数据均包括位于同一位置的标定靶的点云数据,换言之,三维扫描装置在对第一目标空间、第二目标空间分别进行扫描时,标定靶均位于第一目标空间和第二目标空间之间的公共扫描区域的同一个位置。同理,第三点云数据和第四点云数据分别由三维扫描装置均包括位于同一位置的标定靶的点云数据,换言之,三维扫描装置在对第三目标空间、第四目标空间分别进行扫描时,标定靶均位于第三目标空间和第四目标空间之间的公共扫描区域的同一个位置。
第二点云数据与第三点云数据为相同的点云数据,例如为以下情况:三维扫描装置在对第二目标空间、第三目标空间进行扫描时,其位于第二目标空间与第三目标空间之间的公共扫描区域的同一个位置,那么三维扫描装置则以相同的视角进行这两次扫描,且在每一次扫描中三维扫描装置均能同时扫描到第二目标空间和第三目标空间,从而使得第二点云数据与第三点云数据相同。
三维扫描装置实施上述四次扫描过程可按以下方式进行:先将三维扫描装置置于第一目标空间内一个合适的位置,且将标定靶置于第一目标空间和第二目标空间之间的公共扫描区域,这时三维扫描装置扫描后即可扫描得出第一点云数据;之后,将三维扫描装置移动至第二目标空间和第三目标空间之间的公共扫描区域,且保持标定靶的位置不变,这时三维扫描装置扫描后即可扫描得出第二点云数据;之后保持三维扫描装置的位置不变,将标定靶移动至第三目标空间和第四目标空间的公共扫描区域,这时三维扫描装置扫描后即可扫描得出第三点云数据;之后,再将三维扫描装置移动至第四目标空间内的一个合适位置,并保持标定靶的位置不变,这时三维扫描装置即可扫描得出第四点云数据。
步骤S520.从上述第一点云数据、第二点云数据中分别找出第一球心点云、第二球心点云。其中,第一球心点云、所述第二球心点云均包含标定靶内所有标靶球的球心坐标。
步骤S530.根据上述第一球心点云与上述第二球心点云之间的转换关系将第一点云数据与第二点云数据配准到第一坐标系下。
由于第一目标空间和第二目标空间的公共扫描区域的同一个位置设有上述标定靶,因此第一点云数据与第二点云数据之间的配准可采用图1所示实施例提供的三维点云数据配准方法进行配准,即步骤S520、步骤S530分别可以按照上一实施例中步骤S300、步骤S400的具体方式来执行,这里就不再赘述。
步骤540.从第三点云数据、第四点云数据中分别找出第三球心点云、第四球心点云。其中,第三球心点云、第四球心点云均包含标定靶内所有标靶球的球心坐标。
由于第三目标空间和第四目标空间的公共扫描区域的同一个位置也设有该标定靶,因此步骤S540同样可以采用上一实施例提供的三维点云数据配准方法中的步骤S300的具体方式执行,这里就不再赘述。
步骤S550.根据第三球心点云与第四球心点云之间的转换关系,将第四点云数据配准到第三点云数据所在的坐标系下。
该步骤可采用上一实施例中图9提供的方式来执行,这里就不再赘述。
步骤S560.利用第三点云数据与第二点云数据之间的关系,将配准后的第四点云数据配准到第一坐标系下。
由于第三点云数据与第二点云数据相同,因此即可将步骤S550配准后的第四点云数据配准的第一坐标系下。那么,该步骤最终即可将第一点云数据、第二点云数据、第四点云数据均配准至同一坐标系(即第一坐标系)下。
步骤S570.将配准至第一坐标系下的第一点云数据、配准至第一坐标系下的第二点云数据及配准至第一坐标系下的第四点云数据拼接在一起。
该步骤即可完成第一点云数据、第二点云数据、第四点云数据之间的拼接。需要说明的是,由于第二点云数据与第三点云数据相同,因此在第二点云数据和第三点云数据中只需选择一个点云数据即可。
为了更便于说明,接下以公式的方式描述上述拼接方法。
设定步骤S530中将第二点云数据配准到第一点云数据所在的坐标系下,设置第一点云数据为pCloud1、第二点云数据为pCloud2,那么,配准到第一点云数据所在的坐标系下的第二点云数据则为T1*T2*pCloud2,其中T1、T2分别为将第二点云数据配准到第一点云数据所在坐标下的第一变换矩阵、第二变换矩阵。
之后,设置第三点云数据为pCloud3、第四点云数据为pCloud4,那么,步骤S550执行完毕得到的配准到第三点云数据所在坐标系下的第四点云数据为T3*T4**pCloud4。其中T3、T4分别为将第四点云数据配准到第三点云数据所在坐标下的第一变换矩阵、第二变换矩阵。
之后,步骤S560执行完毕后得到的配准到第一坐标系下的第四点云数据为T1*T2*T3*T4*pCloud4。最后,步骤S570执行完毕得到的拼接后的数据为pCloud=pCloud1+T1*T2*pCloud2+T1*T2*T3*T4*pCloud4。
图1、图4至图10为本发明实施例的方法的流程示意图。应该理解的是,虽然图1、图4至图10的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,其可以以其他的顺序执行。而且,图1、图4至图10中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,其执行顺序也不必然是依次进行,而是可以与其他步骤或者其他步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种三维点云数据配准方法,包括:
获取三维扫描装置对目标空间分别以第一视角、第二视角扫描对应得到的第一点云数据、第二点云数据;所述目标空间设有包括多个标靶球的标定靶;所述三维扫描装置分别以第一视角、第二视角扫描时,所述标定靶均放置于所述第一视角和第二视角之间的公共扫描区域的同一个位置;
从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云;所述第一球心点云、所述第二球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到同一坐标系下。
2.根据权利要求1所述的三维点云数据配准方法,其特征在于,还包括:
获取所述目标空间在所述第一视角、所述第二视角下对应的第一二维图像数据、第二二维图像数据;
同时,从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云的步骤包括:
从所述第一二维图像数据、所述第二二维图像数据中根据设定的标靶球判断准则分别检测出所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标;其中,从所述第一二维图像数据中检测的所有所述疑似标靶球、从所述第二二维图像数据中检测的所有所述疑似标靶球均至少包括所述标定靶内的所有标靶球;
根据所有所述第一二维圆心坐标从所述第一点云数据找出包含所述标定靶内所有所述标靶球球心坐标的第一球心点云;
根据所有所述第二二维圆心坐标从所述第二点云数据中找出包含所述标定靶内所有所述标靶球球心坐标的第二球心点云。
3.根据权利要求2所述的三维点云数据配准方法,其特征在于,所述第一二维图像数据包括第一灰度图像和第一深度图像;所述第二二维图像数据包括第二灰度图像和第二深度图像;
同时,从所述第一二维图像数据、所述第二二维图像数据中根据设定的标靶球判断准则分别检测出所有疑似标靶球的第一二维圆心坐标、第二二维圆心坐标的步骤包括:
分别对所述第一灰度图像、所述第二灰度图像进行阈值分割,对应得到第一灰度二值化图像、第二灰度二值化图像;
分别计算所述第一灰度二值化图像、所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标;
判断所述第一灰度二值化图像内疑似标靶球的数量小于所述标定靶内标靶球的总数量时,利用局部阈值分割的方法对所述第一深度图像进行分割得到第一深度二值化图像,并计算所述第一深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第一二维圆心坐标;
判断所述第二灰度二值化图像内疑似标靶球的数量小于所述标靶球的总数量时,利用局部阈值分割的方法对所述第二深度图像进行分割得到第二深度二值化图像,并计算所述第二深度二值化图像内各连通区域的圆形度,以检测出其余疑似标靶球的第二二维圆心坐标。
4.根据权利要求3所述的三维点云数据配准方法,其特征在于,分别计算所述第一灰度二值化图像、所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标的步骤包括:
分别对所述第一灰度二值化图像、所述第二灰度二值化图像进行滤波;
分别计算滤波后的所述第一灰度二值化图像、滤波后的所述第二灰度二值化图像内各连通区域的圆形度,以对应检测出各疑似标靶球的第一二维圆心坐标、第二二维圆心坐标。
5.根据权利要求2所述的三维点云数据配准方法,其特征在于,根据所有所述第一二维圆心坐标从所述第一点云数据找出包含所述标定靶内所有所述标靶球球心坐标的第一球心点云的步骤包括:
根据各所述第一二维圆心坐标从所述第一点云数据找出相应的第一三维球面点云,并拟合出各所述第一三维球面点云的球心;
对所有所述第一三维球面点云的球心进行离群点处理,得到包含所述标定靶内所有所述标靶球球心坐标的第一球心点云。
6.根据权利要求2所述的三维点云数据配准方法,其特征在于,根据所有所述第二二维圆心坐标从所述第二点云数据中找出包含所述标定靶内所有所述标靶球球心坐标的第二球心点云的步骤包括:
根据各所述第二二维圆心坐标从所述第二点云数据找出相应的第二三维球面点云,并拟合出各所述第二三维球面点云的球心;
对所有所述第二三维球面点云的球心进行离群点处理,得到包含所述标定靶内所有所述标靶球球心坐标的第二球心点云。
7.根据权利要求1所述的三维点云数据配准方法,其特征在于,根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到同一坐标系下的步骤包括:
利用配准算法计算将所述第二球心点云配准到所述第一球心点云所在坐标系下的第一变换矩阵;
根据所述第一变换矩阵对所述第二点云数据进行转换;
利用配准算法计算转换后的所述第二点云数据配准到所述第一点云数据所在坐标下的第二变换矩阵;
根据所述第一变换矩阵及所述第二变换矩阵将所述第二点云数据配准到所述第一点云数据所在坐标系下。
8.根据权利要求7所述的三维点云数据配准方法,其特征在于,根据所述第一变换矩阵对所述第二点云数据进行转换的步骤之前还包括:
分别对所述第一点云数据、所述第二点云数据进行抽稀;
同时,根据所述第一变换矩阵对所述第二点云数据进行转换的步骤为:
根据所述第一变换矩阵对抽稀后的所述第二点云数据进行转换,得到转换后的第二抽稀点云数据;
同时,利用配准算法计算转换后的所述第二点云数据配准到所述第一点云数据所在坐标下的第二变换矩阵的步骤为:
利用配准算法计算所述转换后的第二抽稀点云数据配准到抽稀后的所述第一点云数据所在坐标系下的第二变换矩阵。
9.根据权利要求1所述的三维点云数据配准方法,其特征在于,所述标定靶还包括底座及多根连接杆;各所述连接杆的一端安装于所述底座上,且各所述连接杆的另一端分别与相应的所述标靶球连接;其中,所述连接杆的颜色为黑色,所述标靶球的颜色为白色。
10.一种三维点云数据拼接方法,包括:
获取三维扫描装置分别对依次相邻的第一目标空间、第二目标空间、第三目标空间、第四目标空间扫描得到的第一点云数据、第二点云数据、第三点云数据、第四点云数据;其中,所述第一点云数据与所述第二点云数据分别由所述三维扫描装置以不同视角扫描得出,且均包括位于同一位置的标定靶的点云数据,其中,所述标定靶包括多个标靶球;所述第二点云数据和所述第三点云数据为相同的点云;所述第三点云数据和所述第四点云数据分别由所述三维扫描装置以不同视角扫描得出,且均包括位于同一位置的所述标定靶的点云数据;
从所述第一点云数据、所述第二点云数据中分别找出第一球心点云、第二球心点云;所述第一球心点云、所述第二球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第一球心点云与所述第二球心点云之间的转换关系将所述第一点云数据与所述第二点云数据配准到第一坐标系下;
从所述第三点云数据、所述第四点云数据中分别找出第三球心点云、第四球心点云;所述第三球心点云、所述第四球心点云均包含所述标定靶内所有所述标靶球的球心坐标;
根据所述第三球心点云与第四球心点云之间的转换关系,将所述第四点云数据配准到所述第三点云数据所在的坐标系下;
利用所述第三点云数据与所述第二点云数据之间的关系,将配准后的所述第四点云数据配准到所述第一坐标系下;
将配准至所述第一坐标下的第一点云数据、配准至所述第一坐标下的第二点云数据及配准至所述第一坐标下的第四点云数据拼接在一起。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610855932.3A CN106651752B (zh) | 2016-09-27 | 2016-09-27 | 三维点云数据配准方法及拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610855932.3A CN106651752B (zh) | 2016-09-27 | 2016-09-27 | 三维点云数据配准方法及拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106651752A true CN106651752A (zh) | 2017-05-10 |
CN106651752B CN106651752B (zh) | 2020-01-21 |
Family
ID=58854979
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610855932.3A Active CN106651752B (zh) | 2016-09-27 | 2016-09-27 | 三维点云数据配准方法及拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106651752B (zh) |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107452062A (zh) * | 2017-07-25 | 2017-12-08 | 深圳市魔眼科技有限公司 | 三维模型构建方法、装置、移动终端、存储介质及设备 |
CN107861920A (zh) * | 2017-11-27 | 2018-03-30 | 西安电子科技大学 | 点云数据配准方法 |
CN107918753A (zh) * | 2016-10-10 | 2018-04-17 | 腾讯科技(深圳)有限公司 | 点云数据处理方法及装置 |
CN108470323A (zh) * | 2018-03-13 | 2018-08-31 | 京东方科技集团股份有限公司 | 一种图像拼接方法、计算机设备及显示装置 |
CN108876935A (zh) * | 2018-06-06 | 2018-11-23 | 链家网(北京)科技有限公司 | 一种在移动端快速进行房屋三维模型拼接的方法及装置 |
CN109522808A (zh) * | 2018-10-22 | 2019-03-26 | 杭州视在科技有限公司 | 一种停机坪上反光锥筒摆放的自动识别方法 |
CN109636917A (zh) * | 2018-11-02 | 2019-04-16 | 北京微播视界科技有限公司 | 三维模型的生成方法、装置、硬件装置 |
CN109781029A (zh) * | 2019-03-08 | 2019-05-21 | 苏州玻色智能科技有限公司 | 一种产品测量三维数据拼接方法 |
CN109870118A (zh) * | 2018-11-07 | 2019-06-11 | 南京林业大学 | 一种面向绿色植物时序模型的点云采集方法 |
CN110010249A (zh) * | 2019-03-29 | 2019-07-12 | 北京航空航天大学 | 基于视频叠加的增强现实手术导航方法、系统及电子设备 |
CN110110678A (zh) * | 2019-05-13 | 2019-08-09 | 腾讯科技(深圳)有限公司 | 道路边界的确定方法和装置、存储介质及电子装置 |
CN110163900A (zh) * | 2019-05-31 | 2019-08-23 | 北京百度网讯科技有限公司 | 用于调整点云数据方法和装置 |
CN110246167A (zh) * | 2019-06-14 | 2019-09-17 | 北京百度网讯科技有限公司 | 用于处理点云数据的方法和装置 |
CN110619611A (zh) * | 2019-09-09 | 2019-12-27 | 北京迈格威科技有限公司 | 图像校正标定方法、装置、计算机设备和存储介质 |
CN111256592A (zh) * | 2020-03-13 | 2020-06-09 | 易思维(杭州)科技有限公司 | 结构光传感器的外参标定装置及方法 |
CN111325779A (zh) * | 2020-02-07 | 2020-06-23 | 贝壳技术有限公司 | 点云配准方法和装置、电子设备和存储介质 |
CN111539074A (zh) * | 2020-03-31 | 2020-08-14 | 杭州数孪科技有限公司 | 一种工件加工数据处理方法与装置 |
CN111783798A (zh) * | 2020-07-01 | 2020-10-16 | 北京大学深圳研究生院 | 基于显著性特征的模拟残缺点云的遮罩生成方法 |
CN111982018A (zh) * | 2020-09-16 | 2020-11-24 | 广州天域科技有限公司 | 一种跟踪式激光三维扫描仪 |
CN112258494A (zh) * | 2020-10-30 | 2021-01-22 | 北京柏惠维康科技有限公司 | 一种病灶位置确定方法、装置及电子设备 |
CN112292611A (zh) * | 2019-05-24 | 2021-01-29 | 深圳市速腾聚创科技有限公司 | 坐标校正方法、装置、计算设备及计算机存储介质 |
CN112750165A (zh) * | 2019-10-29 | 2021-05-04 | 商汤集团有限公司 | 参数标定方法、智能驾驶方法及其装置、设备和存储介质 |
CN112825190A (zh) * | 2019-11-19 | 2021-05-21 | 广东博智林机器人有限公司 | 一种精度评估方法、系统、电子设备和存储介质 |
CN112884902A (zh) * | 2021-03-17 | 2021-06-01 | 中山大学 | 一种面向点云配准的标靶球位置优化方法 |
CN112880562A (zh) * | 2021-01-19 | 2021-06-01 | 佛山职业技术学院 | 一种机械臂末端位姿误差测量方法及系统 |
CN112967189A (zh) * | 2021-03-16 | 2021-06-15 | 广东电网有限责任公司广州供电局 | 一种点云拼接方法、装置、设备和存储介质 |
CN113532311A (zh) * | 2020-04-21 | 2021-10-22 | 广东博智林机器人有限公司 | 点云拼接方法、装置、设备和存储设备 |
CN113805157A (zh) * | 2021-09-22 | 2021-12-17 | 航天新气象科技有限公司 | 基于标靶的高度测量方法、装置及设备 |
CN113808093A (zh) * | 2021-09-10 | 2021-12-17 | 中铁一局集团第五工程有限公司 | 基于3d激光扫描仪的隧道初支喷射混凝土厚度检测方法 |
CN113927178A (zh) * | 2021-10-21 | 2022-01-14 | 泉州冰点科技有限公司 | 一种去除鞋垫凸点的方法及装置 |
CN115100258A (zh) * | 2022-08-29 | 2022-09-23 | 杭州三坛医疗科技有限公司 | 一种髋关节图像配准方法、装置、设备以及存储介质 |
CN117110319A (zh) * | 2023-10-23 | 2023-11-24 | 汇鼎智联装备科技(江苏)有限公司 | 基于3d成像的球体表面缺陷检测方法和检测系统 |
CN117501322A (zh) * | 2021-04-14 | 2024-02-02 | 利尼芝物流有限责任公司 | 从点云数据生成物理空间的映射 |
CN117579753A (zh) * | 2024-01-16 | 2024-02-20 | 思看科技(杭州)股份有限公司 | 三维扫描方法、装置、计算机设备以及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008573A (zh) * | 2014-06-11 | 2014-08-27 | 清华大学 | 基于三维形态的颗粒分析方法及颗粒放置支架 |
CN105651160A (zh) * | 2014-11-04 | 2016-06-08 | 北京浩宇天地测绘科技发展有限公司 | 一种安装定位装置的三维激光扫描仪球形标靶及其使用方法 |
-
2016
- 2016-09-27 CN CN201610855932.3A patent/CN106651752B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008573A (zh) * | 2014-06-11 | 2014-08-27 | 清华大学 | 基于三维形态的颗粒分析方法及颗粒放置支架 |
CN105651160A (zh) * | 2014-11-04 | 2016-06-08 | 北京浩宇天地测绘科技发展有限公司 | 一种安装定位装置的三维激光扫描仪球形标靶及其使用方法 |
Non-Patent Citations (3)
Title |
---|
DONGHO YUN ET AL.: "Automated registration of multi-view point clouds using sphere targets", 《ADVANCED ENGINEERING INFORMATICS》 * |
MENGMI ZHANG: "Accurate Sphere Marker-Based Registration System of 3D Point Cloud Data in Applications of Shipbuilding Blocks", 《JOURNAL OF INDUSTRIAL AND INTELLIGENT INFORMATION》 * |
曹奇 等: "基于标靶配准的廊道点云数据", 《现代测绘》 * |
Cited By (50)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107918753A (zh) * | 2016-10-10 | 2018-04-17 | 腾讯科技(深圳)有限公司 | 点云数据处理方法及装置 |
US10846874B2 (en) | 2016-10-10 | 2020-11-24 | Tencent Technology (Shenzhen) Company Limited | Method and apparatus for processing point cloud data and storage medium |
CN107452062A (zh) * | 2017-07-25 | 2017-12-08 | 深圳市魔眼科技有限公司 | 三维模型构建方法、装置、移动终端、存储介质及设备 |
CN107452062B (zh) * | 2017-07-25 | 2020-03-06 | 深圳市魔眼科技有限公司 | 三维模型构建方法、装置、移动终端、存储介质及设备 |
CN107861920A (zh) * | 2017-11-27 | 2018-03-30 | 西安电子科技大学 | 点云数据配准方法 |
CN107861920B (zh) * | 2017-11-27 | 2021-11-30 | 西安电子科技大学 | 点云数据配准方法 |
CN108470323A (zh) * | 2018-03-13 | 2018-08-31 | 京东方科技集团股份有限公司 | 一种图像拼接方法、计算机设备及显示装置 |
US11308579B2 (en) | 2018-03-13 | 2022-04-19 | Boe Technology Group Co., Ltd. | Image stitching method, image stitching apparatus, display apparatus, and computer product |
CN108876935A (zh) * | 2018-06-06 | 2018-11-23 | 链家网(北京)科技有限公司 | 一种在移动端快速进行房屋三维模型拼接的方法及装置 |
CN109522808A (zh) * | 2018-10-22 | 2019-03-26 | 杭州视在科技有限公司 | 一种停机坪上反光锥筒摆放的自动识别方法 |
CN109636917A (zh) * | 2018-11-02 | 2019-04-16 | 北京微播视界科技有限公司 | 三维模型的生成方法、装置、硬件装置 |
CN109870118A (zh) * | 2018-11-07 | 2019-06-11 | 南京林业大学 | 一种面向绿色植物时序模型的点云采集方法 |
CN109870118B (zh) * | 2018-11-07 | 2020-09-11 | 南京林业大学 | 一种面向绿色植物时序模型的点云采集方法 |
CN109781029A (zh) * | 2019-03-08 | 2019-05-21 | 苏州玻色智能科技有限公司 | 一种产品测量三维数据拼接方法 |
CN110010249A (zh) * | 2019-03-29 | 2019-07-12 | 北京航空航天大学 | 基于视频叠加的增强现实手术导航方法、系统及电子设备 |
CN110010249B (zh) * | 2019-03-29 | 2021-04-27 | 北京航空航天大学 | 基于视频叠加的增强现实手术导航方法、系统及电子设备 |
CN110110678A (zh) * | 2019-05-13 | 2019-08-09 | 腾讯科技(深圳)有限公司 | 道路边界的确定方法和装置、存储介质及电子装置 |
CN112292611A (zh) * | 2019-05-24 | 2021-01-29 | 深圳市速腾聚创科技有限公司 | 坐标校正方法、装置、计算设备及计算机存储介质 |
CN110163900A (zh) * | 2019-05-31 | 2019-08-23 | 北京百度网讯科技有限公司 | 用于调整点云数据方法和装置 |
CN110163900B (zh) * | 2019-05-31 | 2021-04-20 | 北京百度网讯科技有限公司 | 用于调整点云数据方法和装置 |
CN110246167A (zh) * | 2019-06-14 | 2019-09-17 | 北京百度网讯科技有限公司 | 用于处理点云数据的方法和装置 |
CN110619611A (zh) * | 2019-09-09 | 2019-12-27 | 北京迈格威科技有限公司 | 图像校正标定方法、装置、计算机设备和存储介质 |
CN112750165B (zh) * | 2019-10-29 | 2023-12-15 | 商汤集团有限公司 | 参数标定方法、智能驾驶方法及其装置、设备和存储介质 |
CN112750165A (zh) * | 2019-10-29 | 2021-05-04 | 商汤集团有限公司 | 参数标定方法、智能驾驶方法及其装置、设备和存储介质 |
CN112825190A (zh) * | 2019-11-19 | 2021-05-21 | 广东博智林机器人有限公司 | 一种精度评估方法、系统、电子设备和存储介质 |
CN112825190B (zh) * | 2019-11-19 | 2024-02-09 | 广东博智林机器人有限公司 | 一种精度评估方法、系统、电子设备和存储介质 |
CN111325779B (zh) * | 2020-02-07 | 2020-12-11 | 贝壳找房(北京)科技有限公司 | 点云配准方法和装置、电子设备和存储介质 |
CN111325779A (zh) * | 2020-02-07 | 2020-06-23 | 贝壳技术有限公司 | 点云配准方法和装置、电子设备和存储介质 |
CN111256592A (zh) * | 2020-03-13 | 2020-06-09 | 易思维(杭州)科技有限公司 | 结构光传感器的外参标定装置及方法 |
CN111539074A (zh) * | 2020-03-31 | 2020-08-14 | 杭州数孪科技有限公司 | 一种工件加工数据处理方法与装置 |
CN111539074B (zh) * | 2020-03-31 | 2023-11-14 | 苏州数设科技有限公司 | 一种工件加工数据处理方法与装置 |
CN113532311A (zh) * | 2020-04-21 | 2021-10-22 | 广东博智林机器人有限公司 | 点云拼接方法、装置、设备和存储设备 |
CN111783798A (zh) * | 2020-07-01 | 2020-10-16 | 北京大学深圳研究生院 | 基于显著性特征的模拟残缺点云的遮罩生成方法 |
CN111783798B (zh) * | 2020-07-01 | 2023-09-15 | 北京大学深圳研究生院 | 基于显著性特征的模拟残缺点云的遮罩生成方法 |
CN111982018B (zh) * | 2020-09-16 | 2021-11-30 | 广州天域科技有限公司 | 一种跟踪式激光三维扫描仪 |
CN111982018A (zh) * | 2020-09-16 | 2020-11-24 | 广州天域科技有限公司 | 一种跟踪式激光三维扫描仪 |
CN112258494A (zh) * | 2020-10-30 | 2021-01-22 | 北京柏惠维康科技有限公司 | 一种病灶位置确定方法、装置及电子设备 |
CN112258494B (zh) * | 2020-10-30 | 2021-10-22 | 北京柏惠维康科技有限公司 | 一种病灶位置确定方法、装置及电子设备 |
CN112880562A (zh) * | 2021-01-19 | 2021-06-01 | 佛山职业技术学院 | 一种机械臂末端位姿误差测量方法及系统 |
CN112967189A (zh) * | 2021-03-16 | 2021-06-15 | 广东电网有限责任公司广州供电局 | 一种点云拼接方法、装置、设备和存储介质 |
CN112884902A (zh) * | 2021-03-17 | 2021-06-01 | 中山大学 | 一种面向点云配准的标靶球位置优化方法 |
CN117501322A (zh) * | 2021-04-14 | 2024-02-02 | 利尼芝物流有限责任公司 | 从点云数据生成物理空间的映射 |
CN113808093A (zh) * | 2021-09-10 | 2021-12-17 | 中铁一局集团第五工程有限公司 | 基于3d激光扫描仪的隧道初支喷射混凝土厚度检测方法 |
CN113805157A (zh) * | 2021-09-22 | 2021-12-17 | 航天新气象科技有限公司 | 基于标靶的高度测量方法、装置及设备 |
CN113927178A (zh) * | 2021-10-21 | 2022-01-14 | 泉州冰点科技有限公司 | 一种去除鞋垫凸点的方法及装置 |
CN115100258B (zh) * | 2022-08-29 | 2023-02-07 | 杭州三坛医疗科技有限公司 | 一种髋关节图像配准方法、装置、设备以及存储介质 |
CN115100258A (zh) * | 2022-08-29 | 2022-09-23 | 杭州三坛医疗科技有限公司 | 一种髋关节图像配准方法、装置、设备以及存储介质 |
CN117110319B (zh) * | 2023-10-23 | 2024-01-26 | 汇鼎智联装备科技(江苏)有限公司 | 基于3d成像的球体表面缺陷检测方法和检测系统 |
CN117110319A (zh) * | 2023-10-23 | 2023-11-24 | 汇鼎智联装备科技(江苏)有限公司 | 基于3d成像的球体表面缺陷检测方法和检测系统 |
CN117579753A (zh) * | 2024-01-16 | 2024-02-20 | 思看科技(杭州)股份有限公司 | 三维扫描方法、装置、计算机设备以及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106651752B (zh) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106651752A (zh) | 三维点云数据配准方法及拼接方法 | |
CN108555908B (zh) | 一种基于rgbd相机的堆叠工件姿态识别及拾取方法 | |
CN111486855B (zh) | 一种具有物体导航点的室内二维语义栅格地图构建方法 | |
CN101299270B (zh) | 三维扫描系统中的多个摄像机同步快速标定方法 | |
CN107292927B (zh) | 一种基于双目视觉的对称运动平台位姿测量方法 | |
CN102773862B (zh) | 用于室内移动机器人的快速精确定位系统及其工作方法 | |
CN101839692B (zh) | 单相机测量物体三维位置与姿态的方法 | |
CN106340044B (zh) | 摄像机外参自动标定方法及标定装置 | |
CN104574406B (zh) | 一种360度全景激光与多个视觉系统间的联合标定方法 | |
CN107392929B (zh) | 一种基于人眼视觉模型的智能化目标检测及尺寸测量方法 | |
CN202702247U (zh) | 用于室内移动机器人的快速精确定位系统 | |
CN109544612A (zh) | 基于特征点几何表面描述的点云配准方法 | |
CN109443209A (zh) | 一种基于单应性矩阵的线结构光系统标定方法 | |
CN111028271B (zh) | 基于人体骨架检测的多摄像机人员三维定位跟踪系统 | |
CN103295239A (zh) | 一种基于平面基准影像的激光点云数据的自动配准方法 | |
CN105574812B (zh) | 多角度三维数据配准方法及装置 | |
CN110823252B (zh) | 一种多线激光雷达和单目视觉的自动标定方法 | |
CN113096183B (zh) | 一种基于激光雷达与单目相机的障碍物检测与测量方法 | |
Liang et al. | Automatic registration of terrestrial laser scanning data using precisely located artificial planar targets | |
CN101213440A (zh) | 创建用于检查凹凸图形的主数据的方法 | |
CN112614188A (zh) | 一种基于交比不变性的点阵式标定板及其识别方法 | |
Abdul-Rahman et al. | Innovations in 3D geo information systems | |
CN114140539A (zh) | 一种室内物体的位置获取方法和装置 | |
CN108388854A (zh) | 一种基于改进fast-surf算法的定位方法 | |
CN114137564A (zh) | 一种室内物体自动标识定位方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |