compute_sij.m

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

function [sij]=compute_sij(Ci,Cj)
vsij=eig(Ci,Cj);
sij=[];
for i=1:3
    sijt=vsij(i);
    if imag(sijt)==0
        if det(Ci)*det(sijt*Cj)>0
            sij=[sij;sijt];
        end
    end
end
% if size(sij,1)~=1
%     sij
%     disp('出错！！！');
% end