Discussion :

[ju_cayenne]

Bon j'ai vraiment du mal à me repérer avec des matrices 3D 4D et plus.

J'ai une matrice multidim MAT= 2x2x453, cela correspond aux lon lat time.
Time étant des moyennes mensuelles je souhaiterai extraire des valeurs de ma 3D pour effectuer 4 nouvelles matrices 3D. Je souhaiterais diviser ma matrice saisonnièrement (DJF MAM JJA SON).

La 3D de MAT c'est des valeurs de température tous les mois (attention je commence la série par le mois de mars) donc j'ai M A M J J A S O N D J F M A M J J A .... S O N (453 mois)
j'aimerai 4 matrices (1 par saison) --> M1= 2x2x38 ; M2= 2x2x38 ; M3= 2x2x38 ; M4= 2x2x37 ;

Ainsi la 3D de M1 = M A M M A M ...
La 3D de M2 = J J A J J A ...
La 3D de M3 = S O N S O N ...
La 3D de M4 = D J F D J F ...

Merciii

[le fab]

salut

commences par créer des vecteurs d'index "saison" par rapport à une année commençant par le mois de mars

Code :

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

vPrintemps = [1 1 1 0 0 0 0 0 0 0 0 0];

que tu dupliques sur toute ta période

Code :

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

1
2
vPrintemps = repmat(vPrintemps , 1, ceil(size(M,3)/12));
vPrintemps(453+1:end) = [];

et ensuite tu selectionnes avec ce vecteur d'index

Code :

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

M1 = M(:,:,vPrintemps);

pour l'automatisation, je te conseilles d'adopter un tableau de cellules pour stocker tes matrices saisonnières (M{1} plutôt que M1)

sinon il doit y avoir aussi moyen de déformer ta matrice (reshape) afin de créer une dimension supplémentaires faisant apparaître les trimestres (cf ton post précédent)

Fabien

[ju_cayenne]

Lorsque j'effectue cette commande

Code :

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

>> M1 = OUTPUT_(:,:,vPrintemps);

j'ai l'erreur suivante: :

Code :

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

Subscript indices must either be real positive integers or logicals.

En ce qui concerne la fonction reshape lié à mon dernier post je me suis arraché les cheveux toute une journée. Je suis complètement perdue lorsque j'ai des matrices multidimensionnelles. Cela doit être possible mais lorsque je dois visualiser ce que je veux faire avec 3 ou 4 dimension je suis complètement perdue.

[le fab]

t'as quoi dans vPrintemps ?
normalement que des 0 et des 1
si oui essayes tjrs un vPrintemps = logical(vPrintemps ) ou un vPrintemps = find(vPrintemps )

[ju_cayenne]

C'est parfait ! Merci encore