i = double( imread('fruits2.tif') );
t = i(386:404, 623:673);

subplot(2,2,1); imshow(i,[])
subplot(2,2,2); imshow(t,[])

[xs, ys] = size(i);

tic
i = i - mean(i(:));
t = t - mean(t(:));
I = fft2(i);
T = fft2(t, xs, ys);
c = ifft2( I.*conj(T) );

subplot(2,2,3); imshow(c, [])
subplot(2,2,4); mesh(c)

m = max(c(:));
ind = find( floor(c/m)>0 );
x = mod(ind,xs)
y = ceil(ind/xs)

subplot(2,2,1);
hold on

plot(y,x, 'r*')
hold off
toc