cameracalibrate_sdp.m

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

function [K r3 a]=CameraCalibrate_SDP(C)

s12=compute_sij(C{1},C{2});
s13=compute_sij(C{1},C{3});
count=0;
r3={};
a={};
K={};
vbar={};
for i=1: size(s12,1)
    for j=1:size(s13,1)
        s12_=s12(i);
        s13_=s13(j);
        C_{1}=C{1};
        C_{2}=C{2}*s12_;
        C_{3}=C{3}*s13_;
        vbar_=ComputeV(C_);
        [w]=ComputeW_SDP(C_,vbar_);
       
        if is_positive(w)
            K_=ComputeK_W(w);
%             if validate_K(K_,C_);
%                count=count+1;
%                K{count}=K_;
%             end
            count=count+1;
            K{count}=K_;
            vbar{count}=vbar_;
        end
    end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%recover from
[r3 a]=RecoverExtrinc(C,K,vbar);