solution.m

来自「用MATLAB实现计算机视觉中三维重构所需的基本矩阵的运算」· M 代码 · 共 18 行

%求解Af=0的解f
%X1,X2是两幅图像上的对应点
% n为匹配点的对数
function f=solution(x1,x2,n)
for i=1:n
    A(i,:)=[x2(1,i)*x1(1,i) x2(1,i)*x1(2,i) x2(1,i) x2(2,i)*x1(1,i)...
            x2(2,i)*x1(2,i) x2(2,i) x1(1,i) x1(2,i) 1];
end
rankA=rank(A)

%对A进行SVD分解 V的第九列就是f的最小二乘解
[U D V]=svd(A);
f=V(:,9);
% 将f转换为3×3矩阵
f=reshape(f,[3 3]);
f=f';