CN109035345A - 基于高斯过程回归的tof相机距离校正方法 - Google Patents

基于高斯过程回归的tof相机距离校正方法 Download PDF

Info

Publication number
CN109035345A
CN109035345A CN201810804240.5A CN201810804240A CN109035345A CN 109035345 A CN109035345 A CN 109035345A CN 201810804240 A CN201810804240 A CN 201810804240A CN 109035345 A CN109035345 A CN 109035345A
Authority
CN
China
Prior art keywords
tof camera
distance
target board
gaussian process
central pixel
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.)
Pending
Application number
CN201810804240.5A
Other languages
English (en)
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.)
Qilu University of Technology
Original Assignee
Qilu University of Technology
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 Qilu University of Technology filed Critical Qilu University of Technology
Priority to CN201810804240.5A priority Critical patent/CN109035345A/zh
Publication of CN109035345A publication Critical patent/CN109035345A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T7/85Stereo camera calibration
    • 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/10004Still image; Photographic image

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Measurement Of Optical Distance (AREA)

Abstract

本发明公开了基于高斯过程回归的TOF相机距离校正方法,包括:将TOF相机安装在三角支架上,将TOF相机的镜头平面调整为与镜头前方的目标板平行;目标板上贴有相互垂直的两条反光贴;通过调整相机与目标板相对位置,使得TOF相机的中心像素点定位在两条反光贴的交汇区域的中心位置;用TOF相机获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的待校正距离;用激光测距仪获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的距离真值;基于训练集的待校正距离和训练集的距离真值,建立高斯过程回归模型;基于高斯过程回归模型,对TOF相机的新测距值进行校正。

Description

