CN110763140B - 一种非平行光轴高精度双目测距方法 - Google Patents

一种非平行光轴高精度双目测距方法 Download PDF

Info

Publication number
CN110763140B
CN110763140B CN201910608662.XA CN201910608662A CN110763140B CN 110763140 B CN110763140 B CN 110763140B CN 201910608662 A CN201910608662 A CN 201910608662A CN 110763140 B CN110763140 B CN 110763140B
Authority
CN
China
Prior art keywords
eye camera
camera
optical axis
image
distance
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
Application number
CN201910608662.XA
Other languages
English (en)
Other versions
CN110763140A (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.)
Beijing Aerospace Times Technology Development Co ltd
Original Assignee
Beijing Aerospace Times Technology Development 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 Beijing Aerospace Times Technology Development Co ltd filed Critical Beijing Aerospace Times Technology Development Co ltd
Priority to CN201910608662.XA priority Critical patent/CN110763140B/zh
Publication of CN110763140A publication Critical patent/CN110763140A/zh
Application granted granted Critical
Publication of CN110763140B publication Critical patent/CN110763140B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/02Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
    • G01B11/022Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness by means of tv-camera scanning

Abstract

本发明涉及计算机视觉技术领域,具体涉及一种非平行光轴高精度双目测距方法,首先搭建双目测距装置,并建立与双目测距装置在左右相机光轴非理想平行时配合使用的测距公式;然后标定测距公式中的参数,获得未知量只包含待测场景目标在双目测距装置左右相机图像中的成像位置的测距公式;实际应用时,将直接读取的待测场景目标在左右相机图像中的成像位置数据代入标定后的测距公式即可计算获得待测场景目标与双目测距装置中参考相机的距离。本发明适用于光轴非理想平行状态的双目立体视觉装置,降低了由于安装导致的光轴不平行带来的测距误差,提高了实际工程应用中的测距精度。

Description

