Discussion :

[wahidred]

salut tout le monde
je voudrais lire une image de type JPG avec la fonction imread (suivant la norme JPG) afin de calculer par la suite l'erreur quadratique car je fais de la compression j'ai vu le help de la fonction imread mais j'ai pas compris comment faudrait-il mettre les parametres le prof en tp m'a dit que lorsque vous lisez l'image c-a-dire que vs etes entrain de la decompresser or moi apres l'avoir lu j'ai constaté que l'erreur quadratique =0 ...bizarre d'autant plus que c'est une image de type JPG qui est repondue pour ses pertes et sa compression donc j'ai compris qu'il me manqait quelque chose dans "imread" je ne l'ai pas lu correctement .aidez moi svp merci

[Jean Dumoncel]

Bonjour,

Peux-tu nous montrer ton code qui te permet de lire l'image et de calculer l'erreur quadratique?

[wahidred]

salut ,voici mon code

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
x=multibandread ('D:\S1\compression image\TP_Master_2\aeropport_400_400',...
[400 400 6],'uint8=>uint8',...
   0,'bsq','ieee-le',{'band','direct',[1 2 3 4 5 6]} );
 
% EXTRACTION DE LA BANDE 6 qui va etre mon image de depart
bande6=x(:,:,6);
 
%ENREGISTREMENT DE LIMAGE SOUS FORMAT JPG
imwrite(bande6,'D:\S1\compression image\TP_Master_2\wahid.jpg','jpg')
 
% CALCUL DU TAUX DE COMPRESSION (c'est incomplet c'est le prof qui a dit %cela sans pour autant  savoir comment obtenir  ce taux de compressionn)
fid=fopen('D:\S1\compression image\TP_Master_2\aeropport_400_400','r');
status=fseek(fid,0,'eof');
position=ftell(fid)
 
%DECOMPRESSION LECTURE
jpg=imread('D:\S1\compression image\TP_Master_2\wahid.jpg','jpg');
 
 
% ERREUR QUADRATIQUE MOYENNE
for i=1:400
    for j=1:400
        diff=sum((bande6(i,j)-jpg(i,j)).^2);
    end
end
EQM=diff/(400*400)

EQM affiche 0 ...

[Jean Dumoncel]

Il y a un problème au niveau de ta boucle :

Code :

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

1
2
3
4
5
for i=1:400
    for j=1:400
        diff=sum((bande6(i,j)-jpg(i,j)).^2);
    end
end

Tu écrases diff à chaque itération, donc à la fin diff est égal à la différence au carré des 2 derniers pixels. De plus diff est le nom d'une fonction MATLAB, il faut éviter de l'utiliser pour nommer des variables.
Tu peux faire ceci :

Code :

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

1
2
3
4
5
6
diffe = 0;
for i=1:400
    for j=1:400
        diffe=diffe + (bande6(i,j)-jpg(i,j)).^2;
    end
end

ou directement :

Code :

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

diffe=sum(sum(bande6-jpg).^2));

[Jerome Briot]

ou directement :

Code :

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

diffe=sum(sum(bande6-jpg).^2));

Ou encore en économisant un SUM :

Code :

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

diffe = sum((bande6(:)-jpg(:)).^2);

[wahidred]

merci pour toutes vos reponses

[wahidred]

maintenant je veux savoir la taille de l'image apres l'avoir enregister sous le format jpg par exemple comment dois-je proceder a partir de matlab ?(sans aller au fichier) ,comme le code en haut le prof voulait qu'on utilise les fonction ftell et fseek sans savoir leur interet si vous avez une idée comment on pourait les utiliser pour avoir la taille dites le moi svp

[Jean Dumoncel]

Le format jpeg est un format de fichier binaire avec une syntaxe particulière. A l'intérieur de cet encodage, se situe l'information de la taille (largeur et hauteur) de l'image. Dans un premier temps, explore les octets de ton image en utilisant :

Code :

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

1
2
3
fid=fopen('image.jpg','r');
structure=fread(fid);
fclose(fid);

Chaque case de la variable structure représente un octet du fichier jpeg. En suivant la description du format jpeg, tu devrais pouvoir trouver ou se situe les informations de la taille de l'image.