Discussion :

[Guitoudugte]

Bonjour à tous, fraichement inscrit depuis quelques minutes, je viens poser mon petit problème de matlab. Je précise que j'ai très peu de connaissance et que je galère depuis plusieurs jours. Je me tourne donc vers vous, en espérant trouver une réponse. merci par avance.

Voici une partie du programme:

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
name1='Mask.tif';
fido1=fopen(name1);
entete1=fread(fido1,490,'int8');
mask=fread(fido1,[280 280],'uint8');
mask=mask/255;
 
name='MI=300.add';
fido=fopen(name);
entete=fread(fido,2^16,'int8');
for i=1:20
   entetimage=fread(fido,20,'int8');
   toto=fread(fido,[2*280 280],'uint8');
   tata=toto(1:2:2*280,: )+256.*toto(2:2:2*280,: );
 
   tata=144.87-50/2+tata.*50/4096;
 
   for k=1:280
        for j=1:280 
 
         tata(k,j)=1279.531559/(log(1.767305265/tata(k,j)+8391.636646))-273;
 
        end
   end
end

Pour simplifier l'explication, j'ai pris des sequences d'image par caméra infrarouge (extension "ADD"), que je dois traiter pour obtenir les profils de température. Donc dans un premier temps je lis l'image mask et je l'enregistre dans une matrice. Ensuite je lis les images prises par la caméra que j'enregistre dans une matrice tata. Je transforme les valeurs de la matrice de NN en UI (unité de la caméra). Puis je finis par transformer les valeurs de la matrice tata de UI en T°.

Mon problème se pose sur la dernière ligne :

Code :

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

tata(k,j)=1279.531559/(log(1.767305265/tata(k,j)+8391.636646))-273;

j'obtiens cette erreur:

Code :

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

1
2
3
4
??? Attempted to access tata(1,31); index out of bounds because size(tata)=[280,30].
 
Error in ==> FilmPIR at 23
            tata(k,j)=1279.531559/(log(1.767305265/tata(k,j)+8391.636646))-273;

merci par avance à celui ou celle qui me débloquera.
bien à vous
guillaume

[Jerome Briot]

Le nombre de colonnes de la matrice "tata" est de 30, et tu veux accéder aux colonnes 1 à 280 dans la boucle sur "j", donc cela pose un problème évident d'indexage.

[Guitoudugte]

j'avais compris cela, car dans le workspace il me dis bien que ma matrice
fait <280*30>. Pourtant je pense qu'elle est bien definie non

Code :

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

tata=(256.*toto(1:2:280,:))+256.*toto(2:2:280,:);

car la matrice toto fait <560*280>

Code :

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

  toto=fread(fido,[2*280 280],'uint8');

[Jerome Briot]

Ajoutes la fonction WHOS juste avant la boucle sur "k" pour vérifier la taille des matrices à chaque itération sur "i".

[Guitoudugte]

Il fait 8 passe sur i avant de me sortir l'erreur

Code :

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

1
2
3
4
5
Name      Size                    Bytes  Class
 
  i         1x1                         8  double array
 
Grand total is 1 element using 8 bytes

j'ai ce truc qui s'affiche 8 fois. et i vaut bien 8 dans le workspace

[Jerome Briot]

WHOS devrait trouver plus de variables.
As-tu bien fait ceci :

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
name1='Mask.tif';
fido1=fopen(name1,'r');
entete1=fread(fido1,490,'int8');
mask=fread(fido1,[280 280],'uint8');
mask=mask/255;

name='MI=300.add';
fido=fopen(name,'r');
entete=fread(fido,2^16,'int8');
for i=1:20
   entetimage=fread(fido,20,'int8');
   toto=fread(fido,[2*280 280],'uint8');
   tata=toto(1:2:2*280,: )+256.*toto(2:2:2*280,: );
   
   tata=144.87-50/2+tata.*50/4096;
   whos
   for k=1:280
        for j=1:280 
    
         tata(k,j)=1279.531559/(log(1.767305265/tata(k,j)+8391.636646))-273;
        
        end
   end
end

[Jerome Briot]

Tu peux aussi faire ceci :

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
name1='Mask.tif';
fido1=fopen(name1,'r');
entete1=fread(fido1,490,'int8');
mask=fread(fido1,[280 280],'uint8');
mask=mask/255;