基于高斯过程回归的TOF相机距离校正方法
技术领域
本发明涉及TOF相机技术领域,特别是涉及基于高斯过程回归的TOF相机距离校正方法。
背景技术
基于TOF(飞行时间法,Time-of-flight)技术的三维摄像机是一种新型的、小型化三维成像设备,可同时捕捉动态目标的灰度图像和深度信息,因而轻易解决了机器视觉领域的传统难题,使用这些深度信息不但可以轻松的完成目标图像的分割、标记、识别、跟踪等传统应用,并且能在许多新兴领域发挥重要的作用。
20世纪90年代,首次诞生了第一代的TOF相机,在最近几年,该类型的相机变得越来越流行,并应用在了机器人、无人驾驶车辆、互动娱乐、救援以及医疗和生物等行业。但由于TOF摄像机发展时间比较短,而且尚处在实验研究阶段,并没有进入大规模的生产和实用阶段,因而还存在着比较多的缺点和不足。在性能参数上其主要技术指标还需通过实验做进一步验证,特别是该相机的测距特性及距离误差校正方面。
TOF相机的测距结果受多种因素的影响,如积分时间引起的误差、相机温度变化引起的误差以及系统本身自带的误差等。现今许多学者已对TOF相机的测距误差进行了研究,常用的距离校正方法有:
(1)通过离线测量不同积分时间下的距离误差来统计相机的误差参数,然后采用在线查表的方式对该误差进行校正;
(2)采用正弦函数叠加组合的方式进行距离校正;
(3)高阶多项式插值实现距离校正。
上述校正方法虽可实现距离数据的误差修正,然而查表方式信息量大且不能覆盖所有数据,三角函数叠加组合需要复杂的组合算法而且要对函数周期进行调整,多项式插值存在插值阶次过高的问题。
发明内容
为了解决现有技术的不足,本发明提供了基于高斯过程回归的TOF相机距离校正方法。TOF相机的测距结果受多种因素的影响,如积分时间引起的误差、相机温度变化引起的误差以及系统本身自带的误差等。由于TOF相机一般自带有温度反馈校正,因而该申请不考虑温度对相机的影响,只对积分时间引起的误差及系统本身自带的误差进行分析与校正。
基于高斯过程回归的TOF相机距离校正方法,包括:
将TOF相机安装在三角支架上,将TOF相机的镜头平面调整为与镜头前方的目标板平行;目标板上贴有相互垂直的两条反光贴;通过调整相机与目标板相对位置,使得TOF相机的中心像素点定位在两条反光贴的交汇区域的中心位置;
用TOF相机获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的待校正距离;用激光测距仪获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的距离真值;基于训练集的待校正距离和训练集的距离真值,建立高斯过程回归模型;基于高斯过程回归模型,对TOF相机的新测距值进行校正。
进一步的,用TOF相机获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的待校正距离;用激光测距仪获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的距离真值;基于训练集的待校正距离和训练集的距离真值,建立高斯过程回归模型,具体步骤为:
步骤(1):假设TOF相机包括m个积分时间;在TOF相机中心像素点与目标板之间的设定距离范围内,以设定距离间隔选择若干个观测点;设定TOF相机的当前积分时间;
步骤(2):在设定的积分时间下,当TOF相机位于第一个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;然后,当TOF相机位于第二个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;然后,当TOF相机位于第三个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;以此类推,直至所有观测点的TOF相机中心像素点与目标板之间的距离数据均被读取完毕;设定TOF相机的下一个积分时间,然后,重复执行步骤(2);直至所有积分时间均被选择完毕;进入步骤(3);
用激光测距仪获取每个观测点下,TOF相机的中心像素点与目标板的距离真值;
步骤(3):对每个观测点读取的n次距离数据,并计算每次的距离误差ei,距离误差均值μ以及方差s;
步骤(4):将每个测量点在不同积分时间下的距离误差均值μ按照列向量的形式组成误差观测值矩阵y;
步骤(5):对于高斯过程回归模型
y=f(x)+ε (4)
式中:x为输入向量,其组成元素为每个测量点的距离真值r,f为函数值;
假设噪声得到y的先验分布:
其中,K(X,X)=(kij)为n×n阶对称正定的协方差矩阵,矩阵元素kij=k(xi,xj)用来衡量变量xi和变量xj之间的相关性;In为n维单位矩阵;σn为噪声方差;
将TOF相机在任意距离下的误差预测值设为f*,得误差观测值矩阵y和误差预测值f*的联合先验分布
其中,K(X,x*)=K(x*,X)T,K(X,x*)为测试点x*与训练集的输入X之间的n×1阶协方差矩阵;K(x*,x*)为测试点x*自身的协方差,即步骤(3)求得的方差s;
由公式(6)则计算出误差预测值f*的后验分布为
其中,
为测试点x*对应的误差预测值f*的均值;
为测试点x*对应的误差预测值f*的方差;
协方差函数选择平方指数协方差函数:
其中,∑=diag(l2),l为方差尺度,为信号方差;
参数集合为超参数,通过极大似然法求得;
由于TOF相机的测距误差是关于距离的一元函数,公式(10)简化为:
公式(8)、(9)及(11)组成高斯过程回归模型。
进一步的,所述步骤(3)的每次的距离误差ei,距离误差均值μ以及方差s:
ei=di-r (1)
其中,di为每次TOF相机读取的距离数据,r为由激光测距仪获得的距离真值。
进一步的,基于高斯过程回归模型,对TOF相机的新测距值进行校正的具体步骤为:
针对TOF相机的新测距值D作为输入值x*输入到公式(8)中,通过计算得到
由于从而得到
根据得到的得最终校正后的结果D'为:
与现有技术相比,本发明的有益效果是:
采用高斯过程回归模型对TOF相机的距离数据进行误差校正,解决了采用传统多项式或三角函数模型进行误差修正时模型阶次过高及函数组合复杂的问题。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1是本申请的流程图;
图2是测试场景;
图3误差校正结果;
其中,1.TOF相机,2目标板,3反光贴,4三角支架。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
首先判断TOF相机测距结果是否满足高斯过程回归要求,也即测量结果是否满足高斯分布,方法为:将TOF固定在三角支架上并与前方的目标板平行,以TOF相机中心像素点作为测试对象,获取该像素与目标板不同距离下的测量结果。采用高斯过程回归算法进行误差校正。校正操作:
1、将TOF相机1固定在三角支架4上并与前方的目标板2垂直平行,如图1所示。以TOF相机中心像素点作为标定对象,如图2所示,所分别选择不同积分时间(积分时间即为相机曝光时间)下的测距结果作为校正对象。为确保中心像素点垂直于目标板,目标板中贴有相互垂直的反光贴3。每条反光贴为长方形;通过调整相机角度使得中心像素点位于TOF相机获取图像亮度最强区域的中心时,认为此时中心像素刚好与目标板垂直。为衡量系统本身自带的误差,将距离数据扩展为1000mm~4500mm范围,并以50mm为间隔进行分点采集,每点处连续读取10次距离数据,并计算该点的均值μ、标准差s、及误差e。
e=μ-r (3)
式中,di为每次读取的距离数据,n为测试次数,此处为10,r为距离真值。
2、在不同积分时间下,将距离观测值的均值μ组成观测值矩阵y。对于回归过程考虑如下模型:
y=f(x)+ε (4)
式中:x为输入向量,其组成元素为每个测量点的距离真值r,f为函数值,进一步假设噪声可以得到y的先验分布
将TOF相机在任意距离下的预测值设为f*,可得观测值y和预测值f*的联合先验分布
其中,K(X,X)=(kij)为n×n阶对称正定的协方差矩阵,矩阵元素kij=k(xi,xj)用来衡量变量xi和变量xj之间的相关性;K(X,x*)=K(x*,X)T为测试点x*与训练集的输入X之间的n×1阶协方差矩阵;k(x*,x*)为测试点x*自身的协方差;σn为噪声方差;In为n维单位矩阵。
由上式则可计算出预测值f*的后验分布为
其中,
即为测试点x*对应的预测值f*的均值和方差。
GPR可以选取多种不同的协方差函数,在此选择平方指数协方差函数
式中,∑=diag(l2),l为方差尺度,为信号方差。参数集合即为超参数,一般通过极大似然法求得。
由于TOF相机的测距误差是关于距离的一元函数,公式(10)可简化为
公式(8)、(9)及(11)组成了GPR模型所需的基本函数。
图3展示了采用上述算法进行距离误差修正后的结果,从该图可知,经过校正后距离误差由15mm降到了5mm,可满要求校正要求。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (4)