一种非平行光轴高精度双目测距方法
技术领域
本发明涉及计算机视觉技术领域,具体涉及一种非平行光轴高精度双目测距方法。
背景技术
现有的双目测距技术,一般假设两台摄像机的像平面位于同一平面上、光轴严格平行、距离固定、焦距相同。通过双目视觉原理图和三角测量公式,可根据视差求出测距目标与两台摄像机之间的垂直距离,计算公式如下:
Figure BSA0000185642510000011
其中,Z表示测距目标P与两台摄像机(OL、OR)之间的垂直距离;B表示两台摄像机之间的距离,即基线长度;f表示两台摄像机的焦距;XL和XR分别表示测距目标在左目摄像机和右目摄像机上的像素坐标;k表示像素的实际物理尺寸,|XL-XR|·k即左目摄像机和右目摄像机的视差,见图1。
上述公式的成立是在满足理想模型的情况下才成立的。在实际应用中,左目摄像机和右目摄像机几乎不可能像理想模型中那样严格的光轴平行,且基线长度B和焦距f的真实值与测量值之间存在一定的误差。光轴的平行度、基线长度、焦距的偏差,会显著影响测距的精度,测距误差一般在5%以上。
发明内容
本发明的目的在于提供一种非平行光轴高精度双目测距方法,该测距方法能够在左右相机光轴之间可有小幅夹角时成立,测距精度高。
为了实现上述目的,本发明采用如下技术方案:一种非平行光轴高精度双目测距方法,包括如下步骤:
S1.建立双目测距系统,测距系统包括规格相同的左目摄像机和右目摄像机;
S2.建立与双目测距系统在左、右目摄像机光轴非理想平行时配合使用的测距公式;所述测距公式为:
Figure BSA0000185642510000021
其中:
Figure BSA0000185642510000022
上述公式中:
记待测距目标为P,左目摄像机的光心为OL,右目摄像机的光心为OR,OLOR的长度为B,P到OLOR的垂直距离为Z,垂点为OZ,α为左目摄像机的光轴与测距线POz夹角,β为右目摄像机的光轴与测距线POz夹角;k为摄像机的像素物理尺寸;f为焦距;XL为目标P在左目摄像机中成像点的像素坐标;XR为目标P在右目摄像机中成像点的像素坐标;Q为实测Z值相比真实Z的固定偏差;
S3.选择不同的距离,至少进行3次标定,分别获得标定点的实测距离Zn,标定点在左目图像的像素坐标为Xn R,在右目图像中的坐标为Xn L,n≥3;并将每次获得的Zn、Xn R和Xn L分别带入公式(1),求得M、N、Q,然后将M、N、Q带入公式(1);
S4.读取待测距目标在左、右摄像机中成像的像素坐标,带入由步骤S3获得的公式(1),求得Z。
进一步地,所述光轴非理想平行是指左、右目摄像机光轴之间有小幅度夹角,即夹角小于 15°。
进一步地,在所述步骤S3中,当n=3时,直接可求出M、N、Q;当n>3时,采用
Figure BSA0000185642510000023
从 n中任选3组进行计算,然后求取多次计算结果的平均值。
进一步地,所述测距系统还包括定时指令模块、计算模块和显示模块,定时指令模块分别与左目摄像机和右目摄像机相连,计算模块分别与左目摄像机和右目摄像机相连,定时指令模块与计算模块相连,计算模块与显示模块相连。
进一步地,在所述步骤S3中,每次标定时,触发定时指令模块,分别向左目摄像机和右目摄像机发出采集图像指令,指令中包含时间戳,左目摄像机和右目摄像机收到指令同时拍摄待测距目标的图像,并将图像编号、附带时间戳后转发给计算模块,计算模块通过图像处理获得待测距目标分别在左目图像和右目图像中的像素坐标。
进一步地,所述图像处理包括但不限于以下几种图像处理方法,包括图像匹配、二值化、噪声滤波和图像特征提取。
进一步地,触发定时指令模块时采用手工触发。
本发明提出一种光轴非平行双目测距方法。首先搭建双目测距装置,并建立与双目测距装置在左右相机光轴非理想平行时配合使用的测距公式;然后标定测距公式中的参数,获得未知量只包含待测场景目标在双目测距装置左右相机图像中的成像位置的测距公式;实际应用时,将直接读取的待测场景目标在左右相机图像中的成像位置数据代入标定后的测距公式即可计算获得待测场景目标与双目测距装置中参考相机的距离。本发明适用于光轴非理想平行状态的双目立体视觉装置,降低了由于安装导致的光轴不平行带来的测距误差,提高了实际工程应用中的测距精度。在50m~200m的有效测距范围内,Z值的测量误差在±5‰以内。借助于本发明的上述方法,基于三次自动标定过程的非平行光轴双目测距精度较现有技术提高一个数量级,且不受摄像机光轴平行度、基线长度、焦距、像素物理尺寸等误差的影响,测距精度极高。
附图说明
图1为左、右目摄像机光轴平行时的双目测距示意图;
图2为本发明的双目测距系统;
图3为左、右目摄像机光轴非平行时的双目测距示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
在实际应用中,左目摄像机和右目摄像机几乎不可能像理想模型中那样严格的光轴平行,且基线长度B和焦距f的真实值与测量值之间存在一定的误差。光轴的平行度、基线长度、焦距的偏差,会显著影响测距的精度,测距误差一般在5%以上。
针对上述问题,本发明提供一种非平行光轴高精度双目测距方法,包括以下步骤:
S1.建立双目测距系统,如图2所示,测距系统包括规格相同的左目摄像机和右目摄像机、定时指令模块、计算模块和显示模块,定时指令模块分别与左目摄像机和右目摄像机相连,计算模块分别与左目摄像机和右目摄像机相连,定时指令模块与计算模块相连,计算模块与显示模块相连。其中,左、右目摄像机的光轴处于非理想平行状态,即光轴之间可有小幅夹角(小于15°)。
S2.建立与双目测距系统在左、右目摄像机光轴非理想平行时配合使用的测距公式;
所述测距公式为:
Figure BSA0000185642510000041
其中:
Figure BSA0000185642510000042
上述公式中:
记待测距目标为P,左目摄像机的光心为OL,右目摄像机的光心为OR,OLOR的长度为B,P到OLOR的垂直距离为Z,垂点为OZ,α为左目摄像机的光轴与测距线POz夹角,β为右目摄像机的光轴与测距线POz夹角;k为摄像机的像素物理尺寸;f为焦距;XL为目标P在左目摄像机中成像点的像素坐标;XR为目标P在右目摄像机中成像点的像素坐标;Q为实测Z值相比真实Z的固定偏差。
下面以图3为辅助,进行具体的介绍,记待测距目标为P(有效测距范围50m~200m),两台摄像机的像素物理尺寸为k,焦距为f(50mm~100mm),两台摄像机光心之间的距离为 B(600mm~1000mm),即基线长度为B(600mm~1000mm)。
记左目摄像机的光心为OL,右目摄像机的光心为OR,目标P距离OLOR的垂直距离为Z,垂点为Oz,左目摄像机距离测距线POz的距离为BL;右目摄像机距离测距线POz的距离为BR
记左目摄像机的光轴与测距线POz夹角为α,左目摄像机采集到的目标成像点为PL,PL距离左目相机光轴距离为XL,PL的像素坐标为XL
记右目摄像机的光轴与测距线POz夹角为β,右目摄像机采集到的目标成像点为PR,PR距离右目摄像机光轴距离为XR,PR的像素坐标为XR
记左目摄像机与测距线POz的夹角为γ,记右目摄像机与测距线POz的夹角为θ。
则如图3所示的,各个点、线、角的几何关系如下:
B=BL+BR
Figure BSA0000185642510000051
Figure BSA0000185642510000052
Figure BSA0000185642510000053
Figure BSA0000185642510000054
XL+XR=|XL-XR|·k
其中,|XL-XR|表示测距目标在左目图像和右目图像中的像素坐标差,k为像素的实际物理尺寸。
进一步的,计算得出Z:
Figure BSA0000185642510000055
进一步的,令:
Figure BSA0000185642510000056
进一步的,计算得出:
Figure BSA0000185642510000057
进一步的,记标定时实测Z值时相比真实Z的固定偏差为Q,即摄像机外壳最前端与镜头中心点O之间的距离为Q,则得出就算距离的最终公式为:
Figure BSA0000185642510000058
S3.选择不同的距离,至少进行3次标定,分别获得标定点的实测距离Zn,标定点在左目图像的像素坐标为Xn R,在右目图像中的坐标为Xn L,n≥3;并将每次获得的Zn、Xn R和XN L分别带入公式(1),求得M、N、Q,然后将M、N、Q带入公式(1);
在所述步骤S3中,每次标定时,手工触发定时指令模块,分别向左目摄像机和右目摄像机发出采集图像指令,指令中包含时间戳,左目摄像机和右目摄像机收到指令同时拍摄待测距目标的图像,并将图像编号、附带时间戳后转发给计算模块,计算模块通过图像处理获得待测距目标分别在左目图像和右目图像中的像素坐标。图像处理包括但不限于以下几种图像处理方法,包括图像匹配、二值化、噪声滤波和图像特征提取。
在本实施例中,具体的图像处理步骤为:计算模块对两台摄像机发来的图像进行预处理,根据图像编号、时间戳对两台摄像机的图像进行配对,并进行二值化、噪声滤波处理,然后基于图像特征法提取待测距目标分别在左目图像和右目图像中的像素坐标。
选择三个不同的距离,进行三次标定:
第1次标定,实测距离为Z1=180,标定点在左目图像的像素坐标为X1 L,标定点在右目图像中的坐标为X1 R
第2次标定,实测距离为Z2=190m,标定点在左目图像的像素坐标为X2 L,标定点在右目图像中的坐标为X2 R
第3次标定,实测距离为Z3=210m,标定点在左目图像的像素坐标为X3 L,标定点在右目图像中的坐标为X3 R
进一步的,得出以下3个等式:
Figure BSA0000185642510000061
Figure BSA0000185642510000062
Figure BSA0000185642510000063
进一步的,记三次标定的视差分别为:
X1=|X1 L-X1 R|
X2=|X2 L-X2 R|
X3=|X3 L-X3 R|
进一步的,计算得出测距公式中的三个参数M、N、Q:
Figure BSA0000185642510000071
Figure BSA0000185642510000072
M=(Z1-Q)(X1+N)
进一步的,基于标定得出的参数值M、N、Q,和任一待测距目标在左目图像的像素坐标 XL和在右目图像中的像素坐标XR,根据以下公式计算距离Z:
Figure BSA0000185642510000073
在具体操作时,为了降低误差,可进行多次标定,例如标定n次,n>3,此时可选用
Figure BSA0000185642510000074
的方式,从n次测量中任选3个进行计算,最后求取多次计算结果的均值,得到M、N、Q。
S4.读取待测距目标在左、右摄像机中成像的像素坐标,带入由步骤S3获得的公式(1),求得Z。
触发指令模块,向左目摄像机和右目摄像机发出采集图像指令,指令中包含时间戳。左目摄像机和右目摄像机收到指令同时拍摄待测距目标的图像,并将图像编号、附带时间戳后转发给计算模块。计算模块对两台摄像机发来的图像进行预处理,根据图像编号、时间戳对两台摄像机的图像进行配对,并进行二值化、噪声滤波处理,然后基于图像特征法提取待测距目标分别在左目图像的像素坐标XL和右目图像中的像素坐标XR,代入上述公式计算得出距离Z,输出到显示模块上。
本发明提出一种光轴非平行双目测距方法。首先搭建双目测距装置,并建立与双目测距装置在左右相机光轴非理想平行时配合使用的测距公式;然后标定测距公式中的参数,获得未知量只包含待测场景目标在双目测距装置左右相机图像中的成像位置的测距公式;实际应用时,将直接读取的待测场景目标在左右相机图像中的成像位置数据代入标定后的测距公式即可计算获得待测场景目标与双目测距装置中参考相机的距离。本发明适用于光轴非理想平行状态的双目立体视觉装置,降低了由于安装导致的光轴不平行带来的测距误差,提高了实际工程应用中的测距精度。在50m~200m的有效测距范围内,Z值的测量误差在±5‰以内。借助于本发明的上述方法,基于三次自动标定过程的非平行光轴双目测距精度较现有技术提高一个数量级,且不受摄像机光轴平行度、基线长度、焦距、像素物理尺寸等误差的影响,测距精度极高。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (7)

