Сегментация на основе анализа цвета

clear;

Считывание исходного изображения

L=imread('kids.bmp');
L=imread('im.jpg');
[N M s]=size(L);
L=double(L)./255;
figure, imshow(L);title('Исходное изображение');

Выбор пикселей лица с характерным цветом

[x,y,z]=impixel;
r=median(z(:,1)); std_r=std(z(:,1));
g=median(z(:,2)); std_g=std(z(:,2));
b=median(z(:,3)); std_b=std(z(:,3));

Сегментация на основе анализа цвета

for i=1:N;
disp(i);
for j=1:M;
if (abs(L(i,j,1)-r)<std_r)&(abs(l(i,j,2)-g)<std_g)&(abs(l(i,j,3)-b)L1=ismember(mitka,idx);
L1=1-L1;
figure,imshow(L1);title('После удаления шума (небольших объектов)');

STATS = regionprops(L1,'Eccentricity');
Lmitky = bwlabel(L1,4);
figure,imshow(Lmitky);
imwrite(Lmitky,'Изображение меток.jpg');
figure, imshow(L);
H=hist(Lmitky(,max(max(Lmitky)));
Nomer=find(H>4000);
for r=2:length(Nomer);
MIN_i=N;MIN_j=M;MAX_i=0;MAX_j=0;
for i=1:N;
for j=1:M;
if Lmitky(i,j)==Nomer(r);
if iMAX_i;
MAX_i=i;
end;
if j>MAX_j;
MAX_j=j;
end;
end;
end;
end;
line([MIN_j MAX_j],[MIN_i MIN_i]);hold on;
line([MAX_j MAX_j],[MIN_i MAX_i]);
line([MAX_j MIN_j],[MAX_i MAX_i]);
line([MIN_j MIN_j],[MIN_i MAX_i]);hold off;
end;

(Это Matlab)