Discussion :

[RedShady]

bonjour,

Je cherche à récupérer des valeurs (ici émissions de CO dues aux feux de forêts) seulement pour qq lignes et colonne pour transformer une matrice 720x360x12 en une de 100x160x12

j'ai essayé ça (qui bien sûr ne fonctionne pas)

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
emis_afr=zeros(100,160,12);
 
for i=0:720
    if i>=320 || i<=420
        for j=0:360
            if j>=100 || j<=260
               emis_afr(i,j,1)=((M1(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,2)=((M2(i,j)*pixel(i,j)*10^6*(3600*24*28))/10^9);
               emis_afr(i,j,3)=((M3(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,4)=((M4(i,j)*pixel(i,j)*10^6*(3600*24*30))/10^9);
               emis_afr(i,j,5)=((M5(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,6)=((M6(i,j)*pixel(i,j)*10^6*(3600*24*30))/10^9);
               emis_afr(i,j,7)=((M7(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,8)=((M8(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,9)=((M9(i,j)*pixel(i,j)*10^6*(3600*24*30))/10^9);
               emis_afr(i,j,10)=((M10(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);
               emis_afr(i,j,11)=((M11(i,j)*pixel(i,j)*10^6*(3600*24*30))/10^9);
               emis_afr(i,j,12)=((M12(i,j)*pixel(i,j)*10^6*(3600*24*31))/10^9);   
            end
        end    
    end
end

ce que sont les autres matrices seraient long à expliquer mais le problème est seulement dans la sélection de lignes et colonnes.

merci d'avance.

[shaiHulud]

Bonjour,
Une indication: ta 1ere sous matrices en question s'écrit

Code :

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

M1(320:420,100:260)

La multiplication terme à terme de 2 sous-matrices s'écrit

Code :

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

M1(320:420,100:260) .* pixel(320:420,100:260)

[RedShady]

Merci ça fonctionne