1.基于高斯过程回归的TOF相机距离校正方法,其特征是,包括:
将TOF相机安装在三角支架上,将TOF相机的镜头平面调整为与镜头前方的目标板平行;目标板上贴有相互垂直的两条反光贴;通过调整相机与目标板相对位置,使得TOF相机的中心像素点定位在两条反光贴的交汇区域的中心位置;
用TOF相机获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的待校正距离;用激光测距仪获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的距离真值;基于训练集的待校正距离和训练集的距离真值,建立高斯过程回归模型;基于高斯过程回归模型,对TOF相机的新测距值进行校正。
2.如权利要求1所述的基于高斯过程回归的TOF相机距离校正方法,其特征是,
用TOF相机获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的待校正距离;用激光测距仪获取TOF相机的中心像素点与目标板在不同距离下的测距结果,作为训练集的距离真值;基于训练集的待校正距离和训练集的距离真值,建立高斯过程回归模型,具体步骤为:
步骤(1):假设TOF相机包括m个积分时间;在TOF相机中心像素点与目标板之间的设定距离范围内,以设定距离间隔选择若干个观测点;设定TOF相机的当前积分时间;
步骤(2):在设定的积分时间下,当TOF相机位于第一个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;然后,当TOF相机位于第二个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;然后,当TOF相机位于第三个观测点时,TOF相机读取n次TOF相机中心像素点与目标板之间的距离数据;以此类推,直至所有观测点的TOF相机中心像素点与目标板之间的距离数据均被读取完毕;设定TOF相机的下一个积分时间,然后,重复执行步骤(2);直至所有积分时间均被选择完毕;进入步骤(3);
用激光测距仪获取每个观测点下,TOF相机的中心像素点与目标板的距离真值;
步骤(3):对每个观测点读取的n次距离数据,并计算每次的距离误差ei,距离误差均值μ以及方差s;
步骤(4):将每个测量点在不同积分时间下的距离误差均值μ按照列向量的形式组成误差观测值矩阵y;
步骤(5):对于高斯过程回归模型
y=f(x)+ε (4)
式中:x为输入向量,其组成元素为每个测量点的距离真值r,f为函数值;
假设噪声得到y的先验分布:
其中,K(X,X)=(kij)为n×n阶对称正定的协方差矩阵,矩阵元素kij=k(xi,xj)用来衡量变量xi和变量xj之间的相关性;In为n维单位矩阵;σn为噪声方差;
将TOF相机在任意距离下的误差预测值设为f*,得误差观测值矩阵y和误差预测值f*的联合先验分布
其中,K(X,x*)=K(x*,X)T,K(X,x*)为测试点x*与训练集的输入X之间的n×1阶协方差矩阵;K(x*,x*)为测试点x*自身的协方差,即步骤(3)求得的方差s;
由公式(6)则计算出误差预测值f*的后验分布为
其中,
为测试点x*对应的误差预测值f*的均值;
为测试点x*对应的误差预测值f*的方差;
协方差函数选择平方指数协方差函数:
其中,∑=diag(l2),l为方差尺度,为信号方差;
参数集合为超参数,通过极大似然法求得;
由于TOF相机的测距误差是关于距离的一元函数,公式(10)简化为:
公式(8)、(9)及(11)组成高斯过程回归模型。
3.如权利要求2所述的基于高斯过程回归的TOF相机距离校正方法,其特征是,
所述步骤(3)的每次的距离误差ei,距离误差均值μ以及方差s:
ei=di-r (1)
其中,di为每次TOF相机读取的距离数据,r为由激光测距仪获得的距离真值。
4.如权利要求2所述的基于高斯过程回归的TOF相机距离校正方法,其特征是,
基于高斯过程回归模型,对TOF相机的新测距值进行校正的具体步骤为:
针对TOF相机的新测距值D作为输入值x*输入到公式(8)中,通过计算得到
由于从而得到根据得到的得最终校正后的结果D'为:
CN201810804240.5A 2018-07-20 2018-07-20 基于高斯过程回归的tof相机距离校正方法 Pending CN109035345A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810804240.5A CN109035345A (zh) 2018-07-20 2018-07-20 基于高斯过程回归的tof相机距离校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810804240.5A CN109035345A (zh) 2018-07-20 2018-07-20 基于高斯过程回归的tof相机距离校正方法

