compute_sij.asv

摄像机圆球标定程序（包括SDP,DSDP,IDC,DIDC,ALG,DALG六种算法的实现）

function [sij]=compute_sij(Ci,Cj)

vsij=eig(Ci,Cj);

for i=1:3
    sijt=vsij(i);
    Cit=Ci;
    Cjt=sijt*Cj;

    [Ui Si Vi]=svd(Cit);
    [Uj Sj Vj]=svd(Cjt);
    Tij=Si*Sj;
    if Tij(1,1)>0 &Tij(2,2)>0&T(3)
end