CN105976312A - 基于点特征直方图的点云自动配准方法 - Google Patents
基于点特征直方图的点云自动配准方法 Download PDFInfo
- Publication number
- CN105976312A CN105976312A CN201610371302.9A CN201610371302A CN105976312A CN 105976312 A CN105976312 A CN 105976312A CN 201610371302 A CN201610371302 A CN 201610371302A CN 105976312 A CN105976312 A CN 105976312A
- Authority
- CN
- China
- Prior art keywords
- point
- cloud data
- station
- point cloud
- registered
- 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 68
- 238000005070 sampling Methods 0.000 claims abstract description 116
- 239000011159 matrix material Substances 0.000 claims abstract description 67
- 230000009466 transformation Effects 0.000 claims abstract description 42
- 230000008569 process Effects 0.000 claims description 24
- 230000001133 acceleration Effects 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 13
- 238000002474 experimental method Methods 0.000 description 8
- 239000004575 stone Substances 0.000 description 6
- 238000012795 verification Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 241001577299 Vindex Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 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
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了基于点特征直方图的点云自动配准方法,包括:对被测物体多站采集获取点云数据,确定基准站和待配准站;从两站点云数据中均随机选取m个采样点,使任意三个点都不共线,每两个采样点均满足距离约束原则;对两站进行同名点对匹配,使基准站采样点集和待配准站采样点集形成的同名点对三角网近似全等;获取待配准站点云数据的变换矩阵,计算点特征直方图差异值获取配准误差;筛选出配准误差值最小时对应的变换矩阵;计算两站点云数据的重叠区域点云数据,通过迭代最近点算法获得精细配准矩阵。本发明能够解决对相邻两站点云重叠度较小或重叠区域特征不明显情况下的点云数据配准时容易陷入局部最小解问题,该方法准确性高,整体效率高。
Description
技术领域
本发明属于点云数据配准领域,特别涉及一种基于点特征直方图的点云自动配准方法。
背景技术
三维激光扫描技术早期兴起于工业领域的逆向工程中,主要用于工艺品设计、汽车、飞机以及船舶等工业制造业领域,后来逐渐被应用到测绘行业领域。通常情况下,三维激光扫描工作由于受到被测物体范围较广、扫描设备的测量范围有限、物体遮挡等等因素的影响,需架设多站进行不同角度的扫描,从而获取到被测对象表面的整体点云数据。而根据扫描仪的扫描原理,采集到的三维点云坐标数据并不是出于统一坐标系下,因此需要将不同站、不同时间段或者不同传感器设备采集到的点云数据匹配到统一三维空间直角坐标系底下,这个过程即称之为点云数据配准。点云数据配准主要分为初始配准和精细配准两大模块。
点云初始配准算法主要有以下几大类:基于测量设备的初始配准、基于辅助标识的初始配准、基于纹理特征的初始配准、基于表面几何特征的初始配准、基于RANSAC的初始配准。其中基于测量设备的初始配准指的是依靠测量设备的一些定位信息来获取不同站点云之间的相对移动位置关系,进而利用这些相对关系进行初始配准。这类初始配准对测量设备过于依赖,无法满足大部分实际工程需求。在很多工程应用中,扫描现场摆设适当数量的控制标靶,进而在内业中进行标靶点云的识别。研究者们据此提出了一些基于辅助标识的自动配准方法,但这类方法增加了外业的工作量,并且有时会出现标靶遮挡物体表面的情况出现。此外随着扫描仪的发展,扫描仪逐渐开始搭载相机,因而研究者们提出了一些基于纹理特征的初始配准算法,但由于该类算法需要影像信息,如果影像有变形或者双张影像重叠度不够则无法进行匹配,不具备广泛应用性。除了依靠点云的属性信息进行点云配准工作外,还可以依据点云中潜在的一些几何特征进行点云配准工作,研究人员已经在此方面做出了大量的研究。但这类算法多数只适用于噪声点较少的情况下,很易因噪声点的参与而干扰特征描述子的计算准确度,从而影响配准效果。且点云对象特征不明显的情况下,上述的几种算法很难获取到正确结果。通常情况下,为避免一些错误点或者局外点对某个算法造成影响,应当首先借助某种方法将这些局外点剔除,因而很多研究者提出了一些基于RANSAC的初始配准方法,其中R.B.Rusu等提出了一种基于点特征快速直方图的SAC-IA算法,但该算法针对重叠度较小情况下的点云配准容易造成局部最小解问题,且针对大数据量点云配准该算法无法达到高效实时配准,因此需要进一步改进。
点云精细配准算法大多都是基于ICP(Iterative closest point,迭代最近点)算法或者其变种改进算法进行的。ICP算法由Besl和Mckay以及Chen和Medioni早期提出,该算法为后续基于迭代的配准算法发展奠定了理论和框架基础。基于前人研究这些不同种类的变种ICP算法基础,由于不同种变种算法针对不同传感器和不同环境因素有着各自的优缺点,并没有合适的框架通过比较这些算法针对不同情况提出合适的算法,因此Pomerleau等提出了一种变种ICP算法评价协议标准,此外Pomerleau也提出了一种ICP算法开源库,该库涵盖了一种ICP可重构算法框架,该框架可以融入不同种改进ICP算法,框架内部可以快速比较不同变种算法的性能与准确性。但即便是上述各种ICP算法改进及优化,但针对重叠度较小的点云而言,初始配准后直接进行ICP算法匹配也容易陷入局部最小解问题。
发明内容
本发明的一个目的是提供了一种基于点特征直方图的点云自动配准方法,能够解决现有技术中对相邻两站点云重叠度较小或重叠区域特征不明显情况下的点云数据配准时容易陷入局部最小解的问题,且该方法准确性高,整体效率高。
本发明提供的技术方案为:
一种基于点特征直方图的点云自动配准方法,包括:
对被测物体进行多站采集获取该被测物体表面的整体点云数据,确定基准站和待配准站,其中,基准站和待配准站为相邻的两站;
从所述基准站点云数据和所述待配准站点云数据中均随机选取m个采样点,均使所选取的m个采样点中的任意三个点都不共线,且每两个采样点之间的距离均满足距离约束原则,其中,m>=3;
对所述待配准站采样点集和所述基准站采样点集进行同名点对的对应匹配,使所述基准站采样点集和所述待配准站采样点集形成的同名点对三角网近似全等;
获取所述待配准站点云数据的变换矩阵T,通过计算点特征直方图差异值获取所述待配准站点云数据经过变换矩阵T变换后与所述基准站点云数据的配准误差;
对所述基准站点云数据和所述待配准站点云数据均进行多次采样,从多次采样获取的配准误差中筛选出配准误差值最小时对应的变换矩阵Tmin,进行点云数据的初始配准;
计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据,并通过迭代最近点算法获得精细配准矩阵,进行点云数据的精细配准。
优选的是,所述的基于点特征直方图的点云自动配准方法,从所述基准站点云数据中和所述待配准站点云数据中均随机选取m个采样点之前还包括:
对所述基准站点云数据和所述待配准站点云数据进行预处理;
计算所述基准站点云表面和所述待配准站点云表面的法线;
计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图。
优选的是,所述的基于点特征直方图的点云自动配准方法,从所述基准站点云数据中和所述待配准站点云数据中均随机选取m个采样点具体包括:
对所述基准站点云数据和所述待配准站点云数据进行采样时,当采集一个新样点时,判断所选取的新样点与已有的采样点是否满足距离约束原则;
若不满足,则新样点无效,继续进行采集,若满足,则继续判断所选取的新样点与已有采样点构建的直线的垂直距离是否小于预设阈值∈dist;
若小于,则所采集的新样点与已有的采样点共线,该新样点无效,重新进行采样;
若不小于,则所采集的新样点有效,继续采集下一个样点。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述通过计算点特征直方图差异值获取所述待配准站点云数据经过变换矩阵T变换后与所述基准站点云数据的配准误差包括:
所述待配准站点云数据P经过变换矩阵T变换后为P',在所述基准站点云数据Q中查找所述待配准站点云数据P'中每个点对应的最近邻域点;
根据获取的点Pi对应的最近邻域点,提取出邻域索引和邻域间距离;
判断邻域间距离是否小于预定的邻域点对距离阈值∈;
若不小于,则继续遍历点云数据P'中下一个点进行邻域间距离判断;
若小于,则计算点Pi与其对应的最近邻域点形成的点对之间的快速直方图差异平方值总和error,并统计至总误差值SumError上,同时累计满足预定的邻域点对距离阈值∈的点对个数nr;
判断nr是否大于零,若大于零,返回配准的点特征中误差若不大于零,则返回配准的分值为无穷大。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据包括:
分别计算所述待配准站点云数据P和所述基准站点云数据Q的点云包围盒;
对所述基准站点云数据Q进行空间网格线性划分;
遍历所述待配准站点云数据P中的每个点Pi,并计算每个点Pi对应的网格线性编码值;
在所述基准站点云数据Q的空间网格中查找是否存在与该待配准站点云数据中点Pi的网格线性编码值对应的空间网格;
若不存在,则继续遍历待配准站点云数据P中的点;
若存在,则对所述待配准站点云数据中点Pi的索引进行存储,且标记所述待配准站点云数据中点Pi的网格线性编码值对应所述基准站点云数据Q中的空间网格为重叠网格;
根据存储的所述待配准站点云数据中点Pi的索引,以及所述基准站点云数据Q中重叠的空间网格,得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述通过迭代最近点算法获得精细配准矩阵包括下列步骤:
步骤一、得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据后,根据获取的重叠区域点云数据执行迭代最近点算法,得到配准矩阵Ti;
步骤二、将配准矩阵Ti应用至经过初始配准后的所述待配准站点云数据P'上,即P′i=P′i-1*Ti;
步骤三、循环执行步骤一至步骤二,达到迭代循环次数N,其中,N>=2;
步骤四、计算所述待配准站点云数据P'经过配准矩阵Ti配准后与所述基准站点云数据的配准误差ε;
步骤五、判断配准误差ε是否小于等于点云采样分辨率的0.5倍范围内,
若不满足,则使迭代循环次数N加1;
若满足,则输出精细配准矩阵:T=T1*T2*...Tn。
优选的是,所述的基于点特征直方图的点云自动配准方法,在所述待配准站点云数据的初始配准阶段和精细配准阶段,均采用开放运算语言并行加速模块进行加速。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述待配准站点云数据初始配准阶段的开放运算语言并行加速模块包括:点云邻域搜索模块、点特征直方图计算模块以及初始配准迭代过程模块。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述待配准站点云数据初始配准阶段的开放运算语言并行加速过程包括:
中央处理器将去噪后的所述基准站点云数据和所述待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,点云邻域并行搜索和计算所述基准站点云表面和所述待配准站点云表面的法线同时执行;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征简易直方图;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图;
并行执行初始配准迭代过程;
图形处理器将初始配准迭代过程中获取的变换矩阵及配准误差发送至中央处理器中,筛选出配准误差最小时对应的变换矩阵Tmin。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述待配准站点云数据精细配准阶段的开放运算语言并行加速模块包括:每次迭代的邻域搜索模块、最小二乘估计变换矩阵模块以及点云矩阵变换模块。
优选的是,所述的基于点特征直方图的点云自动配准方法,所述待配准站点云数据精细配准阶段的开放运算语言并行加速过程包括:
中央处理器将初始配准后的基准站点云数据和待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,进行邻域并行搜索;
对变换矩阵变量G、h并行填充,并发送矩阵G、h至中央处理器中,进行楚列斯基分解法求解变换矩阵T_i;
中央处理器发送变换矩阵T_i至图形处理器中进行点云矩阵变换并行计算,同时累乘精细配准矩阵
其中,
本发明至少包括以下有益效果:首先在点云数据初始配准阶段采用了基于点特征直方图的采样一致性初始配准算法(SAC-IA),并对该算法的不足之处,从算法准确度和效率角度提出了改进,算法在点云数据初始配准阶段提出了两项改进,一是在算法中加入了采样点几何形状约束以保证采样点不共线,二是通过计算点特征直方图差异代替了欧式距离差异改进了配准误差评定方法,从而保证了算法的准确性,避免了局部最小解问题;附加了同名点对三角网是否近似全等判断条件,避免了大部分不符合条件点对的迭代计算,从而提高了算法的整体效率。在点云数据精细配准阶段通过计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据,并通过迭代最近点算法获得精细配准矩阵,从而进行所述待配准站点云数据的精细配准,避免了现有算法不能应对点云重叠度较小的点云数据配准情况。在初始配准阶段和精细配准阶段,均采用开放运算语言(Open ComputingLanguage,简称OpenCL)并行加速算法进行加速,从而在整体上很大程度地提高了算法效率,使得算法足以完成海量点云实时性配准工作。
附图说明
图1为本发明所述的基于点特征直方图的点云自动配准方法的流程示意图;
图2为本发明所述的基于点特征直方图的点云自动配准方法的初始配准阶段的流程示意图;
图3为本发明所述的基于点特征直方图的点云自动配准方法的精细配准阶段的流程示意图;
图4为点对三角网近似全等示意图;
图5为点云数据空间网格划分效果图;
图6为基于重叠区域的精细配准算法迭代效果图;
图7为点云数据初始配准并行加速流程示意图;
图8为点云数据精细配准并行加速流程示意图;
图9(a)为对故宫云龙大石雕扫描的两站去噪后原始数据视图;
图9(b)为对故宫云龙大石雕扫描的两站原始点云数据物理叠合效果图;
图9(c)为对故宫云龙大石雕扫描的两站点云数据自动配准后的点云叠加效果图;
图10(a)为对一高校扫描的两站去燥后原始数据视图;
图10(b)为对一高校扫描的两站原始点云数据物理叠合效果图;
图10(c)为对一高校扫描的两站点云数据自动配准后的点云叠加效果图;
图11(a)为对一块地形扫描的两站去燥后原始数据视图;
图11(b)为对一块地形扫描的两站原始点云数据物理叠合效果图;
图11(c)为对一块地形扫描的两站点云数据自动配准后的点云叠加效果图;
图12(a)为对一棵树扫描的两站去燥后原始数据视图;
图12(b)为对一棵树扫描的两站原始点云数据物理叠合效果图;
图12(c)为对一棵树扫描的两站点云数据自动配准后的点云叠加效果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
为使本发明技术方案的优点更加清楚,下面结合附图和实施例对本发明作详细说明。
如图1,图2所示,本发明实施例提供的基于点特征直方图的点云自动配准方法,包括下列步骤:
S01、对被测物体进行多站采集获取该被测物体表面的整体点云数据,确定基准站和待配准站,其中,基准站和待配准站为相邻的两站;
需要说明的是:对被测物体进行多站采集,一般将第一站作为基准站,但是具体将哪一站作为基准站,本发明并不做具体限定。待配准站和基准站为相邻的两站,如果将第一站作为基准站,一般将相邻的第二站作为待配准站,当将第二站和第一站进行点云数据配准后,第二站的点云数据和第一站的点云数据位于统一的三维空间直角坐标系中,这时将统一后的第一站和第二站作为基准站,将和第二站相邻的第三站作为待配准站。最后将采集的所有站点云数据都配准至统一的三维空间直角坐标系中。
S02、从所述基准站点云数据和所述待配准站点云数据中均随机选取m个采样点,均使所选取的m个采样点中的任意三个点都不共线,且每两个采样点之间的距离均满足距离约束原则,其中,m>=3;
其中,如图2所示,在所述步骤S02之前,还包括步骤:
A、对所述基准站点云数据和所述待配准站点云数据进行预处理;其中,对所述基准站点云数据和所述待配准站点云数据进行预处理,主要是对点云数据进行去噪处理。
B、计算所述基准站点云表面和所述待配准站点云表面的法线。
C、计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图。
其中,步骤S02的过程具体为:
S021、对所述基准站点云数据和所述待配准站点云数据进行采样时,当采集一个新样点时,判断所选取的新样点与已有的采样点是否满足距离约束原则;
S022、若不满足,则新样点无效,继续进行采集,
S023、若满足,则继续判断所选取的新样点与已有采样点构建的直线的垂直距离是否小于预设阈值∈dist;
S024、若小于,则所采集的新样点与已有的采样点共线,该新样点无效,重新进行采样;
S025、若不小于,则所采集的新样点有效,继续采集下一个样点。
例如:在进行SAC-IA算法的采样点时,依次在待配准站点云中依据距离约束原则采样点,当采样到第二个点时,由已有的两个采样点P1、P2构建成一条直线L;继续在待配准站点云中采集点,每采样一个新的样点Pnew,就开始首先判断与已有样点的距离是否满足距离约束原则,进而再判断该点到构建的直线L的垂直距离是否小于预设阈值∈dist,如果小于则代表与已有样点共线,该样点无效,需重新筛选,否则该样点有效,继续采样下一个点,直到采集样点的个数达到要求个数。
其中,预设阈值∈dist,一般设置为点云采样分辨率值的若干倍,例如3到5倍之间。
S03、对所述待配准站采样点集和所述基准站采样点集进行同名点对的对应匹配,使所述基准站采样点集和所述待配准站采样点集形成的同名点对三角网近似全等;
需要说明的是:对所述待配准站采样点集和所述基准站采样点集进行同名点对的对应匹配,是利用SAC-IA算法寻找同名点对,同名点对搜索完成后,分别依照顺序将基准站点云数据Q,待配准站点云数据P中各自的m个采样点每三个点构成一个三角形,然后各自组成一个三角网形式。依次判断所有组两站点对的对应三角形是否近似全等,如果都近似全等则继续SAC-IA算法的矩阵求解步骤,否则SAC-IA算法执行返回到采样点步骤。
例如:当采样点m=3时,判断单个三角形是否近似全等,m>3时则判断三角网是否近似全等。两站点云对应的采样点组成的三角网如图4所示。
S04、估计所述待配准站点云数据的变换矩阵T,通过计算点特征直方图差异值获取所述待配准站点云数据经过变换矩阵T变换后与所述基准站点云数据的配准误差;
其中,获取配准误差的具体过程为:
S041、所述待配准站点云数据P经过变换矩阵T变换后为P',在所述基准站点云数据Q中查找所述待配准站点云数据P'中每个点对应的最近邻域点;
S042、根据获取的点Pi对应的最近邻域点,提取出邻域索引index和邻域间距离dist;
S043、判断邻域间距离dist是否小于预定的邻域点对距离阈值∈;
S044、若不小于,则继续遍历点云数据P'中下一个点进行邻域间距离判断;
S045、若小于,则计算点Pi与其对应的最近邻域点形成的点对之间的快速直方图差异平方值总和error,并统计至总误差值SumError上,即:SumError=SumError+error,同时累计满足预定的邻域点对距离阈值∈的点对个数nr,即:nr=nr+1;
S046、判断nr是否大于零,若大于零,返回配准的点特征中误差若不大于零,则返回配准的分值为无穷大。
S05、循环执行步骤S02至S04,对所述基准站点云数据和所述待配准站点云数据进行多次采样,从多次采样获取的配准误差中筛选出配准误差值最小时对应的变换矩阵Tmin,进行所述待配准站点云数据的初始配准;
其中,循环执行步骤S02至S04,也就是循环执行一个迭代过程,用来获取配准误差,提高了配准的准确性。
S06、计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据,并通过迭代最近点算法获得精细配准矩阵,进行所述待配准站点云数据的精细配准。
其中,如图3所示,所述计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据包括下列步骤:
步骤一、分别计算所述待配准站点云数据P和所述基准站点云数据Q的点云包围盒;
步骤二、对所述基准站点云数据Q进行空间网格线性划分;
步骤三、遍历所述待配准站点云数据P中的每个点Pi,并计算每个点Pi对应的网格线性编码值;
步骤四、在所述基准站点云数据Q的空间网格中查找是否存在与该待配准站点云数据中点Pi的网格线性编码值对应的空间网格;
步骤五、若不存在,则执行步骤三继续遍历待配准站点云数据P中的点;
步骤六、若存在,则对所述待配准站点云数据中点Pi的索引进行存储,存储至vIndex_P中,且标记所述待配准站点云数据中点Pi的网格线性编码值对应的空间网格为重叠网格;
步骤七、根据存储的所述待配准站点云数据中点Pi的索引,以及所述基准站点云数据Q中重叠的空间网格,得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据,将重叠区域点云数据存储于P_overlap、Q_overlap;
具体的点云空间网格划分效果如图5(a)(b)所示,图5(b)为最终的重叠部分网格筛选结果。
其中,所述通过迭代最近点算法获得精细配准矩阵包括下列步骤:
步骤一、得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据后,根据获取的重叠区域点云数据执行迭代最近点算法,得到配准矩阵Ti;
步骤二、将配准矩阵Ti应用至经过初始配准后的所述待配准站点云数据P'上,即P′i=P′i-1*Ti;
步骤三、循环执行步骤一至步骤二,达到迭代循环次数N,其中,N>=2;
步骤四、计算所述待配准站点云数据P'经过配准矩阵Ti配准后与所述基准站点云数据的配准误差ε;
步骤五、判断配准误差ε是否小于等于点云采样分辨率的0.5倍范围内,
若不满足,则使迭代循环次数N加1;
若满足,则输出精细配准矩阵:T=T1*T2*...Tn。
其中,算法的迭代收敛条件是采用迭代次数收敛,一般情况下迭代最大次数N设置为2或者3即可。具体的算法迭代效果如图6所示。
进一步的,在所述待配准站点云数据的初始配准阶段和精细配准阶段,均采用开放运算语言并行加速模块进行加速。
其中,开放运算语言为OpenCL。
其中,所述待配准站点云数据初始配准阶段的开放运算语言并行加速模块包括:点云邻域搜索模块、点特征直方图计算模块以及初始配准迭代过程模块。
如图7所示,所述待配准站点云数据初始配准阶段的开放运算语言并行加速过程包括:
中央处理器将去噪后的所述基准站点云数据和所述待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,点云邻域并行搜索和计算所述基准站点云表面和所述待配准站点云表面的法线同时执行;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征简易直方图;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图;
并行执行初始配准迭代过程;
图形处理器将初始配准迭代过程中获取的变换矩阵及配准误差发送至中央处理器中,筛选出配准误差最小时对应的变换矩阵Tmin。
其中,中央处理器为CPU,图形处理器为GPU。
其中,所述待配准站点云数据精细配准阶段的开放运算语言并行加速模块包括:每次迭代的邻域搜索模块、最小二乘估计变换矩阵模块以及点云矩阵变换模块。
如图8所示,所述待配准站点云数据精细配准阶段的开放运算语言并行加速过程包括:
中央处理器将初始配准后的基准站点云数据和待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,进行邻域并行搜索;
对变换矩阵变量G、h并行填充,并发送矩阵G、h至中央处理器中,进行楚列斯基分解法求解变换矩阵T_i;
中央处理器发送变换矩阵T_i至图形处理器中进行点云矩阵变换并行计算,同时累乘精细配准矩阵
其中,
为了验证本发明实施例的广泛适用性,本发明实施例准备了三组实验,分别对应不同种类的点云数据。针对配准的应用,实际数据主要分为三大类:线特征明显点云、面特征明显点云以及无明显特征点云。其中线特征明显点云指的是点云的上线性特征比较明显、线性轮廓较为突出,如:石雕文物等等;面特征明显点云指的是点云中有很多突出的平面特征,如:城市建筑物等等;无明显点特征点云指的是点云中没有任何突出的线、面等特征,如:地形、树木等等。因此本发明实施例在实际实验前准备了以上三种种类数据,实际采集点云三维激光扫描设备为RIEGL VZ1000,具体数据分别为:故宫云龙大石雕点云数据、某高校科研楼点云数据、某块地形数据、单棵树木点云。针对每类数据,均采取两站扫描,两站间点云重叠度控制在15%-20%左右。各种不同类型点云数据配准实验如下:
(1)线特征明显点云自动配准实验
故宫云龙大石雕扫描两站数据数量分别为:2,637,262、3,779,241,点云采样分辨率为1.0mm,通过实际算法验证,本发明实施例的算法以及并行加速后的算法均能完成两站点云自动配准过程,如图9(a)(b)(c)(d)所示。图(a)透明区域分别代表两站点云重叠区域。
并行加速前后,平均运行时间对比如下表:
表格1故宫云龙大石雕点云配准运行时间消耗对比
其中,需要说明的是,本发明实施例原始算法是在初始配准阶段和精细配准阶段没有采用并行加速时的算法。
(2)面特征明显点云自动配准实验
某高校扫描两站数据数量分别为:36,074,395、34,973,429,点云采样分辨率为5.0mm,通过实际算法验证,本发明实施例原始算法以及并行加速后算法均能完成两站点云自动配准过程,如图10(a)(b)(c)所示,图(a)透明区域分别代表两站点云主要重叠区域。
算法并行前后,平均运行时间对比如下表:
表格2某高校科研楼点云配准运行时间消耗对比
(3)无明显特征点云自动配准实验
1、某块地形扫描两站数据数量分别为:1,507,694、1,154,906,点云采样分辨率为10.0mm,通过实际算法验证,本发明实施例原始算法以及并行加速后算法均能完成两站点云自动配准过程,如图11(a)(b)(c)(d)所示,(a)图透明区域即代表两站大致重叠区域。
算法并行前后,平均运行时间对比如下表:
表格3某块地形点云配准运行时间消耗对比
2、某单棵树木扫描两站数据数量分别为:1,507,694、1,154,906,点云采样分辨率为6.0mm,通过实际算法验证,本发明实施例算法以及并行加速后算法均能完成两站点云自动配准过程,如图12(a)(b)(c)(d)所示。
表格4某单棵树木点云配准运行时间消耗对比
通过上述配准实验可知,本发明实施例的算法能够很好应用于不同类型的点云数据,具有较广的适用性,且每次实验的配准中误差均在点云采样分辨率的0.5倍以下,满足配准中误差要求。且通过不同类型点云配准运行时间消耗的对比可发现,通过并行加速优化后,点云自动配准算法效率提高了大约8倍左右。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (11)
1.一种基于点特征直方图的点云自动配准方法,其特征在于,包括:
对被测物体进行多站采集获取该被测物体表面的整体点云数据,确定基准站和待配准站,其中,基准站和待配准站为相邻的两站;
从所述基准站点云数据和所述待配准站点云数据中均随机选取m个采样点,均使所选取的m个采样点中的任意三个点都不共线,且每两个采样点之间的距离均满足距离约束原则,其中,m>=3;
对所述待配准站采样点集和所述基准站采样点集进行同名点对的对应匹配,使所述基准站采样点集和所述待配准站采样点集形成的同名点对三角网近似全等;
获取所述待配准站点云数据的变换矩阵T,通过计算点特征直方图差异值获取所述待配准站点云数据经过变换矩阵T变换后与所述基准站点云数据的配准误差;
对所述基准站点云数据和所述待配准站点云数据均进行多次采样,从多次采样获取的配准误差中筛选出配准误差值最小时对应的变换矩阵Tmin,进行点云数据的初始配准;
计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据,并通过迭代最近点算法获得精细配准矩阵,进行点云数据的精细配准。
2.如权利要求1所述的基于点特征直方图的点云自动配准方法,其特征在于,从所述基准站点云数据中和所述待配准站点云数据中均随机选取m个采样点之前还包括:
对所述基准站点云数据和所述待配准站点云数据进行预处理;
计算所述基准站点云表面和所述待配准站点云表面的法线;
计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图。
3.如权利要求1所述的基于点特征直方图的点云自动配准方法,其特征在于,从所述基准站点云数据中和所述待配准站点云数据中均随机选取m个采样点具体包括:
对所述基准站点云数据和所述待配准站点云数据进行采样时,当采集一个新样点时,判断所选取的新样点与已有的采样点是否满足距离约束原则;
若不满足,则新样点无效,继续进行采集,若满足,则继续判断所选取的新样点与已有采样点构建的直线的垂直距离是否小于预设阈值∈dist;
若小于,则所采集的新样点与已有的采样点共线,该新样点无效,重新进行采样;
若不小于,则所采集的新样点有效,继续采集下一个样点。
4.如权利要求2所述的基于点特征直方图的点云自动配准方法,其特征在于,所述通过计算点特征直方图差异值获取所述待配准站点云数据经过变换矩阵T变换后与所述基准站点云数据的配准误差包括:
所述待配准站点云数据P经过变换矩阵T变换后为P',在所述基准站点云数据Q中查找所述待配准站点云数据P'中每个点对应的最近邻域点;
根据获取的点Pi对应的最近邻域点,提取出邻域索引和邻域间距离;
判断邻域间距离是否小于预定的邻域点对距离阈值∈;
若不小于,则继续遍历点云数据P'中下一个点进行邻域间距离判断;
若小于,则计算点Pi与其对应的最近邻域点形成的点对之间的快速直方图差异平方值总和error,并统计至总误差值SumError上,同时累计满足预定的邻域点对距离阈值∈的点对个数nr;
判断nr是否大于零,若大于零,返回配准的点特征中误差若不大于零,则返回配准的分值为无穷大。
5.如权利要求1所述的基于点特征直方图的点云自动配准方法,其特征在于,所述计算所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据包括:
分别计算所述待配准站点云数据P和所述基准站点云数据Q的点云包围盒;
对所述基准站点云数据Q进行空间网格线性划分;
遍历所述待配准站点云数据P中的每个点Pi,并计算每个点Pi对应的网格线性编码值;
在所述基准站点云数据Q的空间网格中查找是否存在与该待配准站点云数据中点Pi的网格线性编码值对应的空间网格;
若不存在,则继续遍历待配准站点云数据P中的点;
若存在,则对所述待配准站点云数据中点Pi的索引进行存储,且标记所述待配准站点云数据中点Pi的网格线性编码值对应所述基准站点云数据Q中的空间网格为重叠网格;
根据存储的所述待配准站点云数据中点Pi的索引,以及所述基准站点云数据Q中重叠的空间网格,得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据。
6.如权利要求5所述的基于点特征直方图的点云自动配准方法,其特征在于,所述通过迭代最近点算法获得精细配准矩阵包括下列步骤:
步骤一、得到所述基准站点云数据和所述待配准站点云数据的重叠区域点云数据后,根据获取的重叠区域点云数据执行迭代最近点算法,得到配准矩阵Ti;
步骤二、将配准矩阵Ti应用至经过初始配准后的所述待配准站点云数据P'上,即P'i=P'i-1*Ti;
步骤三、循环执行步骤一至步骤二,达到迭代循环次数N,其中,N>=2;
步骤四、计算所述待配准站点云数据P'经过配准矩阵Ti配准后与所述基准站点云数据的配准误差ε;
步骤五、判断配准误差ε是否小于等于点云采样分辨率的0.5倍范围内,
若不满足,则使迭代循环次数N加1;
若满足,则输出精细配准矩阵:T=T1*T2*...Tn。
7.如权利要求4所述的基于点特征直方图的点云自动配准方法,其特征在于,在所述待配准站点云数据的初始配准阶段和精细配准阶段,均采用开放运算语言并行加速模块进行加速。
8.如权利要求7所述的基于点特征直方图的点云自动配准方法,其特征在于,所述待配准站点云数据初始配准阶段的开放运算语言并行加速模块包括:点云邻域搜索模块、点特征直方图计算模块以及初始配准迭代过程模块。
9.如权利要求8所述的基于点特征直方图的点云自动配准方法,其特征在于,所述待配准站点云数据初始配准阶段的开放运算语言并行加速过程包括:
中央处理器将去噪后的所述基准站点云数据和所述待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,点云邻域并行搜索和计算所述基准站点云表面和所述待配准站点云表面的法线同时执行;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征简易直方图;
并行计算所述基准站点云数据和所述待配准站点云数据的点特征快速直方图;
并行执行初始配准迭代过程;
图形处理器将初始配准迭代过程中获取的变换矩阵及配准误差发送至中央处理器中,筛选出配准误差最小时对应的变换矩阵Tmin。
10.如权利要求7所述的基于点特征直方图的点云自动配准方法,其特征在于,所述待配准站点云数据精细配准阶段的开放运算语言并行加速模块包括:每次迭代的邻域搜索模块、最小二乘估计变换矩阵模块以及点云矩阵变换模块。
11.如权利要求10所述的基于点特征直方图的点云自动配准方法,其特征在于,所述待配准站点云数据精细配准阶段的开放运算语言并行加速过程包括:
中央处理器将初始配准后的基准站点云数据和待配准站点云数据发送至图形处理器并存储至全局内存中;
中央处理器创建基准站点云邻域搜索KD-tree并发送至图形处理器中,进行邻域并行搜索;
对变换矩阵变量G、h并行填充,并发送矩阵G、h至中央处理器中,进行楚列斯基分解法求解变换矩阵T_i;
中央处理器发送变换矩阵T_i至图形处理器中进行点云矩阵变换并行计算,同时累乘精细配准矩阵
其中,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610371302.9A CN105976312B (zh) | 2016-05-30 | 2016-05-30 | 基于点特征直方图的点云自动配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610371302.9A CN105976312B (zh) | 2016-05-30 | 2016-05-30 | 基于点特征直方图的点云自动配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105976312A true CN105976312A (zh) | 2016-09-28 |
CN105976312B CN105976312B (zh) | 2019-03-01 |
Family
ID=57009818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610371302.9A Expired - Fee Related CN105976312B (zh) | 2016-05-30 | 2016-05-30 | 基于点特征直方图的点云自动配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105976312B (zh) |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108470358A (zh) * | 2018-02-13 | 2018-08-31 | 南京旅游职业学院 | 一种基于二阶震荡人工蜂群算法的点云配准方法 |
CN108648167A (zh) * | 2018-03-06 | 2018-10-12 | 深圳市菲森科技有限公司 | 一种口内扫描的三维扫描方法 |
CN109345620A (zh) * | 2018-08-13 | 2019-02-15 | 浙江大学 | 融合快速点特征直方图的改进icp待测物体点云拼接方法 |
CN109859254A (zh) * | 2019-02-28 | 2019-06-07 | 北京百度网讯科技有限公司 | 用于发送信息的方法和装置 |
CN110136181A (zh) * | 2019-05-17 | 2019-08-16 | 百度在线网络技术(北京)有限公司 | 用于生成信息的方法和装置 |
WO2019161558A1 (en) * | 2018-02-26 | 2019-08-29 | Intel Corporation | Method and system of point cloud registration for image processing |
CN110322492A (zh) * | 2019-07-03 | 2019-10-11 | 西北工业大学 | 一种基于全局优化的空间目标三维点云配准方法 |
CN110675507A (zh) * | 2019-08-27 | 2020-01-10 | 腾讯科技(深圳)有限公司 | 一种三维模型处理方法、装置和存储介质 |
CN110874849A (zh) * | 2019-11-08 | 2020-03-10 | 安徽大学 | 一种基于局部变换一致的非刚性点集配准方法 |
CN112085793A (zh) * | 2020-09-04 | 2020-12-15 | 上海理工大学 | 一种基于组合透镜组的三维成像扫描系统及点云配准方法 |
CN112446844A (zh) * | 2020-11-27 | 2021-03-05 | 广东电网有限责任公司肇庆供电局 | 一种基于点云特征提取与配准融合方法 |
CN112650790A (zh) * | 2020-12-28 | 2021-04-13 | 北京字节跳动网络技术有限公司 | 目标点云平面确定方法,装置,电子设备及存储介质 |
CN112950682A (zh) * | 2021-02-07 | 2021-06-11 | 南京航空航天大学 | 一种基于小波变换与尺度表征的跨尺度点云配准算法 |
CN113379807A (zh) * | 2021-05-27 | 2021-09-10 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种移动平台点云数据的配准方法和系统 |
CN113496505A (zh) * | 2020-04-03 | 2021-10-12 | 广州极飞科技股份有限公司 | 图像配准方法、装置、多光谱相机、无人设备及存储介质 |
CN114067555A (zh) * | 2020-08-05 | 2022-02-18 | 北京万集科技股份有限公司 | 多基站数据的配准方法、装置、服务器和可读存储介质 |
CN115147551A (zh) * | 2022-07-18 | 2022-10-04 | 南昌工程学院 | 远程地面激光扫描技术提取监测对象表面同名特征点的方法 |
CN115283172A (zh) * | 2022-07-22 | 2022-11-04 | 颖态智能技术(上海)有限公司 | 一种基于点云处理的机器人自动喷涂方法 |
CN115482195A (zh) * | 2022-08-03 | 2022-12-16 | 西南交通大学 | 一种基于三维点云的列车部件变形检测方法 |
US11535400B2 (en) | 2020-05-09 | 2022-12-27 | Nanjing University Of Aeronautics And Astronautics | Fairing skin repair method based on measured wing data |
CN116188546A (zh) * | 2022-12-30 | 2023-05-30 | 北京斯年智驾科技有限公司 | 激光点云配准方法、装置、电子设备及存储介质 |
CN116246121A (zh) * | 2023-05-12 | 2023-06-09 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
CN116468762A (zh) * | 2023-04-17 | 2023-07-21 | 南京邮电大学 | 一种变电站点云配准方法、装置、电子设备及存储介质 |
CN116721137A (zh) * | 2023-08-08 | 2023-09-08 | 北京爱康宜诚医疗器材有限公司 | 配准方法及装置、存储介质及电子设备 |
CN116777961A (zh) * | 2023-06-26 | 2023-09-19 | 安徽开源路桥有限责任公司 | 一种利用kd树搜索的并行化点云配准的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646156A (zh) * | 2013-12-30 | 2014-03-19 | 北京建筑大学 | 一种基于球标靶探测的激光点云数据自动配准方法 |
CN103679741A (zh) * | 2013-12-30 | 2014-03-26 | 北京建筑大学 | 一种基于三维线特征的激光点云数据自动配准方法 |
US20150104105A1 (en) * | 2013-10-14 | 2015-04-16 | Fu Tai Hua Industry (Shenzhen) Co., Ltd. | Computing device and method for jointing point clouds |
-
2016
- 2016-05-30 CN CN201610371302.9A patent/CN105976312B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150104105A1 (en) * | 2013-10-14 | 2015-04-16 | Fu Tai Hua Industry (Shenzhen) Co., Ltd. | Computing device and method for jointing point clouds |
CN103646156A (zh) * | 2013-12-30 | 2014-03-19 | 北京建筑大学 | 一种基于球标靶探测的激光点云数据自动配准方法 |
CN103679741A (zh) * | 2013-12-30 | 2014-03-26 | 北京建筑大学 | 一种基于三维线特征的激光点云数据自动配准方法 |
Non-Patent Citations (1)
Title |
---|
沈萦华: "基于法向特征直方图的点云配准算法", 《光学精密工程》 * |
Cited By (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108470358A (zh) * | 2018-02-13 | 2018-08-31 | 南京旅游职业学院 | 一种基于二阶震荡人工蜂群算法的点云配准方法 |
CN108470358B (zh) * | 2018-02-13 | 2021-11-30 | 南京大学 | 一种基于二阶震荡人工蜂群算法的点云配准方法 |
WO2019161558A1 (en) * | 2018-02-26 | 2019-08-29 | Intel Corporation | Method and system of point cloud registration for image processing |
US11423508B2 (en) | 2018-02-26 | 2022-08-23 | Intel Corporation | Method and system of point cloud registration for image processing |
CN108648167A (zh) * | 2018-03-06 | 2018-10-12 | 深圳市菲森科技有限公司 | 一种口内扫描的三维扫描方法 |
CN108648167B (zh) * | 2018-03-06 | 2021-10-01 | 深圳市菲森科技有限公司 | 一种口内扫描的三维扫描方法 |
CN109345620B (zh) * | 2018-08-13 | 2022-06-24 | 浙江大学 | 融合快速点特征直方图的改进icp待测物体点云拼接方法 |
CN109345620A (zh) * | 2018-08-13 | 2019-02-15 | 浙江大学 | 融合快速点特征直方图的改进icp待测物体点云拼接方法 |
CN109859254A (zh) * | 2019-02-28 | 2019-06-07 | 北京百度网讯科技有限公司 | 用于发送信息的方法和装置 |
CN110136181B (zh) * | 2019-05-17 | 2021-08-20 | 百度在线网络技术(北京)有限公司 | 用于生成信息的方法和装置 |
CN110136181A (zh) * | 2019-05-17 | 2019-08-16 | 百度在线网络技术(北京)有限公司 | 用于生成信息的方法和装置 |
CN110322492A (zh) * | 2019-07-03 | 2019-10-11 | 西北工业大学 | 一种基于全局优化的空间目标三维点云配准方法 |
CN110322492B (zh) * | 2019-07-03 | 2022-06-07 | 西北工业大学 | 一种基于全局优化的空间三维点云配准方法 |
CN110675507A (zh) * | 2019-08-27 | 2020-01-10 | 腾讯科技(深圳)有限公司 | 一种三维模型处理方法、装置和存储介质 |
CN110874849B (zh) * | 2019-11-08 | 2023-04-18 | 安徽大学 | 一种基于局部变换一致的非刚性点集配准方法 |
CN110874849A (zh) * | 2019-11-08 | 2020-03-10 | 安徽大学 | 一种基于局部变换一致的非刚性点集配准方法 |
CN113496505B (zh) * | 2020-04-03 | 2022-11-08 | 广州极飞科技股份有限公司 | 图像配准方法、装置、多光谱相机、无人设备及存储介质 |
CN113496505A (zh) * | 2020-04-03 | 2021-10-12 | 广州极飞科技股份有限公司 | 图像配准方法、装置、多光谱相机、无人设备及存储介质 |
US11535400B2 (en) | 2020-05-09 | 2022-12-27 | Nanjing University Of Aeronautics And Astronautics | Fairing skin repair method based on measured wing data |
CN114067555A (zh) * | 2020-08-05 | 2022-02-18 | 北京万集科技股份有限公司 | 多基站数据的配准方法、装置、服务器和可读存储介质 |
CN112085793A (zh) * | 2020-09-04 | 2020-12-15 | 上海理工大学 | 一种基于组合透镜组的三维成像扫描系统及点云配准方法 |
CN112085793B (zh) * | 2020-09-04 | 2022-07-05 | 上海理工大学 | 一种基于组合透镜组的三维成像扫描系统及点云配准方法 |
CN112446844A (zh) * | 2020-11-27 | 2021-03-05 | 广东电网有限责任公司肇庆供电局 | 一种基于点云特征提取与配准融合方法 |
CN112446844B (zh) * | 2020-11-27 | 2024-01-23 | 广东电网有限责任公司肇庆供电局 | 一种基于点云特征提取与配准融合方法 |
CN112650790A (zh) * | 2020-12-28 | 2021-04-13 | 北京字节跳动网络技术有限公司 | 目标点云平面确定方法,装置,电子设备及存储介质 |
CN112650790B (zh) * | 2020-12-28 | 2023-10-13 | 北京字节跳动网络技术有限公司 | 目标点云平面确定方法,装置,电子设备及存储介质 |
CN112950682A (zh) * | 2021-02-07 | 2021-06-11 | 南京航空航天大学 | 一种基于小波变换与尺度表征的跨尺度点云配准算法 |
CN113379807A (zh) * | 2021-05-27 | 2021-09-10 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种移动平台点云数据的配准方法和系统 |
CN115147551A (zh) * | 2022-07-18 | 2022-10-04 | 南昌工程学院 | 远程地面激光扫描技术提取监测对象表面同名特征点的方法 |
CN115283172A (zh) * | 2022-07-22 | 2022-11-04 | 颖态智能技术(上海)有限公司 | 一种基于点云处理的机器人自动喷涂方法 |
CN115283172B (zh) * | 2022-07-22 | 2024-05-17 | 颖态智能技术(上海)有限公司 | 一种基于点云处理的机器人自动喷涂方法 |
CN115482195B (zh) * | 2022-08-03 | 2023-06-20 | 西南交通大学 | 一种基于三维点云的列车部件变形检测方法 |
CN115482195A (zh) * | 2022-08-03 | 2022-12-16 | 西南交通大学 | 一种基于三维点云的列车部件变形检测方法 |
CN116188546B (zh) * | 2022-12-30 | 2023-07-14 | 北京斯年智驾科技有限公司 | 激光点云配准方法、装置、电子设备及存储介质 |
CN116188546A (zh) * | 2022-12-30 | 2023-05-30 | 北京斯年智驾科技有限公司 | 激光点云配准方法、装置、电子设备及存储介质 |
CN116468762A (zh) * | 2023-04-17 | 2023-07-21 | 南京邮电大学 | 一种变电站点云配准方法、装置、电子设备及存储介质 |
CN116246121B (zh) * | 2023-05-12 | 2023-08-11 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
CN116246121A (zh) * | 2023-05-12 | 2023-06-09 | 山东科技大学 | 一种基于改进最远点采样的点云数据处理方法 |
CN116777961A (zh) * | 2023-06-26 | 2023-09-19 | 安徽开源路桥有限责任公司 | 一种利用kd树搜索的并行化点云配准的方法 |
CN116721137A (zh) * | 2023-08-08 | 2023-09-08 | 北京爱康宜诚医疗器材有限公司 | 配准方法及装置、存储介质及电子设备 |
CN116721137B (zh) * | 2023-08-08 | 2023-10-27 | 北京爱康宜诚医疗器材有限公司 | 配准方法及装置、存储介质及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN105976312B (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105976312B (zh) | 基于点特征直方图的点云自动配准方法 | |
CN109410256B (zh) | 基于互信息的点云与影像自动高精度配准方法 | |
CN110443836B (zh) | 一种基于平面特征的点云数据自动配准方法及装置 | |
CN108846888B (zh) | 一种古木建筑构件精细尺寸信息自动化提取方法 | |
CN104331699B (zh) | 一种三维点云平面化快速搜索比对的方法 | |
CN106023298B (zh) | 基于局部泊松曲面重建的点云刚性配准方法 | |
CN109919984A (zh) | 一种基于局部特征描述子的点云自动配准方法 | |
CN109186608B (zh) | 一种面向重定位的稀疏化三维点云地图生成方法 | |
Ge et al. | Surface-based matching of 3D point clouds with variable coordinates in source and target system | |
CN116721228B (zh) | 一种基于低密度点云的建筑物高程提取方法及系统 | |
CN115222884A (zh) | 一种基于人工智能的空间对象分析及建模优化方法 | |
CN114463396B (zh) | 一种利用平面形状和拓扑图投票的点云配准方法 | |
Farella et al. | Quality features for the integration of terrestrial and UAV images | |
CN117994484A (zh) | 一种大型建筑预制构件三维图像拼接装置及方法 | |
CN112231848B (zh) | 一种构建车辆喷涂模型的方法及系统 | |
Withers et al. | Modelling scene change for large-scale long term laser localisation | |
CN117706577A (zh) | 基于激光雷达三维点云算法的船舶尺寸测量方法 | |
CN113870326A (zh) | 一种基于图像和三维点云配准的结构损伤映射、量化及可视化方法 | |
CN117765039A (zh) | 点云粗配准方法、装置及设备 | |
CN113628170A (zh) | 一种基于深度学习的激光线提取方法及系统 | |
Dos Santos et al. | Building boundary extraction from LiDAR data using a local estimated parameter for alpha shape algorithm | |
CN117315518A (zh) | 一种增强现实目标初始注册方法及系统 | |
JP6761388B2 (ja) | 推定装置及びプログラム | |
CN116778266A (zh) | 多尺度邻域扩散的遥感点云投影图像处理方法 | |
Bethmann et al. | Object-based semi-global multi-image matching |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190301 |