Выделение контуров на основе дескрипторов Фурье

% global bin_img
clear all; clc;
%наши изображения
Img1 = imread('triangle.bmp');
Img2 = imread('star.bmp');
Img3 = imread('ellipse.bmp');
Img4 = imread('triangle2.bmp');
Img5 = imread('ellipse.bmp');
%диалог выбора изображения
CellStr = {'Треугольник'; 'Звезда'; 'Круг'};
[picture] = listdlg('ListString', CellStr,...
'SelectionMode','single','Name', 'Окно выбора', 'ListSize', [150 50], 'PromptString', 'Выберите изображение', 'fus', 5);
%переключение между изображениями с целью выбора необходимого для работы
disp(['Вы выбрали: ' CellStr{picture}]);
switch picture
   case 1
        Img = Img1;
   case 2
        Img = Img2;
   case 3
        Img = Img3;
end
%Выделим одну цветовую составляющую
Img=Img(:,:,1);
Img=logical(Img);
[N M]=size(Img);
figure, imshow(Img);
title ('Исходное изображение');
 
%Определение начальной точки обработки
     for i1=1:N;
         for j1=1:M;
            if Img(i1,j1)==1;
               p_i=i1;p_j=j1;               
               break
            end;            
         end;
            if Img(i1,j1)==1;               
               break
            end; 
     end;      
% Координаты начальной точки
x(1)=p_i;
y(1)=p_j;
 
contour=bwtraceboundary(Img, [p_i p_j], 'W');
hold on
plot(contour(:,2),contour(:,1),'g','LineWidth',2);
 
%Формирование двумерного массива данных
s=[x',y'];
 
[np,nc]=size(s);
if nc~=2
    error('!!!!!!!!!')
end
if np/2~=round(np/2);
    s(end+1,:)=s(end,:);
    np=np+1;
end
 
% Центройд
 
Img=bwlabel(Img);
STATS=regionprops(Img,'centroid');
 
% Прямое преобразование Фурье
z=fft(STATS.Centroid,64); % вынести в качестве параметра алгоритма
 
fdesc = abs(z);
fdesc = fdesc(2:32)/fdesc(1);

STATS=regionprops(Img,'centroid');
 
for k=1:length(contour)
cx = STATS.Centroid(1);
cy = STATS.Centroid(2);
 
r=sqrt((cx-k)^2+(cy-k)^2);
 
end