一种单目立体视觉系统三维光路分析方法

1.本发明属于机器视觉测量领域，涉及一种单目立体视觉系统三维光路分析方法，通过光路分析实现三维层面上的反射镜有效几何参数以及单目立体视觉系统成像参数的计算。


背景技术：

2.三维视觉测量技术是机器视觉技术的一种，该技术通过图像信息和相机成像模型实现目标三维信息的测量，具有非接触、三维、实时、高精度以及空间全域的测量优势，被广泛应用于智能制造、智慧交通、航空航天、海洋船舶、能源动力等领域。由四个反射镜和一个真实相机组成的单目立体视觉系统有效的解决了传统单目视觉方法测量精度差的问题，同时也克服了多目视觉方法使用多相机造成的成本高、同步性差、占据空间大的劣势。通过光路分析可确定单目立体视觉系统反射镜有效几何参数以及单目立体视觉系统成像参数，对于设计合理的单目立体视觉系统具有重要意义，对该方面的研究颇具价值。
3.华南农业大学的赵祚喜等发明的专利号为cn 109099838 a“一种基于折反射的单目立体视觉系统三维信息测量方法”提出一种单目立体视觉测量方法，根据待测物体大小建立了单目立体视觉系统，该单目立体视觉系统包括若干矩形反射镜和一个单目相机，但该方法仅对二维光路进行分析计算，且该系统中反射镜要求互不遮挡且结构不对称，测量过程复杂且反射镜面得不到充分利用。此外，2006年，天津大学的闫丽在chinese journal of sensors and actuators期刊上发表了题为“单目立体视觉传感器的优化设计及精度分析”的论文，论文采用一台ccd相机和两组对称反射镜组建单目立体视觉系统，建立了透视成像数学模型，分析了传感器的几何结构参数，但其仅在二维层面建立单目立体视觉系统的数学模型并进行分析计算，同样未涉及在三维层面的光路分析和计算。


技术实现要素：

