equal.m

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

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