validate_k.m

来自「摄像机圆球标定程序（包括SDP,DSDP,IDC,DIDC,ALG,DALG六种」· M 代码 · 共 19 行

function b=validate_K(K,C)

one=[];
b=0;
for i=1:3
    Q{i}=K'*C{i}*K;
    D_=eig(Q{i});
    if abs(D_(1)/D_(2)-1)<1e-6 || abs(D_(1)/D_(3)-1)<1e-6
        one(i)=D_(1);
    elseif abs(D_(2)/D_(3)-1)<1e-6
        one(i)=D_(2);
    else
        return;
    end
end
if abs(one(1)/one(2)-1)<1e-6&&abs(one(2)/one(3)-1)<1e-6
    b=1;
end