4.本发明要解决的技术难题是克服现有技术的缺陷，发明了一种单目立体视觉系统空间光路分析方法，目的是通过对由四个反射镜(即两个外镜与两个内镜)和一个真实相机组成的单目立体视觉系统进行空间光路分析，进而更加完备的计算反射镜的有效几何参数和单目立体视觉系统的成像参数。首先，使用四个反射镜和一个真实相机构建单目立体视觉系统；其次，选择合适的原点建立世界坐标系并初步确定本发明所采用的单目立体视觉系统的结构参数；然后，建立单目立体视觉系统的多镜面反射模型；最后，根据计算的反射光线和反射点求解出反射镜的有效几何参数与单目立体视觉系统的成像参数。
5.本发明所采用的技术方案是发明了单目立体视觉系统空间光路分析方法。其特征为：首先，通过使用四个反射镜和一个真实相机创建单目立体视觉系统；其次，在单目立体视觉系统上建立坐标系并定义结构参数；在此基础上，构建单目立体视觉系统多镜面反射模型来计算各反射镜的反射光线和反射点；最后，根据反射光线和反射点求解各反射镜的有效几何参数和单目立体视觉系统的成像参数。通过本发明的三维光路分析方法，可以设
计最佳的反射镜形状，更好地描述单目立体视觉系统成像参数，使得搭建的单目立体视觉系统结构更加紧凑、更能满足测量系统成像要求。
6.一种单目立体视觉系统三维光路分析方法，步骤如下：
7.(1)单目立体视觉系统搭建
8.单目立体视觉系统由反射单元和一个真实相机1组成，反射单元置于真实相机1的前方，由两对、共四面呈对称布置的高精度反射镜组成。借助反射镜对光路的多次反射，真实相机1的靶面被分成对称的右靶面2和左靶面3，右靶面2和左靶面3分别对应的目标成像的两个视角，即虚拟右相机4与虚拟左相机5，这样单目立体视觉系统就等效为双目立体视觉系统。
9.(2)坐标系建立与结构参数定义
10.以两个内镜8交线的中点为原点，建立世界坐标系o
‑
xyz7。世界坐标系o
‑
xyz7的z轴与真实相机1的光轴共线。以真实相机1的光心为原点建立相机坐标系o
c
‑
x
c
y
c
z
c
6，内镜8、外镜9呈对称分布，且均垂直于世界坐标系的xoz平面，令内镜8与xoy平面之间的夹角锐角为α，外镜9与xoy平面之间的夹角锐角为β，内镜8交点到相机光心之间的距离为d'，世界坐标系o
‑
xyz7原点到外镜9与x轴交点之间的距离为l，θ为真实相机1视场角的一半。
11.(3)多镜面反射模型
12.分析像方到物方的三维光路，令某一反射镜法向量为n＝(n
x
,n
y
,n
z
)
t
，入射光线的方向向量及其上一点分别为l和l0，那么入射光线或反射光线和反射镜的交点p表示为：
[0013][0014]
其中，d是从世界坐标系o
‑
xyz7原点到该反射镜的距离。从上式可得(l0+gl)
t
·
n+d＝0，整理得到经过s个反射镜反射后得到的矢量可表示为：
[0015][0016]
其中，n
k
为第k个反射镜的单位法向量，满足d
k
为世界坐标系o
‑
xyz 7原点到第k个反射镜的距离。
[0017]
令物方平面π＝o
c
x
c
z
c
，变量左上标中的符号“r”和“l”分别表示真实相机1的右靶面2和左靶面3上的成像视图。然后，对于虚拟右相机4，确定从真实相机1的右靶面2到视场10的空间光路。自真实相机1右靶面2上点出发并穿过o
c
点的四条光线，经过内镜8和外镜9的反射后，反射光线与π处相交于则经过单目立体视觉系统多镜面作用后的反射光线和反射点可表示为：
[0018][0019]
其中，为入射到第k个反射镜的第i条光线的方向向量；为光线上一点，为第i条入射光线与第k个反射镜的交点，亦为第(k+1)个反射镜第i条入射光线上的点，即
r
n
k
为第k个反射镜的单位法向量，单位法向量可以结构参数计算得到；为第k个反射镜的映射矩阵；为第(k+1)个反射镜的第i条入射光线的方向向量，即为第k个反射镜的第i条反射光线；为外镜9第i条反射光线上的一点；为外镜9的第i条反射光线与物方平面π的交点；n
π
为空间平面π的单位法向量；d
π
为世界坐标系原点到空间平面π的距离；为空间平面π第i条入射光线的方向向量。
[0020]
(4)单目立体视觉系统反射镜有效几何参数与成像参数计算
[0021]
第k个反射镜的第n条入射光线围成区域的几何属性为：
[0022][0023]
其中，s为入射光线在反射镜上围成区域的面积，c为入射光线在反射镜上围成区域的周长。
[0024]
四条光线与各反射镜交点围成等腰梯形，该等腰梯形即为反射镜的有效区域，等腰梯形各边分别为并且对于内镜8，该区域的有效几何参数可以通过以下表达式计算：
[0025][0026]
其中，mirrinl为内镜8的有效长度；ain为光线与内镜8相交形成的等腰梯形的面积；infh为光线与内镜8相交形成的等腰梯形的底边长；inbh为光线与内镜8相交形成的等腰梯形的顶边长；inc为光线与内镜8相交形成的等腰梯形的周长。
[0027]
同理，外镜9的有效几何参数可表示为：
[0028][0029]
其中，mirrexl外镜9的有效长度；aex光线与外镜9相交形成的等腰梯形所占的面积；exfh光线与外镜9相交形成的等腰梯形的底边长；exbh光线与外镜9相交形成的等腰梯形的顶边长；exc光线与外镜9相交形成的等腰梯形的周长。
[0030]
左靶面3和右靶面2到物方的三维光路对称，同理，根据上述步骤计算左靶面3即虚拟左相机5经过反射镜作用后的反射光线和反射点。
[0031]
对于单目立体视觉系统，两外镜9与反射光线即和的交点和共同围成了单目立体视觉系统的视场10，两外镜9反射光线的交点可表示为：
[0032][0033]
其中，“∩”表示两反射向量相交的符号，则单目立体视觉系统成像参数可表示为：
[0034][0035]
其中，dof为景深；fovl是在y
c
轴方向上的最大视场；bl为基线长度；fovh和fovhm分别为x
c
轴方向上的最小视场和最大视场；f为焦距；o
c
为真实相机1的光心；dis是最大视场处的物距。
[0036]
本发明的有益之处为提出了一种单目立体视觉系统三维光路分析方法，相比于其他二维光路分析方法，本发明在三维层面对空间光路进行分析，可以更加全面的求解反射镜的有效几何参数和单目立体视觉系统的成像参数，这些结果可为单目立体视觉系统的设
计提供更好的依据，从而使单目立体视觉系统结构更加紧凑、更能满足测量系统成像要求。
附图说明
[0037]
图1为单目立体视觉系统三维光路分析方法原理图。其中，1
‑
真实相机，2
‑
右靶面，3
‑
左靶面，4
‑
虚拟右相机，5
‑
虚拟左相机，6
‑
真实相机坐标系o
c
‑
x
c
y
c
z
c
，7
‑
世界坐标系o
‑
xyz，8
‑
内镜，9
‑
外镜，10
‑
视场。
[0038]
图2为单目立体视觉系统三维光路分析方法流程图。
具体实施方式
[0039]
以下结合技术方案和附图1、2详细叙述本发明的具体实施方式。附图1为单目立体视觉系统三维光路分析方法原理图。实施例选用实验室搭建的单目立体视觉系统，包括真实相机1、内镜8、外镜9、视场10。采用附图2所示测量流程求解反射镜的有效几何参数及单目立体视觉系统的成像参数。首先，利用四个反射镜和一个真实相机1构建单目立体视觉系统；其次，选择合适的原点建立世界坐标系并确定本发明所采用的单目立体视觉系统的结构参数；然后，建立单目立体视觉系统多镜面反射模型；最后，根据计算的反射光线和反射镜反射点求解反射镜的有效几何参数与单目立体视觉系统的成像参数。以下对具体实施进行详细说明。
[0040]
(1)搭建单目立体视觉系统
[0041]
使用四个反射镜(两个内镜8和两个外镜9)和一个真实相机1创建单目立体视觉系统，真实相机1的采集帧频为60帧，相机配备15mm的定焦镜头。
[0042]
(2)坐标系建立与结构参数定义
[0043]
在单目立体视觉系统上建立真实相机坐标系o
c
‑
x
c
y
c
z
c
6，在两内镜8交线的中心建立世界坐标系o
‑
xyz7。对结构参数进行赋值，令内镜8与xoy平面之间的夹角(锐角)α＝45
°
，外镜9与xoy平面之间的夹角(锐角)β＝60
°
，内镜8交点到相机光心之间的距离d'＝40mm；水平方向内镜8交点到(外镜9与x轴)交点之间的距离l＝100mm，相机视场角的一半θ＝10
°
。
[0044]
(3)建立多镜面反射模型
[0045]
以虚拟右相机(4)为例，分析从右靶面2到视场10方向上的光路。令右靶面2上以虚拟右相机(4)为例，分析从右靶面2到视场10方向上的光路。令右靶面2上同时令光心点o
c
＝(0,0,
‑
40)
t
。o
c
与连接构成的光线入射到内镜8上，根据结构参数计算得到内镜8的单位法向量n1＝(
‑
0.73,0,
‑
0.68)
t
，在此基础上，根据公式(1)、公式(2)和公式(3)可得到内镜8上等腰梯形的顶点坐标(3)可得到内镜8上等腰梯形的顶点坐标进一步计算得到内镜8的光线与外镜9相交等腰梯形的各个顶点和对于外镜9，由结构参数可计算其法向量n2＝(
‑
0.91,0,0.41)
t
，然后根据公式(1)、公式(2)和公式(3)计算得到由外镜9出发指向视场10方向的四条光线和公式(3)计算得到由外镜9出发指向视场10方向的四条光线和最后，计算可得到视场10上各个点的坐标，即
[0046]
(4)单目立体视觉系统反射镜几何参数与成像参数计算
[0047]
在步骤(3)的基础上进一步计算。四条光线与各反射镜交点围成等腰梯形，该等腰梯形即为反射镜的有效区域。对于内镜8，根据公式(4)和公式(5)可得到内镜8的几何参数为：内镜8的有效长度mirrinl＝50.43mm；光线与内镜8相交形成等腰梯形的面积ain＝1.56
×
103mm2；等腰梯形的底边长infh＝51.82mm；等腰梯形的顶边长inbh＝10.30mm；等腰梯形的周长inc＝181.28mm。
[0048]
同理，由公式(6)可计算外镜9的有效几何参数为：外镜9的有效长度mirrexl＝58.55mm；光线所占外镜9面积aex＝1.82
×
103mm2；光线与外镜9相交等腰梯形的底边长exfh＝55.38mm；光线与外镜9相交等腰梯形的顶边长exbh＝10.30mm；光线与外镜9相交等腰梯形的周长exc＝185.32mm。
[0049]
由公式(7)和公式(8)可求得单目立体视觉系统成像参数为：景深dof＝190.90mm、在y
c
轴方向上的最大视场fovl＝101.80mm、基线长度bl＝497.15mm、x
c
轴方向上的最小视场fovh＝10.30mm、x
c
轴方向上的最大视场fovhm＝51.82mm、最大视场处的物距dis＝443.99mm。
[0050]
本发明在三维层面对单目立体视觉系统光路进行分析和计算，从而得出单目立体视觉系统成像参数与结构参数间的关系，并计算出各反射镜的有效几何参数与单目立体视觉系统的成像参数。相较于现有二维光路分析方法，本发明可得到的参数更多，对整个单目立体视觉系统表征的更全面，能够为单目立体视觉系统设计提供更好的依据，使单目立体视觉系统结构更加紧凑、更能满足被测物的测量要求。