calculer le barycentre d'une région avec plusieurs pixel de valeur différente

Bonjour à tous !

Je post ce message car j'ai besoin de votre aide.

Rentrons dans le vif du sujet : j'ai une région (qui s'assimile à une particule) qui contient plusieurs pixel. Chaque pixel a un système de coordonée i et j ainsi que la valeur de son intensité (en niveau de gris). Selon la région dans laquelle je me trouve, j'ai un nombre de pixel qui peut varier.

L'objectif de mon travail est de calculer le barycentre de chacune de ces régions en fonction, donc, de l'intensité des pixels qui les composent.

J'ai toutes les données nécéssaires pour calculer le barycentre d'une région c'est à dire : les coordonnées ainsi que les valeurs de chaque pixel. En fait j'ai fait le plus compliqué..

Le résultat que je cherche est donc sous la forme : i=fonction(i,j,v) j=fonction(i,j,v)

avec i, j et v des vecteurs (qui peuvent également avoir une taille 1x1 si ma région est composé d'uniquement 1 pixel) contenant les coordonées ainsi que les valeurs de chacun des pixel.

J'espere que vous pourrez m'aider et vous remercie d'avance. Amicalement.

[Guillaume]

Quick summary of the above post: The OP is asking how to calculate the centre of mass of a pixel region. He's already found the pixels location and intensity.

Answer 1

[Guillaume]

The barycentre is simply the mean of the point coordinates weighted by their intensity thus:

[i, j] = barycentre(i, j, v)
   %the (:) is just to make sure the inputs to mean are vectors with the same shape. it may not be needed
   i = sum(i(:) .* v(:)) / sum(v(:));  
   j = sum(j(:) .* v(:)) / sum(v(:));  
end

Le barycentre is simplement la moyenne des coordinates ponderees par l'intensite

Answer 2

[Image Analyst]

See my "Image Segmentation Tutorial" http://www.mathworks.com/matlabcentral/fileexchange/?term=authorid%3A31862

[Image Analyst]

You can ask regionprops() for an unweighted 'Centroid', or a centroid weighted by the gray levels with 'WeightedCentroid'.

See additional example attached. I do both in that example.

[Hidir ABAY]

Thanks a lot I will use this function to compare it with the others.

Answer 3

[Nut]

Bonjour,

excuse-moi pour mon francais, mais ils sont beaucoup d'ans que je ne parle cette langue.

Je pense que la fonction "regionprops" peux t'aider, il fait part de le Image Processing Toolbox. Observe le "Step 3":

http://it.mathworks.com/help/images/examples/measuring-regions-in-grayscale-images.html

[Hidir ABAY]

Bonjour, merci beaucoup cela me sera utile.Bonne journée

Answer 4

[Hidir ABAY]

I thank you for your answers and especially for your responsiveness =)

Have all a nice day.