Discussion :

[Potté75]

Bonjour à tous

voila, je travaille actuellement sur des matrices 4D qui ont les dimensions suivantes:

Code :

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

1
2
3
4
 
size(POP)
ans =
   300   200     5    12

donc dans cet exemple 300 et 200 correspondent aux longitudes et latitudes (on a donc une grille), "5" correspond a 5 catégories différentes de populations et "12" correspond aux 12 mois de l'année qui m'intéresse
et enfin les valeurs qui remplissent toute cette matrice sont des fractions de concentrations de population

or je souhaite étudier une 6eme catégorie, qui est les zones sans populations
donc en gros je veux faire : 1-somme(des concentrations des 5 autres catégories)
ce que je souhaiterais faire c'est donc faire la somme de toutes les catégories a chaque point de grille (longitude, latitude)

mais je n'arrive pas au résultat que je souhaiterais! en effet, je sais que mon résultat doit etre de dimension 300x200x1x12
puisque je n'aurais plus qu'une categorie (celle sans population) mais je veux garder le meme nombre de longitude, latitude et mois

mais quoi que je tente sur matlab, je n'arrive pas a ces dimensions

Code :

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

SANS_POP=1-nansum(POP(:,:,1;5,:));

ici je pensais faire la somme sur l'ensemble des categories, mais pourtant le resultat est de dimension 1x200x5x2
quoi que je fasse matlab fait aussi un calcul sur les longitude! je ne comprends pas pourquoi

Si l'un d'entre vous peut m'éclairer, je vous remercie beaucoup
Bonne journée

[Potté75]

le seul truc que j'ai trouvé et qui semble correspondre a ce que je souhaite, c'est de permuter la variable POP

Code :

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

1
2
3
4
5
POP_perm=permute(POP,[3 4 1 2]);
size(POP_perm)
ans =
 
     5    12   300   200

puis de refaire la ligne de code

Code :

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

1
2
3
4
5
SANS_POP=1-nansum(POP_perm(1:5,:,:,:));
size(SANS_POP)
ans =
 
     1    12   300   200

mais bon, je ne suis pas vraiment convaincue de ce resultat, je ne vois pas pourquoi la permutation permettrait d'arranger quoi que ce soit...
Des idées??
Merci a vous

[Jean Dumoncel]

Bonjour,

la réponse est dans la documentation de nansum :

Code :

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

y = nansum(X,DIM)

Il suffit de faire la somme selon la 3ème dimension :

Code :

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

SANS_POP = nansum(POP,3);

[Potté75]

youps.. Tellement simple vu comme ca!
Merci beaucoup