Discussion :

[Shargat]

Bonjour,

Aujourd'hui j'essaie de constuire l'edge oriented histogram (EOH) d'une image. Pour cela je voudrais appliquer les formules de calcul de magnitude de gradient et de direction ci-dessous :

m( x,y ) = sqrt( Gx( x, y )² + Gy( x, y )² )
teta( x,y ) = arctan( Gy( x,y ) / Gx( x,y ) )

Il me faut donc les gradients horizontaux et verticaux Gx et Gy à chaque pixel p(x, y).
Et c'est là que je bloque. Je sais que la fonction Canny fait cela en interne et sort en output une matrice contenant les contours. Mais je ne sais pas comment lire cette matrice pour récupérer ces informations.
Comment suis-je censé m'y prendre ?

Merci d'avance et désolé si cette question vous semble absurde, mais je n'ai pas trouvé de doc sur comment récupérer ces infos.

[math_lab]

L'aide d'OpenCV laisse sous entendre que leur opérateur Canny ne fait quasiment qu'appeler l’opérateur Sobel. Passes directement par Sobel et tu auras tes Gx et Gy. Apres, pour retrouver Canny a partir de ça, il faut faire un seuillage sur la magnitude du gradient (si gradient fort, alors il y a un edge).

[Shargat]

Ah d'accord très bien. Je vais regarder du côté de Sobel dans ce cas-là.

Merci!

[Shargat]

Bonjour,

J'ai un souci pour accéder aux gradients récupérés via Sobel.
Voici l'appel à Sobel et la récupération des gradients

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
/// Generate grad_x and grad_y
    Mat grad_x, grad_y;
    Mat abs_grad_x, abs_grad_y;
 
    /// Gradient X
    Sobel( src_gray, grad_x, ddepth, 1, 0, 3, scale, delta, BORDER_DEFAULT );
    convertScaleAbs( grad_x, abs_grad_x );
 
 
    /// Gradient Y
    Sobel( src_gray, grad_y, ddepth, 0, 1, 3, scale, delta, BORDER_DEFAULT );
    convertScaleAbs( grad_y, abs_grad_y );

Le code suivant m'affiche de jolis gradients positifs et négatifs comme on les aime

Code :

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

cout << grad_x;

En revanchhe, si j'essaie d'y accéder pixel par pixel, les nombres affichés sont différents et tous compris entre 0 et 255 :

Code :

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

1
2
3
4
5
6
7
8
 
 for( int i = 0; i < grad_x.rows; i++ )
    {
        for ( int j = 0; j < grad_x.cols; j++ )
        {
            cout << "Gx(x, y) = " << (float)(grad_x.data[i*grad_x.step + j]) << endl;
        }
    }

Je sais que c'est un problème d'accès aux données et de typage. Mais je ne parviens pas à trouver le bon typage. J'ai aussi tenté via la méthode at(), mais quel type mettre dans le template at<T> ? Il ne s'agit pas là d'une matrice à plusieurs channels donc je suis un peu perturbé ( en tant normal, je travaille sur 3 channels de couleurs, donc j'utilise <Vec3b> )

[math_lab]

Le type est probablement cv::Vec<float,2> ou quelque chose comme ça. Sinon, si tu as 3 dimensions (hauteur, largeur, channels), il me semble que l’opérateur at a une surcharge qui prend 3 entiers en paramètre (et dans ce cas la, le type est simplement float). Mais dans ton cas, je pense que le cv::Vec est le plus "logique" vu que ton image contient bien des vecteurs 2D.

[Shargat]

Merci pour vos réponses.

J'ai tenté ces deux typages :

Code :

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

1
2
3
4
 
            cout << "Gx = " << (float)(((Vec<float, 2>)(grad_x.at< Vec<float, 2> >(i, j)))[0]) << endl;
            cout << "Gx = " << (float)(grad_x.at<float>(i, j)) << endl;
            cout << "Gx = " << (float)(((Vec2f)(grad_x.at<Vec2f>(i, j)))[0]) << endl;

et ça ne convient pas. J'obtiens même des "not a number" (NaN).
Pour info, voici les valeurs des attributs du headers de grad_x. A ce sujet, je suis confus sur la valeur de depth. Estc-e que ça signifie que la valeur est codée sur 3 bits ou sur 3 octets ( d'après ce que j'avais noté, cela correspond à 3 bits )

channels = 1
step = 800
depth = 3