1.一种非平行光轴高精度双目测距方法,其特征在于,包括如下步骤:
S1.建立双目测距系统,测距系统包括规格相同的左目摄像机和右目摄像机;
S2.建立与双目测距系统在左、右目摄像机光轴非理想平行时配合使用的测距公式;所述测距公式为:
Figure FSB0000192350410000011
其中:
Figure FSB0000192350410000012
Figure FSB0000192350410000013
上述公式中:
记待测距目标为P,左目摄像机的光心为OL,右目摄像机的光心为OR,OLOR的长度为B,P到OLOR的垂直距离为Z,垂点为OZ,α为左目摄像机的光轴与测距线POz夹角,β为右目摄像机的光轴与测距线POz夹角;k为摄像机的像素物理尺寸;f为焦距;XL为目标P在左目摄像机中成像点的像素坐标;XR为目标P在右目摄像机中成像点的像素坐标;Q为实测Z值相比真实Z的固定偏差;左目摄像机采集到的目标成像点为PL,PL距离左目相机光轴距离为XL,PL的像素坐标为XL;右目摄像机采集到的目标成像点为PR,PR距离右目摄像机光轴距离为XR,PR的像素坐标为XR,|XL-XR|表示测距目标在左目图像和右目图像中的像素坐标差;
S3.选择不同的距离,至少进行3次标定,分别获得标定点的实测距离Zn,标定点在左目图像的像素坐标为Xn R,在右目图像中的坐标为Xn L,n≥3;并将每次获得的Zn、Xn R和Xn L分别带入公式(1),求得M、N、Q,然后将M、N、Q带入公式(1);
S4.读取待测距目标在左、右摄像机中成像的像素坐标,带入由步骤S3获得的公式(1),求得Z。
2.根据权利要求1所述的非平行光轴高精度双目测距方法,其特征在于,所述光轴非理想平行是指左、右目摄像机光轴之间有小幅度夹角,即夹角小于15°。
3.根据权利要求1所述的非平行光轴高精度双目测距方法,其特征在于,在所述步骤S3中,当n=3时,直接可求出M、N、Q;当n>3时,采用
Figure FSB0000192350410000014
从n中任选3组进行计算,然后求取多次计算结果的平均值。
4.根据权利要求1所述的非平行光轴高精度双目测距方法,其特征在于,所述测距系统还包括定时指令模块、计算模块和显示模块,定时指令模块分别与左目摄像机和右目摄像机相连,计算模块分别与左目摄像机和右目摄像机相连,定时指令模块与计算模块相连,计算模块与显示模块相连。
5.根据权利要求4所述的非平行光轴高精度双目测距方法,其特征在于,在所述步骤S3中,每次标定时,触发定时指令模块,分别向左目摄像机和右目摄像机发出采集图像指令,指令中包含时间戳,左目摄像机和右目摄像机收到指令同时拍摄待测距目标的图像,并将图像编号、附带时间戳后转发给计算模块,计算模块通过图像处理获得待测距目标分别在左目图像和右目图像中的像素坐标。
6.根据权利要求5所述的非平行光轴高精度双目测距方法,其特征在于,所述图像处理包括但不限于以下几种图像处理方法,包括图像匹配、二值化、噪声滤波和图像特征提取。
7.根据权利要求5所述的非平行光轴高精度双目测距方法,其特征在于,触发定时指令模块时采用手工触发。
CN201910608662.XA 2019-07-08 2019-07-08 一种非平行光轴高精度双目测距方法 Active CN110763140B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910608662.XA CN110763140B (zh) 2019-07-08 2019-07-08 一种非平行光轴高精度双目测距方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910608662.XA CN110763140B (zh) 2019-07-08 2019-07-08 一种非平行光轴高精度双目测距方法

