correlation / image

bonsoir,

je dispose d'une image et de sa copie translatée dans une direction donnée. je cherche en utilisant la fonction d'intercorrelation à déterminer ce vecteur de translation. quelqu'un saurait il comment faire? je ne vois pas la relation avec l'intercorrelation.

merci

Il me semble que tu cherches plus à comprendre la méthode qu'à l'implémenter en MATLAB (pour l'instant). Je déplace donc dans le bon forum ;)

Sinon, si tu as l'Image Processing Toolbox, tu peux jeter un oeil à ceci : Image registration

Bonjour,

Ton image de correlation devrait te donner un pic de correlation à l'endroit ou l'objet dans l'image s'est déplacé. En le localisant et en considérant d'autre part le centre de ton image d'origine, tu devrais obtenir ton vecteur de translation.

J'espère avoir pu t'aider...
a+

bonsoir,
merci pour votre aide.
j'ai donc implementé le script de test suivant :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
A=[ 1   2   3   0;
    4   5   6   0;
    7   8   9   0;
    0   0   0   0];
B=[ 0   0   0   0;
    1   2   3   0;
    4   5   6   0;
    7   8   9   0];
 
x=0;y=0;z=0;
 
C=zeros(6);
D=zeros(6);
 
C(2:5,2:5)=A;
D(2:5,2:5)=B;
 
for i=2:5
    for j=2:5
        x=x+(C(i,j)*D(i+1,j+1));
        y=y+(C(i,j)*C(i,j));
        z=z+(D(i,j)*D(i,j));
    end
end
 
t=x/sqrt(y*z);

j'obtiens donc bien un coef de 1 pour un vecteur de translation de (1,0)

Retour à Matlab avec deux questions:
y a-t-il moyen de se passer du "for dans ce code:

Code:

1
2
3
4
5
6
7
for i=2:5
    for j=2:5
        x=x+(C(i,j)*D(i+1,j+1));
        y=y+(C(i,j)*C(i,j));
        z=z+(D(i,j)*D(i,j));
    end
end

tout comme on s'est passé du for ici:

Code:

1
2
i=1:26866;
vect_bouc1(i)=vect_bouc1(i)-sum(vect_bouc1)/26866;

Citation:

Envoyé par einstein84

y a-t-il moyen de se passer du "for dans ce code:
Code:

1 2 x = C(2:end-1,2:end-1).*D(3:end,3:end); x = sum(x(:));
Idem pour les deux autres...

Sinon, as-tu résolu le problème de la compréhension de la méthode ?
Si oui, je déplacerai cette discussion dans le forum MATLAB ;)

merci pour le code. je le modifierais ulterieurement.

Citation:

Sinon, as-tu résolu le problème de la compréhension de la méthode ?
Si oui, je déplacerai cette discussion dans le forum MATLAB ;)

oui.
à present j'ai 7 images de 96x96 bits decalées les unes desautres de 1 bit, je dois trouver la trajectoire pour passer d'une image à une autre. j'ai fait ce code:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
% texture:192x192 | camera souris:12x12

% trajectoire: 0,0 ; 1,0 ; 2,1 ; 3,1 ; 3,2 ; 2,3 ; 2,4

% creation de la texture

texture=round(rand(192,192)*255);
% colormap(gray);
% image(texture);

i1=texture(1:96,1:96);
i2=texture(2:97,1:96);
i3=texture(3:98,2:97);
i4=texture(4:99,2:97);
i5=texture(4:99,3:98);
i6=texture(3:98,4:99);
i7=texture(3:98,5:100);

j1=zeros(98);   j1(2:97,2:97)=i1;
j2=zeros(98);   j2(2:97,2:97)=i2;
j3=zeros(98);   j3(2:97,2:97)=i3;
j4=zeros(98);   j4(2:97,2:97)=i4;
j5=zeros(98);   j5(2:97,2:97)=i5;
j6=zeros(98);   j6(2:97,2:97)=i6;
j7=zeros(98);   j7(2:97,2:97)=i7;

