Discussion :

[abel413]

Bonjour à tous,

Disons que je dispose d'une matrice de grande dimension qui est constituée de 0 et de groupements de nombres 1, 2, 3...par exemple avec 1 seulement :

Code :

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

A=[0 0 0 0 0 ; 0 0 1 1 0  ; 0 0 1 1 0  ; 0 0 0 0 0 ; 0 0 0 0 0]

J'aurais besoin de calculer les coordonnées des différents 1 et de retourner un vecteur colonne de ces coordonnées. Ici il est évident étant donnée la place des 1 que les coordonnées respectives attendues sont :

Code :

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

1
2
3
4
5
C=[2 3;
2 4;
3 3;
3 4;
]

J'imagine qu'il faut se servir de find mais find(A==1) ne donne que la position lue dans l'ordre d’apparition.

Une fois ces coordonnées obtenues j'aurais aussi besoin de prendre le barycentre. Comme je travaille avec des images binaires (les 1 font en fait référence aux pixels d'un domaine blanc de l'image), j'ai pu lire qu'il suffisait de faire la moyenne des coordonnées de C. Mais n'y a t-il pas une manière plus exacte de calculer un barycentre ?

[Jerome Briot]

Pour répondre à ta première question (même si je n'en vois pas l’intérêt) :

Code :

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

1
2
3
[i,j] = find(A.');
 
C = [j i]

Pour la seconde question, puisque tous les points ont la même valeur, il suffit de moyenner leurs coordonnées.

Si tu possèdes l'Image Processing Toolbox, tu peux le vérifier à l'aide de :

Code :

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

regionprops(A, 'centroid')

[abel413]

Effectivement, merci Dut c'est ce que je voulais faire, et j'ai également moyenné. En gros j'ai un fond noir et des cellules blanches pour lesquelles j'ai placé les centroïdes ainsi calculé.

C'est un peu hors sujet, mais j'ai utilisé la fonction regionprops dont tu parles comme suit :

Code :

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

STATS = regionprops(label_image, 'Eccentricity','MajorAxisLength','MinorAxisLength');

Pour chaque cellule j'obtiens les axes majeurs et mineurs qui sont apparemment des points mais j'ignore comment les représenter. J'ai entendu parler de quiver :

Code :

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

quiver(STATS(1,1).MinorAxisLength,STATS(1,1).MajorAxisLength)

Mais cela ne représente rien sur la première cellule par exemple.

[Jerome Briot]

Mais cela ne représente rien sur la première cellule par exemple.

As-tu au moins lu la documentation de la fonction quiver ?

[abel413]

Oui, j'ai même essayé avec les centroïdes de regionprops (mais j'obtiens les mêmes manuellement) :

Code :

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

hold on, quiver(STATS(1,1).Centroid(1),STATS(1,1).Centroid(2),STATS(1,1).MinorAxisLength,STATS(1,1).MajorAxisLength)

Ceci ne trace pas deux axes mais une sorte de vecteur partant du centre de la première cellule orienté vers le bas. Peut-être est-ce normal ? En fait je m'attendais à deux axes orthogonaux comme ceux d'une ellipse par exemple.

[Jerome Briot]

Il te manque une information :

Code :

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

STATS = regionprops(label_image, 'Eccentricity','MajorAxisLength','MinorAxisLength','Orientation');