Publications (2)

Publication Number Publication Date
CN110763140A CN110763140A (zh) 2020-02-07
CN110763140B true CN110763140B (zh) 2021-04-16

Family

ID=69329316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910608662.XA Active CN110763140B (zh) 2019-07-08 2019-07-08 一种非平行光轴高精度双目测距方法

Country Status (1)

Country Link
CN (1) CN110763140B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111814659B (zh) * 2020-07-07 2024-03-29 杭州海康威视数字技术股份有限公司 一种活体检测方法、和系统
CN112785647A (zh) * 2021-01-27 2021-05-11 北京宇航时代科技发展有限公司 一种三目立体图像检测方法和系统
CN113848964A (zh) * 2021-09-08 2021-12-28 金华市浙工大创新联合研究院 一种非平行光轴双目测距方法
CN115797459B (zh) * 2022-08-29 2024-02-13 南京航空航天大学 一种任意焦距组合的双目视觉系统测距方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
CN105091849A (zh) * 2014-05-05 2015-11-25 南京理工大学 一种光轴非平行双目测距方法
WO2017080451A1 (en) * 2015-11-11 2017-05-18 Zhejiang Dahua Technology Co., Ltd. Methods and systems for binocular stereo vision
CN107147891A (zh) * 2017-05-17 2017-09-08 浙江大学 光轴可调节式三目深度获取摄像机
CN109813278A (zh) * 2019-02-26 2019-05-28 北京中科慧眼科技有限公司 测距模型修正方法、测距方法、装置以及自动驾驶系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
CN105091849A (zh) * 2014-05-05 2015-11-25 南京理工大学 一种光轴非平行双目测距方法
WO2017080451A1 (en) * 2015-11-11 2017-05-18 Zhejiang Dahua Technology Co., Ltd. Methods and systems for binocular stereo vision
CN107147891A (zh) * 2017-05-17 2017-09-08 浙江大学 光轴可调节式三目深度获取摄像机
CN109813278A (zh) * 2019-02-26 2019-05-28 北京中科慧眼科技有限公司 测距模型修正方法、测距方法、装置以及自动驾驶系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
低照度双目立体显著目标距离测定方法与实现;万一龙 等;《红外与激光工程》;20150331;第44卷(第3期);1053-1060 *

