fcmc.m

来自「Advanced Human Mood detection technique 」· M 代码 · 共 61 行

function dist=test(i)
%a=imageread();
a=rgb2gray(imread('iii.jpg'));
%a = imread('32.bmp');
a=imresize(a,[160,128]);
a=im2double(a);
[r c d]=size(a)
data=reshape(a,r*c,1);
[cen u]= fcm(data,5);
l=max(u);
for i=1:r*c
    if u(1,i)==l(i)
        p1(i)=255;
    else
        p1(i)=0;
    end
    if u(2,i)==l(i)
        p2(i)=255;
    else
        p2(i)=0;
    end
    if u(3,i)==l(i)
        p3(i)=255;
    else
        p3(i)=0;
    end
     if u(4,i)==l(i)
        p4(i)=255;
    else
        p4(i)=0;
     end
     if u(5,i)==l(i)
        p5(i)=255;
    else
        p5(i)=0;
     end
end

da =reshape(p1,r,c);
if da(1,1) ==255
    ret=da;
end

da2=reshape(p2,r,c);
if da2(1,1) ==255
    ret=da2;
end
da3=reshape(p3,r,c);
if da3(1,1) ==255
    ret=da3;
end
da4=reshape(p4,r,c);
if da4(1,1) ==255
    ret=da4;
end
da5=reshape(p5,r,c);
if da5(1,1) ==255
    ret=da5;
end
figure;imshow(ret);
dist=cut(ret);