equal.m

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

function [b]=equal(a,b,e)
if nargin<3
    e=1e-6;
end
if size(a)==size(b)
   for i=1:size(a,1)
       for j=1:size(a,2)
           if abs(a(i,j)/b(i,j)-1)<e
               b(i,j)=1;
           else
               b(i,j)=0;
           end
       end
   end
end