Discussion :

[myriamulrik]

bonjour ,

quelqu'un saurait il implementer l'algorithme kmean en matlab ?

j'aivu qu'il existait une fonction kmeans dans matlab mais je n'arrive pas a l'utiliser pour segmenter une image.

http://www.mathworks.com/help/toolbox/stats/kmeans.html

en effet je pars d'une image 512x512.

comment puis-je utiliser la fonction kmeans alors pour en déduire mon image segmentée en k cluster et l'afficher ?


merci !

[Jean Dumoncel]

Bonjour,

un exemple de segmentation utilisant les kmeans si tu travailles sur les couleurs :
Color-Based Segmentation Using K-Means Clustering

[myriamulrik]

non je travaille sur des images de cerveau IRM en noir et blanc

[Jean Dumoncel]

plutôt en niveau de gris, non?

le principe reste le même. Qu'as-tu essayé de faire et ou bloques-tu?

[myriamulrik]

oui pardon en niveaux de gris

jai repris le code du lien que vous m'avez envoyé mais la redimension des resultats via la commande reshape ne semble pas fonctionner :


j'ai mon image X à l'origine de taille 512 x 512 en niveau de gris.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
nrows = size(X,1);
ncols = size(X,2);
X = reshape(X,nrows*ncols,2);
 
nColors = 3;
% repeat the clustering 3 times to avoid local minima
 
[cluster_idx cluster_center] = kmeans(ab,nColors,'distance','sqEuclidean', .
..
                                      'Replicates',3);
 
 
pixel_labels = reshape(cluster_idx,nrows,ncols);
imshow(pixel_labels,[]), title('image labeled by cluster index');

j'obtiens une erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
??? Error using ==> reshape
To RESHAPE the number of elements must not change.
 
Error in ==> script7 at 13
X = reshape(X,nrows*ncols,2);

pourriez vous m'aider ?
merci beaucoup

[Jean Dumoncel]

C'est normal puisque tu as une image en niveau de gris (donc une seule composante par pixel).
Change comme ceci :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

X = reshape(X,nrows*ncols,1);

[1wallla]

bonjour, j'ai presque la même problème, mon question est concernant le code donné par le modérateur: Mr. pour quoi le nombre de répétition de l'algorithme k-means est égal à 3? l'algorithme dit qu'on doit répéter les itérations jusqu’à "pas de changement dans les clusters". n'est ce pas?
et merci d'avance

[Invité]

Bonjour 1walla,

Tu as bien compris le principe, le paramètre 'replicate' comme l'indique la documentation :

'replicates'

Number of times to repeat the clustering, each with a new set of initial cluster centroid positions. kmeans returns the solution with the lowest value for sumd. You can supply 'replicates' implicitly by supplying a 3D array as the value for the 'start' parameter.

sert à exécuter plusieurs fois l'algorithme kmeans à partir de centroïds de départ différents. On garde à la fin la solution donnant les distances minimales (sortie sumd).