CN115205354B - 基于ransac和icp点云配准的相控阵激光雷达成像方法 - Google Patents
基于ransac和icp点云配准的相控阵激光雷达成像方法 Download PDFInfo
- Publication number
- CN115205354B CN115205354B CN202210727991.8A CN202210727991A CN115205354B CN 115205354 B CN115205354 B CN 115205354B CN 202210727991 A CN202210727991 A CN 202210727991A CN 115205354 B CN115205354 B CN 115205354B
- Authority
- CN
- China
- Prior art keywords
- point
- point cloud
- matrix
- registration
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 62
- 230000009466 transformation Effects 0.000 claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 27
- 230000008569 process Effects 0.000 claims abstract description 18
- 238000012216 screening Methods 0.000 claims abstract description 5
- 239000013598 vector Substances 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 9
- 238000013519 translation Methods 0.000 claims description 9
- 230000004927 fusion Effects 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 239000003102 growth factor Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000000007 visual effect Effects 0.000 abstract description 9
- 238000009616 inductively coupled plasma Methods 0.000 abstract description 4
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008570 general process Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- 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
-
- 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/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及相控阵激光雷达快速成像技术,尤其涉及一种基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,用于解决实际采集三维点云数据过程中单一视角下场景范围小、三维目标信息不完整的问题。根据相控阵激光雷达扫描过程中不同视角下的点云存在目标重合区域的情况,先提取点云的几何特征,寻找重合部分,利用RANSAC算法筛除误匹配点信息,从而完成不同视角的粗配准;随后根据粗配准得到的变换矩阵,进一步利用ICP算法通过不断迭代继续优化变换矩阵得到最佳变换矩阵,完成点云配准,最终将不同视角下的点云通过该变换矩阵进行拼接,获得场景更大、信息更完整的三维点云数据。
Description
技术领域
本发明涉及相控阵激光雷达快速成像技术,尤其涉及一种基于随机采样一致性(Random Sample Consensus,RANSAC)和迭代最近点(Iterative Closed Point,ICP)点云配准的相控阵激光雷达成像方法。
背景技术
扫描式激光雷达是目前获取三维点云的主要设备之一,目前可以分为两大类:机械扫描式激光雷达和相控阵扫描式激光雷达。机械扫描式激光雷达利用电机转动带动探测器旋转,从而扫描目标整体,这种方法受电机机械惯性的影响,无法满足雷达成像中的实时性要求,且功耗大。相控阵扫描式激光雷达借鉴了微波雷达成像技术中利用相控阵天线实现快速扫描的方法,成功解决了机械式扫描激光雷达中实时性低的问题,可实现非机械式多角度快速扫描,且体积更小,更易集成化。
在利用相控阵激光雷达采集三维点云数据的过程中,受到雷达视角限制、目标自身遮挡以及其他物体遮挡目标等因素的影响,采集到的三维点云目标一般不完整,不能很好的表征整个三维场景/目标。为了解决这一问题,需要从不同视角下采集三维点云数据,将采集到的数据通过变换矩阵关联到同一坐标系下,该关联过程称为点云配准,然后对配准后的点云进行拼接,从而获得场景更大、信息更完整的三维点云数据,该过程称为点云视场拼接,点云视场拼接中最为关键的环节为点云配准。
发明内容
针对上述现有技术中存在的问题,本发明提出一种基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,用于解决实际场景中相控阵激光雷达在采集三维点云数据过程中单一视角下场景范围小、三维目标信息不完整的问题。根据相控阵激光雷达扫描过程中不同视角下的点云存在目标重合区域的情况,先提取点云的几何特征,寻找重合部分,利用RANSAC算法筛除误匹配点信息,从而完成不同视角的粗配准;随后根据粗配准得到的变换矩阵,进一步利用ICP算法通过不断迭代继续优化变换矩阵得到最佳变换矩阵,完成点云配准,最终将不同视角下的点云通过该变换矩阵进行拼接和融合,获得场景更大、信息更完整的三维点云数据。
本发明采用的技术方案为:一种基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,该方法包括以下步骤:
S1:相控阵激光雷达获取点云。
在本发明中,通过相控阵激光雷达采集点云。与光学相机不同,相控阵激光雷达通过向目标发射激光光束后收集目标反射的激光信息来获取点云。在单一视角下,相控阵激光雷达采集的数据视场有限、信息量不足,因此需要将在多个视角下采集到的点云数据进行配准拼接,以得到场景更大、信息更丰富的点云数据。在本发明中利用相控阵激光雷达首先获取两帧点云P1和P2,P1、P2之间的重叠范围大于60%。
S2:点云体素化下采样。
对于分别包含有N和M个点的三维点云P1和P2,由于原始点云P1、P2中点数较多,如果直接在原始点云上进行配准计算的话计算量繁重。针对这一问题,将P1、P2进行体素化下采样以减少后续配准算法的计算量。以P1体素化下采样为例,具体步骤如下:
S2.1为P1创建大小为l×w×h、体积为v的三维包围框,设置体素的边长为L,所创建的三维包围框根据体素的边长L被划分为Q个体素:
S2.2体素的大小和数量是体素化下采样的重点,同一点云数据若体素边长小、数量多,则后续数据处理耗时长,但能够表征的点云信息丰富;若体素边长大、数量少,则后续数据处理耗时短,但丢失的点云信息多。为了得到最佳的体素边长和数量,根据点云密度γ和点云下采样率ρ估计体素的边长L′:
S2.4用每个体素中所有点云的重心点(xm,ym,zm)表示当前体素:
其中n表示体素中点云的个数,i表示第i个点,(xi,yi,zi)表示第i个点的坐标值。
S3:计算点云法向量。
点云体素化下采样后,可以根据点的法向量之间的差异找到两个点云中相似区域,这就需要提取点云的法向量特征。以P1为例,点云法向量的计算过程具体如下:
S3.1采用KD-Tree搜索得到点云P1中某一点p1,i的K邻域点集p1,j,对其构建协方差矩阵C:
计算得出协方差矩阵C的特征值e和特征向量λ,协方差矩阵的最小特征值emin对应的特征向量λi即为法向量f。
S4:计算点云的快速点特征直方图描述子(Fast Point Feature Histograms,FPFH)特征。
为了更加准确的寻找两个点云的相似区域,需要进一步提取点云的局部几何特征:快速点特征直方图描述子(Fast Point Feature Histograms,FPFH)。点云的FPFH特征是在简化点特征直方图描述子(Simplified Point Feature Histogram,SPFH)特征上进一步优化的一种特征描述符,用于描述点云的局部几何特征。以P1为例,点云的FPFH特征计算步骤具体如下:
S4.1利用KD-Tree搜索点云P1中某一点p1,i的K邻域,得到K近邻点集O,通过参数化p1,i与点集O中任意一点在空间上的差别和用一个计算得到的33维的向量直方图来描述K近邻点的空间几何特征。
假设p1,i与其K近邻点集O中的一点O1及对应的法向量f1、f2在局部坐标系vuw的空间位置和空间关系为:
U=f1 (7)
V=(p1,i-O1)·U (8)
W=U·V (9)
α=V·f2 (10)
θ=arctan(W·f1,U·f2) (12)
V、U、W表示局部坐标系的坐标轴方向上的单位向量,α表示f2与局部坐标系V轴的夹角,β表示O1向p1,i方向的单位向量与U轴作内积的夹角,θ表示f2在WU平面投影与U轴的夹角。将点p1,i和其K近邻点集O中的点组成的所有点对求出(α,β,θ),将(α,β,θ)按照11个区间划分,得到3×11=33维的SPFH,根据p1,i及其K近邻点的SPFH特征计算FPFH:
其中SPFH(p1,i)表示点云P1中的一点p1,i的SPFH特征,SPFH(Ok)表示p1,i的K近邻点Ok的SPFH特征,ωk表示p1,i与K近邻点Ok之间的距离。公式(13)最终得到33维的FPFH特征。
S5:RANSAC筛除误匹配信息。
在S4中得到的FPFH特征可以很好的描述点云之间的局部几何关系,但是由于噪声和算法自身性能的影响,计算得到的FPFH特征中包含一些错误的特征点会影响后续配准算法的性能,所以需要利用RANSAC算法筛除这些错误的特征点。
假设原始点云P1中有一点p1,i,待配准点云P2中有一对应点p2,i,在经过S4步骤得到FPFH特征后,通过利用RANSAC算法筛除误匹配的FPFH特征得到最佳的粗配准变换矩阵。利用RANSAC算法筛除FPFH特征中误匹配点的流程为:
S5.1在原始点云P1中随机选择不共线的4个点p1,i,其中i=0,1,2,3,在待配准点云P2中寻找与这4个点的FPFH特征最相似的点p2,i,将其作为这四个点的对应点,对应点对的两个点p1,i、p2,i满足以下关系:
p1,i=Rp2,i+T (14)
S5.2一般来说,正确的对应点应使下式L为0:
由于数据中含有噪声、对应点对之间的空间关系{R|T}变换矩阵求解不准确及计算误差等因素的影响,L通常情况下不为0,为了得到最佳的{R|T}变换矩阵,一般通过优化算法循环迭代公式(15),当公式(15)的L为最小值时,此时{R|T}为最佳变换矩阵。
{R|T}变换矩阵通过单元四元数法求解。单元四元数由四矢量a0,a1,a2,a3构成:
a=a0+a1I+a2J+a3L=[a0,a1,a2,a3]T (16)
其中a0≥0且a0 2+a1 2+a2 2+a3 2=1,I、J、L为虚部单位。通过单元四元数求解(Closed-form solution of absolute orientation using unit quaternions,BKP Horn;AMethod for Registration of 3-D Shapes,PJ Besl)的旋转矩阵R和平移矩阵T分别为公式(15)、(16):
T=[a4,a5,a6]t (18)
S6:ICP点云精细配准。
ICP点云精细配准是在S5的基础上进一步精细化变换矩阵{R|T},其基本流程为:
S6.1搜索最近点以确定对应点对:在原始点云P1中选取一个点p1,i,计算待配准点云P2中所有点与p1,i的距离,得到与p1,i距离最近的点p2,i,这一组距离最近的点构成一组对应点对对集合P1中所有的点都在集合Y中搜索距离最近的对应点。
S6.2求解配准变换:求解配准变换即计算点云在不同坐标系中产生的变换矩阵,本文所说的点云数据相对原始数据不会发生形变,只存在平移和旋转。变换矩阵即采用单位四元素法求解的3×3旋转矩阵和3×1平移矩阵。
S6.3应用变换:对三维点云Q根据刚体变换矩阵按照公式(14)进行变换。
S6.4重复迭代:重复迭代以上三个步骤,设置ICP收敛阈值TH,通常,TH=0.02。当L<TH时得到最优刚体变换矩阵,此时得到最佳配准效果,迭代停止。
S7:点云融合。
点云融合是在点云配准完成后将两个点云融合为一个点云,融合具体步骤为:
S7.1将待配准点云P2与点云配准计算得到的变换矩阵{R|T}相乘,然后将其和原始点云P1拼接得到配准后的点云P3:
P3=[P1,P2·{R|T}] (23)
S7.2计算点云P3中两个相邻点之间的距离d,设置距离阈值dth,通常,dth=0.02;
S7.3判断d和dth的关系,若d<dth,说明这两个相邻点重复,去除掉其中的一点;若d≥dth,说明两个相邻点没有重复。
本发明具有以下有益效果:通过本发明,根据相控阵激光雷达扫描成像过程中在单一视场视角小、扫描得到的三维场景信息有限的问题,通过将激光雷达扫描过程中各个视角下的小视场配准拼接成一个大的点云场景,通过RANSAC算法进行点云粗配准,通过ICP算法在粗配准的前提下进一步进行点云精细配准。本发明对促进大视场下相控阵激光雷达快速成像向工程化、使用化方向发展有重要作用。
附图说明
图1是本发明的处理流程图;
图2基于RANSAC的点云粗配准流程图;
图3是基于ICP的点云精细配准流程图;
图4是最小二乘法和RANSAC算法拟合数据效果;
图5是点云Q中某一点的K近邻点O集合中的两点O1、O2的局部坐标系;
图6是点云配准Simu Data数据集示例;
图7在Simu Data上RANSAC算法的实验结果;
图8在Simu Data上ICP算法的实验结果;
图9在Actu Data上RANSAC算法的实验结果;
图10在Actu Data上ICP算法的实验结果;
具体实施方式
下面结合附图对本发明进行进一步说明:
图1为本发明总处理流程。
本发明所述一种基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,该方法包括以下步骤:
S1:相控阵激光雷达获取点云;
S2:点云体素化下采样;
S3:计算点云法向量;
S4:计算点云的快速点特征直方图描述子(Fast Point Feature Histograms,FPFH)特征;
S5:RANSAC筛除误匹配信息;
S6:ICP点云精细配准;
S7:点云融合。
下面以仿真数据来检验本发明的有益效果。
本发明在两个数据集上测试算法效果,数据集中包含的目标以及各个目标信息如表1所示。第一个数据集是通过仿真生成大场景点云,在本发明中称之为Simu Data,通过将Simu Data中的一个大场景点云切分为若干个小场景点云来模拟相控阵激光雷达的扫描过程;第二个数据集是通过Realsense L515在实际场景中获取军事目标模型的深度图像,将深度图转换为点云数据,在本发明中称之为Actu Data。对于Simu Data,本发明将切分后小场景点云绕Z轴旋转度,使之模拟相控阵激光雷达的扫描过程,然后将小场景点云进行配准拼接;对于Actu Data,本发明将整个点云场景绕Z轴旋转度,然后对其进行配准。旋转角度服从高斯分布(5.73,0.29)。图6展示了切分后的用于视场拼接的两个仿真场景。
表1 Simu Data场景信息
表2 Actu Data场景信息
在本发明中,RANSAC算法的超参数设定为:点云搜索半径为10cm;对应近邻点数为30;点云之间的距离阈值为2cm;随机采样点数为4。ICP算法超参数设定为:点云之间的距离阈值为0.2cm。为验证算法的抗噪声性能,本文向p1,...,i、p′1,...i、G和G‘中分别添加不同的高斯噪声(Gauss Noise,GN),GN满足正态分布(μ,σ),其中μ=0,在[0,0.01]中等间隔选取15个值作为标准差σ。
图7和图9分别展示了基于RANSAC算法的点云粗配准技术在Simu Data和ActuData上的配准效果,从图中可以看出,随着高斯噪声标准差σ的增大,点云配准的角度均方误差逐渐增大,点云配准效果逐渐变差。在高斯噪声的标准差小于0.005时,噪声对RANSAC算法的影响很小,可以得到较好的变换矩阵。图8和图10分别展示了基于ICP算法的点云精细配准技术在Simu Data和ActuData上的配准效果,从图中可以看出,随着高斯噪声标准差σ的增大,点云配准的角度均方误差逐渐增大,点云配准效果逐渐变差。从图7至图10可知,点云粗配准的角度均方误差大于点云精细配准的均方误差,说明点云精细配准可以在点云粗配准的前提下进一步优化变化矩阵,在点云精细配准配准后,角度均方误差保持了一个较小的值,证明了本发明方法的配准拼接能力。
Claims (3)
1.一种基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,其特征在于,该方法包括以下步骤:
S1:相控阵激光雷达获取点云;
S2:点云体素化下采样
对于分别包含有N和M个点的三维点云P1和P2,P1体素化下采样具体步骤如下:
S2.1为P1创建大小为l×w×h、体积为v的三维包围框,设置体素的边长为L,所创建的三维包围框根据体素的边长L被划分为Q个体素:
S2.2为了得到最佳的体素边长和数量,根据点云密度γ和点云下采样率ρ估计体素的边长L′:
S2.4用每个体素中所有点云的重心点(xm,ym,zm)表示当前体素:
其中n表示体素中点云的个数,i表示第i个点,(xi,yi,zi)表示第i个点的坐标值;
S3:计算点云法向量
P1点云法向量的计算过程具体如下:
S3.1采用KD-Tree搜索得到点云P1中某一点p1,i的K邻域点集p1,j,对其构建协方差矩阵C:
计算得出协方差矩阵C的特征值e和特征向量λ,协方差矩阵的最小特征值emin对应的特征向量λi即为法向量f;
S4:计算点云的快速点特征直方图描述子特征
P1点云的快速点特征直方图描述子特征计算步骤具体如下:
S4.1利用KD-Tree搜索点云P1中某一点p1,i的K邻域,得到K近邻点集O,通过参数化p1,i与点集O中任意一点在空间上的差别和用一个计算得到的33维的向量直方图来描述K近邻点的空间几何特征;
假设p1,i与其K近邻点集O中的一点O1及对应的法向量f1、f2在局部坐标系vuw的空间位置和空间关系为:
U=f1 (7)
V=(p1,i-O1)·U (8)
W=U·V (9)
α=V·f2 (10)
θ=arctan(W·f1,U·f2) (12)
V、U、W表示局部坐标系的坐标轴方向上的单位向量,α表示f2与局部坐标系V轴的夹角,β表示O1向p1,i方向的单位向量与U轴作内积的夹角,θ表示f2在WU平面投影与U轴的夹角;将点p1,i和其K近邻点集O中的点组成的所有点对求出(α,β,θ),将(α,β,θ)按照11个区间划分,得到3×11=33维的SPFH,根据p1,i及其K近邻点的SPFH特征计算FPFH:
其中SPFH(p1,i)表示点云P1中的一点p1,i的SPFH特征,SPFH(Ok)表示p1,i的K近邻点Ok的SPFH特征,ωk表示p1,i与K近邻点Ok之间的距离;公式(13)最终得到33维的FPFH特征;
S5:RANSAC筛除误匹配信息
利用RANSAC算法筛除FPFH特征中误匹配点的流程为:
S5.1在原始点云P1中随机选择不共线的4个点p1,i,其中i=0,1,2,3,在待配准点云P2中寻找与这4个点的FPFH特征最相似的点p2,i,将其作为这四个点的对应点,对应点对的两个点p1,i、p2,i满足以下关系:
p1,i=Rp2,i+T (14)
S5.2正确的对应点应使下式L为0:
由于数据中含有噪声、对应点对之间的空间关系{R|T}变换矩阵求解不准确及计算误差等因素的影响,L通常情况下不为0,为了得到最佳的{R|T}变换矩阵,通过优化算法循环迭代公式(15),当公式(15)的L为最小值时,此时{R|T}为最佳变换矩阵;
{R|T}变换矩阵通过单元四元数法求解,单元四元数由四矢量a0,a1,a2,a3构成:
a=a0+a1I+a2J+a3L=[a0,a1,a2,a3]T (16)
其中a0≥0且a0 2+a1 2+a2 2+a3 2=1,I、J、L为虚部单位;通过单元四元数求解的旋转矩阵R和平移矩阵T分别为公式(15)、(16):
T=[a4,a5,a6]t (18)
S6:ICP点云精细配准
ICP点云精细配准是在S5的基础上进一步精细化变换矩阵{R|T},其基本流程为:
S6.1搜索最近点以确定对应点对:在原始点云P1中选取一个点p1,i,计算待配准点云P2中所有点与p1,i的距离,得到与p1,i距离最近的点p2,i,这一组距离最近的点构成一组对应点对对集合P1中所有的点都在集合Y中搜索距离最近的对应点;
S6.2求解配准变换:求解配准变换即计算点云在不同坐标系中产生的变换矩阵,变换矩阵即采用单位四元素法求解的3×3旋转矩阵和3×1平移矩阵;
S6.3应用变换:对三维点云Q根据刚体变换矩阵按照公式(14)进行变换;
S6.4重复迭代:重复迭代以上三个步骤,设置ICP收敛阈值TH,当L<TH时得到最优刚体变换矩阵,此时得到最佳配准效果,迭代停止;
S7:点云融合
点云融合是在点云配准完成后将两个点云融合为一个点云,融合具体步骤为:
S7.1将待配准点云P2与点云配准计算得到的变换矩阵{R|T}相乘,然后将其和原始点云P1拼接得到配准后的点云P3:
P3=[P1,P2·{R|T}] (23)
S7.2计算点云P3中两个相邻点之间的距离d,设置距离阈值dth;
S7.3判断d和dth的关系,若d<dth,说明这两个相邻点重复,去除掉其中的一点;若d≥dth,说明两个相邻点没有重复。
2.一种根据权利要求1所述基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,其特征在于:S6.4中,ICP收敛阈值TH=0.02。
3.一种根据权利要求1所述基于RANSAC和ICP点云配准的相控阵激光雷达成像方法,其特征在于:S7.2中,距离阈值dth=0.02。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210727991.8A CN115205354B (zh) | 2022-06-23 | 2022-06-23 | 基于ransac和icp点云配准的相控阵激光雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210727991.8A CN115205354B (zh) | 2022-06-23 | 2022-06-23 | 基于ransac和icp点云配准的相控阵激光雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115205354A CN115205354A (zh) | 2022-10-18 |
CN115205354B true CN115205354B (zh) | 2023-04-07 |
Family
ID=83579099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210727991.8A Active CN115205354B (zh) | 2022-06-23 | 2022-06-23 | 基于ransac和icp点云配准的相控阵激光雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115205354B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117689536B (zh) * | 2024-02-01 | 2024-05-10 | 浙江华是科技股份有限公司 | 激光雷达拼接配准方法、***、装置及计算机存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104463953A (zh) * | 2014-11-11 | 2015-03-25 | 西北工业大学 | 基于惯性测量单元与rgb-d传感器的三维重建方法 |
CN107038717A (zh) * | 2017-04-14 | 2017-08-11 | 东南大学 | 一种基于立体栅格自动分析3d点云配准误差的方法 |
CN107123161A (zh) * | 2017-06-14 | 2017-09-01 | 西南交通大学 | 一种基于narf和fpfh的接触网零全网三维重建方法 |
CN108133458A (zh) * | 2018-01-17 | 2018-06-08 | 视缘(上海)智能科技有限公司 | 一种基于目标物体空间点云特征的自动拼接方法 |
CN111009005A (zh) * | 2019-11-27 | 2020-04-14 | 天津大学 | 几何信息与光度信息相结合的场景分类点云粗配准方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6348093B2 (ja) * | 2015-11-06 | 2018-06-27 | ファナック株式会社 | 入力データから検出対象物の像を検出する画像処理装置および方法 |
CN113240603A (zh) * | 2020-12-28 | 2021-08-10 | 武汉纺织大学 | 基于点云数据的三维人体测量***和测量方法 |
CN113048871B (zh) * | 2021-03-11 | 2022-04-05 | 中国人民解放军国防科技大学 | 一种基于谐波分量的电容位移检测非线性实时校准方法 |
CN114359042A (zh) * | 2021-11-30 | 2022-04-15 | 深圳市纵维立方科技有限公司 | 点云拼接方法和装置、三维扫描仪及电子设备 |
CN114170279A (zh) * | 2021-11-30 | 2022-03-11 | 哈尔滨工程大学 | 一种基于激光扫描的点云配准方法 |
-
2022
- 2022-06-23 CN CN202210727991.8A patent/CN115205354B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104463953A (zh) * | 2014-11-11 | 2015-03-25 | 西北工业大学 | 基于惯性测量单元与rgb-d传感器的三维重建方法 |
CN107038717A (zh) * | 2017-04-14 | 2017-08-11 | 东南大学 | 一种基于立体栅格自动分析3d点云配准误差的方法 |
CN107123161A (zh) * | 2017-06-14 | 2017-09-01 | 西南交通大学 | 一种基于narf和fpfh的接触网零全网三维重建方法 |
CN108133458A (zh) * | 2018-01-17 | 2018-06-08 | 视缘(上海)智能科技有限公司 | 一种基于目标物体空间点云特征的自动拼接方法 |
CN111009005A (zh) * | 2019-11-27 | 2020-04-14 | 天津大学 | 几何信息与光度信息相结合的场景分类点云粗配准方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115205354A (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107705328B (zh) | 用于3d对准算法的平衡探针位置选择 | |
Lenssen et al. | Deep iterative surface normal estimation | |
CN104376552B (zh) | 一种3d模型与二维图像的虚实配准方法 | |
Digne et al. | Scale space meshing of raw data point sets | |
Pan et al. | Iterative global similarity points: A robust coarse-to-fine integration solution for pairwise 3d point cloud registration | |
CN111340862B (zh) | 一种基于多特征融合的点云配准方法、装置及存储介质 | |
CN113628263A (zh) | 一种基于局部曲率及其近邻特征的点云配准方法 | |
CN114359042A (zh) | 点云拼接方法和装置、三维扫描仪及电子设备 | |
Kroemer et al. | Point cloud completion using extrusions | |
Sveier et al. | Object detection in point clouds using conformal geometric algebra | |
CN116452644A (zh) | 基于特征描述子的三维点云配准方法、装置及存储介质 | |
Wang et al. | Density-invariant registration of multiple scans for aircraft measurement | |
CN113012063A (zh) | 一种动态点云修复方法、装置及计算机设备 | |
CN115205354B (zh) | 基于ransac和icp点云配准的相控阵激光雷达成像方法 | |
Fua et al. | Reconstructing surfaces from unstructured 3d points | |
CN114332125A (zh) | 点云重建方法、装置、电子设备和存储介质 | |
CN113838069B (zh) | 基于平面度约束的点云分割方法和*** | |
Ventura et al. | P1ac: Revisiting absolute pose from a single affine correspondence | |
CN110969650B (zh) | 一种基于中心投影的强度图像与纹理序列间的配准方法 | |
Sandström et al. | Learning online multi-sensor depth fusion | |
Liu et al. | A dense visual SLAM method in dynamic scenes | |
CN113111741A (zh) | 一种基于三维特征点的装配状态识别方法 | |
Peralta et al. | A Preliminary Study on Low Overlapping Unorganized Point Clouds Registration Using Hough Voting | |
Lian et al. | 3D-SIFT Point Cloud Registration Method Integrating Curvature Information | |
Avidar et al. | Point cloud registration using a viewpoint dictionary |
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 |