algorithmekmeans.m

来自「matlab programming for clustering pam , 」· M 代码 · 共 70 行

function [moy,classes]= AlgorithmeKMEANS(ima,k)


tic %Debut du timer

% Initialisation de l'image
ima=double(ima);
copy=ima;         % copie de l'image
ima=ima(:);       % convertir en vecteur
mi=min(ima);     
ima=ima-mi+1;    

s=length(ima);


m=max(ima)+1;
h=zeros(1,m);
hc=zeros(1,m);

for i=1:s
  if(ima(i)>0) h(ima(i))=h(ima(i))+1;end;
end
ind=find(h);
hl=length(ind);

% Choix arbitraire des centres

moy=(1:k)*m/(k+1);

% Debut du clustering

while(true)
  
  oldmoy=moy;
  % Classification courante  
 
  for i=1:hl
      c=abs(ind(i)-moy);
      cc=find(c==min(c));
      hc(ind(i))=cc(1);
  end
  
  %recalculer les moyennes 
  
  for i=1:k, 
      a=find(hc==i);
      moy(i)=sum(a.*h(a))/sum(h(a));
  end
  
  if(moy==oldmoy) break;end;  %S'il ya plus de changements on sort de la boucle
  
end

% attribution des num閞os de clusters pour chaque pixels (classes c'est la matrice en r閟ultat)
s=size(copy);
classes=zeros(s);
for i=1:s(1),
for j=1:s(2),
  c=abs(copy(i,j)-moy);
  a=find(c==min(c));  
  classes(i,j)=a(1);
end
end

moy=moy+mi-1;   % reconstruire les moyennes avec leurs ordres r閑ls

DureeTraitements=toc % Donner la dur閑 du traitement

imagesc (classes); figure(gcf)