CN115375713A - 地面点云的分割方法、装置及计算机可读存储介质 - Google Patents

地面点云的分割方法、装置及计算机可读存储介质 Download PDF

Info

Publication number
CN115375713A
CN115375713A CN202211311144.XA CN202211311144A CN115375713A CN 115375713 A CN115375713 A CN 115375713A CN 202211311144 A CN202211311144 A CN 202211311144A CN 115375713 A CN115375713 A CN 115375713A
Authority
CN
China
Prior art keywords
subspace
area
ground
region
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
Application number
CN202211311144.XA
Other languages
English (en)
Other versions
CN115375713B (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.)
Hangzhou Huacheng Software Technology Co Ltd
Original Assignee
Hangzhou Huacheng Software Technology Co Ltd
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 Hangzhou Huacheng Software Technology Co Ltd filed Critical Hangzhou Huacheng Software Technology Co Ltd
Priority to CN202211311144.XA priority Critical patent/CN115375713B/zh
Publication of CN115375713A publication Critical patent/CN115375713A/zh
Application granted granted Critical
Publication of CN115375713B publication Critical patent/CN115375713B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20104Interactive definition of region of interest [ROI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Processing (AREA)

Abstract

本申请公开了一种地面点云的分割方法、装置及计算机可读存储介质,该分割方法包括:获取当前帧的目标点云;对目标点云所处的空间区域进行切割,得到多个第一子空间区域;分别根据每个第一子空间区域包括的子点云,确定每个第一子空间区域对应的预估地面;分别根据每个第一子空间区域相邻的至少一个第一子空间区域对应的预估地面,对每个第一子空间区域对应的预估地面进行调整,得到每个第一子空间区域对应的目标地面;分别根据每个第一子空间区域对应的目标地面,从每个第一子空间区域包括的子点云中分割出地面点云。本申请所提供的分割方法能够提高地面点云的分割精度。

Description

地面点云的分割方法、装置及计算机可读存储介质
技术领域
本申请属于点云分割技术领域,特别是涉及一种地面点云的分割方法、装置及计算机可读存储介质。
背景技术
移动机器人平台如无人驾驶汽车、物流机器人、服务机器人以及家庭扫拖机器人等,主要通过多维传感系统(例如激光雷达、深度相机、摄像头等)来获取周边的环境数据。然后障碍物感知算法对这些环境数据进行综合分析计算后,经决策控制模块发出指令使移动机器人平台实现全自主运行。由此可见,障碍物感知是实现移动机器人平台智能化的前提条件之一。
但是目前在利用激光雷达、深度相机等传感设备采集的点云检测障碍物时,很容易将地面误检为障碍物,所以点云地面分割算法对障碍物感知系统来说至关重要。
发明内容
本申请提供一种地面点云的分割方法、装置及计算机可读存储介质,能够提高地面点云的分割精度。
本申请实施例第一方面提供一种地面点云的分割方法,所述方法包括:获取当前帧的目标点云;对所述目标点云所处的空间区域进行切割,得到多个第一子空间区域;获取前一历史帧对应的多个第二子空间区域以及每个所述第二子空间区域对应的历史地面;在多个所述第二子空间区域中,分别确定每个所述第一子空间区域匹配的所述第二子空间区域;分别根据每个所述第一子空间区域匹配的所述第二子空间区域对应的所述历史地面,对每个所述第一子空间区域包括的子点云进行筛选,得到每个所述第一子空间区域对应的候选点云;分别根据每个所述第一子空间区域对应的所述候选点云,确定每个所述第一子空间区域对应的预估地面;分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述预估地面,对每个所述第一子空间区域对应的所述预估地面进行调整,得到每个所述第一子空间区域对应的目标地面;分别根据每个所述第一子空间区域对应的所述目标地面,从每个所述第一子空间区域包括的所述子点云中分割出地面点云。
本申请实施例第二方面提供一种分割装置,所述分割装置包括处理器、存储器以及通信电路,所述处理器分别耦接所述存储器、所述通信电路,所述存储器中存储有程序数据,所述处理器通过执行所述存储器内的所述程序数据以实现上述方法中的步骤。
本申请实施例第三方面提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序能够被处理器执行以实现上述方法中的步骤。
有益效果是:本申请的分割方法首先将目标点云所处的空间区域进行分割,得到多个第一子空间区域,然后根据各个第一子空间区域所包括的子点云,初步得到每个第一子空间区域的预估地面,然后考虑到各个第一子空间区域之间的相邻特性,利用每个第一子空间区域相邻的至少一个第一子空间区域对应的预估地面,对每个第一子空间区域对应的预估地面进行调整,得到每个第一子空间区域对应的目标地面,使得每个第一子空间区域对应的目标地面更加贴近真实地面平滑光顺的特性,而后将每个第一子空间区域对应的目标地面分别作为每个第一子空间区域各自的真实地面,根据每个第一子空间区域各自的真实地面,从每个第一子空间区域所包括的子点云中分割出地面点云。其中由于是根据第一子空间区域对应的目标地面,从第一子空间区域所包括的子点云中分割出地面点云,而目标地面更加贴近真实地面平滑光顺的特性,因此可以提高最终分割出的地面点云的精度。同时本申请在确定每个第一子空间区域对应的预估地面时,还利用第一子空间区域匹配的第二子空间区域的历史地面作为先验地面,可以大大减少计算量,提高算法速度。
附图说明
为了更清楚地说明本申请实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1是本申请地面点云的分割方法一实施方式的流程示意图;
图2是本申请地面点云的分割方法另一实施方式的部分流程示意图;
图3是本申请地面点云的分割方法一实施方式中步骤S160的流程示意图;
图4是本申请多个第一子空间区域在真实地面上的投影示意图;
图5是本申请地面点云的分割方法一实施方式中步骤S130的流程示意图;
图6是本申请地面点云的分割方法一实施方式中步骤S140的流程示意图;
图7是本申请地面点云的分割方法一实施方式中步骤S150的流程示意图;
图8是本申请分割系统一实施方式的结构示意图;
图9是本申请分割装置一实施方式的结构示意图;
图10是本申请分割装置另一实施方式的结构示意图;
图11是本申请计算机可读存储介质一实施方式的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本申请的一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。
首先需要说明的是,本申请障碍点云的分割方法由分割装置执行,该分割装置可以是手机、电脑或者机器人等任一项具有算法处理能力的装置,在此不做限制。
参阅图1,在本申请一实施方式中,地面点云的分割方法包括:
S110:获取当前帧的目标点云。
目标传感器对所处周围环境进行数据采集,得到连续的多帧点云数据,其中,当前帧的目标点云为当前时刻待处理帧对应的点云数据。
其中,目标传感器包括但不限于深度相机、单/多线激光雷达、TOF相机(飞行时间相机)以及单/双目结构光设备。
其中,当前帧的目标点云包括多个点,每个点对应有一个三维坐标值。
S120:对目标点云所处的空间区域进行切割,得到多个第一子空间区域。
S130:分别根据每个第一子空间区域包括的子点云,确定每个第一子空间区域对应的预估地面。
具体地,分别对每个第一子空间区域包括的子点云进行处理,确定每个第一子空间区域中的预估地面,该预估地面是对第一子空间区域所包括真实地面的一个预估平面。
S140:分别根据每个第一子空间区域相邻的至少一个第一子空间区域对应的预估地面,对每个第一子空间区域对应的预估地面进行调整,得到每个第一子空间区域对应的目标地面。
具体地,为了使各个第一子空间区域对应的预估地面更加贴近真实地面平滑光顺的特性,本实施方式考虑到各个第一子空间区域之间的相邻特性,对于任意第一子空间区域A而言,根据其周围相邻的至少一个第一子空间区域B对应的预估地面,对第一子空间区域A对应的预估平面进行调整,最终得到第一子空间区域A的目标地面。
需要说明的是,不管是预估地面,还是目标地面,其都是用平面方程进行表示。
其中,对于任意第一子空间区域A而言,可以是根据其相邻的部分第一子空间区域B对应的预估平面对第一子空间区域A的预估地面进行调整,也可以是根据其相邻的全部第一子空间区域B对应的预估平面对第一子空间区域A的预估地面进行调整。
S150:分别根据每个第一子空间区域对应的目标地面,从每个第一子空间区域包括的子点云中分割出地面点云。
具体地,在得到每个第一子空间区域对应的目标地面后,分别将每个第一子空间区域对应的目标地面,确定为每个第一子空间区域各自所包含的真实地面,进而根据每个第一子空间区域各自的真实地面,从每个第一子空间区域所包括的子点云中分割出地面点云。
可以理解的是,所有第一子空间区域对应的地面点云构成了目标点云中的所有地面点云。
同时在分割出各个第一子空间区域中的地面点云后,可以进一步分割出各个第一子空间区域中的障碍物点云,例如,分别将各个子点云中不在地面点云的点所构成的点云确定为障碍物点云,从而可以用于后续的障碍物提取。
从上述内容可以看出,上述方法获得的目标地面更加贴近真实地面平滑光顺的特性,因此根据第一子空间区域对应的目标地面,从第一子空间区域所包括的子点云中分割出地面点云,可以提高最终分割出的出地面点云的精度。
其中,目标传感器直接输出的点云数据是传感器坐标系的点云数据,即目标传感器输出的每个点对应的坐标是该点在传感器坐标系下的坐标,在本实施方式中,在步骤S110中,当目标点云中的点对应的坐标是该点在传感器坐标系下的坐标时,将点的坐标从传感器坐标系转换到世界坐标系。具体可以通过传感器坐标系到世界坐标系的转换矩阵将点的坐标从传感器坐标系转换到世界坐标系。
需要说明的是,本申请并非限定步骤S110一定要将点的坐标从传感器坐标系转换到世界坐标系,在其他实施方式中,步骤S110中目标点云中点的坐标也可以是点在传感器坐标系下的坐标,具体可参见下文。
可以理解的是,在世界坐标系的z轴与地面垂直,即在理想状态下,真实地面对应的平面方程为z=0。
在本实施方式中,为了减少数据处理量,参阅图2,在步骤S120对目标点云所处的空间区域进行分割之前,还包括:
S160:分别判断目标点云中的各个点是否处于感兴趣空间区域。
具体地,感兴趣空间区域也称为ROI(region of interest)空间区域,在后续处理过程中,只需要对感兴趣空间区域中的点云进行处理。
S170:响应于点不处于感兴趣空间区域,将点从目标点云中去除。
在经过步骤S160后,目标点云中的点均处于感兴趣空间区域。
S180:对目标点云中的点进行滤波处理。
具体地,为了进一步减少计算量,在只保留感兴趣空间区域中的点后,还会对剩余的点进行滤波处理,从而以滤波处理后剩余的点进行后续步骤。
其中,滤波处理可以是任何一种类型的滤波处理,例如滤波处理为体素下采样。
需要说明的是,在其他实施方式中,在步骤S170将不处于感兴趣空间区域中的点从目标点云中去除后,可以不执行步骤S180。
或者在其他实施方式中,可以不执行步骤S160和步骤S170,而直接对目标点云进行滤波处理。
在一应用场景中,如图3所示,分别对目标点云中的每个点都执行如下步骤,以判断其是否处于感兴趣空间区域:
S161:判断点对应的坐标中的第一分量、第二分量、第三分量是否分别处于点对应的第一范围、第二范围、第三范围。
若判断结果为是,则执行步骤S162,否则执行步骤S163。
S162:确定点处于感兴趣空间区域。
S163:确定点不处于感兴趣空间区域。
具体地,第三分量与真实地面垂直的方向一致,也就是说,当点的坐标是该点在世界坐标系下的坐标时,第三分量是点坐标中的z坐标值,相应地,第一分量是点坐标中的x坐标值,第二分量是点坐标中的y坐标值。
其中考虑理想状态下,地面对应的平面方程为z=0,因此为了不遗漏第一子空间区域中的地面点云,第三范围的最小值小于零,最大值大于零。
同时在该应用场景中,考虑到目标传感器的标定误差、测距特性(随着点与目标传感器的距离增大,点在z轴上的噪声增大),设置点与采集目标点云的目标传感器的距离越大,点对应的第三范围越宽。
例如,如果点满足以下条件,则确定该点处于感兴趣区域:
xi∈[xmin,xmax],yi∈[ymin,ymax],zi∈[zmin×max(1.0,rα),zmax×max(1.0,rα)],其中,xi表示点的x坐标值,yi表示点的y坐标值,zi表示点的z坐标值,[xmin,xmax]为上述的第一范围,[ymin,ymax]为上述的第二范围,[zmin×max(1.0,rα),zmax×max(1.0,rα)],为上述的第三范围,其中,xmin、xmax、ymin、ymax、zmin以及zmax均为预设阈值,可以根据实际需求设置,其中,zmin小于零,zmax大于零,同时r为点到目标传感器中心的距离值,α与目标传感器的采集精度相关,当目标传感器为深度相机时,α通常取为1.08,通常而言,α的取值在1与2之间,且目标传感器的采集精度越高,α越接近1。
也就是说,不同点对应的第一范围相同,对应的第二范围相同,但是对应的第三范围不同。
可以理解的是,当步骤S110中点的坐标是该点在传感器坐标系下的坐标,且传感器坐标系的y轴垂直于真实地面时,第三分量是点坐标中的y坐标值,相应地,第一分量是点坐标中的x坐标值,第二分量是点坐标中的z坐标值,相应地,此时如果点满足以下条件,则确定该点处于感兴趣区域:
xi∈[xmin,xmax],yi∈[ymin×max(1.0,rα),ymax×max(1.0,rα)],zi∈[zmin,zmax]。
在经过上述步骤后,完成了对目标点云的数据预处理,接下来介绍步骤S120对目标点云所处的空间区域进行切割的过程:
在经过上述步骤后,目标点云所处的空间区域为感兴趣空间区域,且感兴趣空间区域在真实地面上的投影为矩形。
本实施方式在切割时,先沿着世界坐标系的y轴方向将目标点云所处的空间区域平均划分为m份,然后再沿着x轴方向将目标点云所处的空间区域平均划分为n份,从而得到m×n个第一子空间区域,此时m×n个第一子空间区域(图4中用标号101表示)在真实地面上的正投影如图4所示。在一应用场景中,m≤n。也就是说,此时目标点云所处的空间区域在真实地面上的垂直投影为第一矩形,多个第一子空间区域在真实地面上的垂直投影均为第二矩形,多个第二矩形呈阵列排布且多个第二矩形构成第一矩形。
需要说明的是,在其他实施方式中,也可以采用其他的切割方式对目标点云所处的空间区域进行切割,例如,先沿着与真实地面垂直但与世界坐标系的y轴方向相交的第一切割方向将目标点云所处的空间区域划分若干份,然后再沿着与真实地面垂直但与世界坐标系的x轴方向相交的第二切割方向将目标点云所处的空间区域划分若干份。此时多个第一子空间区域在真实地面上的正投影中,有的是平行四边形,有的可能是三角形,有的还可能是其他异形,需要说明的是,具有这些形状正投影的多个第一子空间区域所能实现的功能与具有矩形正投影的多个第一子空间区域所能实现的功能等同。
总而言之,本申请对具体的切割方式不做限制。
参阅图5,在本实施方式中,步骤S130分别对每个第一子空间区域都执行如下步骤,以确定第一子空间区域对应的预估地面的过程:
S131:判断当前帧是否是初始帧。
若判断结果为是,则进入步骤S132,若判断结果为否,则进入步骤S133。
具体地,初始帧指的是目标传感器输出的连续多帧点云数据中的第一帧。
S132:分别对每个第一子空间区域包括的子点云进行平面拟合,得到每个第一子空间区域对应的预估地面。
具体地,可以采用任意的平面拟合算法包括但不限于随机采样一致性(RandomSample Consensus,RANSAC)、奇异值分解(Singular Value Decomposition,SVD)以及布面模拟地面滤波(Cloth Simulation Filtering,CSF)等算法分别对每个第一子空间区域包括的子点云进行平面拟合,从而得到每个第一子空间区域对应的预估地面。
S133:将每个第一子空间区域及各自对应的预估地面分别保存为对于下一帧而言的第二子空间区域以及第二子空间区域对应的历史地面。
具体地,在得到每个第一子空间区域对应的预估地面后,将第一子空间区域、第一子空间区域对应的预估地面对应保存为对于下一帧而言的第二子空间区域以及第二子空间区域对应的历史地面。
S134:获取保存的前一历史帧对应的多个第二子空间区域以及多个历史地面。
具体地,如果当前帧不是初始帧,则获取上一帧得到的各个子空间区域以及各个子空间区域对应的预估地面,也就是,保存的前一历史帧对应的多个第二子空间区域以及多个历史地面。
S135:在多个第二子空间区域中,分别确定各个第一子空间区域匹配的第二子空间区域。
具体地,对于当前帧的每个第一子空间区域而言,在前一历史帧的多个子空间区域中,查找最匹配的子空间区域,也就是,在多个第二子空间区域中,确定每个第一子空间区域匹配的第二子空间区域。
其中,多个第二子空间区域中与第一子空间区域匹配的第二子空间区域指的是,与第一子空间区域是同一空间区域的可能性最大的第二子空间区域。
在一应用场景中,在多个第二子空间区域中,确定第一子空间区域匹配的第二子空间区域的步骤,包括:分别确定第一子空间区域与每个第二子空间区域的重叠度;将最大重叠度对应的第二子空间区域,确定为与第一子空间区域匹配。
具体地,第一子空间区域与第二子空间区域的重叠度越大,表明第一子空间区域与第二子空间区域是同一子空间区域的可能性越大。
其中,可以将第一子空间区域与第二子空间区域的交并比确定为两者的重叠度,也可以将第一子空间区域与第二子空间区域重叠部分的体积作为两者的重叠度。具体确定第一子空间区域与第二子空间区域的重叠度的过程不做具体限制。
同时确定第一子空间区域与第二子空间区域的交并比的过程可以是,确定第二子空间区域在真实地面上的正投影与第一子空间区域在真实地面上的正投影的交并比。
其中,考虑到存在最大重叠度有可能为零的情况,本应用场景在最大重叠度不为零时,将最大重叠度对应的第二子空间区域,确定为与第一子空间区域匹配,但是当最大重叠度为零时,将距离第一子空间区域最近的第二子空间区域,确定为第一子空间区域对应的第二子空间区域。
其中,可以将第一子空间区域中心到第二子空间区域中心的距离,确定为第一子空间区域与第二子空间区域的距离,也可以是将第一子空间区域、第二子空间区域相对设置的两表面的垂直距离,确定为第一子空间区域与第二子空间区域的距离。
S136:分别根据每个第一子空间区域包括的子点云、每个第一子空间区域对应的第二子空间区域对应的历史地面,确定每个第一子空间区域对应的预估地面。
具体地,在找到每个第一子空间区域对应的第二子空间后,分别对每个第一子空间区域执行如下步骤:
S1361:根据第一子空间区域包括的子点云、该第一子空间区域对应的第二子空间区域对应的历史地面,确定该第一子空间区域的预估地面。
在一应用场景中,S1361具体包括:分别从第一子空间区域包括的子点云中,筛选出到第一子空间区域对应的历史地面的距离不超过第二阈值的点,得到第一子空间区域对应的候选点云;对候选点云进行平面拟合,得到第一子空间区域对应的预估地面。
具体地,利用第一子空间区域对应的第二子空间区域对应的历史地面,对第一子空间区域所包括的子点云进行筛选,得到候选点云,然后对候选点云进行平面拟合,得到第一子空间区域对应的预估地面。
其中,筛选的标准是:如果点到历史地面的距离小于第二阈值,则保留,否则删除。其中,第二阈值可以根据实际需求设置,在此不做限制。例如,根据对地面点云分割精度的要求,设置第二阈值,该过程可以是:对地面点云分割精度的要求越高,设置第二阈值越小,对地面点云分割精度的要求越低,设置第二阈值越大。
此时以第二子空间区域对应的历史地面作为第一子空间区域的先验地面,对第一子空间区域所包括的子点云进行筛选,可以大大减少计算量,提高算法速度。
S137:将多个第一子空间区域、多个第一子空间区域各自对应的预估地面分别作为第二子空间区域、历史地面(相对于下一帧)进行保存。
具体地,在得到当前帧对应的多个第一子空间区域、多个第一子空间区域各自对应的预估地面后进行保存,以供下一帧进行使用。
参阅图6,在本实施方式中,步骤S140包括:
S141:分别根据每个第一子空间区域对应的预估地面,确定每个第一子空间区域对应的第一法向量,其中,第一子空间区域对应的第一法向量与第一子空间区域对应的预估地面垂直。
具体地,在得到第一子空间区域对应的预估平面后,根据预估平面的平面方程可以确定与预估平面垂直的第一法向量,从而每个第一子空间区域都有一个第一法向量。
S142:分别根据每个第一子空间区域相邻的至少一个第一子空间区域对应的第一法向量,对每个第一子空间区域对应的第一法向量进行调整,得到每个第一子空间区域对应的第二法向量。
具体地,对于第一子空间区域A而言,根据其相邻的至少一个第一子空间区域B对应的第一法向量,对第一子空间区域A的第一法向量进行调整,从而得到每个第一子空间区域A对应的第二法向量。
在一应用场景中,步骤S142具体包括:分别将每个第一子空间区域对应的第一法向量与第一子空间区域各自相邻的至少一个第一子空间区域对应的第一法向量按照第一法向量各自对应的权重进行加权平均计算,得到第一子空间区域对应的第二法向量;其中,第一法向量对应的第一子空间区域与采集目标点云的目标传感器的距离越大,第一法向量对应的权重越小。
具体地,假设第一子空间区域A对应的第一法向量为(xA,yA,zA),且该第一法向量对应的权重为λA,与第一子空间区域A相邻的至少一个第一子空间区域B对应的第一法向量分别为(xB1,yB1,zB1)、(xB2,yB2,zB2)……以及(xBn,yBn,zBn),这些第一法向量对应的权重分别为λB1、λB2……λBn,其中,n为至少一个第一子空间区域B的数量,其中在对(xA,yA,zA)进行调整后,第一子空间区域A对应的第二法向量为(x’,y’,z’),其中,x’、y’、z’的计算公式如下:
x’=(λA×xAB1×xB1B2×xB2+……+λBn×xBn)/(n+1);
y’=(λA×yAB1×yB1B2×yB2+……+λBn×yBn)/(n+1);
z’=(λA×zAB1×zB1B2×zB2+……+λBn×zBn)/(n+1)。
其中,考虑到目标传感器在采集点云数据时,点与目标传感器的距离越大,目标传感器采集到的点的信息的准确率越低,因此为了提高算法的准确率,设置第一子空间区域与采集目标点云的目标传感器的距离越大,第一子空间区域对应的第一法向量对应的权重越小。
例如可以按照如下公式确定第一法向量对应的权重:
β=kRcenter ,其中,β为第一法向量对应的权重,k为预设的比例系数,可根据实际需求设定,Rcenter为第一法向量对应的第一子空间区域的中心到目标传感器中心的距离值,α与目标传感器的采集精度相关,其中关于α已经在上文进行了具体介绍,在此不做赘述。
需要说明的是,本申请对确定第一法向量对应的权重的具体过程不做限制,具体可以根据实际需求进行设定,只要保证第一子空间区域与采集目标点云的目标传感器的距离越大,第一子空间区域对应的第一法向量对应的权重越小即可。
在另一应用场景中,也可以按照如下公式确定第一子空间区域A对应的第二法向量为(x’,y’,z’),其中,x’、y’、z’的计算公式如下:
x’=(xA+xB1+xB2+……+xBn)/(n+1);
y’=(yA+yB1+yB2+……+yBn)/(n+1);
z’=(zA+zB1+zB2+……+zBn)/(n+1)。
在本实施方式中,上述步骤中与第一子空间区域A相邻的至少一个第一子空间区域B为第一子空间区域A的八邻域内所有的第一子空间区域。
例如在图4中,对于第一子空间区域A而言,其相邻的至少一个第一子空间区域包括:第一子空间区域B1、第一子空间区域B2、第一子空间区域B3、第一子空间区域B4、第一子空间区域B5、第一子空间区域B26、第一子空间区域B7以及第一子空间区域B8。
需要说明的是,在其他实施方式中,与第一子空间区域A相邻的至少一个第一子空间区域B可以是第一子空间区域A的四邻域或者对角邻域内所有的第一子空间区域。
S143:分别根据每个第一子空间区域对应的第二法向量,确定每个第一子空间区域对应的目标地面,其中,每个第一子空间区域对应的目标地面与第一子空间区域对应的第二法向量垂直。
具体地,对于任意第一子空间区域而言,在得到其对应的第二法向量后,将与第二法向量垂直的平面作为第一子空间区域对应的目标地面。
在上述实施方式中,通过分别对每个第一子空间区域对应的第一法向量进行调整,进而实现直接对第一子空间区域对应的预估地面进行调整,但是本申请并不限制于此,在其他实施方式中,也可以分别根据每个第一子空间区域相邻的至少一个第一子空间区域对应的预估地面,直接对每个第一子空间区域对应的预估地面进行调整,例如,对于任意第一子空间区域A而言,在其相邻的至少第一个第一子空间区域B对应的预估地面C和该第一子空间区域A对应的预估地面C中,确定任意两个预估地面C之间的垂直距离,然后分别将每个预估地面C对应的垂直距离进行相加,而得到每个预估地面C对应的和值,然后将最小和值对应的预估地面C作为第一子空间区域A对应的目标地面。
参阅图7,在本实施方式中,步骤S150包括:
S151:分别确定每个第一子空间区域包括的子点云中的每个点到第一子空间区域各自对应的目标地面的距离。
S152:响应于点对应的距离小于或者等于点所在的第一子空间区域对应的第一阈值,确定点是第一子空间区域中的地面点。
具体地,对于任意第一子空间区域A而言,分别确定其所包括的子点云中的各个点到第一子空间区域A对应的目标地面的距离,如果点对应的距离小于或者等于该第一子空间区域A对应的第一阈值,则确定该点是第一子空间区域A中的地面点。
在一应用场景中,第一子空间区域不同,对应的第一阈值相同,也就是说,第一阈值始终是一个固定值。
在另一应用场景中,考虑到目标传感器的特性,即点距离目标传感器的距离越大,目标传感器输出的该点的信息的精度越低,因此设置第一子空间区域与采集目标点云的目标传感器的距离越大,第一子空间区域对应的第一阈值越大。
例如,按照如下公式确定第一子空间区域对应的第一阈值;
h=kh×max(1.0,rcenter α),其中,h为第一子空间区域对应的第一阈值,kh为预设的分割阈值参数,可以根据实际需求设定,rcenter为第一子空间区域的中心到目标传感器中心的距离值,α与目标传感器的采集精度相关。
需要说明的是,本申请对确定第一子空间区域对应的第一阈值的具体过程不做限制,只要保证第一子空间区域与采集目标点云的目标传感器的距离越大,第一子空间区域对应的第一阈值越大即可。
下面结合图8,对本实施方式一具体实施例中地面点云的分割方法进行详细的说明:
目标传感器10用于对周围环境进行采集,输出多帧连续的点云数据,而分割装置20在接收到目标传感器10输出的多帧连续的点云数据后,对每帧点云数据进行处理,其中将分割装置20当前时刻需要处理的点云数据定义为当前帧的目标点云,分割装置20对目标点云的处理过程包括:
S1:将目标点云中每个点的坐标由传感器坐标系转换到世界坐标系,具体可以采用传感器坐标系与世界坐标系之间的转换矩阵进行转换。
S2:设置感兴趣(ROI)空间区域,并分别判断目标点云中的每个点是否处于感兴趣空间区域,如果点处于感兴趣空间区域,则在目标点云中保留该点,否则去除该点,该具体过程包括:如果点的坐标满足以下条件,则确定该点处于感兴趣空间区域,否则确定该点不处于感兴趣空间区域:
xi∈[xmin,xmax],yi∈[ymin,ymax],zi∈[zmin×max(1.0,rα),zmax×max(1.0,rα)],其中,xmin、xmax、ymin、ymax、zmin以及zmax均为预设阈值,可以根据实际需求设置,同时zmin小于零,zmax大于零,同时r为点到目标传感器中心的距离值,α与目标传感器的采集精度相关,当目标传感器为深度相机时,α通常取为1.08。
S3:考虑到后续算法的执行效率与硬件计算资源,将目标点云进行滤波处理,其中滤波处理的方式可以是体素下采样处理等方式,本申请对滤波处理的具体过程不做限制。
S4:将目标点云所处的空间区域进行切割,得到m×n个第一子空间区域,其中,m×n个第一子空间区域在真实地面上的投影均为第二矩形,且多个第二矩形呈阵列排布,同时多个第二矩形构成目标点云所处的空间区域在真实地面上为第一矩形的正投影。
S5:判断获取的当前点帧是否为第一帧,若是,执行步骤S6;反之,则执行步骤S7。
S6:分别对每个第一子空间区域包括的子点云进行平面拟合,得到每个第一子空间区域对应的预估地面。其中,此处可选用任意的平面拟合算法,包括但不限于随机采样一致性(Random SampleConsensus,RANSAC)、奇异值分解(Singular Value Decomposition,SVD)以及布面模拟地面滤波(Cloth Simulation Filtering,CSF)等。之后,执行步骤S8。
S7:利用当前帧的上一帧所得到的第二子空间区域对应的预估地面作为当前帧各个第一子空间区域的地面先验信息,利用该先验信息加速当前帧各个第一子空间区域的平面拟合,最后获取当前帧各个第一子空间区域对应的预估地面,该具体过程包括:
S701:将上一帧的各个第二子空间区域与当前帧各个第一子空间区域进行匹配,分别计算上一帧各个第二子空间区域与当前帧各个第一子空间区域的重叠度,该重叠度具体可以是交并比(Intersection Over Union,IOU)。
S702:对于当前帧的任意第一子空间区域A而言,在上一帧的各个第二子空间区域中查找与其交并比最大的第二子空间区域C,并将第二子空间区域C对应的预估地面作为第一子空间区域A的先验地面,但是如果与第一子空间区域A对应的最大交并比为0,则将上一帧中与第一子空间区域A最相邻的第二子空间区域(中心与第一子空间区域A的距离最近)对应的预估地面作为第一子空间区域A的先验地面。
S703:对于当前帧的任意第一子空间区域A而言,计算其所包括的子点云中的每个点到第一子空间区域A的先验地面的距离,并将距离满足要求的点作为第一子空间区域A的候选点。
S704:分别对当前帧中每个第一子空间区域对应的候选点云进行平面拟合,最终得到每个第一子空间区域对应的预估地面。然后执行步骤S8。
S8:为了使各个第一子空间区域对应的最终地面更贴近真实的地面特性(平滑光顺),此处考虑各个第一子空间区域之间的相邻特性,同时考虑到目标传感器10的特性,对各个第一子空间区域分配距离权重参数β,β的计算公式为β=krcenter ,其中,k为预设的比例系数,可根据实际需求设定,rcenter为第一子空间区域的中心到目标传感器中心的距离值,α与目标传感器的采集精度相关。
接着,对于任意第一子空间区域A而言,将与其预估地面垂直的第一法向量和其八邻域内所有第一子空间区域B对应的预估地面所对应的第一法向量按照各个第一子空间区域对应的权重进行加权求平均值,得到第一子空间区域A对应的第二法向量。然后根据第一子空间区域A对应的第二法向量,确定第一子空间区域对应的目标地面,其中,目标地面与第二法向量垂直。
S9:分别根据每个第一子空间区域对应的目标地面,从每个第一子空间区域包括的子点云中分割出地面点云,分割的标准是:如果第一子空间区域对应的子点云中的点到第一子空间区域对应的目标地面的距离小于或者等于第一子空间区域对应的第一阈值,则确定该点是第一子空间区域中的地面点。其中,第一子空间区域对应的第一阈值的计算公式为:h=kh×max(1.0,rcenter α),其中,h为第一子空间区域对应的第一阈值,kh为预设的分割阈值参数,可以根据实际需求设定,rcenter为第一子空间区域的中心到目标传感器中心的距离值,α与目标传感器的采集精度相关。
参阅图9,图9是本申请分割装置一实施方式的结构示意图。该分割装置200包括处理器210、存储器220以及通信电路230,处理器210分别耦接存储器220、通信电路230,存储器220中存储有程序数据,处理器210通过执行存储器220内的程序数据以实现上述任一项实施方式方法中的步骤,其中详细的步骤可参见上述实施方式,在此不再赘述。
其中,分割装置200可以是电脑、手机、机器人等任一项具有算法处理能力的装置,在此不做限制。
参阅图10,图10是本申请分割装置另一实施方式的结构示意图。该分割装置300包括获取模块310、第一分割模块320、预估模块330、调整模块340以及第二分割模块350。
获取模块310用于获取当前帧的目标点云。
第一分割模块320与获取模块310连接,用于对目标点云所处的空间区域进行切割,得到多个第一子空间区域。
预估模块330与第一分割模块320连接,用于获取前一历史帧对应的多个第二子空间区域以及每个第二子空间区域对应的历史地面;在多个第二子空间区域中,分别确定每个第一子空间区域匹配的第二子空间区域;分别根据每个第一子空间区域匹配的第二子空间区域对应的历史地面,对每个第一子空间区域包括的子点云进行筛选,得到每个第一子空间区域对应的候选点云;分别根据每个第一子空间区域对应的候选点云,确定每个第一子空间区域对应的预估地面。
调整模块340与预估模块330连接,用于分别根据每个第一子空间区域相邻的至少一个第一子空间区域对应的预估地面,对每个第一子空间区域对应的预估地面进行调整,得到每个第一子空间区域对应的目标地面。
第二分割模块350与调整模块340连接,用于分别根据每个第一子空间区域对应的目标地面,从每个第一子空间区域包括的子点云中分割出地面点云。
其中,分割装置300可以是电脑、手机、机器人等任一项具有算法处理能力的装置,在此不做限制。且该分割装置300在工作时执行上述任一项实施方式中地面点云的分割方法中的步骤,其中详细的步骤可参见上述实施方式,在此不再赘述。
参阅图11,图11是本申请计算机可读存储介质一实施方式的结构示意图。该计算机可读存储介质400存储有计算机程序410,计算机程序410能够被处理器执行以实现上述任一项方法中的步骤。
其中,计算机可读存储介质400具体可以为U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等可以存储计算机程序410的装置,或者也可以为存储有该计算机程序410的服务器,该服务器可将存储的计算机程序410发送给其他设备运行,或者也可以自运行该存储的计算机程序410。
以上所述仅为本申请的实施例,并非因此限制本申请的专利范围,凡是利用本申请说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本申请的专利保护范围内。

Claims (18)

1.一种地面点云的分割方法,其特征在于,所述方法包括:
获取当前帧的目标点云;
对所述目标点云所处的空间区域进行切割,得到多个第一子空间区域;
获取前一历史帧对应的多个第二子空间区域以及每个所述第二子空间区域对应的历史地面;
在多个所述第二子空间区域中,分别确定每个所述第一子空间区域匹配的所述第二子空间区域;
分别根据每个所述第一子空间区域匹配的所述第二子空间区域对应的所述历史地面,对每个所述第一子空间区域包括的子点云进行筛选,得到每个所述第一子空间区域对应的候选点云;
分别根据每个所述第一子空间区域对应的所述候选点云,确定每个所述第一子空间区域对应的预估地面;
分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述预估地面,对每个所述第一子空间区域对应的所述预估地面进行调整,得到每个所述第一子空间区域对应的目标地面;
分别根据每个所述第一子空间区域对应的所述目标地面,从每个所述第一子空间区域包括的所述子点云中分割出地面点云。
2.根据权利要求1所述的方法,其特征在于,所述分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述预估地面,对每个所述第一子空间区域对应的所述预估地面进行调整,得到每个所述第一子空间区域对应的目标地面的步骤,包括:
分别根据每个所述第一子空间区域对应的所述预估地面,确定每个所述第一子空间区域对应的第一法向量,其中,所述第一子空间区域对应的所述第一法向量与所述第一子空间区域对应的所述预估地面垂直;
分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述第一法向量,对每个所述第一子空间区域对应的所述第一法向量进行调整,得到每个所述第一子空间区域对应的第二法向量;
分别根据每个所述第一子空间区域对应的所述第二法向量,确定每个所述第一子空间区域对应的所述目标地面,其中,每个所述第一子空间区域对应的所述目标地面与所述第一子空间区域对应的所述第二法向量垂直。
3.根据权利要求2所述的方法,其特征在于,所述分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述第一法向量,对每个所述第一子空间区域对应的所述第一法向量进行调整,得到每个所述第一子空间区域对应的第二法向量的步骤,包括:
分别将每个所述第一子空间区域对应的所述第一法向量与所述第一子空间区域各自相邻的至少一个所述第一子空间区域对应的所述第一法向量按照所述第一法向量各自对应的权重进行加权平均计算,得到所述第一子空间区域对应的所述第二法向量;
其中,所述第一法向量对应的所述第一子空间区域与采集所述目标点云的目标传感器的距离越大,所述第一法向量对应的所述权重越小。
4.根据权利要求3所述的方法,其特征在于,所述方法还包括:
按照如下公式确定所述第一法向量对应的所述权重:
β=kRcenter ,其中,β为所述第一法向量对应的所述权重,k为预设的比例系数,Rcenter为所述第一法向量对应的所述第一子空间区域的中心到采集所述目标点云的目标传感器中心的距离值,α与所述目标传感器的采集精度相关。
5.根据权利要求2所述的方法,其特征在于,所述目标点云所处的空间区域在真实地面上的垂直投影为第一矩形,多个所述第一子空间区域在所述真实地面上的垂直投影均为第二矩形,多个所述第二矩形呈阵列排布且多个所述第二矩形构成所述第一矩形;
所述分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述第一法向量,对每个所述第一子空间区域对应的所述第一法向量进行调整,得到每个所述第一子空间区域对应的第二法向量的步骤,包括:
分别根据每个所述第一子空间区域的八邻域内所有的所述第一子空间区域对应的所述第一法向量,对每个所述第一子空间区域对应的所述第一法向量进行调整,得到每个所述第一子空间区域对应的所述第二法向量。
6.根据权利要求1所述的方法,其特征在于,所述分别根据每个所述第一子空间区域对应的所述目标地面,从每个所述第一子空间区域包括的所述子点云中分割出地面点云的步骤,包括:
分别确定每个所述第一子空间区域包括的所述子点云中的每个点到所述第一子空间区域各自对应的所述目标地面的距离;
响应于所述点对应的所述距离小于或者等于所述点所在的所述第一子空间区域对应的第一阈值,确定所述点是所述第一子空间区域中的地面点。
7.根据权利要求6所述的方法,其特征在于,所述第一子空间区域与采集所述目标点云的目标传感器的距离越大,所述第一子空间区域对应的第一阈值越大。
8.根据权利要求7所述的方法,其特征在于,所述方法还包括:
按照如下公式确定所述第一子空间区域对应的所述第一阈值;
h=kh×max(1.0,rcenter α),其中,h为所述第一子空间区域对应的所述第一阈值,kh为预设的分割阈值参数,rcenter为所述第一子空间区域的中心到采集所述目标点云的目标传感器中心的距离值,α与所述目标传感器的采集精度相关。
9.根据权利要求1所述的方法,其特征在于,在所述获取前一历史帧对应的多个第二子空间区域以及每个所述第二子空间区域对应的历史地面之前,还包括:
判断所述当前帧是否是初始帧;
响应于所述当前帧不是初始帧,执行所述获取前一历史帧对应的多个第二子空间区域以及每个所述第二子空间区域对应的历史地面的步骤;
否则,分别对每个所述第一子空间区域包括的所述子点云进行平面拟合,得到每个所述第一子空间区域对应的所述预估地面,然后执行所述分别根据每个所述第一子空间区域相邻的至少一个所述第一子空间区域对应的所述预估地面,对每个所述第一子空间区域对应的所述预估地面进行调整,得到每个所述第一子空间区域对应的目标地面的步骤。
10.根据权利要求1所述的方法,其特征在于,所述分别根据每个所述第一子空间区域匹配的所述第二子空间区域对应的所述历史地面,对每个所述第一子空间区域包括的子点云进行筛选,得到每个所述第一子空间区域对应的候选点云的步骤,包括:
分别从每个所述第一子空间区域包括的所述子点云中,筛选出到所述第一子空间区域各自对应的所述历史地面的距离不超过第二阈值的点,得到每个所述第一子空间区域对应的所述候选点云。
11.根据权利要求9所述的方法,其特征在于,所述在多个所述第二子空间区域中,分别确定每个所述第一子空间区域匹配的所述第二子空间区域的步骤,包括:
分别确定每个所述第二子空间区域与每个所述第一子空间区域的重叠度;
在每个所述第一子空间区域对应的所述重叠度中,确定最大重叠度;
分别将每个所述第一子空间区域对应的所述最大重叠度对应的所述第二子空间区域,确定为所述第一子空间区域各自匹配的所述第二子空间区域。
12.根据权利要求11所述的方法,其特征在于,所述方法还包括:
响应于所述第一子空间区域对应的所述最大重叠度不为零,将所述第一子空间区域对应的所述最大重叠度对应的所述第二子空间区域,确定为所述第一子空间区域匹配的所述第二子空间区域;
否则,将距离所述第一子空间区域最近的所述第二子空间区域,确定为所述第一子空间区域匹配的所述第二子空间区域。
13.根据权利要求1所述的方法,其特征在于,在所述对所述目标点云所处的空间区域进行切割,得到多个第一子空间区域之前,还包括:
分别判断所述目标点云中的各个点是否处于感兴趣空间区域;
响应于所述点不处于所述感兴趣空间区域,将所述点从所述目标点云中去除。
14.根据权利要求13所述的方法,其特征在于,所述方法还包括:
响应于所述点对应的坐标中的第一分量、第二分量、第三分量分别处于所述点对应的第一范围、第二范围、第三范围,确定所述点处于所述感兴趣空间区域;
其中,所述第三分量与垂直于真实地面的方向一致,且所述第三范围的最小值小于零,所述第三范围的最大值大于零,且所述点与采集所述目标点云的目标传感器的距离越大,所述点对应的所述第三范围越宽。
15.根据权利要求14所述的方法,其特征在于,所述方法还包括:
按照如下公式确定所述点对应的所述第三范围的最小值z1
z1=zmin×max(1.0,rα);
以及,按照如下公式确定所述点对应的所述第三范围的最大值z2
z2=zmax×max(1.0,rα);
其中,zmin为第一预设阈值、zmax为第二预设阈值,r为所述点到采集所述目标点云的目标传感器中心的距离值,α与所述目标传感器的采集精度相关。
16.根据权利要求13所述的方法,其特征在于,在所述响应于所述点不处于所述感兴趣空间,将所述点从所述目标点云中去除之后,还包括:
对所述目标点云中的点进行滤波处理。
17.一种分割装置,其特征在于,所述分割装置包括处理器、存储器以及通信电路,所述处理器分别耦接所述存储器、所述通信电路,所述存储器中存储有程序数据,所述处理器通过执行所述存储器内的所述程序数据以实现如权利要求1-16任一项所述方法中的步骤。
18.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序能够被处理器执行以实现如权利要求1-16任一项所述方法中的步骤。
CN202211311144.XA 2022-10-25 2022-10-25 地面点云的分割方法、装置及计算机可读存储介质 Active CN115375713B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211311144.XA CN115375713B (zh) 2022-10-25 2022-10-25 地面点云的分割方法、装置及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211311144.XA CN115375713B (zh) 2022-10-25 2022-10-25 地面点云的分割方法、装置及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN115375713A true CN115375713A (zh) 2022-11-22
CN115375713B CN115375713B (zh) 2023-03-24

Family

ID=84074212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211311144.XA Active CN115375713B (zh) 2022-10-25 2022-10-25 地面点云的分割方法、装置及计算机可读存储介质

Country Status (1)

Country Link
CN (1) CN115375713B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100209013A1 (en) * 2009-02-13 2010-08-19 Harris Corporation Registration of 3d point cloud data to 2d electro-optical image data
CN112558035A (zh) * 2019-09-24 2021-03-26 北京百度网讯科技有限公司 用于估计地面的方法和装置
CN112630798A (zh) * 2019-09-24 2021-04-09 北京百度网讯科技有限公司 用于估计地面的方法和装置
CN112767456A (zh) * 2021-01-18 2021-05-07 南京理工大学 一种三维激光点云快速重定位方法
US20210148709A1 (en) * 2019-11-14 2021-05-20 Nio Usa, Inc. Real-time ground surface segmentation algorithm for sparse point clouds
CN113835099A (zh) * 2021-02-01 2021-12-24 贵州京邦达供应链科技有限公司 点云地图更新方法及装置、存储介质、电子设备
CN113870435A (zh) * 2021-09-28 2021-12-31 浙江华睿科技股份有限公司 一种点云分割方法、装置、电子设备及存储介质
CN114200477A (zh) * 2021-12-13 2022-03-18 上海无线电设备研究所 一种激光三维成像雷达地面目标点云数据处理方法
CN114663526A (zh) * 2022-03-17 2022-06-24 深圳市优必选科技股份有限公司 障碍物检测方法、装置、机器人及计算机可读存储介质
CN114820657A (zh) * 2022-04-27 2022-07-29 北京主线科技有限公司 地面点云分割方法、系统、地面建模方法及介质
CN114981840A (zh) * 2020-11-11 2022-08-30 深圳元戎启行科技有限公司 基于点云数据的地面分割方法、装置和计算机设备
CN115170772A (zh) * 2022-05-18 2022-10-11 上海大学 一种基于点云地图交互式可通过性分析的三维路径规划方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100209013A1 (en) * 2009-02-13 2010-08-19 Harris Corporation Registration of 3d point cloud data to 2d electro-optical image data
CN112558035A (zh) * 2019-09-24 2021-03-26 北京百度网讯科技有限公司 用于估计地面的方法和装置
CN112630798A (zh) * 2019-09-24 2021-04-09 北京百度网讯科技有限公司 用于估计地面的方法和装置
US20210148709A1 (en) * 2019-11-14 2021-05-20 Nio Usa, Inc. Real-time ground surface segmentation algorithm for sparse point clouds
CN114981840A (zh) * 2020-11-11 2022-08-30 深圳元戎启行科技有限公司 基于点云数据的地面分割方法、装置和计算机设备
CN112767456A (zh) * 2021-01-18 2021-05-07 南京理工大学 一种三维激光点云快速重定位方法
CN113835099A (zh) * 2021-02-01 2021-12-24 贵州京邦达供应链科技有限公司 点云地图更新方法及装置、存储介质、电子设备
CN113870435A (zh) * 2021-09-28 2021-12-31 浙江华睿科技股份有限公司 一种点云分割方法、装置、电子设备及存储介质
CN114200477A (zh) * 2021-12-13 2022-03-18 上海无线电设备研究所 一种激光三维成像雷达地面目标点云数据处理方法
CN114663526A (zh) * 2022-03-17 2022-06-24 深圳市优必选科技股份有限公司 障碍物检测方法、装置、机器人及计算机可读存储介质
CN114820657A (zh) * 2022-04-27 2022-07-29 北京主线科技有限公司 地面点云分割方法、系统、地面建模方法及介质
CN115170772A (zh) * 2022-05-18 2022-10-11 上海大学 一种基于点云地图交互式可通过性分析的三维路径规划方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHIXIN LENG 等: "An Improved Fast Ground Segmentation Algorithm for 3D Point Cloud", 《2020 CHINESE CONTROL AND DECISION CONFERENCE》 *
张凯 等: "基于自适应阈值的三维激光点云地面分割算法研究", 《汽车工程》 *
游天 等: "一种面向室内导航的通行区域模型及其自动提取算法", 《武汉大学学报·信息科学版》 *
葛科迪 等: "面向复杂曲面加工的虚拟夹具辅助机器人示教编程研究", 《中国机械学报》 *

Also Published As

Publication number Publication date
CN115375713B (zh) 2023-03-24

Similar Documents

Publication Publication Date Title
CN108985171B (zh) 运动状态估计方法和运动状态估计装置
US8199977B2 (en) System and method for extraction of features from a 3-D point cloud
EP1329850B1 (en) Apparatus, program and method for detecting both stationary objects and moving objects in an image
Kang et al. Automatic targetless camera–lidar calibration by aligning edge with gaussian mixture model
CN109446886B (zh) 基于无人车的障碍物检测方法、装置、设备以及存储介质
US10762704B2 (en) Method for establishing a deformable 3D model of an element, and associated system
Weon et al. Object Recognition based interpolation with 3d lidar and vision for autonomous driving of an intelligent vehicle
CN111582054B (zh) 点云数据处理方法及装置、障碍物检测方法及装置
CN112651944B (zh) 基于cad模型的3c部件高精度六维位姿估计方法及系统
CN112097732A (zh) 一种基于双目相机的三维测距方法、系统、设备及可读存储介质
CN113015924A (zh) 基于来自一个或多个位置传感器的测量样本表征对象的装置和方法
CN113052907B (zh) 一种动态环境移动机器人的定位方法
CN111709988A (zh) 一种物体的特征信息的确定方法、装置、电子设备及存储介质
CN111311651A (zh) 点云配准方法和装置
Wallner et al. Real-time map refinement by fusing sonar and active stereo-vision
CN115496923A (zh) 一种基于不确定性感知的多模态融合目标检测方法及装置
CN114241448A (zh) 障碍物航向角的获取方法、装置、电子设备及车辆
US10223803B2 (en) Method for characterising a scene by computing 3D orientation
Hamzah et al. Stereo matching algorithm based on illumination control to improve the accuracy
CN115375713B (zh) 地面点云的分割方法、装置及计算机可读存储介质
CN116921932A (zh) 焊接轨迹识别方法、装置、设备及存储介质
Kokovkina et al. The algorithm of EKF-SLAM using laser scanning system and fisheye camera
CN116740160A (zh) 一种复杂交通场景中的毫秒级多平面实时提取方法及装置
CN116385997A (zh) 一种车载障碍物精确感知方法、系统及存储介质
CN113409446B (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