Discussion :

[Pegasus42]

Bonjour,

J'ai segmenté, à partir de la méthode watershed, une image de cellules.

Pourriez-vous me dire les fonctions permettant de compter le nombre de parties segmentées que j'ai sur mon image afin de connaitre le nombre de cellules?

J'aimerai également trouver pour chacune d'elles leur aire:

Code :

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

STATS = regionprops(mon_image, 'area');

et classer ensuite ces données dans un tableau disant "j'ai tant de cellules de X pixels, tant de cellules de Y pixels etc.".

J'ai vu que l'on pouvait créer un tableau grâce à la fonction "ones".

Mes questions sont donc:
- comment compter le nombre de parties segmentés sur une image?
- la solution du tableau dynamique est-elle la meilleure? (mettre une nouvelle ligne à chaque fois que nous avons une nouvelle aire de cellules, et mettre +1 dans la case cellule dans la ligne appropriée si cette aire existe déjà). Dans ce cas quelles sont les fonctions que je pourrai utiliser?

J'espère que ce n'est pas trop embrouillé...

Merci,

Pegasus.

[Jean Dumoncel]

Bonjour,

- comment compter le nombre de parties segmentés sur une image?

C'est un algorithme d'étiquetage en composantes connexes que tu cherches. Comme tu parles de regionprops, je suppose que tu as accès à la toolbox traitement d'image, donc tu peux utiliser bwlabel :

Code :

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

help bwlabel

- la solution du tableau dynamique est-elle la meilleure? (mettre une nouvelle ligne à chaque fois que nous avons une nouvelle aire de cellules, et mettre +1 dans la case cellule dans la ligne appropriée si cette aire existe déjà). Dans ce cas quelles sont les fonctions que je pourrai utiliser?

Je ne suis pas sur d'avoir compris ta question, mais j'ai l'impression que tu cherches à faire un histogramme :

Code :

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

help hist

[Pegasus42]

Merci pour ta réponse! J'ai réussi à extraire l'information qui m'intéressait de chaque cellules:

Voici les lignes que j'écris (bw9 étant mon image segmentée):

Code :

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

1
2
3
4
5
6
7
 
bw9 = bwlabel(L,4);
stats = regionprops(bw9,'Area')
 
for n = 1:length(stats)
  disp(sprintf('Cellule numéro %d, Aire = %g',n,stats(n).Area))
end

Donc à chaque fois j'ai une cellule associée à une aire. J'aimerai pouvoir construire et afficher un tableau qui classe le nombre de cellules suivant leur aire.

Encore merci pour ton aide,

Pegasus

[Jerome Briot]

Essaie ceci :

Code :

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

[A,idx] = sort([stats.Area]);

[Pegasus42]

D'accord, la fonction que tu m'a fournie permet de trier mon tableau.
Mais si j'ai par exemple le tableau suivant:
[24 24 24 26 27 28 28 28]

comment pourrai-je créer le tableau suivant (j'ai 3x la valeur 24, 1x la valeur 26 etc.):

[3 1 1 3
24 26 27 28] ?

Cela nécessite-il une boucle?

Merci,

Pegasus

[Jean Dumoncel]

Par exemple :

Code :

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

1
2
A = [24 24 24 26 27 28 28 28];
resultat = [hist(A,unique(A));unique(A)];