gradient d'image (Prewitt et Sobel)

[etibat]

Bonsoir, suite à un exo en informatique, je devais codé les gradient de sobel et de prewitt mais je n'obtiens pas un bon résultat.

voila mes lignes de codes qui doivent poser problème

for (int a=1; a<hauteurImagePgm-1;++a)
for(int b=1;b<largeurImagePgm;++b)
{
tmp[a][b]=(-PixelsPgm[a-1][b-1]-k*PixelsPgm[a][b-1]-PixelsPgm[a+1][b-1]+PixelsPgm[a-1][b+1]+k*PixelsPgm[a][b+1]+PixelsPgm[a+1][b+1]);
tmp1[a][b]=(PixelsPgm[a-1][b-1]+k*PixelsPgm[a-1][b]+PixelsPgm[a-1][b+1]-PixelsPgm[a+1][b-1]-k*PixelsPgm[a+1][b]-PixelsPgm[a+1][b+1]);
tmp[a][b]=sqrt(tmp[a][b]*tmp[a][b]+tmp1[a][b]*tmp1[a][b]);
/*if (tmp[a][b]<0)
tmp[a][b]=-tmp[a][b];
if(tmp[a][b]>255)
tmp[a][b]=255;*/
}

tmp et tmp1 sont des tableaux temporaires, k peut prendre la valeur de 1 ou de 2 selon si c'est prewitt et sobel.
PixelsPgm sont des unsigned char tout comme les tmp

je ne sais vraiment pas pourquoi cela ne marche pas

merci de votre aide
-----

[invite4492c379]

Hello,

for (int a=1; a<hauteurImagePgm-1;++a)
  for(int b=1;b<largeurImagePgm;++b)
  {
    tmp[a][b]=(-PixelsPgm[a-1][b-1]-k*PixelsPgm[a][b-1]-PixelsPgm[a+1][b-1]+PixelsPgm[a-1][b+1]+k*PixelsPgm[a][b+1]+PixelsPgm[a+1][b+1]);
    tmp1[a][b]=(PixelsPgm[a-1][b-1]+k*PixelsPgm[a-1][b]+PixelsPgm[a-1][b+1]-PixelsPgm[a+1][b-1]-k*PixelsPgm[a+1][b]-PixelsPgm[a+1][b+1]);
    tmp[a][b]=sqrt(tmp[a][b]*tmp[a][b]+tmp1[a][b]*tmp1[a][b]);
    /*if (tmp[a][b]<0)
      tmp[a][b]=-tmp[a][b];
      if(tmp[a][b]>255)
      tmp[a][b]=255;*/
  }

Ça ne marche pas signifie-t-il :

* ça compile pas
* ça compile mais l'exécution est en erreur
* ça s'exécute mais le résultat est faux


Remarque sans savoir exactement ce qui se passe dans le morceau de code je me demande pourquoi on a b<largeurImagePgm et pas b<largeurImagePgm-1 ? (simple remarque visuelle).

[etibat]

le résultat obtenu n'est pas le bon, au lieu d'avoir une image plutôt noire avec les contours blanc, j'ai une image "poivre et sel"
en effet il y a une erreur sur j'ai bien oublié le -1 j'ai corrigé mais ça ne change rien.

Merci

[invite4492c379]

Pour autant que je puisse en juger le morceau de code à l'air correct, tu utilises bien les deux matrices


Pourquoi penses-tu que le problème vient de ce morceau de code ? As-tu essayé en mode debug pour voir les valeurs intermédiares ? As-tu généré les contours de ton image originelle avec un autre soft pour être certain du résultat ?

[A voir en vidéo sur Futura]

[etibat]

je ne sais pas trop utilisé le mode débug...

[invite4492c379]

Pourquoi penses-tu que le problème vient de ce morceau de code ?
(...)
As-tu généré les contours de ton image originelle avec un autre soft pour être certain du résultat ?

Et les autres questions ?

[etibat]

voila j'ai modifié une ligne pour :
tmp[a][b]=((tmp[a][b]^2)+(tmp1[a][b]^2))^(1/2);

et ceci a nettement améliorer mon image de sorti,
par contre je ne comprend pas la différence entre les deux. Si quelqu'un peut expliquer

[invite4492c379]

Tu programmes en quel langage ? (ce qui aurait dû être ma première question en fait ...)

[invite4492c379]

Pour info en C, il y a une fonction puissance qui est pow, ^ est un XOR bit à bit ...

[lou_ibmix_xi]

Fait attention au format des pixels, si ils sont encodés en "char" (de -128 à +127) ou "unsigned char" (de 0 à 255), tes calculs (même intermédiaire) dessus risque de dépasser la dynamique et donc de produire des resultats faux. Un moyen simple de palier à ce problème et de passer en "float" juste le temps des calculs, puis revenir au format originel avec éventuellement une normalisation pour recadrer la dynamique de sortie.