Discussion :

[hanane78]

Bonjour tout le monde


je veux implémenté l'algorithme de harris avec l’algorithme suivant :
Calculer du gradient (I) suivant x Ix;
Calculer du gradient (I) suivant y Iy;
Pour chaque pixel la matrice :

Code :

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

1
2
3
4
 
   |Ix²    Ix*Iy|
M=|              | 
   |Ix*Iy    Iy²|

R=Determinant (M)-k(Trace(M)²);
comment je peux determiner le gradient de chaque pixel de mon image suivant x et y

pour que je puisse determine la matrice M

Merci D'avence

[Jerome Briot]

Tu avais trouvé tout seul la réponse : Gradient d'une image

Quel est le problème avec GRADIENT ?
MATLAB renvoie-t-il un message d'erreur ? Si oui, lequel ?

[hanane78]

Merci pour votre reponse
j'ai fait ce test

Code :

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

1
2
3
4
5
 
 
I= imread('1a000.pgm');
 
[Fx,Fy]= gradient(I);

voila l'erreur qui me donner

??? Error using ==> rdivide
Integers can only be combined with integers of the same class, or scalar doubles.

Error in ==> gradient at 70
g(2:n-1, : ) = (f(3 : n, : )-f(1 : n-2, : ))./h( : ,ones(p,1));

PS: mon image est a niveau de gris

[Jerome Briot]

Essaie ceci :

Code :

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

[Fx,Fy]= gradient(double(I));

[hanane78]

oui tres bien ça marche ,

maintenant comment je peut determiner la matrice M

Code :

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

1
2
3
4
5
 
 
   |Fx²    Fx*Fy|
M=|              | 
   |Fx*Fy    Fy²|

ainsi que le trace de M

Trace(M)²

j'ai fait cette instruction je pense que sa marche ?

Code :

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

1
2
 
M=[Fx.^2 Fx.*Fy;Fy.*Fx Fy.^2];

Merci d'avence

[Jerome Briot]

Tu vas rencontrer des problèmes de dimension de matrice si tu fais ainsi car Fx et Fy sont des matrices pas des scalaires.

Je ne connais pas bien la méthode mais je pense qu'il faut que tu calcules M pour chaque pixel de l'image dans une (ou deux) boucle(s) FOR-END.

Si je me souviens bien, la trace est simplement la somme des termes diagonaux d'une matrice (à vérifier).

[Edit] C'est d'ailleurs ce qui est expliqué dans le lien qui as dû te servir de référence, je pense : point d'interet harris

[hanane78]

oui je dois calculer M pour chaque pixel

[happy_1177]

hanane vous pouvez me dire comment tu as fait pr tracer M.

[hanane78]

voila mon code pour determiner l'equation M

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
dx = [-1 0 1; -1 0 1; -1 0 1];  
dy = dx';
 
Ix = conv2(InImage, dx, 'same');      
Iy = conv2(InImage, dy, 'same');    
 
W = gausswin(Taillefenetre); 
Ix2 = imfilter(Ix.^2, W, 'symmetric');
Iy2 = imfilter(Iy.^2, W, 'symmetric');
Ixy = imfilter(Ix.*Iy, W, 'symmetric');
 
kh = 0.04;
M= (Ix2.*Iy2 - Ixy.^2) - kh*(Ix2 + Iy2).^2;

[happy_1177]

merci hanane , tu as 1e idé sur la proramation dynamique ds le domaine de traiment d'image