coefficient de correlation entre des images bizarre

Version imprimable

bonjour,
j'ai 5 images que je dois comparer à 1.
j'affiche les images pour verifier qu'elles ne sont effectivement pas identique.
je fais mon code matlab pour calculer la norme de chaque vecteur, le produit scalaire de chaque vecteur par le premier.
par definition le coef de correlation est le cotion du produit scalaires par le produit des normes et là je trouve un coef identique pour toutes les images meme si leurs normes et leurs produits scalaires ont differents.
je n'arrive pas à tomber sur la couile:

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
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
%***********************************
%INIT
%***********************************
importfile signalbsb1.txt;
importfile signalbsb2.txt;
importfile signalbsb3.txt;
importfile signalbsb4.txt;
importfile signalbsba.txt;
importfile signalbsbb.txt;
 
signalbsb1=signalbsb1/100000;
signalbsb2=signalbsb2/100000;
signalbsb3=signalbsb3/100000;
signalbsb4=signalbsb4/100000;
signalbsba=signalbsba/100000;
signalbsbb=signalbsbb/100000;
 
norm_bord1=0;
norm_bord2=0;
norm_bord3=0;
norm_bord4=0;
norm_borda=0;
norm_bordb=0;
norm_bouc1=0;
norm_bouc2=0;
norm_bouc3=0;
norm_bouc4=0;
norm_bouca=0;
norm_boucb=0;
%***********************************
%***********************************
 
 
 
 
%***********************************
%SCRIPT
%***********************************
 
%***********************************
%MISE EN FORME DES DONNEES
%***********************************
for i=1:1152
    image1(i,:)=signalbsb1((i-1)*445+1:i*445);
    image2(i,:)=signalbsb2((i-1)*445+1:i*445);
    image3(i,:)=signalbsb3((i-1)*445+1:i*445);
    image4(i,:)=signalbsb4((i-1)*445+1:i*445);
    imagea(i,:)=signalbsba((i-1)*445+1:i*445);
    imageb(i,:)=signalbsbb((i-1)*445+1:i*445);
end
 
mat_bord1=image1(600:800,90:355);
mat_bord2=image2(600:800,90:355);
mat_bord3=image3(600:800,90:355);
mat_bord4=image4(600:800,90:355);
mat_borda=imagea(600:800,90:355);
mat_bordb=imageb(600:800,90:355);
 
mat_bouc1=image1(900:1000,90:355);
mat_bouc2=image2(900:1000,90:355);
mat_bouc3=image3(900:1000,90:355);
mat_bouc4=image4(900:1000,90:355);
mat_bouca=imagea(900:1000,90:355);
mat_boucb=imageb(900:1000,90:355);
 
for i=1:266
    vect_bord1((i-1)*201+1:i*201,1)=mat_bord1(:,i);
    vect_bord2((i-1)*201+1:i*201,1)=mat_bord2(:,i);
    vect_bord3((i-1)*201+1:i*201,1)=mat_bord3(:,i);
    vect_bord4((i-1)*201+1:i*201,1)=mat_bord4(:,i);
    vect_borda((i-1)*201+1:i*201,1)=mat_borda(:,i);
    vect_bordb((i-1)*201+1:i*201,1)=mat_bordb(:,i);
 
 
    vect_bouc1((i-1)*101+1:i*101,1)=mat_bouc1(:,i);
    vect_bouc2((i-1)*101+1:i*101,1)=mat_bouc2(:,i);
    vect_bouc3((i-1)*101+1:i*101,1)=mat_bouc3(:,i);
    vect_bouc4((i-1)*101+1:i*101,1)=mat_bouc4(:,i);
    vect_bouca((i-1)*101+1:i*101,1)=mat_bouca(:,i);
    vect_boucb((i-1)*101+1:i*101,1)=mat_boucb(:,i);
end
%***********************************
%***********************************
 
%***********************************
%NORME
%***********************************
for i=1:53466
    norm_bord1=norm_bord1+sqrt(vect_bord1(i).*vect_bord1(i));
    norm_bord2=norm_bord2+sqrt(vect_bord2(i).*vect_bord2(i));
    norm_bord3=norm_bord3+sqrt(vect_bord3(i).*vect_bord3(i));
    norm_bord4=norm_bord4+sqrt(vect_bord4(i).*vect_bord4(i));
    norm_borda=norm_borda+sqrt(vect_borda(i).*vect_borda(i));
    norm_bordb=norm_bordb+sqrt(vect_bordb(i).*vect_bordb(i));
