Discussion :

[tasse2café]

Bonjour,

J'ai besoin d'un petit coup de main s'il vous plaît...
Je travail sur OpenCV depuis peu, et en me fixant sur une image, j'arrive à obtenir l'histogramme correspondant.
A partir de là je désire connaître la valeur moyenne et l'écart-type des pixels de l'image. (J'ai également un logiciel qui le fait automatiquement, ce qui me permet de f=vérifier mes résultats).

La boucle que j'utilise est la suivante.


for( int i = 0; i < hist_size; i++ )
{
mean += *cvGetHistValue_1D(hist,i) * i;
}
mean /= nb_lignes*nb_colonnes;

for( int i = 0; i < hist_size; i++ )
{
variance += pow((*cvGetHistValue_1D(hist,i) *i - mean),2);
}
variance /= (nb_lignes*nb_colonnes);


Si le résultat pour la moyenne est exact, ce n'est pas le cas pour la variance...
Est-ce que quelqu'un peut m'aider sur le je ne sais quoi que je ne vois pas...

D'avance merci,

[snowpy]

variance += pow((*cvGetHistValue_1D(hist,i) *i - mean),2);

à quoi te serve les deux parties en gras ?

[tasse2café]

à quoi te serve les deux parties en gras ?

Je cherche à calculer la variance :
- x(i) - moy(x) : (valeur de l'histogramme pour l'indice i) - (valeur du pixel moyen de l'image)
- je met cela au carré
- je fait cela pour i variant de 0 à 255 et je somme le tout
- je divise le résultat par le nombre d'éléments pris en considération, soit le nombre de pixels de mon image.

De la variance je déduis l'écart-type en calculant sa racine carrée.


Le "i" est donc l'indice de l'histogramme.
Je pondère la valeur de l'histogramme (*cvGetHistValue_1D(hist,i) = nombre de pixel à "i") par i.
D'où le " (*cvGetHistValue_1D(hist,i) "multiplié par" i.

Merci pour le temps consacré au passage.

[tasse2café]

Coucou, j'ai compris...
... et poste et clos le sujet au cas où quelqu'un passe et effectue la même erreur.


variance += *cvGetHistValue_1D(hist,i) * (pow( (i-mean),2) );

est à utiliser au lieu de

variance += pow( (*cvGetHistValue_1D(hist,i)*i - mean),2) ;


Ben, oui, ça paraît évident après coup...
Pour un nombre de pixels *cvGetHistValue_1D(hist,i) à la valeur i, il fallait faire *cvGetHistValue_1D(hist,i) fois (i - moy)^2 (pour accumuler les différences au carré),
et non ce que j'ai fait.

En espérant que ce soit assez clair... Merci pour la réponse en tout cas, ça m'a fait réfléchir sur la bonne piste.

[snowpy]

j'étais pas sur mais c'est ce que j'avais remarqué d'ou ma première question =)