Also Published As

Publication number Publication date
CN110763140A (zh) 2020-02-07

Similar Documents

Publication Publication Date Title
CN110763140B (zh) 一种非平行光轴高精度双目测距方法
CN109767476B (zh) 一种自动对焦双目摄像头标定及深度计算方法
CN109859272B (zh) 一种自动对焦双目摄像头标定方法及装置
CN102494609B (zh) 一种基于激光探针阵列的三维摄影方法及装置
CN111243002A (zh) 应用于高精度三维测量的单目激光散斑投影系统标定及深度估计方法
KR101706093B1 (ko) 3차원 좌표 추출 시스템 및 그 방법
CN109579695B (zh) 一种基于异构立体视觉的零件测量方法
CN109443245B (zh) 一种基于单应性矩阵的多线结构光视觉测量方法
CN106709955B (zh) 基于双目立体视觉的空间坐标系标定系统和方法
WO2011125937A1 (ja) キャリブレーションデータ選択装置、選択方法、及び選択プログラム、並びに三次元位置測定装置
CN107941153B (zh) 一种激光测距优化标定的视觉系统
CN105717511B (zh) 基于线束激光器和普通摄像头芯片的多点测距方法
CN102221331A (zh) 一种基于不对称双目立体视觉技术的测量方法
CN111854622B (zh) 一种大视场光学动态变形测量方法
CN112556594A (zh) 一种融合红外信息的应变场和温度场耦合测量方法及系统
CN110849269A (zh) 一种田间玉米棒几何尺寸测量系统及方法
CN113255643A (zh) 一种应用于位移监测的机器视觉识别算法
CN110243338B (zh) 一种实时动态双目测距方法及系统
CN112907647B (zh) 一种基于固定式单目摄像机的三维空间尺寸测量方法
CN103487036B (zh) 基于空间前方交会和链码形似性的同名标记点匹配方法
CN111721194A (zh) 一种多激光线快速检测方法
CN109682312B (zh) 一种基于摄像头测量长度的方法及装置
CN111829435A (zh) 一种多双目相机与线激光协同检测方法
CN111768448A (zh) 一种基于多相机检测的空间坐标系标定方法
CN112950727B (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