end
 
for i=1:26866
    norm_bouc1=norm_bouc1+sqrt(vect_bouc1(i).*vect_bouc1(i));
    norm_bouc2=norm_bouc2+sqrt(vect_bouc2(i).*vect_bouc2(i));
    norm_bouc3=norm_bouc3+sqrt(vect_bouc3(i).*vect_bouc3(i));
    norm_bouc4=norm_bouc4+sqrt(vect_bouc4(i).*vect_bouc4(i));
    norm_bouca=norm_bouca+sqrt(vect_bouca(i).*vect_bouca(i));
    norm_boucb=norm_boucb+sqrt(vect_boucb(i).*vect_boucb(i));
end
%***********************************
%***********************************
 
%***********************************
%PRODUIT SCALAIRE
%***********************************
for i=1:53466
    ps_bord2=norm_bord2+(vect_bord1(i).*vect_bord2(i));
    ps_bord3=norm_bord3+(vect_bord1(i).*vect_bord3(i));
    ps_bord4=norm_bord4+(vect_bord1(i).*vect_bord4(i));
    ps_borda=norm_borda+(vect_bord1(i).*vect_borda(i));
    ps_bordb=norm_bordb+(vect_bord1(i).*vect_bordb(i));
end
 
for i=1:26866
    ps_bouc2=norm_bouc2+(vect_bouc1(i).*vect_bouc2(i));
    ps_bouc3=norm_bouc3+(vect_bouc1(i).*vect_bouc3(i));
    ps_bouc4=norm_bouc4+(vect_bouc1(i).*vect_bouc4(i));
    ps_bouca=norm_bouca+(vect_bouc1(i).*vect_bouca(i));
    ps_boucb=norm_boucb+(vect_bouc1(i).*vect_boucb(i));
end
%***********************************
%***********************************
 
%***********************************
%COEFFICIENT DE CORRELATION
%***********************************
% la bouteille 1 est suppoée bonne et donc prise comme reference %
cor_bord2=ps_bord2/(norm_bord2*norm_bord1);
cor_bord3=ps_bord3/(norm_bord3*norm_bord1);
cor_bord4=ps_bord4/(norm_bord4*norm_bord1);
cor_borda=ps_borda/(norm_borda*norm_bord1);
cor_bordb=ps_bordb/(norm_bordb*norm_bord1);
 
 
cor_bouc2=ps_bouc2/(norm_bouc2*norm_bouc1);
cor_bouc3=ps_bouc3/(norm_bouc3*norm_bouc1);
cor_bouc4=ps_bouc4/(norm_bouc4*norm_bouc1);
cor_bouca=ps_bouca/(norm_bouca*norm_bouc1);
cor_boucb=ps_boucb/(norm_boucb*norm_bouc1);
%***********************************
%***********************************
 
%***********************************
%AFFICHAGE
%***********************************
% colormap(gray);
% image(mat_bord1*256);
% image(mat_bord2*256);
% image(mat_bord3*256);
% image(mat_bord4*256);
% image(mat_borda*256);
% image(mat_bordb*256);
% 
% image(mat_bouc1*256);
% image(mat_bouc2*256);
% image(mat_bouc3*256);
% image(mat_bouc4*256);
% image(mat_bouca*256);
% image(mat_boucb*256);
%***********************************
%***********************************
 
%***********************************
%***********************************

merci

Il me semble que tu peux grandement simplifier ton code...
Code:

1 2 3 for i=1:1152 image1(i,:)=signalbsb1((i-1)*445+1:i*445); end
devient

Code:

image1 = reshape(signalbsb1,445,[]).';

et
Code:

1 2 3 for i=1:53466 norm_bord1=norm_bord1+sqrt(vect_bord1(i).*vect_bord1(i)); end
devient :
Code:

1 2 3 i=1:53466; norm_bord1=sum(sqrt(vect_bord1(i).*vect_bord1(i)));

en effet:

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
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
%***********************************
%INIT
%***********************************
importfile signalbsb1.txt;
importfile signalbsb2.txt;
importfile signalbsb3.txt;
importfile signalbsb4.txt;
importfile signalbsba.txt;
importfile signalbsbb.txt;
 
norm_bord1=0;
norm_bord2=0;
norm_bord3=0;
norm_bord4=0;
norm_borda=0;
norm_bordb=0;
norm_bouc1=0;
norm_bouc2=0;
norm_bouc3=0;
norm_bouc4=0;
norm_bouca=0;
norm_boucb=0;
%***********************************
%***********************************
 
 
 
 
%***********************************
%SCRIPT
%***********************************
 
