ngjordanweiss2002.m

Normalized Cut Image Segmentation Code

%本程序为Ng,Jordan,Weiss(2002)为例进行的程序
clear;
close all;
load fisheriris;

W=SimilarMatrix(meas,1);
D=diag(sum(W));
L=D-W;
Lsym=(D^(-1/2))*L*(D^(-1/2));
[ev ed]=eig(Lsym);
d=diag(ed);
V=ev(:,1:3);
for i=1:150
    for j=1:3
        U(i,j)=V(i,j)/sqrt(sum(V(i,:).^2));
    end
end
idx=kmeans(U,3);
figure;plot(idx);
result=[length(find(idx==1)),length(find(idx==2)),length(find(idx==3))]