Discussion :

[moskaper]

Bonjour,
Je veux calculer le centre de masse de plusieurs forme dans une image j'ai tapé le code suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
xm=0;
for i=1:n
    for j=1:m
 
        xm=j*filter(i,j)+xm;
    end
end
xm=(1/(m*n))*xm;
 
 
ym=0;
for i=1:n
    for j=1:m
        ym=i*filter(i,j)+ym;
    end
end
ym=(1/(m*n))*ym;

Mais c'est juste pour une forme dans une image et j'ai utilisé le centroid mais ça marche que pour des formes fermées.
Merci.

[mihaispr]

Salut!

Je veux calculer le centre de masse de plusieurs forme dans une image

Peux-tu nous donner 3-4 examples des divers formes pour que tu veux calculer le centre de masse. Donc il faut que tu nous donne l'image pour que tu veux calculer le centre de masse pour qu'on a des chances de t'aider.


J'attends ton reponse!

Amicalment,

Michel

[mihaispr]

A c'est la matrice de ton image que tu lis avec imread.

Code :

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

1
2
 
>>doc imread

Une solution:

Code :

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

1
2
3
4
5
6
7
 
C=cellfun(@(n) 1:n, num2cell(size(A)),'uniformoutput',0);
[C{:}]=ndgrid(C{:});
C=cellfun(@(x) x(:), C,'uniformoutput',0);
C=[C{:}];
 
CenterOfMass=A(:).'*C/sum(A(:),'double'),

Aussi regarde la documenation pour ndgrid et cellfun:

Code :

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

1
2
3
4
 
>>doc ndgrid
>>doc cellfun
>>doc num2cell

Aussi:

Code :

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

1
2
3
 
>>doc bwlabel
>>doc regionprops

Si tu reussi postes ici la solution entiere ca va aidera les autres qui vont se confronter le meme probleme que toi!


Ou si tu as encore des problemes continue cette discussion!

J'espere que ca t'aides!

Amicalment,

Michel

[Jerome Briot]

Code :

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

1
2
3
4
5
 
C=cellfun(@(n) 1:n, num2cell(size(A)),'uniformoutput',0);
[C{:}]=ndgrid(C{:});
C=cellfun(@(x) x(:), C,'uniformoutput',0);
C=[C{:}];

... tout ça pour simplement générer les permutations des indices de chaque pixel ... on devrait pourvoir faire plus simple :

Code :

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

1
2
[r,c] = meshgrid(1:size(A,2),1:size(A,1));
C = [c(:) r(:)]

[moskaper]

merci michel et dut, mais j'y arrive toujours pas avec votre code je trouve la meme chose qu'avec le centroid je trouve qu' un point au centre de l'image.
J'aimerais trouver le centre de toutes les grosses tache blanche de l'image sur un fond noir.
Merci encore je vous envoie mon image ci-joind.

[Jerome Briot]

Code :

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

1
2
3
 
>>doc bwlabel
>>doc regionprops

[jarlaxle]

si tu veux analyser les différentes tâches, il faut d'abord les isoler,
bwlabel te donnera une valeur pour chacun de tes objets
tu peux effectuer des prétraitements à l'aide de morphologie mathématique pour éviter d'avoir trop d'objet.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
%BW image binaire, n connectivité : 4 ou 8
[L, num] = bwlabel(BW, n)
 
%L image labelisée, num nombre d'objets
 
%initialisation d'un tableau stockant tes centroïdes
centroid=zeros(num,2);
 
for indice=1:num
%pour chaque objet identifiée, on fait la moyenne des coordonnées trouvées %pour calculer le centroide.
centroid(indice,1:2)=mean(find(bwlabel==indice));
 
end

[mihaispr]

J'ai regarde ton code jarlaxle et a partir de ton code j'ai fais ca:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
>> img=imread('center_of_mass.jpg');imshow(img)
>> BW=im2bw(img);imshow(BW);
>> %BW image binaire, n connectivité : 4 ou 8
[L, num] = bwlabel(BW, n)
 
