Discussion :

[b4sols]

bonjour, je dois calculer l'erreur entre 2 images noir et blanc et j'ai fait la fonction suivante pour calculer mse et nmse :

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
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
function [MSE,NMSE] = proj18_Erreurs(picture1,picture2)
 
% affichage des deux images ‡ l'aide du toolbox Image
 
figure; subplot(1,2,1) ; imshow(picture1);
subplot(1,2,2); imshow(picture2);
 
%lire les images dans les variables
 
picture1 = imread(picture1);
picture2 = imread(picture2);
 
% calcul de la taille des images 
 
[Lignes1,Colonnes1] = size(picture1);
[Lignes2,Colonnes2] = size(picture2);
 
%calcul de la somme des diffÈrences
 
if size(picture1) == size(picture2);
 
    test = 0;
    test2 = 0;
 
    for i = 1 : Lignes1,
 
        for j = 1 : Colonnes1,
 
            test = test + double(picture1(i,j)^2);  
            test2 = test2 + double((picture1(i,j) - picture2(i,j))^2); 
 
        end
 
    end
 
%calcul des des 2 erreurs
 
MSE = test2/double(Lignes1*Colonnes1)
NMSE = test2/test
 
% Si les images ne sont pas comparable (dimension diffÈrentes), afficher un message d'erreur
 
else 'erreur images'
 
end

les calculs de mse et nmse vous semble-t-il correct???
merci de vos réponses

[Jerome Briot]

Que signifie MSE en NMSE ?

Sinon, Matlab est optimisé pour le calcul matriciel.
Les boucles FOR-END doivent être évitées (si possible) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
    test = 0;
    test2 = 0;
 
    for i = 1 : Lignes1,
 
        for j = 1 : Colonnes1,
 
            test = test + double(picture1(i,j)^2);  
            test2 = test2 + double((picture1(i,j) - picture2(i,j))^2); 
 
        end
 
    end

devient :

Code :

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

1
2
3
 
test=sum(double(picture1(:)).^2);
test2=sum((double(picture1(:))-double(picture2(:))).^2);

Note : DOUBLE peut être enlevé dans les versions récentes de Matlab.