%***********************************
%MISE EN FORME DES DONNEES
%***********************************
image1 = reshape(signalbsb1,445,[]).';
image2 = reshape(signalbsb2,445,[]).';
image3 = reshape(signalbsb3,445,[]).';
image4 = reshape(signalbsb4,445,[]).';
imagea = reshape(signalbsba,445,[]).';
imageb = reshape(signalbsbb,445,[]).';
 
mat_bord1=image1(600:800,90:355);
mat_bord2=image2(600:800,90:355);
mat_bord3=image3(600:800,90:355);
mat_bord4=image4(600:800,90:355);
mat_borda=imagea(600:800,90:355);
mat_bordb=imageb(600:800,90:355);
 
mat_bouc1=image1(900:1100,90:355);
mat_bouc2=image2(900:1100,90:355);
mat_bouc3=image3(900:1100,90:355);
mat_bouc4=image4(900:1100,90:355);
mat_bouca=imagea(900:1100,90:355);
mat_boucb=imageb(900:1100,90:355);
 
for i=1:266
    vect_bord1((i-1)*201+1:i*201,1)=mat_bord1(:,i);
    vect_bord2((i-1)*201+1:i*201,1)=mat_bord2(:,i);
    vect_bord3((i-1)*201+1:i*201,1)=mat_bord3(:,i);
    vect_bord4((i-1)*201+1:i*201,1)=mat_bord4(:,i);
    vect_borda((i-1)*201+1:i*201,1)=mat_borda(:,i);
    vect_bordb((i-1)*201+1:i*201,1)=mat_bordb(:,i);
 
 
    vect_bouc1((i-1)*201+1:i*201,1)=mat_bouc1(:,i);
    vect_bouc2((i-1)*201+1:i*201,1)=mat_bouc2(:,i);
    vect_bouc3((i-1)*201+1:i*201,1)=mat_bouc3(:,i);
    vect_bouc4((i-1)*201+1:i*201,1)=mat_bouc4(:,i);
    vect_bouca((i-1)*201+1:i*201,1)=mat_bouca(:,i);
    vect_boucb((i-1)*201+1:i*201,1)=mat_boucb(:,i);
end
%***********************************
%***********************************
 
%***********************************
%NORME
%***********************************
norm_bord1=sqrt(sum(vect_bord1.*vect_bord1));
norm_bord2=sqrt(sum(vect_bord2.*vect_bord2));
norm_bord3=sqrt(sum(vect_bord3.*vect_bord3));
norm_bord4=sqrt(sum(vect_bord4.*vect_bord4));
norm_borda=sqrt(sum(vect_borda.*vect_borda));
norm_bordb=sqrt(sum(vect_bordb.*vect_bordb));
 
norm_bouc1=sqrt(sum(vect_bouc1.*vect_bouc1));
norm_bouc2=sqrt(sum(vect_bouc2.*vect_bouc2));
norm_bouc3=sqrt(sum(vect_bouc3.*vect_bouc3));
norm_bouc4=sqrt(sum(vect_bouc4.*vect_bouc4));
norm_bouca=sqrt(sum(vect_bouca.*vect_bouca));
norm_boucb=sqrt(sum(vect_boucb.*vect_boucb));
%***********************************
%***********************************
 
%***********************************
%PRODUIT SCALAIRE
%***********************************
ps_bord2=norm_bord2+(vect_bord1.*vect_bord2);
ps_bord3=norm_bord3+(vect_bord1.*vect_bord3);
ps_bord4=norm_bord4+(vect_bord1.*vect_bord4);
ps_borda=norm_borda+(vect_bord1.*vect_borda);
ps_bordb=norm_bordb+(vect_bord1.*vect_bordb);
 
ps_bouc2=norm_bouc2+(vect_bouc1.*vect_bouc2);
ps_bouc3=norm_bouc3+(vect_bouc1.*vect_bouc3);
ps_bouc4=norm_bouc4+(vect_bouc1.*vect_bouc4);
ps_bouca=norm_bouca+(vect_bouc1.*vect_bouca);
ps_boucb=norm_boucb+(vect_bouc1.*vect_boucb);
%***********************************
%***********************************
 
