find_a.m

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

function a=find_a(d)

a=[];
d=sort(d);
% if equal(d(1),d(2),1e-2)
if d(1)*d(2)>0
    a=d(3)*2/(d(1)+d(2));
% elseif equal(d(2),d(3),1e-2)
elseif d(2)*d(3)>0
    a=d(1)*2/(d(2)+d(3));
end
if size(a,1)==0
    disp('无法重建！');
end