G = rgb2gray(imread('motyle.png'));
GM = rgb2gray(imread('motyl3.png'));

% oznacenie samostatnych regionov
[L, num] = bwlabel(not(imbinarize(G)));
[LM, num] = bwlabel(not(imbinarize(GM)));

% struct s priznakmi
s = regionprops(L,'MajorAxisLength', 'MinorAxisLength');
sm = regionprops(LM,'MajorAxisLength', 'MinorAxisLength');

% priznaky
f = sm.MajorAxisLength / sm.MinorAxisLength;
fs = [s.MajorAxisLength]./[s.MinorAxisLength];

[whatever,idx] = min(abs(f-fs));

% nakreslenie a oznacenie
pos = c(idx).Centroid;
imshow(G);
hold on;
plot(pos(1),pos(2), 'r*')
hold off