[kmeans] Seuillage

[math_infomatique]

salut

je veux appliquer k-means pour le seuillage

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13

I= imread('img.jpg');
X=I;
% k-means clustering
k=2; % nombre de groupes 
nb=3; % nbre de réitération
[m n]=size(X);
X=reshape(X,m*n,1);
[clusterIdx clusterCtr]=kmeans(X,k,'distance', 'sqEuclidean', 'Replicates', nb);
% l'affichage
clusterlabel=reshape(clusterIdx, m, n);
imshow(clusterlabel, []);

mais j'ai obtenu cette erreur

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
??? Error using ==> minus
Integers can only be combined with integers of the same
class, or scalar doubles.

Error in ==> kmeans>onlineUpdate at 576
            Del(:,i) = (m(i) ./ (m(i) + sgn)) .* sum((X -
            C(repmat(i,n,1),:)).^2, 2);

Error in ==> kmeans at 341
            converged = onlineUpdate();

Error in ==> tester1 at 12
[clusterIdx clusterCtr]=kmeans(X,k,'distance',
'sqEuclidean', 'Replicates', nb);

[Winjerome]

Bonjour,

Ta matrice obtenue avec IMREAD est de type UINT8, et c'est ce qui pose problème.
Il te faut la convertir en double au préalable.

[math_infomatique]

Envoyé par Winjerome

Bonjour,

Ta matrice obtenue avec IMREAD est de type UINT8, et c'est ce qui pose problème.
Il te faut la convertir en double au préalable.

ouii ça marche merci


mais j'ai remarqué un truc avec les méthodes prédéfinies de Matlab, elles ne donnent pas de bons résultats!!
j'ai voulu faire un seuillage automatique alors j'ai testé Otsu par graythresh et k-means par kmeans
le résultat est mauvais

peux tu m expliquer pourquoi ?

[magelan]

Envoyé par math_infomatique

mais j'ai remarqué un truc avec les méthodes prédéfinies de Matlab, elles ne donnent pas de bons résultats!!
j'ai voulu faire un seuillage automatique alors j'ai testé Otsu par graythresh et k-means par kmeans
le résultat est mauvais

C'est le problème de ce genre de méthode : il n'existe pas de méthode universelles. Le résultat dépend de tes images et de ce que tu cherches à l'intérieur. Sans cela, il sera difficile de t'indiquer une méthode adéquate...

[math_infomatique]

Envoyé par magelan

C'est le problème de ce genre de méthode : il n'existe pas de méthode universelles. Le résultat dépend de tes images et de ce que tu cherches à l'intérieur. Sans cela, il sera difficile de t'indiquer une méthode adéquate...

mmm je vois
alors peut être je dois les implémenter et voir si ça peut me donner un meilleur résultat !!

je vais commencer par écrire un code pour Otsu et un pour Kmeans

[math_infomatique]

voilà j'ai écrit le code de Otsu
mais toujours le même résultat que la méthode de graythresh

qu 'en pensez vous du code ??

sinon je suis entrain de tester les différentes méthodes possibles du seuillage et j'essaie de faire une comparaison entre elles !!!

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30

% *** Otsu sans méthode prédefinie 
I=rgb2gray(I);
h= imhist(I);
[n m]=size(I);
ntotal=n*m;
sum=0;
for i=1:256
    sum=sum+((i-1)*h(i));
end
sumB=0; wb=0; wf=0; varMax=0; threshold=0;
for t=1:256
    wb=wb+h(t);
    wf=ntotal-wb;
    sumB=sumB+((t-1)*h(t));
    mb=sumB/wb; 
    mf=(sum-sumB)/wf;  % moyenne
    var=wb*wf*(mb-mf)*(mb-mf);
    if var>varMax
        varMax=var; 
        threshold=t-1;
    end
end

%binarisation
I1 = zeros(size(I));
I1(I<threshold) = 0;
I1(I>threshold) = 255;
imshow(I1);

[magelan]

Envoyé par math_infomatique

voilà j'ai écrit le code de Otsu
mais toujours le même résultat que la méthode de graythresh

Si tu as codé exactement la même méthode, tu obtiendras exactement le même résultat puisque le seuil est calculé en fonction des valeurs de ton image...

On ne pourra pas t'aider beaucoup plus si tu ne nous expliques pas le contexte, ce que tu cherches et les résultats que tu as obtenus avec tes tests (avec des images par exemple).

[math_infomatique]

Envoyé par magelan

Si tu as codé exactement la même méthode, tu obtiendras exactement le même résultat puisque le seuil est calculé en fonction des valeurs de ton image...

On ne pourra pas t'aider beaucoup plus si tu ne nous expliques pas le contexte, ce que tu cherches et les résultats que tu as obtenus avec tes tests (avec des images par exemple).

oui t as raison !!
je m excuse

alors je travaille sur des images échographiques pour détecter la tumeur

voilà les images avant et apres application seuillage ostu

[bertry]

Salut,
Quelle différence entre ce que tu obtiens et ce que tu souhaite obtenir?
La tumeur est-elle la forme arrondie que l'on distingue en bas de l'image ( le truc qui semble pousser à l'intérieur de la cavité en noir ) ?

[math_infomatique]

Envoyé par bertry

Salut,
Quelle différence entre ce que tu obtiens et ce que tu souhaite obtenir?
La tumeur est-elle la forme arrondie que l'on distingue en bas de l'image ( le truc qui semble pousser à l'intérieur de la cavité en noir ) ?

oui c'est la forme arrondie qui représente la tumeur !!
oui peut être le résultat n'est pas aussi mauvais..
j'ai donné le mauvais exemple

bon voilà deux autres cas où j'ai appliqué les deux méthode Otsu et kmeans :

le premier pour un cancer du colon, le résultat n'est pas mauvais
mais pour le 2eme cas cancer du sein le résultat est mauvais !!

donc j'ai remarqué que ces méthodes ne donnent pas forcément un bon résultat..

[azertyuio]

j'ai rencontré le même problème.. le seuillage otsu ne donne pas un résultat pour mes images

regardes ici ça peut t'aider

[math_infomatique]

oui ça marche
mais c'est une autre méthode ?! entropique


donc conclusion
chaque image a ses propres spécifications et les méthodes de seuillages ne sont pas universelles pour toutes les images...