%***********************************
%COEFFICIENT DE CORRELATION
%***********************************
% la bouteille 1 est suppoée bonne et donc prise comme reference %
cor_bord2=ps_bord2/(norm_bord2*norm_bord1);
cor_bord3=ps_bord3/(norm_bord3*norm_bord1);
cor_bord4=ps_bord4/(norm_bord4*norm_bord1);
cor_borda=ps_borda/(norm_borda*norm_bord1);
cor_bordb=ps_bordb/(norm_bordb*norm_bord1);
 
cor_bouc2=ps_bouc2/(norm_bouc2*norm_bouc1);
cor_bouc3=ps_bouc3/(norm_bouc3*norm_bouc1);
cor_bouc4=ps_bouc4/(norm_bouc4*norm_bouc1);
cor_bouca=ps_bouca/(norm_bouca*norm_bouc1);
cor_boucb=ps_boucb/(norm_boucb*norm_bouc1);
%***********************************
%***********************************
 
%***********************************
%AFFICHAGE
%***********************************
colormap(gray);
image(mat_bord1*256);
% image(mat_bord2*256);
% image(mat_bord3*256);
% image(mat_bord4*256);
% image(mat_borda*256);
% image(mat_bordb*256);
% 
% image(mat_bouc1*256);
% image(mat_bouc2*256);
% image(mat_bouc3*256);
% image(mat_bouc4*256);
% image(mat_bouca*256);
% image(mat_boucb*256);
%***********************************
%***********************************
 
%***********************************
%***********************************

l'execution est aussi plus rapide.
par contre mon probleme persiste.*
voici le code épuré ne gardant que la partie me posant probleme:

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
mat_bord1=image1(600:800,90:355);
mat_bord2=image2(600:800,90:355);
mat_bord3=image3(600:800,90:355);
mat_bord4=image4(600:800,90:355);
mat_borda=imagea(600:800,90:355);
mat_bordb=imageb(600:800,90:355);
 
for i=1:266
    vect_bord1((i-1)*201+1:i*201,1)=mat_bord1(:,i);
    vect_bord2((i-1)*201+1:i*201,1)=mat_bord2(:,i);
    vect_bord3((i-1)*201+1:i*201,1)=mat_bord3(:,i);
    vect_bord4((i-1)*201+1:i*201,1)=mat_bord4(:,i);
    vect_borda((i-1)*201+1:i*201,1)=mat_borda(:,i);
    vect_bordb((i-1)*201+1:i*201,1)=mat_bordb(:,i);
end
%***********************************
%***********************************
 
%***********************************
%NORME
%***********************************
norm_bord1=sqrt(sum(vect_bord1.*vect_bord1));
norm_bord2=sqrt(sum(vect_bord2.*vect_bord2));
norm_bord3=sqrt(sum(vect_bord3.*vect_bord3));
norm_bord4=sqrt(sum(vect_bord4.*vect_bord4));
norm_borda=sqrt(sum(vect_borda.*vect_borda));
norm_bordb=sqrt(sum(vect_bordb.*vect_bordb));
%***********************************
%***********************************
 
%***********************************
%PRODUIT SCALAIRE
%***********************************
ps_bord2=norm_bord2+(vect_bord1.*vect_bord2);
ps_bord3=norm_bord3+(vect_bord1.*vect_bord3);
ps_bord4=norm_bord4+(vect_bord1.*vect_bord4);
ps_borda=norm_borda+(vect_bord1.*vect_borda);
ps_bordb=norm_bordb+(vect_bord1.*vect_bordb);
%***********************************
%***********************************
 
%***********************************
%COEFFICIENT DE CORRELATION
%***********************************
% la bouteille 1 est suppoée bonne et donc prise comme reference %
cor_bord2=ps_bord2/(norm_bord2*norm_bord1);
cor_bord3=ps_bord3/(norm_bord3*norm_bord1);
cor_bord4=ps_bord4/(norm_bord4*norm_bord1);
cor_borda=ps_borda/(norm_borda*norm_bord1);
cor_bordb=ps_bordb/(norm_bordb*norm_bord1);
%***********************************
%***********************************
 
%***********************************
%AFFICHAGE
%***********************************
colormap(gray);
image(mat_bord1*256);
% image(mat_bord2*256);
% image(mat_bord3*256);
% image(mat_bord4*256);
% image(mat_borda*256);
% image(mat_bordb*256);
%***********************************
%***********************************
 
%***********************************
%***********************************

10/12/2008, 00h26
ABN84

Resolu:
le problème résidait dans le fait que j'exprimais le coefficient de corrélation directement sur les vecteurs sans les avoir centré auparavant.