x1=0;y1=0;z1=0;
x2=0;y2=0;z2=0;
x3=0;y3=0;z3=0;
x4=0;y4=0;z4=0;
x5=0;y5=0;z5=0;
x6=0;y6=0;z6=0;

for i=-1:1
    for j=-1:1
        %%%%%%%%%%%%%%%
        for k=2:97
            for l=2:97
                x1=x1+(j1(k,l)*j2(k+i,l+j));
                y1=y1+(j1(k,l)*j1(k,l));
                z1=z1+(j2(k,l)*j2(k,l));
                
                x2=x2+(j2(k,l)*j3(k+i,l+j));
                y2=y2+(j2(k,l)*j2(k,l));
                z2=z2+(j3(k,l)*j3(k,l));
                
                x3=x3+(j3(k,l)*j4(k+i,l+j));
                y3=y3+(j3(k,l)*j3(k,l));
                z3=z3+(j4(k,l)*j4(k,l));
                
                x4=x4+(j4(k,l)*j5(k+i,l+j));
                y4=y4+(j4(k,l)*j4(k,l));
                z4=z4+(j5(k,l)*j5(k,l));
                
                x5=x5+(j5(k,l)*j6(k+i,l+j));
                y5=y5+(j5(k,l)*j5(k,l));
                z5=z5+(j6(k,l)*j6(k,l));
                
                x6=x6+(j6(k,l)*j7(k+i,l+j));
                y6=y6+(j6(k,l)*j6(k,l));
                z6=z6+(j7(k,l)*j7(k,l));
                
            end
        end

        t12(i+2,j+2)=x1/sqrt(y1*z1);
        t23(i+2,j+2)=x2/sqrt(y2*z2);
        t34(i+2,j+2)=x3/sqrt(y3*z3);
        t45(i+2,j+2)=x4/sqrt(y4*z4);
        t56(i+2,j+2)=x5/sqrt(y5*z5);
        t67(i+2,j+2)=x6/sqrt(y6*z6);
        %%%%%%%%%%%%%%%
    end
end

t12
t23
t34
t45
t56
t67

%v1=[2-position x du max de t12 dans t12 , 2-position y du max de t12] dans t12

à priori le resultat est correct. comment pourrais-je faire ce qui est ecrit en rouge?
merci

j'ai trouvé: find(A==x)

la seconde facon de faire:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
% texture:192x192 | camera souris:96x96
 
% trajectoire: 0,0 ; 1,0 ; 2,1 ; 3,1 ; 3,2 ; 2,3 ; 2,4
% **************************************************** %
 
% creation de la texture 
 
texture=round(rand(192,192)*255);
% colormap(gray);
% image(texture);
 
