Discussion :

[totoc1001]

Bonjour a tous.

Je viens de trouver un code qui permet de mesurer la valeur d'un déplacement en x et en y en faisant une corélation de 2 images.
Ce code fonctionne parfaitement, je l'ai testé, mais mon problème c'est que ce code n'est pas commenté, et j'ai énormément de mal a le comprendre.

Est ce que vous pourriez m'aider a le traduire en langage plus explicite? au moins les grandes ligne? Parce que la je suis paumé!!!

MErci beaucoup.

Voici mon code:

La fonction principale

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
 function [tx,ty] = register(imas1,imas2,t)
  corr = phase_correlation(imas1,imas2);
  if (max(max(abs(corr)))>t)
      max(max(corr));
      [k,l] = find(corr == max(max(corr)));
      tx = round(size(corr,1)/2-k+1);
      ty = round(size(corr,2)/2-l+1);
  else
      tx = 0;
      ty = 0;
  end

Cette fonction fait appel a la fonction phase_correlation dont voici le code ci-dessous

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
function ima = PhaseCorrelation(ima1,ima2)
 
s1 = size(ima1);
s2 = size(ima2);
 
s = max(s1,s2);
sx = s(1);
sy = s(2);
image1 = zeros(sx,sy);
image1(sx/2-s1(1)/2+1:sx/2+s1(1)/2,sy/2+1-s1(2)/2:sy/2+s1(2)/2) = ima1;
H1 = hamming(sx);
H2 = hamming(sy);
H = H1*H2';
image1 = image1.*H;
image2 = zeros(sx,sy);
image2(sx/2-s2(1)/2+1:sx/2+s2(1)/2,sy/2+1-s2(2)/2:sy/2+s2(2)/2) = ima2;
image2 = image2.*H;
f1 = fft2(fftshift(image1));
f1 = f1./abs(f1);
f2 = fft2(fftshift(image2));
f2 = f2./abs(f2);
corr = fftshift(ifft2(f1.*conj(f2)));
ima = abs(corr);

Encore Merci beaucoup
Thomas

[progfou]

Je ne commenterais pas le code, mais quelles sont les lignes qui te posent problème ?
As-tu trouver l'utilité de certaines ?

[totoc1001]

Bonjour

Apres avoir compilé et recompilé ce programme, j'ai pigé comment il fonctionne dans son ensemble, mais il ya a tjs des subtilités que je n'arrive pas a comprendre.

Par exemple, dans la premièr partie, que signifie:

Code :

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

1
2
3
 
 if (max(max(abs(corr)))>t)
      [k,l] = find(corr == max(max(corr)));

Quel est l'interet d'appeller 2 fois la valeur max?

et dans la seconde partie que signifient les lignes:

Code :

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

1
2
3
 
image1(sx/2-s1(1)/2+1:sx/2+s1(1)/2,sy/2+1-s1(2)/2:sy/2+s1(2)/2) = ima1;
image2(sx/2-s2(1)/2+1:sx/2+s2(1)/2,sy/2+1-s2(2)/2:sy/2+s2(2)/2) = ima2;

Je vais essayer de le commenter parfaitement d'ici la fin de ce week end et si poster ce que je pense sur le forum.
Si quelqu'un pouvais jetter un oeil pour etre sur que je n'ai pas fait de contre sens, ca serait vmt sympa.

Merci beaucoup en tout cas.
Thomas

[larimoise]

Salut,

Code :

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

1
2
 
if (max(max(abs(corr)))>t) [k,l] = find(corr == max(max(corr)));

max(max(A)) permet d'avoir le max d'une matrice A.

Code :

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

1
2
image1(sx/2-s1(1)/2+1:sx/2+s1(1)/2,sy/2+1-s1(2)/2:sy/2+s1(2)/2) = ima1; 
image2(sx/2-s2(1)/2+1:sx/2+s2(1)/2,sy/2+1-s2(2)/2:sy/2+s2(2)/2) = ima2;

ce code permet d'avoir une partie de la matrice.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
>> A =
 
     1     2     3
     4     5     6
     7     8     9
    10    11    12
 
>> A(1:4,2:3)
 
ans =
 
     2     3
     5     6
     8     9
    11    12

A+

[rostomus]

Bonjour,

la fonction max(A) calcule les maxima par colonne si A est une matrice 2D, donc elle envoie un vecteur ligne qui contient les max de chaque colonne, et pour determiner le max de A il faut encore calculer le max de ce vecteur.
mais aussi tu peux faire max(A( : )) au lieu de max(max(A))

image1(sx/2-s1(1)/2+1:sx/2+s1(1)/2,sy/2+1-s1(2)/2:sy/2+s1(2)/2) = ima1;
image2(sx/2-s2(1)/2+1:sx/2+s2(1)/2,sy/2+1-s2(2)/2:sy/2+s2(2)/2) = ima2;

c'est pour copier la matrice ima1 exactement au milieu de la matrice image1
idem pour ima2

[totoc1001]

OK, merci beaucoup.
Ca me rend bien service

A bientot
Thomas

[totoc1001]

Bonjour a tous,

Désolé de revenir sur ce sujet apres si longtemps, mais je n'ai aps eu trop le temps de bosser dessus depuis un petit moment.
Enfin bref, j'ai encore quelques questions sur ce code.
J'ai passé pas mal de temps a chercher a comprendre comment fonctionnait ce programme, et j'avoue etre un peu perdu.
JE comprend le code, mais je ne comprend pas trop le principe de la correlation de phase. J'ai trouvé des dizaines de documents associés a ce sujet, mais aucun appliqué a l'image processing et au dépalcement d'image. Je ne comprend pas comment fait ce programme pour obtenir un déplacement.

Est ce que quelqu'un pourrait éclairer ma lenterne?

Merci beaucoup d'avance.
Thomas

[hanane78]

ce que je connais c'est comment calculer la correlation entre les composantes dans un espaces quelconque:
cor(Xk,Xd)=cov(Xk,Xd)/(ecarttype(Xk)ecarttype(Xd))
cov(Xk,Xd) : matrice de co-variance

[progfou]

hanane78: Tu vas être le nouveau champion du forum toi...
Ce n'est même pas la question posée !

totoc1001: et "phase correlation" te donnes des pistes .

[totoc1001]

OK, lol.Merci du renseignement
Mais comme je le dis dans mon message, j'ai déja cherché sur google, et j'ai trouvé des dizaines de documents sur la corrélation de phase, mais ts expliqué par des formules mathématiques et jamais associé a un déplacement d'un pixel par rapport a un autre. C'est pourquoi j'ai demandé de l'aide, parce que les math et moi...
Je vais continuer a investiguer de ce coté la ;-)

Thx
Thomas