基本矩阵FF的代数推导
l′:(P′C)×(P′P+x)=[e′]×P′P+xl′:(P′C)×(P′P+x)=[e′]×P′P+x
F=[e′]×P′P+F=[e′]×P′P+其他表示
P=K[I|0],P′=K′[R|t]P=K[I|0],P′=K′[R|t]
则,P+=[K−10T],C=[01]则,P+=[K−10T],C=[01]
e=P[−RTt1]=KRTt, e′=P′[01]=K′te=P[−RTt1]=KRTt, e′=P′[01]=K′t
F=[e′]×K′RK−1=K′−T[t]×RK−1=K′−TR[RTt]×K−1=K′−TRKT[e]×F=[e′]×K′RK−1=K′−T[t]×RK−1=K′−TR[RTt]×K−1=K′−TRKT[e]×因为[e′]×[e′]×的秩是2,射影FF的秩是2
基本矩阵的性质
(1)
(2) l′=Fx, l=FTx′l′=Fx, l=FTx′
(3) 对极点:Fe=0,FTe′=0Fe=0,FTe′=0
(4) FF是秩2,自由度是7的齐次矩阵对极线单应,是对极线,k是不过对极点e的任何直线
l′=F[k]×l, l=FT[k′]×l′l′=F[k]×l, l=FT[k′]×l′
由特殊运动产生的基本矩阵
纯平移,F=[e′]×K′RK−1=[e′]×F=[e′]×K′RK−1=[e′]×
x=(x,y,1)T, Zx=PX=K[I|0](X,Y,Z,1)T,因为等号右边的第三个数是Zx=(x,y,1)T, Zx=PX=K[I|0](X,Y,Z,1)T,因为等号右边的第三个数是Z
(X,Y,Z)T=ZK−1X, Zx′=P′X=K[I|t](X,Y,Z,1)T=Zx+Kt/Z(X,Y,Z)T=ZK−1X, Zx′=P′X=K[I|t](X,Y,Z,1)T=Zx+Kt/Z
可以看出,相同的tt,越大,变化越小。可参考由火车窗口向外看的情况一般运动,可以先旋转,在针对K不同可以增加变换,使得转换为上面的情况
P=K[I|0], P′=K′[R|t]P=K[I|0], P′=K′[R|t]
射影变换H=K′RK−1H=K′RK−1(不改变摄像机中心),将第一幅图P↦K′[R|0]P↦K′[R|0],之后有x′T[e′]xHx=0, F=x′T[e′]xHxx′T[e′]xHx=0, F=x′T[e′]xHx,得x′=K′RK−1x+K′t/Zx′=K′RK−1x+K′t/Z
恢复摄像机矩阵
HH表示3维射影变换的一个4×4矩阵,那么的基本矩阵是相同的。也就是说由基本矩阵来确定摄像机矩阵最好的结果也要相差一个右乘3D射影变换
非零矩阵FF是的基本矩阵⟺P′TFP⟺P′TFP是反对称矩阵
基本矩阵FF,可以选择,P′=[[e′]×F+e′vT|λe′]P′=[[e′]×F+e′vT|λe′]。其中,vv是任何3维矢量,是一个正标量
本质矩阵
- 归一化坐标:x^=K−1X=[R|t]Xx^=K−1X=[R|t]X
- P=[I|0], P′=[R|t], E=[t]×R=K′TFKP=[I|0], P′=[R|t], E=[t]×R=K′TFK
EE的自由度是5,3+3-1
是本质矩阵⟺⟺它的奇异值中有2个相等而第三个是0
已知本质矩阵E=Udiag(1,1,0)VTE=Udiag(1,1,0)VT和第一个摄像机矩阵P=[I|0]P=[I|0],那么P′P′有4中可能
P′=[UWVT|u3], [UWVT|−u3], [UWTVT|u3], [UWTVT|−u3]P′=[UWVT|u3], [UWVT|−u3], [UWTVT|u3], [UWTVT|−u3]
W=⎡⎣⎢010−100001⎤⎦⎥W=[0−10100001]
只要用一个重构点XX作测试,验证它是否两个摄像机前面就可以确定唯一一个
应用
求出,进而求出本质矩阵EE,进而可以求出。这样2个相机的相对位置和旋转就可知