cameracalibrate_alg.asv

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

function [K r3 a]=CameraCalibrate_ALG(C)

s12=compute_sij(C{1},C{2});
s13=compute_sij(C{1},C{3});
count=0;
r3={};
a={};
K={};
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_;
        
         [w A pp]=ComputeW_D(C_,1,1);
%        vbar=ComputeV(C_);
%        [w A pp]=ComputeW(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=ComputeV(C_);
            vbar{count}=ComputeV(C_);
        end
    end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%recover from
[r3 a]=RecoverExtrinc(C,K_,vbar);