对极几何和基本矩阵
1.对极几何
本质上,两幅视图之间的对极几何是图像平面与以基线(基线是连接两摄像机中心的直线)为轴的平面束的交的几何。
- 对极点是连接两摄像机中心的直线(基线)与像平面的交点。等价地,对极点是在一幅视图中另一个摄像机中心的像。它也是基线(平移)方向的消影点。
- 对极平面是一张包含基线的平面。存在着对极平面的一个单参数簇(束)。
- 对极线是对极平面与图像平面的交线。所有对极线相交于对极点。一张对极平面与左或右像平面相交于对极线,并定义了对极线之间的对应。
- 平行于图像平面的运动。对于平行于图像平面作平移以及旋转轴垂直图像平面的特殊运动情形,基线与图像平面的交点位于无穷远。因此,对极点是无穷远点而对极线是平行线。
2.基本矩阵
F
F
F
结论1: 基本矩阵
F
F
F可以记为
F
=
[
e
′
]
×
H
π
F = [e'] _\times H_\pi
F=[e′]×Hπ, 其中
H
π
=
P
′
P
+
H_\pi=P'P^+
Hπ=P′P+是从一幅图像到另一幅图像通过任意平面
π
π
π的转移映射。 而且, 因为
[
e
′
]
×
[e'] _\times
[e′]×秩为2 和
H
π
H_\pi
Hπ秩为3 ,所以
F
F
F是秩2 的矩阵。
假设摄像机矩阵是一个巳标定的双眼装置且世界原点在第一个摄像机上:
F
=
[
e
′
]
×
K
′
R
K
−
1
=
K
′
−
T
[
t
]
×
R
K
−
1
=
K
′
−
T
R
[
R
T
t
]
×
K
−
1
=
K
′
−
T
R
K
T
[
e
]
×
F=[e']_ \times K'RK^{-1}=K'^{-T}[t]_ \times RK^{-1}=K'^{-T}R[R^Tt]_ \times K^{-1}=K'^{-T}RK^T[e]_ \times
F=[e′]×K′RK−1=K′−T[t]×RK−1=K′−TR[RTt]×K−1=K′−TRKT[e]×
2.1对应条件
结论2: 对两幅图像中任何一对对应点
x
↔
x
′
x\leftrightarrow x'
x↔x′,基本矩阵都满足条件:
x
′
T
F
x
=
0
x'^TFx=0
x′TFx=0
F F F可以从两个摄像机矩阵 P P P、 P ′ P' P′计算出来,即 F F F可在相差一个整体尺度因子的意义下由摄像机唯一确定。
2.2基本矩阵的性质
定义1: 假设两幅图像由中心不重合的摄像机获得,则基本矩阵
F
F
F为对所有的对应点
x
↔
x
′
x\leftrightarrow x'
x↔x′, 都满足:
x
′
T
F
x
=
0
x'^TFx=0
x′TFx=0
的秩2的唯一的 3 × 3 3\times 3 3×3齐次矩阵。
- F F F是秩2 、自由度7的齐次矩阵。
- 点对应: 如果 x x x和 x ′ x' x′’是对应的图像点,那么 x ′ T F x = 0 x'^TFx=0 x′TFx=0
- 对极线:
- l ′ = F x l'=Fx l′=Fx是对应于 x x x的对极线。
- l = F T x ′ l=F^Tx' l=FTx′是对应于 x ′ x' x′的对极线。
- 对极点:
- F e = 0 Fe=0 Fe=0
- F T e ′ = 0 F^Te'=0 FTe′=0
- 由摄像机矩阵 P P P 、 P ′ P' P′进行计算:
- 一般摄像机: F = [ e ′ ] × P ′ P + F=[e']_\times P'P^+ F=[e′]×P′P+,其中 P + P^+ P+是 P P P的伪逆, e ′ = P ′ C 且 P C = 0 e'=P'C且PC=0 e′=P′C且PC=0
- 规范摄像机:
P
=
[
I
∣
0
]
P=[I|0]
P=[I∣0],
P
′
=
[
M
∣
m
]
P'=[M|m]
P′=[M∣m]
F = [ e ′ ] × M = M − T [ e ] × F=[e']_\times M=M^{-T}[e]_\times F=[e′]×M=M−T[e]×,其中 e ′ = m e'=m e′=m且 e = M − 1 m e=M^{-1}m e=M−1m - 非无穷远摄像机:
P
=
K
[
I
∣
0
]
P=K[I|0]
P=K[I∣0],
P
′
=
K
′
[
R
∣
t
]
P'=K'[R|t]
P′=K′[R∣t]
F = K ′ − T [ t ] × R K − 1 = [ K ′ t ] × K ′ R K − 1 = K ′ − T R K T [ K R T t ] × F=K'^{-T}[t]_ \times RK^{-1}=[K't]_ \times K'RK^{-1}=K'^{-T}RK^T[KR^Tt]_ \times F=K′−T[t]×RK−1=[K′t]×K′RK−1=K′−TRKT[KRTt]× - 转置:如果 F F F是摄像机对作 ( P , P ′ ) (P,P' ) (P,P′)的基本矩阵,则 F T F^T FT是其反序对 ( P ′ , P ) (P',P) (P′,P)的基本矩阵。
2.3对极线单应
结论3: 假设
l
l
l和
l
′
l '
l′是对应的对极线,且
k
k
k是不过对极点
e
e
e的任何直线, 则
l
l
l和
l
′
l '
l′间的关系是
l
′
=
F
[
k
]
×
l
l' =F[ k ]_\times l
l′=F[k]×l,对称地有
l
=
F
T
[
k
′
]
×
l
′
l=F^T[ k' ]_\times l'
l=FT[k′]×l′。
3.由特殊运动产生的基本矩阵
3.1纯平移
假设摄像机的运动是没有旋转并且内参数不变的纯平移,可以假设两个摄像机矩阵分别是
P
=
K
[
I
∣
0
]
P=K[I|0]
P=K[I∣0]和
P
′
=
K
[
I
∣
t
]
P'=K[I|t]
P′=K[I∣t],则:
F
=
[
e
′
]
×
K
K
−
1
=
[
e
′
]
×
F=[e']_\times KK^{-1}=[e']_\times
F=[e′]×KK−1=[e′]×
从像点
x
x
x到像点
x
′
x'
x′的映射是:
x
′
=
x
+
K
t
/
Z
x'=x+Kt/Z
x′=x+Kt/Z
其中 Z Z Z是点 X \mathbf{X} X的深度(在第一个摄像机的主轴上测量从 X \mathbf{X} X到摄像机中心的距离)。
- 注意在纯平移情形, F = [ e ′ ] × F=[e']_\times F=[e′]×是反对称的而且只有2个自由度,它们对应于对极点的位置。
- x x x, x ′ x' x′与 e = e ′ e = e' e=e′是共线的(假设把两幅图像叠起来) 。这种共线性质称为自对极,并且对于一般运动不成立。
3.2一般运动
给定两个任意摄像机,我们可以把第一幅图像的摄像机旋转使它与第二个摄像机平行。该旋转可以通过对第一幅图像施加一个射影变换来仿真。针对两幅图像的标定矩阵的差别还可对第一幅图像施加进一步的矫正。这两次矫正的结果等价于对第一幅图像作射影变换
H
H
H。假设这些矫正已经完成,那么这两个摄像机之间的有效关系是一个纯平移。
假设两个摄像机矩阵分别是
P
=
K
[
I
∣
0
]
P=K[I|0]
P=K[I∣0]和
P
′
=
K
′
[
R
∣
t
]
P'=K'[R|t]
P′=K′[R∣t],需要的射影变换是
H
=
K
′
R
K
−
1
=
H
∞
H=K'RK^{-1}=H_\infty
H=K′RK−1=H∞,其中
H
∞
H_\infty
H∞是无穷单应,且
F
=
[
e
′
]
×
H
∞
F=[e']_\times H_\infty
F=[e′]×H∞,则图像点
x
x
x到图像点
x
′
x'
x′的映射是:
x
′
=
K
′
R
K
−
1
x
+
K
′
t
/
Z
x'=K'RK^{-1}x+K't/Z
x′=K′RK−1x+K′t/Z
3.3纯平面运动
在纯平面运动时旋转轴与平移方向正交。正交性给这种运动增加了一个约束,如果 K ′ = K K '= K K′=K,则 F F F 的对称部分在这种平面运动情形时的秩为2。使它的自由度个数从一般运动的7个减少到纯平面运动的6 个。
4.基本矩阵的几何表示
基本矩阵的对称和反对称部分分别是
F
s
=
(
F
+
F
T
)
/
2
F
a
=
(
F
−
F
T
)
/
2
F_s=(F+F^T)/2 \;\;\;\;\;\;F_a=(F-F^T)/2
Fs=(F+FT)/2Fa=(F−FT)/2
满足: F = F s + F a F=F_s+F_a F=Fs+Fa
4.1对称部分
矩阵
F
s
F_s
Fs是对称的并且一般地秩为3 ,它有5 个自由度并等同于一条点二次曲线,称为Steiner 二次曲线。
4.2反对称部分
矩阵
F
a
F_a
Fa是反对称的并可以记为
F
a
=
[
x
a
]
×
F_a= [ x_a] _\times
Fa=[xa]×,其中
x
a
x_a
xa是
F
a
F_a
Fa的零矢量。反对称部分有2 个自由度且得同于点
x
a
x_a
xa。
4.3对极线对应
由一个单应关联的两个直线束,它们对应直线的交点的轨迹是一条二次曲线。这两个束都是对极线束, 一个过
e
e
e而另一个过
e
′
e'
e′。这些对极线由一个1D单应相关联。其交点的轨迹是二次曲线
F
s
F_s
Fs
4.4纯平面运动
平面运动情形,其中
F
s
F_s
Fs的秩为2。此时,Steiner 二次曲线是退化的,它等价于两条不重合直线。
F
s
=
l
h
l
s
T
+
l
s
l
h
T
F_s=l_hl_s^T+l_sl_h^T
Fs=lhlsT+lslhT
F
F
F可以记为:
F
=
[
e
′
]
×
[
l
s
]
×
[
e
]
×
F=[e']_\times [l_s]_\times [e]_\times
F=[e′]×[ls]×[e]×
两个对极点的分别有2 个自由度和直线 l s l_s ls有2 个自由度。
5.恢复摄像机矩阵
5.1射影不变性和规范摄像机
结论4: 如果
H
H
H是表示3维射影变换的一个4 X 4 矩阵,那么对应于摄像机矩阵
(
P
,
P
′
)
(P,P')
(P,P′)和
(
P
H
,
P
′
H
)
( PH , P'H)
(PH,P′H)的基本矩阵是相同的。
5.2摄像机矩阵的规范形式
结论5: 对应于摄像机矩阵对
P
=
[
I
∣
0
]
P =[I|0]
P=[I∣0]和
P
′
=
[
M
∣
m
]
P'= [M |m]
P′=[M∣m]的基本矩阵等于
[
m
]
×
M
[ m ]_\times M
[m]×M。
5.3给定
F
F
F后的摄像机射影多义性
定理1: 令
F
F
F为基本矩阵而
(
P
,
P
′
)
(P,P')
(P,P′)和
(
P
~
,
P
~
′
)
(\tilde{P},\tilde{P}')
(P~,P~′)都是与基本矩阵
F
F
F对应的两组摄像机矩阵对,则存在一个非奇异的4 X 4矩阵
H
H
H使得
P
~
=
P
H
\tilde{P} = PH
P~=PH和
P
~
′
=
P
′
H
\tilde{P}' = P'H
P~′=P′H。
引理1: 若秩2的矩阵
F
F
F可以分解为两个不同的形式:
F
=
[
a
]
×
A
F=[a]_\times A
F=[a]×A和
F
=
[
a
~
]
×
A
~
F=[\tilde{a} ]_\times \tilde{A}
F=[a~]×A~则对某非零常数
k
k
k和3维矢量
v
v
v有
a
~
=
k
a
\tilde{a}=ka
a~=ka和
A
~
=
k
−
1
(
A
+
a
v
T
)
\tilde{A}=k^{-1}(A+av^T)
A~=k−1(A+avT)。
5.4给定
F
F
F求规范摄像机对
结论6: 一个非零矩阵
F
F
F是对应于一对摄像机矩阵
P
P
P和
P
′
P'
P′的基本矩阵的充要条件是
P
′
T
F
P
P'^TFP
P′TFP是反对称矩阵。
结论7: 令
F
F
F是基本矩阵,
S
S
S是任意反对称矩阵。 定义摄像机矩阵对为:
P
=
[
I
∣
0
]
和
P
′
=
[
S
F
∣
e
′
]
P=[I|0]和P'=[SF|e']
P=[I∣0]和P′=[SF∣e′]
其中
e
′
e'
e′是满足
e
′
T
F
=
0
e'^T F = 0
e′TF=0的对极点,并假定所定义的
P
′
P '
P′是有效的摄像机矩阵(秩为3) ,则
F
F
F是对应子
(
P
,
P
′
)
(P , P' )
(P,P′)的基本矩阵。
结论8: 基本矩阵 F F F的摄像机矩阵可以选择为 P = [ I ∣ 0 ] P=[I|0] P=[I∣0]和 P ′ = [ [ e ′ ] × F ∣ e ′ ] P'=[[e']_ \times F|e'] P′=[[e′]×F∣e′]
结论9: 对应于基本矩阵
F
F
F,一对规范形式的摄像机矩阵的一般公式是:
P
=
[
I
∣
0
]
P
′
=
[
[
e
′
]
×
F
+
e
′
v
T
∣
λ
e
′
]
P=[I|0]\;\;\;\;\;P'=[[e']_ \times F+e'v^T|\lambda e']
P=[I∣0]P′=[[e′]×F+e′vT∣λe′]
其中 v v v是任何3维矢量, λ \lambda λ是一个正标量。
6.本质矩阵
本质矩阵是归一化图像坐标下的基本矩阵的特殊形式。基本矩阵可以看作是本质矩阵的推广,基本矩阵把要求标定摄像机这样一个(非本质的)假设去掉了。与基本矩阵相比较,本质矩阵有较少的自由度但增加了一些性质。
归一化坐标
- 如果知道标定矩阵 K K K,用它的逆矩阵作用于点 x x x得到点 x ^ = K − 1 x \hat{x}=K^{-1}x x^=K−1x。则 x ^ = [ R ∣ t ] X \hat{x}=[R|t]X x^=[R∣t]X,其中 x ^ \hat{x} x^是图像上的点在归一化坐标下的表示。
- 摄像机矩阵
K
−
1
P
=
[
R
∣
t
]
K^{-1}P=[R|t]
K−1P=[R∣t]称为归一化摄像机矩阵。
定理2: 用归一化图像坐标表示对应点 x ↔ x ′ x\leftrightarrow x' x↔x′时,本质矩阵的定义方程是:
x ^ ′ T E x ^ = 0 \hat{x}'^TE\hat{x}=0 x^′TEx^=0
基本矩阵和本质矩阵之间的关系是:
E
=
K
′
T
F
K
E=K'^TFK
E=K′TFK
6.1本质矩阵的性质
结论10: 一个3X3矩阵是本质矩阵的充要条件是它的奇异值中有两个相等而第三个是0。
6.2由本质矩阵恢复摄像机矩阵
结论11: 着
E
E
E的SVD分解为
U
d
i
a
g
(
1
,
1
,
0
)
V
T
Udiag (1, 1, 0)V^T
Udiag(1,1,0)VT,
E
=
S
R
E = SR
E=SR有如下两种可能的分解 :
S
=
U
Z
U
T
R
=
U
W
V
T
或
者
U
W
T
V
T
S=UZU^T\;\;\;\;\;\;R=UWV^T或者UW^TV^T
S=UZUTR=UWVT或者UWTVT
结论12: 己知本质矩阵
U
d
i
a
g
(
1
,
1
,
0
)
V
T
Udiag (1, 1, 0)V^T
Udiag(1,1,0)VT和第一个摄像机矩阵
P
=
[
I
∣
0
]
P = [ I|0 ]
P=[I∣0], 那么第二个摄像机矩阵
P
′
P'
P′有下列四种可能的选择:
P
′
=
[
U
W
V
T
∣
u
3
]
;
[
U
W
V
T
∣
−
u
3
]
;
[
U
W
T
V
T
∣
u
3
]
;
[
U
W
T
V
T
∣
−
u
3
]
;
P'=[UWV^T|u_3];\;\;\;[UWV^T|-u_3];\;\;\;[UW^TV^T|u_3];\;\;\;[UW^TV^T|-u_3];\;\;\;
P′=[UWVT∣u3];[UWVT∣−u3];[UWTVT∣u3];[UWTVT∣−u3];