%L image labelisée, num nombre d'objets
 
%initialisation d'un tableau stockant tes centroïdes
centroid=zeros(num,2);
 
for indice=1:num
%pour chaque objet identifiée, on fait la moyenne des coordonnées trouvées %pour calculer le centroide.
centroid(indice,1:2)=mean(find(bwlabel==indice));
 
end

Mais ca ne marche toujours

Pourquoi?

Il me donne cette erreur:

Code :

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

1
2
 
??? Undefined function or variable 'n'.

How to define variable n?

n=6?

6 regions de l'image center_of_mass que j'ai sauve dans le dossier work?

Si je fais l'initialisation ca me donne toujour une erreur:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
n=6;img=imread('center_of_mass.jpg');imshow(img)
BW=im2bw(img);imshow(BW);
%BW image binaire, n connectivit? : 4 ou 8
[L, num] = bwlabel(BW, n)
%L image labelis?e, num nombre d'objets
%initialisation d'un tableau stockant tes centro?des
centroid=zeros(num,2);
for indice=1:num
%pour chaque objet identifi?e, on fait la moyenne des coordonn?es trouv?es %pour calculer le centroide.
centroid(indice,1:2)=mean(find(bwlabel==indice));
end

Erreur:

Code :

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

1
2
3
4
5
 
??? Function bwlabel1 expected its second input argument, CONN, to be either 4 or 8.
 
Error in ==> bwlabel at 65
[sr,er,sc,labels,i,j] = bwlabel1(BW,mode);

J'attends vos reponses s'il vous plait!

[jarlaxle]

n représente la connectivité, c'est à dire qu'on regarde les n plus proches voisins autour du point P.

on 2 connectivités possibles :
n=4 ,
0 1 0
1 P 1
0 1 0

n=8
1 1 1
1 P 1
1 1 1

donc n=6 te renvoie toujours une erreur, n vaut 4 ou 8 dans la fonction bwlabel.

[mihaispr]

Ah oui je me souviens maintenant j'ai oublie il y a des pixels tetra et octo-connectes. J'ai regarde un lien sur google.

Bon. Ca fait tres bien le calcul mais toujours une erreur

Code :

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

1
2
3
4
 
??? Error using ==> iptchecknargin
Function BWLABEL expected at least 1 input argument
but was called instead with 0 input arguments.

Error in ==> bwlabel at 47
iptchecknargin(1,2,nargin,mfilename);


J'ai remplace dans le code n=6 par n=4 et puis j'ai teste aussi avec n=8. Ca donne bien le resultat.

Mais une petite erreur encore persiste

Merci beacoup pour ton reponse!

[moskaper]

Merci jarlaxle et mihaispr pour votre aide.
Je trouve la meme erreur que toi je sait pas c'est a cause de quoi , je vais essayer de trouver la solution et si j' y arrive j'envoie la reponse .
encore merci .

??? Error using ==> iptchecknargin at 57
Function BWLABEL expected at least 1 input argument
but was called instead with 0 input arguments.

Error in ==> bwlabel at 47
iptchecknargin(1,2,nargin,mfilename);

Error in ==> detection_visage at 208
Centroid(indice,1:2)=mean(find(bwlabel==indice));

[Caro-Line]

A mon avis il doit juste s'agir d'une faute d'inattention.
BWLABL est une fonction, on ne peut l'utiliser comme variable.
D'ailleurs dans le code de jarlaxle elle est utilisée dès la 1ère ligne.
Par contre, n'y connaissant rien, je ne sais pas quelle variable il voulait réellement mettre. BW peut-être ?
Enfin vous vous devriez savoir normalement

[mihaispr]

OUi je remplacais comment tu as dit.


Il faisait confusion entre le nom de variable et le nom de la fonction.

Merci pour ton reponse caro! C'etait une faute de manque d'attention de ma part!

[jarlaxle]

désolé pour l'erreur, en fait ça sera l'image nommé "L" (labellisée) que tu vas utiliser dans cette méthode

[moskaper]

Merci a vous tous j'ai enfin compris .