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