name='MI=300.add';
fido=fopen(name,'r');
entete=fread(fido,2^16,'int8');
for i=1:20
   entetimage=fread(fido,20,'int8');
   [toto,count]=fread(fido,[2*280 280],'uint8');
   if count~=(2*280*280)
      error('Moins de 2*280*280 octets ont été lu')
   end
   tata=toto(1:2:2*280,: )+256.*toto(2:2:2*280,: );
   
   tata=144.87-50/2+tata.*50/4096;

   for k=1:280
        for j=1:280 
    
         tata(k,j)=1279.531559/(log(1.767305265/tata(k,j)+8391.636646))-273;
        
        end
   end
end

[Guitoudugte]

hum, non désolé j'avais mis que Whos('i').

je viens de le faire uniquement avec Whos voici le resultat:

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
Name             Size                    Bytes  Class
 
  entete       65536x1                    524288  double array
  entete1        490x1                      3920  double array
  entetimage      20x1                       160  double array
  fido             1x1                         8  double array
  fido1            1x1                         8  double array
  i                1x1                         8  double array
  j                1x1                         8  double array
  k                1x1                         8  double array
  mask           280x280                  627200  double array
  name             1x10                       20  char array
  name1            1x8                        16  char array
  tata           280x30                    67200  double array
  totalimage       1x1                         8  double array
  toto           560x30                   134400  double array
 
Grand total is 169670 elements using 1357252 bytes

[Guitoudugte]

je viens d'essayer ta manip, apparement pas de surprise

Code :

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

1
2
3
4
5
??? Error using ==> FilmPIR
Moins de 2*280*280 octets on été lu
 
Error in ==> FilmPIR at 17
      error('Moins de 2*280*280 octets on été lu')

[Jerome Briot]

Le problème vient donc de FREAD. Es-tu sur que la fin du fichier n'est pas atteinte ? Es-tu sur pour la taille de l'en-tête ? 2^16 c'est un peu beaucoup non ?

Essaies avec ceci :

Code :

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

1
2
3
   if count~=(2*280*280)
      error(sprintf('Moins de 2*280*280 octets ont été lu à l'itération %d',i))
   end

[Guitoudugte]

aucune idée pour la taille de l'entête, car c'est un programme que j'ai recuperé (qui a priori donnait des résultats avant, mais plus maintenant), et j'en aurais vraiment besion.

voici le dernier relevé avec le dernier changement de 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
Name             Size                    Bytes  Class
 
  count            1x1                         8  double array
  entete       65536x1                    524288  double array
  entete1        490x1                      3920  double array
  entetimage      20x1                       160  double array
  fido             1x1                         8  double array
  fido1            1x1                         8  double array
  i                1x1                         8  double array
  j                1x1                         8  double array
  k                1x1                         8  double array
  mask           280x280                  627200  double array
  name             1x10                       20  char array
  name1            1x8                        16  char array
  tata           280x280                  627200  double array
  totalimage       1x1                         8  double array
  toto           560x280                 1254400  double array
 
Grand total is 379671 elements using 3037260 bytes
 
??? Error using ==> FilmPIR
Moins de 2*280*280 octets ont été lu à litération 8
 
Error in ==> FilmPIR at 17
      error(sprintf('Moins de 2*280*280 octets ont été lu à litération %d',i))

désolé de t'harceler comme ca, mais en tout cas merci pour ton aide

[Jerome Briot]

Si je calcul bien, la taille de ton fichier dois être (au minimum) de <entete+20*entetimage+20*toto> soit en octets :

Code :

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

1
2
3
4
2^16+20*20+20*2*280*280
ans =
 
     3201936

Pour connaitre la taille de ton fichier (tu peux aussi utiliser ton OS) :

Code :

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

1
2
3
4
fid=fopen('MI=300.add','r');
fseek(fid,0,'eof');
taille=ftell(fid);
fclose(fid);

Alors ?

[Guitoudugte]

taille = 1.1796e+006 octets

[Jerome Briot]

Donc ton fichier ne peut pas contenir toutes les données <entete+20*entetimage+20*toto> puisque sa taille est inférieur à 3201936 octets !

Désolé, mais la seule solution est de revoir ton format de lecture

[Guitoudugte]

oki je vais voir tout ca, en tout cas merci pour ton aide. @ bientot

[Jerome Briot]

Si tu es sur que ton fichier contient bien 20 groupes successifs de données, tu peux calculer la taille de l'entête comme ceci :

Code :

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

taille_entete=taille_fichier-20*20+20*2*280*280

Ensuite tu remplaces le 2^16 dans ton code par cette valeur :

Code :

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

entete=fread(fido,taille_entete,'int8');

Je reconnais que c'est de la bidouille