Publications (1)

Publication Number Publication Date
CN109035345A true CN109035345A (zh) 2018-12-18

Family

ID=64644752

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810804240.5A Pending CN109035345A (zh) 2018-07-20 2018-07-20 基于高斯过程回归的tof相机距离校正方法

Country Status (1)

Country Link
CN (1) CN109035345A (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109816735A (zh) * 2019-01-24 2019-05-28 哈工大机器人(合肥)国际创新研究院 一种快速标定及校正方法及其tof相机
CN110058212A (zh) * 2019-05-15 2019-07-26 上海炬佑智能科技有限公司 标靶和tof相机标定集成系统
CN110196422A (zh) * 2019-06-26 2019-09-03 Oppo广东移动通信有限公司 激光测距的测试方法、测试装置及移动终端
CN111562562A (zh) * 2020-04-28 2020-08-21 重庆市天实精工科技有限公司 基于tof的3d成像模块校准方法
CN112415492A (zh) * 2019-08-22 2021-02-26 Enscape有限公司 具有tof功能的相机模块的检查装置
CN112965048A (zh) * 2021-02-09 2021-06-15 深圳市灵明光子科技有限公司 激光测距误差校正方法、装置、电子设备和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011094366A1 (en) * 2010-01-27 2011-08-04 Intersil Americas Inc. Gesture recognition with principal component anaysis
CN107270890A (zh) * 2017-04-25 2017-10-20 深圳大学 一种无人机上tof距离传感器的测距方法及系统
CN108197380A (zh) * 2017-12-29 2018-06-22 南京林业大学 基于偏最小二乘的高斯回归软测量建模方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011094366A1 (en) * 2010-01-27 2011-08-04 Intersil Americas Inc. Gesture recognition with principal component anaysis
CN107270890A (zh) * 2017-04-25 2017-10-20 深圳大学 一种无人机上tof距离传感器的测距方法及系统
CN108197380A (zh) * 2017-12-29 2018-06-22 南京林业大学 基于偏最小二乘的高斯回归软测量建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张慧 等: "四足机器人地形识别与路径规划算法", 《机器人》 *
张慧: "四足机器人环境感知、识别与领航员跟随算法研究", 《中国博士学位论文全文数据库 信息科技辑(月刊)》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109816735A (zh) * 2019-01-24 2019-05-28 哈工大机器人(合肥)国际创新研究院 一种快速标定及校正方法及其tof相机
CN109816735B (zh) * 2019-01-24 2022-10-21 哈工大机器人(合肥)国际创新研究院 一种快速标定及校正方法及其tof相机
CN110058212A (zh) * 2019-05-15 2019-07-26 上海炬佑智能科技有限公司 标靶和tof相机标定集成系统
CN110196422A (zh) * 2019-06-26 2019-09-03 Oppo广东移动通信有限公司 激光测距的测试方法、测试装置及移动终端
CN112415492A (zh) * 2019-08-22 2021-02-26 Enscape有限公司 具有tof功能的相机模块的检查装置
CN111562562A (zh) * 2020-04-28 2020-08-21 重庆市天实精工科技有限公司 基于tof的3d成像模块校准方法
CN111562562B (zh) * 2020-04-28 2023-04-14 重庆市天实精工科技有限公司 基于tof的3d成像模块校准方法
CN112965048A (zh) * 2021-02-09 2021-06-15 深圳市灵明光子科技有限公司 激光测距误差校正方法、装置、电子设备和存储介质
WO2022170933A1 (zh) * 2021-02-09 2022-08-18 深圳市灵明光子科技有限公司 激光测距误差校正方法、装置、电子设备和存储介质

Similar Documents

Publication Publication Date Title
CN109035345A (zh) 基于高斯过程回归的tof相机距离校正方法
CN107229930B (zh) 一种指针式仪表数值智能识别方法
CN105069743B (zh) 探测器拼接实时图像配准的方法
CN106548462B (zh) 基于薄板样条插值的非线性sar图像几何校正方法
CN104730802B (zh) 光轴夹角的校准、对焦方法和系统和双摄像头设备
CN109598762A (zh) 一种高精度双目相机标定方法
WO2018196303A1 (zh) 多方位投影的投影仪标定方法及装置
CN104104937B (zh) 图像处理装置以及图像处理方法
CN107726975B (zh) 一种基于视觉拼接测量的误差分析方法
CN111210468A (zh) 一种图像深度信息获取方法及装置
CN107680137B (zh) 一种基于球体投影公切线的多相机标定及参数优化方法
CN109827502A (zh) 一种标定点图像补偿的线结构光视觉传感器高精度标定方法
CN105378794A (zh) 3d拍摄装置、用于建立3d图像的方法和用于构建3d拍摄装置的方法
CN104766309A (zh) 一种平面特征点导航定位方法与装置
CN108921864A (zh) 一种光条中心提取方法及装置
CN110223355B (zh) 一种基于双重极线约束的特征标志点匹配方法
CN107292925A (zh) 基于Kinect深度相机测量方法
CN109523595A (zh) 一种建筑工程直线棱角间距视觉测量方法
CN108470356A (zh) 一种基于双目视觉的目标对象快速测距方法
CN109974618A (zh) 多传感器视觉测量系统的全局标定方法
CN109799073A (zh) 一种光学畸变测量装置及方法、图像处理系统、电子设备和显示设备
CN109186347A (zh) 一种枪用瞄准镜零位走动量的检测装置及方法
CN109727291A (zh) 一种变焦距摄像机的高精度在线标定方法
CN109003312A (zh) 一种基于非线性优化的相机标定方法
CN110414101B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20181218