i1=texture(1:96,1:96);
i2=texture(2:97,1:96);
i3=texture(3:98,2:97);
i4=texture(4:99,2:97);
i5=texture(4:99,3:98);
i6=texture(3:98,4:99);
i7=texture(3:98,5:100);
 
 
j1=zeros(98);   j1(2:97,2:97)=i1;
j2=zeros(98);   j2(2:97,2:97)=i2;
j3=zeros(98);   j3(2:97,2:97)=i3;
j4=zeros(98);   j4(2:97,2:97)=i4;
j5=zeros(98);   j5(2:97,2:97)=i5;
j6=zeros(98);   j6(2:97,2:97)=i6;
j7=zeros(98);   j7(2:97,2:97)=i7;
 
 
x1=0;y1=0;z1=0;
x2=0;y2=0;z2=0;
x3=0;y3=0;z3=0;
x4=0;y4=0;z4=0;
x5=0;y5=0;z5=0;
x6=0;y6=0;z6=0;
 
 
for i=-1:1
    for j=-1:1
        %%%%%%%%%%%%%%%
 
        x1=j1(2:end-1,2:end-1).*j2(2+i:end-1+i,2+j:end-1+j);
        x1=sum(x1(:));
        y1=j1(2:end-1,2:end-1).*j1(2+i:end-1+i,2+j:end-1+j);
        y1=sum(y1(:));
        z1=j2(2:end-1,2:end-1).*j2(2+i:end-1+i,2+j:end-1+j);
        z1=sum(z1(:));
 
        x2=j2(2:end-1,2:end-1).*j3(2+i:end-1+i,2+j:end-1+j);
        x2=sum(x2(:));
        y2=j2(2:end-1,2:end-1).*j2(2+i:end-1+i,2+j:end-1+j);
        y2=sum(y2(:));
        z2=j3(2:end-1,2:end-1).*j3(2+i:end-1+i,2+j:end-1+j);
        z2=sum(z2(:));
 
        x3=j3(2:end-1,2:end-1).*j4(2+i:end-1+i,2+j:end-1+j);
        x3=sum(x3(:));
        y3=j3(2:end-1,2:end-1).*j3(2+i:end-1+i,2+j:end-1+j);
        y3=sum(y3(:));
        z3=j4(2:end-1,2:end-1).*j4(2+i:end-1+i,2+j:end-1+j);
        z3=sum(z3(:));
 
        x4=j4(2:end-1,2:end-1).*j5(2+i:end-1+i,2+j:end-1+j);
        x4=sum(x4(:));
        y4=j4(2:end-1,2:end-1).*j4(2+i:end-1+i,2+j:end-1+j);
        y4=sum(y4(:));
        z4=j5(2:end-1,2:end-1).*j5(2+i:end-1+i,2+j:end-1+j);
        z4=sum(z4(:));
 
        x5=j5(2:end-1,2:end-1).*j6(2+i:end-1+i,2+j:end-1+j);
        x5=sum(x5(:));
        y5=j5(2:end-1,2:end-1).*j5(2+i:end-1+i,2+j:end-1+j);
        y5=sum(y5(:));
        z5=j6(2:end-1,2:end-1).*j6(2+i:end-1+i,2+j:end-1+j);
        z5=sum(z5(:));
 
        x6=j6(2:end-1,2:end-1).*j7(2+i:end-1+i,2+j:end-1+j);
        x6=sum(x6(:));
        y6=j6(2:end-1,2:end-1).*j6(2+i:end-1+i,2+j:end-1+j);
        y6=sum(y6(:));
        z6=j7(2:end-1,2:end-1).*j7(2+i:end-1+i,2+j:end-1+j);
        z6=sum(z6(:));
 
        t12(i+2,j+2)=x1/sqrt(y1*z1);
        t23(i+2,j+2)=x2/sqrt(y2*z2);
        t34(i+2,j+2)=x3/sqrt(y3*z3);
        t45(i+2,j+2)=x4/sqrt(y4*z4);
        t56(i+2,j+2)=x5/sqrt(y5*z5);
        t67(i+2,j+2)=x6/sqrt(y6*z6);
        %%%%%%%%%%%%%%%
    end
end
 
 
[a12,b12]=find(t12==max(max(t12)));
[a23,b23]=find(t23==max(max(t23)));
[a34,b34]=find(t34==max(max(t34)));
[a45,b45]=find(t45==max(max(t45)));
[a56,b56]=find(t56==max(max(t56)));
[a67,b67]=find(t67==max(max(t67)));
 
 
v12=[2-a12,2-b12]
v23=[2-a23,2-b23]
v34=[2-a34,2-b34]
v45=[2-a45,2-b45]
v56=[2-a56,2-b56]
v67=[2-a67,2-b67]
 
% on reptrouve les points de la trajectoire citée au depart
 
pt=[0 0;v12;v12+v23;v12+v23+v34;v12+v23+v34+v45;v12+v23+v34+v45+v56;v12+v23+v34+v45+v56+v67]

la forme que tu m'as indiqué ne donne pas la meme valeur des x y z. et c'est autant plus bizarre qu'on trouve des valeurs de t superieures à 1

bonjour,
nouveau probleme plus delicat: imaginons que je compare deux images translatées et dont l'une a tourné d'un certain angle. comment on fait dans ce cas pour trouver le vecteur translation?
merci