fndc.m

来自「这是一个关于简单matlab motino trakcing 算法的工具箱,十分」· M 代码 · 共 52 行

function [ey,ex,minval] = fndc(IL,IS,nop) %[Y-Pos, X-Pos, Min-Value] = fndc(Larger Matrix, Smaller Matrix, #points to be found)

FL1 = double(IL(:,:,1));
FS1 = double(IS(:,:,1));
FL2 = double(IL(:,:,2));
FS2 = double(IS(:,:,2));
FL3 = double(IL(:,:,3));
FS3 = double(IS(:,:,3));

sf1=size(FL1);
sf2=size(FS1);

yf1=sf1(1);xf1=sf1(2);
yf2=sf2(1);xf2=sf2(2);

dy = yf1-yf2;
dx = xf1-xf2;
for ey = 1:dy;
    for ex = 1:dx;
        FF1 = FL1(ey:yf2+ey-1,ex:xf2+ex-1);
        I1 = (FF1-FS1).^2;
        k1 = (sum(sum(I1)));
        FF2 = FL2(ey:yf2+ey-1,ex:xf2+ex-1);
        I2 = (FF2-FS2).^2;
        k2 = (sum(sum(I2)));
        FF3 = FL3(ey:yf2+ey-1,ex:xf2+ex-1);
        I3 = (FF3-FS3).^2;
        k3 = (sum(sum(I3)));
        k = k1+k2+k3;
        avge(ey,ex) = k/((yf2)*(xf2));
    end
end
clear ey;
clear ex;

% for i=1:nop
% minval(i) = min(min(avge));
% [ey(i),ex(i)]=find(avge == minval(i));
% avge(ey(i),ex(i)) = 1e9;
% end

minval = min(min(avge));
[eyy,exx]=find(avge == minval);

ex = exx(1);
ey = eyy(1);

%j = ey;
%i = ex;

end