CN105574905B - 一种三维激光点云数据的二维图像化表述方法 - Google Patents
一种三维激光点云数据的二维图像化表述方法 Download PDFInfo
- Publication number
- CN105574905B CN105574905B CN201510938881.6A CN201510938881A CN105574905B CN 105574905 B CN105574905 B CN 105574905B CN 201510938881 A CN201510938881 A CN 201510938881A CN 105574905 B CN105574905 B CN 105574905B
- Authority
- CN
- China
- Prior art keywords
- point cloud
- laser
- gray
- point
- dimensional
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
一种三维激光点云数据的二维图像化表述方法,属于影像技术领域。主要由移动机器人、激光扫描器和计算机系统组成;激光扫描器和计算机系统均安装在移动机器人平台上,激光扫描器和计算机系统之间通过数据线进行通信;激光扫描器采集场景的三维点云数据信息,计算机系统将三维激光点云数据转化为二维灰度图像。本发明为基于激光点云数据的场景分割与场景理解提供了有力的支持,能用于大范围数据,适应多种扫描方式,提高了图像的景物边界区分度以及细节质量,纹理更加清晰。
Description
技术领域
本发明属于机器人自主环境感知技术及激光扫描成像领域,涉及到将移动机器人系统采集的三维激光点云数据转化为二维灰度图像以完成对室内外场景的表述。
背景技术
视觉是智能机器人与智能系统进行环境感知的重要手段之一。由于视觉传感器是基于被动感知机理,因此在成像过程中容易受到亮度及光照强度变化等环境因素的影响。此外单目视觉传感器采集到的图像信息缺少空间深度信息,不利于区分具有不同深度信息的景物。
三维激光数据的获取不受亮度及光照强度变化等环境因素的影响,可提供精准的三维空间测量信息。但利用三维点云数据对大范围环境进行表述时,观测视点的选取会直接影响环境感知效果;此外三维激光扫描点的分布是由实际场景特性所决定,三维点的分布无规律可循。考虑二维图像表述符合人类视觉感知习惯且具有像素排列格式规范等优点,因此将三维点云数据转化为景物边缘显著、纹理清晰且有较强场景层次的二维灰度图像,将有助于基于该图像模型的场景认知与理解。目前,各领域研究中使用三维激光点云数据生成的灰度图像模型主要有深度图和方位角图两种。
传统的深度图是通过把激光传感器或其他测距传感器点云数据的深度信息直接映射为灰度值,从而得到灰度图像,参考文献(Anguelov D,Dulong C,Filip D,etal.Google street view:Capturing the world at street level[J].Computer,2010(6):32-38.)。其优点是成像原理简单、计算耗时较少,能够区分出不同深度层次的景物。但其在图像细节和纹理的描述上很不理想,尤其是在描述大范围场景时,对距离的变化很不敏感。生成的图像往往会弱化景物的细节和近距离景物的边界信息,从而给之后的图像处理和理解造成困难。并且由于传统深度图景深方向取决于激光的扫描方向,使得深度图的效果会受到激光扫描方向的影响。
方位角图是另一种常用的由三维激光点云数据生成的灰度图像。此方法是DavideScaramuzza等学者首次提出的,参考文献(Scaramuzza D,Harati A,SiegwartR.Extrinsic self calibration of a camera and a 3d laser range finder fromnatural scenes[C]//Intelligent Robots and Systems,2007.IROS 2007.IEEE/RSJInternational Conference on.IEEE,2007:4164-4169)。方位角图的原理是将每个激光点对应的方位角映射为对应像素的灰度值。其中方位角定义为视点(激光传感器扫描时位置)到当前激光点位置的向量与当前激光点到其右上方邻近激光点的向量的夹角。参考文献(Zhuang Y,Lin X,Hu H,et al.Using Scale Coordination and Semantic Informationfor Robust 3-D Object Recognition by a Service Robot[J].Sensors Journal,IEEE,2015,15(1):37-47.)应用方位角图模型将定点扫描的三维点云数据转换为方位角图,并在此基础上对场景进行识别与理解。从以上文献中可知,方位角图中景物的细节和纹理相对清晰。但图像灰度值完全由方位角决定,不能清楚的区分不同深度层次的景物,并且图像质量受视点选择的影响较大。场景中扫描位置(视点)唯一,限制了实际应用中移动机器人平台的连续工作性能。计算的方位角相对视点产生渐变的趋势,使得整幅图像在横向具有灰度值渐变。
在实际应用中,大范围的场景一般都是采用激光连续扫描的方式采集点云数据,每一帧激光数据都有自己的扫描中心(局部视点),即使通过优化方法计算出相对较优的全局视点,仍无法有效消除基于视点的构图方法(如方位角图)的局限性。
发明内容
为解决上述问题,本发明提出了一种全新的灰度图像生成方法,完全不依赖于场景视点并且不使用角度计算,这也是本发明与方位角图的本质区别。以最优深度代替传统深度图的深度概念,消除了扫描方向对图像质量的影响。不同于传统的深度图、方位角图的单变量映射到灰度值的做法,本方法通过计算点云对应的最优深度值和方位向量长度,经加权后映射为灰度值,从而生成图像。基于本方法生成的二维灰度图像既具有清晰的细节和纹理特征,又利用深度分布信息强化了图像的空间层次感。本发明可以用于连续扫描的大范围数据,也可以用于定点扫描数据,并且对扫描方式和扫描位置没有特别的要求,有着良好的适应性。
本发明的具体技术方案如下:
一种三维激光点云数据的二维图像化表述方法,主要由移动机器人、激光扫描器和计算机系统组成;激光扫描器和计算机系统均安装在移动机器人平台上,激光扫描器和计算机系统之间通过数据线进行通信;激光扫描器采集图像的三维点云数据信息,计算机系统将三维激光点云数据转化为二维灰度图像;计算机系统软件要求采集的三维激光点云数据可以区分属于每一帧激光线的激光数据点,以便能够使三维激光点云数据以二维数组的形式存储,并且使每个激光点和灰度图的像素点一一对应;将获得的激光点云数据存入二维数组中,使得二维数组的每一列激光点属于同一帧激光线,并在二维数组中截取所需大小的场景,用截取场景的三维点云数据构建灰度图;
具体将三维激光点云数据转化为二维灰度图像计算方法如下:
(1)最优深度对应灰度值的计算
最优深度方向即为能使场景不同层次景物间区分度最大的深度方向;由于场景点云分布的主方向反映了激光点的统计学方向,而激光点的方向即为场景内景物的层次方向,所以垂直于场景点云分布主方向时,最能清楚的区分景物的层次,即最优深度方向垂直于场景点云分布的主方向,本发明将场景点云投影到水平面,运用统计学原理计算场景点云分布的主方向,由此获得场景的最优深度信息:
式(1)为点云在水平面投影的协方差矩阵,其中(x,y)为点云在水平面投影的坐标,cov代表协方差运算;通过计算协方差矩阵的特征向量V1、V2即可估计出场景的点云分布主方向和最优深度方向;设|V1|>|V2|,则V1为场景点云分布的主方向;由点云在水平面上的投影区域的主方向向量V1可求出过原点且平行于V1的直线l1的斜率K,进而确定直线l1的方程:
K*X+Y=0 (2)
由点到直线的距离公式:
其中D为位于i行j列的激光点Pi,j的投影(x,y)到直线l1的距离,这个距离便是点云中激光点Pi,j的最优深度Di,j;将最优深度Di,j映射至0-255的灰度范围,得到点云中激光点Pi,j对应于最优深度Di,j的灰度值
其中Dmin为点云最优深度值的最小值,Dmax为点云最优深度值的最大值;
(2)方位向量长度的计算
对于激光点云中位于第i行j列列的激光点Pi,j,取其上方相邻的激光点Pi-1,j,设点Pi,j到点Pi-1,j的向量为Vd;向量Vd的长度|Vd|即为点Pi,j的方位向量长度Li,j;方位向量长度能够反应场景中景物的纹理变化,对景物边缘的表现效果尤其明显,所以为了保证较好的效果,在映射为灰度值时不能采用单纯的线性映射;
(3)方位向量长度灰度数值均衡化
求出激光点数量为I*J的点云中每个点所对应的方位向量长度后,在映射为灰度值时,为了保证得到的灰度图具有高对比度、多变的灰度色调和明显的景物边界,对方位向量长度采取灰度数值均衡化;
对I*J个激光点对应的方位向量长度Li,j按数值大小升序排列,把0-255的灰度值范围视作256个灰度值区间,求出每个灰度值区间所对应的向量长度数据量N:
通过N与升序排列的方位向量长度Li,j,即可确定每个灰度值区间的端点长度值;从而得到点云中激光点Pi,j对应于方位向量长度Li,j的灰度值
其中L0,L1,L2…,L255为256个灰度值区间的端点长度值;
(4)权重动态调整
将点云中激光点Pi,j对应于最优深度Di,j的灰度值与对应于方位向量长度Li,j的灰度值加权平均得到激光点Pi,j在灰度图中对应的灰度值PGi,j:
其中PGi,j为灰度图中第i行j列像素的灰度值;a与b为权值,且a+b=10;此时所生成灰度图的灰度值矩阵的方差δ2为:
其中为灰度图中像素灰度均值;当方差最大时,则意味着灰度图具有高对比度和多变的灰度色调;取使方差最大的权值a与b生成灰度图,此时得到的灰度图是最优的;利用传统求最优解的方法求解式(8),得到使δ2为最大值的a与b,记做与
(5)生成灰度图
求出使δ2为最大值的权值与后,遍历点云数据中的I*J个激光点,求出每个激光点对应的灰度值:
PGi,j即为大小为I*J的灰度图中第i行j列像素的灰度值。
本发明灰度图的像素点与点云中的激光点一一对应,有效的提供了信息融合的方法,为基于激光点云数据的场景分割与场景理解提供了有力的支持。与传统的深度图与方位角图构图方法相比,本发明能用于大范围数据,适应多种扫描方式,提高了图像的景物边界区分度以及细节质量,纹理更加清晰。
附图说明
图1为移动机器人平台激光数据采集过程示意图
图2为最优深度计算方式示意图
图3为方位向量长度计算方式示意图
图4为一幅场景对应的三维激光点云数据
图5为图4对应的本发明的灰度图
图6为图4对应的深度图
图7为图4对应的方位角图
图8(a)为图5的超像素分割结果
图8(b)为图6的超像素分割结果
图8(c)为图7的超像素分割结果
具体实施方式
本发明的具体实施方式包括两个方面,一是三维激光点云数据的获取,二是根据三维激光点云数据生成灰度图。
本发明使用的三维激光点云数据可以是由地面移动机器人平台或空中移动机器人平台采集的,也可以是由其他任意采集方式获得的三维激光点云数据。本发明对场地的选取无过多的要求,室内室外环境均可。对场景的扫描方式选用定点扫描方式或连续激光扫描方式均可。要求采集的三维激光点云数据可以区分属于每一帧激光线的激光数据点,以便能够使三维激光点云数据以二维数组的形式存储,并且使每个激光点和灰度图的像素点一一对应。将获得的激光点云数据存入二维数组中,使得二维数组的每一列激光点属于同一帧激光线,并在二维数组中截取所需大小的场景,用截取场景的三维点云数据构建灰度图。
为了验证本方法的有效性,利用场景的三维激光点云数据(如图4所示)来进行灰度图的构建。场景点云的大小为400*300(数据形式为400行300列的二维数组),即每幅场景中有300帧激光线,每一帧激光线包含400个激光数据点。将点云数据中所有点的Z坐标置零,即把点云所有点投影到水平面上。计算得到点云投影的协方差矩阵:
求协方差矩阵的特征向量,得到点云投影的两个方向向量:
V1=(0.998163,0.060591)
V2=(-0.060591,0.998163)
根据投影区域的方向向量V1可求出直线l1的斜率K,进而确定直线l1的方程:
16.4738x+y=0
计算位于i行j列的激光点的投影(x,y)到直线l1的距离D:
这个距离便是点云中激光点Pi,j的最优深度Di,j。将最优深度Di,j映射至0-255的灰度范围,得到点云中激光点Pi,j对应于最优深度Di,j的灰度值
对于激光点云中位于第i行j列的激光点Pi,j,取其上方相邻的激光点Pi-1,j,设点Pi,j到点Pi-1,j的向量为Vd(如图3所示)。向量Vd的长度|Vd|即为点Pi,j的方位向量长度Li,j。对120000个激光点对应的方位向量长度Li,j按数值大小升序排列,把0-255的灰度值范围视作256个灰度值区间,求出每个灰度值区间所对应的向量长度数据量N:
通过N与升序排列的方位向量长度Li,j,确定每个灰度值区间的端点向量长度值。即可得到点云中激光点Pi,j对应于方位向量长度Li,j的灰度值
根据灰度图的灰度值矩阵方差公式:
其中
a+b=10
解得,当使δ2为最大值时
则灰度图中第i行j列像素的灰度值:
得到本发明的灰度图如图4所示。
为了直观的证明本发明较传统方法的优势,生成图3所示的三维点云数据对应的深度图(如图6所示)、方位角图(如图7所示)以及本发明的灰度图(如图5所示),可见传统深度图在纹理信息上丢失严重,而方位角图的景物边界区分度十分模糊。分别对图5、图6、图7进行超像素分割比较(分割数量为500块),效果如图8所示。计算超像素分割后,超像素块中的总错分像素点在总像素点中的比例(错误率),结果如下表所示:
本发明灰度图 | 深度图 | 方位角图 | |
错误率 | 2.13% | 4.41% | 4.25% |
可见本发明有着更好的分割效果,与传统的深度图与方位角图构图方法相比,提高了图像的景物边界区分度以及细节质量,图像纹理更加清晰。能为后续的图像分割、物体识别及场景理解等进一步研究提供有效的支持。
Claims (1)
1.一种三维激光点云数据的二维图像化表述方法,其特征在于,完成该方法的系统主要由移动机器人、激光扫描器和计算机系统组成;激光扫描器和计算机系统均安装在移动机器人平台上,激光扫描器和计算机系统之间通过数据线进行通信;激光扫描器采集图像的三维点云数据信息,计算机系统将三维激光点云数据转化为二维灰度图像;计算机系统软件要求采集的三维激光点云数据可以区分属于每一帧激光线的激光数据点,以便能够使三维激光点云数据以二维数组的形式存储,并且使每个激光数据点和灰度图的像素点一一对应;将获得的激光点云数据存入二维数组中,使得二维数组的每一列激光数据点属于同一帧激光线,并在二维数组中截取所需大小的场景,用截取场景的三维点云数据构建灰度图;
具体将三维激光点云数据转化为二维灰度图像计算方法如下:
(1)最优深度对应灰度值的计算
将场景点云投影到水平面,运用统计学原理计算场景点云分布的主方向,由此获得场景的最优深度信息:
<mrow>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>cov</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cov</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>cov</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>cov</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)为点云在水平面投影的协方差矩阵,其中(x,y)为点云在水平面投影的坐标,cov代表协方差运算;通过计算协方差矩阵的特征向量V1、V2即可估计出场景的点云分布主方向和最优深度方向;设|V1|>|V2|,则V1为场景点云分布的主方向;由点云在水平面上的投影区域的主方向向量V1可求出过原点且平行于V1的直线l1的斜率K,进而确定直线l1的方程:
K*X+Y=0 (2)
由点到直线的距离公式:
<mrow>
<mi>D</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mrow>
<mi>K</mi>
<mo>*</mo>
<mi>x</mi>
<mo>+</mo>
<mi>y</mi>
</mrow>
<mo>|</mo>
</mrow>
<msqrt>
<mrow>
<msup>
<mi>K</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中D为位于i行j列的激光点Pi,j的投影(x,y)到直线l1的距离,这个距离便是点云中激光点Pi,j的最优深度Di,j;将最优深度Di,j映射至0-255的灰度范围,得到点云中激光点Pi,j对应于最优深度Di,j的灰度值
<mrow>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>min</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mn>255</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Dmin为点云最优深度值的最小值,Dmax为点云最优深度值的最大值;
(2)方位向量长度的计算
对于激光点云中位于第i行j列的激光点Pi,j,取其上方相邻的激光点Pi-1,j,设点Pi,j到点Pi-1,j的向量为Vd;向量Vd的长度|Vd|即为点Pi,j的方位向量长度Li,j;方位向量长度能够反应场景中景物的纹理变化,对景物边缘的表现效果尤其明显,所以为了保证较好的效果,在映射为灰度值时不能采用单纯的线性映射;
(3)方位向量长度灰度数值均衡化
求出激光点数量为I*J的点云中每个点所对应的方位向量长度后,在映射为灰度值时,为了保证得到的灰度图具有高对比度、多变的灰度色调和明显的景物边界,对方位向量长度采取灰度数值均衡化;
对I*J个激光点对应的方位向量长度Li,j按数值大小升序排列,把0-255的灰度值范围视作256个灰度值区间,求出每个灰度值区间所对应的向量长度数据量N:
<mrow>
<mi>N</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>I</mi>
<mo>*</mo>
<mi>J</mi>
</mrow>
<mn>256</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
通过N与升序排列的方位向量长度Li,j,即可确定每个灰度值区间的端点长度值;从而得到点云中激光点Pi,j对应于方位向量长度Li,j的灰度值
<mrow>
<msub>
<mi>P</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mo>{</mo>
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo><</mo>
<mo>=</mo>
<msub>
<mi>L</mi>
<mn>0</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>n</mi>
</mtd>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo><</mo>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo><</mo>
<mo>=</mo>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3...</mn>
<mo>,</mo>
<mn>255</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中L0,L1,L2…,L255为256个灰度值区间的端点长度值;
(4)权重动态调整
将点云中激光点Pi,j对应于最优深度Di,j的灰度值与对应于方位向量长度Li,j的灰度值加权平均得到激光点Pi,j在灰度图中对应的灰度值PGi,j:
<mrow>
<msub>
<mi>PG</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>aP</mi>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>+</mo>
<msub>
<mi>bP</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
</mrow>
<mrow>
<mi>a</mi>
<mo>+</mo>
<mi>b</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中PGi,j为灰度图中第i行j列像素的灰度值;a与b为权值,且a+b=10;此时所生成灰度图的灰度值矩阵的方差δ2为:
<mrow>
<msup>
<mi>&delta;</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>PG</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mrow>
<mi>P</mi>
<mi>G</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>I</mi>
<mo>*</mo>
<mi>J</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中为灰度图中像素灰度均值;当方差最大时,则意味着灰度图具有高对比度和多变的灰度色调;取使方差最大的权值a与b生成灰度图,此时得到的灰度图是最优的;利用传统求最优解的方法求解式(8),得到使δ2为最大值的a与b,记做与
(5)生成灰度图
求出使δ2为最大值的权值与后,遍历点云数据中的I*J个激光点,求出每个激光点对应的灰度值:
<mrow>
<msub>
<mi>PG</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
<mo>+</mo>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
<msub>
<mi>P</mi>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
</msub>
</mrow>
<mrow>
<mover>
<mi>a</mi>
<mo>~</mo>
</mover>
<mo>+</mo>
<mover>
<mi>b</mi>
<mo>~</mo>
</mover>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
PGi,j即为大小为I*J的灰度图中第i行j列像素的灰度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510938881.6A CN105574905B (zh) | 2015-12-15 | 2015-12-15 | 一种三维激光点云数据的二维图像化表述方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510938881.6A CN105574905B (zh) | 2015-12-15 | 2015-12-15 | 一种三维激光点云数据的二维图像化表述方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105574905A CN105574905A (zh) | 2016-05-11 |
CN105574905B true CN105574905B (zh) | 2018-01-16 |
Family
ID=55884992
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510938881.6A Active CN105574905B (zh) | 2015-12-15 | 2015-12-15 | 一种三维激光点云数据的二维图像化表述方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105574905B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018039871A1 (zh) * | 2016-08-29 | 2018-03-08 | 北京清影机器视觉技术有限公司 | 三维视觉测量数据的处理方法和装置 |
WO2018122335A1 (en) * | 2016-12-30 | 2018-07-05 | Robert Bosch Gmbh | Mobile robotic device that processes unstructured data of indoor environments to segment rooms in a facility to improve movement of the device through the facility |
CN111174788B (zh) * | 2018-11-13 | 2023-05-02 | 北京京东乾石科技有限公司 | 一种室内二维建图方法和装置 |
CN110163064B (zh) * | 2018-11-30 | 2022-04-05 | 腾讯科技(深圳)有限公司 | 一种道路标志物的识别方法、装置及存储介质 |
CN111327906B (zh) * | 2018-12-13 | 2022-08-09 | 华为技术有限公司 | 点云编解码方法和编解码器 |
CN110136130A (zh) * | 2019-05-23 | 2019-08-16 | 北京阿丘机器人科技有限公司 | 一种检测产品缺陷的方法及装置 |
CN111220967B (zh) * | 2020-01-02 | 2021-12-10 | 小狗电器互联网科技(北京)股份有限公司 | 一种激光雷达数据有效性的检测方法与装置 |
CN115412721B (zh) * | 2021-05-26 | 2024-05-28 | 荣耀终端有限公司 | 一种点云二维规则化平面投影方法及装置 |
CN113538547A (zh) * | 2021-06-03 | 2021-10-22 | 苏州小蜂视觉科技有限公司 | 一种3d线激光传感器的深度处理方法及点胶设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4474487B2 (ja) * | 2009-04-13 | 2010-06-02 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
CN104268933A (zh) * | 2014-09-11 | 2015-01-07 | 大连理工大学 | 一种车载二维激光运动中三维环境扫描成像方法 |
-
2015
- 2015-12-15 CN CN201510938881.6A patent/CN105574905B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4474487B2 (ja) * | 2009-04-13 | 2010-06-02 | キヤノン株式会社 | 画像処理装置及び画像処理方法 |
CN104268933A (zh) * | 2014-09-11 | 2015-01-07 | 大连理工大学 | 一种车载二维激光运动中三维环境扫描成像方法 |
Non-Patent Citations (2)
Title |
---|
Using Scale Coordination and Semantic Information for Robust 3-D Object Recognition by a Service Robot;Yan Zhuang等;《IEEE SENSORS JOURNAL》;20150131;第15卷(第1期);第37-47页 * |
机载三维激光点云数据处理技术研究;孙兴华 等;《科技资讯》;20150623(第14期);第1-2页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105574905A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105574905B (zh) | 一种三维激光点云数据的二维图像化表述方法 | |
CN110264567B (zh) | 一种基于标记点的实时三维建模方法 | |
CN109074083B (zh) | 移动控制方法、移动机器人及计算机存储介质 | |
CN110097553B (zh) | 基于即时定位建图与三维语义分割的语义建图系统 | |
CN113436260B (zh) | 基于多传感器紧耦合的移动机器人位姿估计方法和系统 | |
CN110853075B (zh) | 一种基于稠密点云与合成视图的视觉跟踪定位方法 | |
CN109903327A (zh) | 一种稀疏点云的目标物尺寸测量方法 | |
CN106709947A (zh) | 一种基于rgbd相机的三维人体快速建模系统 | |
CN106774296A (zh) | 一种基于激光雷达和ccd摄像机信息融合的障碍检测方法 | |
CN108230242B (zh) | 一种从全景激光点云到视频流的转换方法 | |
CN103632167B (zh) | 类地重力场环境下单目视觉空间识别方法 | |
CN205451195U (zh) | 一种基于多摄像机的实时三维点云重建系统 | |
CN106875437A (zh) | 一种面向rgbd三维重建的关键帧提取方法 | |
CN101777129A (zh) | 一种基于特征检测的图像匹配方法 | |
CN113050074B (zh) | 无人驾驶环境感知中相机与激光雷达标定系统及标定方法 | |
CN107514745A (zh) | 一种智能空调立体视觉定位的方法及系统 | |
CN111292369B (zh) | 激光雷达的伪点云数据生成方法 | |
CN113112588A (zh) | 基于rgb-d深度相机重建的地下管井三维可视化方法 | |
CN108362205A (zh) | 基于条纹投影的空间测距方法 | |
CN114298151A (zh) | 一种基于点云数据与图像数据融合的3d目标检测方法 | |
CN106709432B (zh) | 基于双目立体视觉的人头检测计数方法 | |
CN113487631B (zh) | 基于lego-loam的可调式大角度探测感知及控制方法 | |
CN114463357A (zh) | 一种在重介质选煤中实时确定介质堆动态信息的方法 | |
CN112581505B (zh) | 一种简易的激光雷达点云与光学影像自动配准方法 | |
CN108734148A (zh) | 一种基于云计算的公共场合图像信息采集无人机控制系统 |
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 |