Discussion :

[j_ere_my]

Bonjour à tous,

voilà, j'ai commencé une thèse où je dois utiliser matlab pour traiter mes données.
Cependant, mes connaissances datant du deug, autant dire que je découvre matlab.

Je commence à m'y habituer, je réfléchis en matrice, vecteur ; mais là, j'ai pas bcp de temps pour réaliser une petit papier et je bute sur un truc tout con ...

g un fichier .txt avec ma matrice,

dont, une colonne pour la pression atmosphérique et une pour des concentrations chimique
cependant les pressions sont mélangées car j'ai filtré des données d'une année

j'aimerais tracer un profil vertical de la concentration moyennée pour chaque pression ...

si qq'un peut m'aider, merci bcp.

[Jerome Briot]

Attache le fichier txt (en zip) et dis nous comment les valeurs sont "mélangées" ?

Sinon, tu pourrais préciser ce que tu entends par "profil vertical" ?

[j_ere_my]

en fait, j'ai mes données trier par journée (en fait par vols d'avions ==> donc par journée).
Et donc, j'ai classé ces données sur un fichier .txt pour chaque saison, ce qui signifie que dans mon fichier txt j'ai mes données tiées par jours de vols, et comme pour chak vols, j'ai les memes pressions atmospherique, c'est tout mélangé.

Je sais pas si je m'exprime bien

En tout, cas, je cherche à tracer un plot(x,y) avec y en pression et x en concentration pour avoir un profil vertical de mes concentrations (du fait que la pression définit l'atitude). mais je veux un profil moyenné, cad, que pour chak pression, j'ai plusieurs données de concentrations (du à plusieurs vols) que je veux moyenné...

[ricololo]

Si je comprends bien, ta première colonne correspond aux pressions atmosphériques, la 2eme à la concentration d'une espèce (correspondant aux niveaux de pression) pour le vol 1, la 3eme à la concentration de cette espèce pour le vol 2,... et tu veux une concentration moyenné par niveaux.
C'est ça?
eh bien, tu n'as qu'à former une nouvelle matrice ou la premiere colonne serait tes niveuax de pression et la deuxième serait constituée de :

Code :

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

mean(tamatrice(:,2:size(tamatrice,2)),2)

[j_ere_my]

en fait, j'ai seulement 2 colonnes car tous les vols sont réunis à la suite sur une seule colonne. C'est pour ca que j'ai mes pressions mélangées ...

sinon je peux faire un scatplot avec toutes les valeurs de ma concentration en espece mais j'aurais préféré moyenné tout ca sur un plot normal, ce que ca je sais faire

[ricololo]

Tu ferais mieux d'envoyer un fichier, ce serait plus clair!

Tu as plusieurs blocs de 2 colonnes (chaque bloc correspondant à un vol) contenant dans la 1ere tes niveaux de pression et dans la 2ème tes concentrations. Quelques choses de ce genre (avec P la pression et C la concentration):

P(niveau1,vol1) C(niveau1,vol1)
P(niveau2,vol1) C(niveau2,vol1)
...
P(dernierNiveau,vol1) C(dernierNiveau,vol1)
P(niveau1,vol2) C(niveau1,vol2)
P(niveau2,vol2) C(niveau2,vol2)
...
P(dernierNiveau,vol2) C(dernierNiveau,vol2)
...

As tu le même nombre de niveau par vol?
Tu ne travaillerai pas par hasard sur des données du projet MOSAIC? (ou MOZAIC je ne sais plus)

[j_ere_my]

effectivement, je travaille sur MOZAIC

non je n'ai pas le meme nombre de niveau car, logiquement ils dépendent du vol, c'est pour cette raison que je bute d'ailleurs...

mais ce que tu as mis resume bien ma matrice

merci de m'aider

[ricololo]

Faudrait que tu repères le premier niveau de chaque vol pour pouvoir distinguer tes blocs.
Ensuite, je suppose que ta grille d'altitude est différente d'un vol à l'autre donc interpole tes concentrations sur une grille commune (fonction matlab polyfit et polyval). A toi de déterminer les limites physico chimiques de cette interpolation!
Tu as donc maintenant le même nombre de niveau pour tous tes vols. Ils ne te restent qu'à les moyenner!
Envoie moi ton fichier, je pourrai t'aider plus efficacement

[j_ere_my]

je vais essayer cette méthode, ca devrait aller ...

si ca "résiste", je reposterai. Mais avec cette fonction d'interpolation, que je connais maintenant, je devrais pouvoir y arriver.

Merci beaucoup.

[j_ere_my]

bon en fait, je me galere quand meme

et le fichier est bcp trop lourd

mais ce que tu décrivais précedemment resumait, en effet, ma matrice mais j'ai du mal à utiliser les fonctions polyfit et polyval sans exemples.

[Jerome Briot]

et le fichier est bcp trop lourd

Même sous forme d'archive zip ?

[j_ere_my]

ah si biensur


a decompresser avec gunzip (linux)

la pression est colonne 5
la concentration est colonne 7

[ricololo]

Salut,
merci pour ton fichier, c'est plus clair. Si je ne m'abuse, tu as dans ce fichier 4 vols:
1er vol: ligne 1 à 364
2eme vol: ligne 365 à 662
3eme vol: ligne 663 à 1038
4eme vol: ligne 1039 à 1355

Question:
tes pressions me paraissent bizarre. Pour le 1er vol, elles diminuent de 25029 à 19713 puis augmentent jusque 73976! c'est quoi l'unité?

pour faire ton interpolation, utilise la fonction matlab interp1. Ci-dessous un exemple d'utilisation en notant 'alt' la grille de ta mesure (ta colonne 5), 'conc' tes concentrations (ta colonne 7) et 'z' ta grille qui sera commune à tous tes vols

Code :

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

1
2
3
 
z=[0:15]' % grille d'altitude  commune en km 
y=interp1(alt,conc,z,'linear'); % pour une interpolation lineaire, il y en a d'autres

Le vecteur y obtenu correspon aux concentrations interpolées sur la grille z. Tu fais ça pour chaque vol donc tu as 4 vecteurs ''y'' ou une matrice (de dimension 16*4 pour les 16 niveaux de mon exemple et tes 4 vols) puis tu calcule les moyennes ligne par ligne et tu auras tes concentrations moyennes pour chaque niveau.

[j_ere_my]

merci pour ta réponse, je vais essayer au plus vite.

Quant aux pressions ...
si tu preferes, 73976 correspond à 739 hPa, donc en effet, elles augmentent quand l'avion commence son atterissage. Et quand elles diminuent au début, c'est tout simplement que l'avion prend de l'atitude pour arriver à son altitude de croisière entre 9 et 12 km.
De plus, si je n'arrive pas à "environ"1000hPa, cad, vers le sol, c'est dû au fait que j'ai ifltré mes données en ne récuperant que celles où les concentrations de mon espèce chimiques sont présentes, cad, quand l'instrument marchait...

voilà, pour ta question, sinon ben merci pour ton aide.

[j_ere_my]

petite remarque,

chaque vol ne va pas avoir les memes hauteurs de niveau interpolées puisque chaque vol n'a pas les meme niveaux quantitativement et qualitativement, donc je ne px pas faire une moyenne de mes vols sur chaque niveau interpolé.

j'ai pas encore essayé ta méthode, mais j'ai réfléchis à ca, et ca m'interpelle ...

[ricololo]

Je ne comprends pas ce que tu veux faire. Et je suis pas sûr que tu comprennes bien le sens de 'interpoler'

Tes grilles de pression pour tes 4 vols sont différentes. Exemple de grille de pression avec 2 vols pour simplifier:
vol 1 vol 2
25000 21600
24995 21590
... ...
19000 20000
Tes mesures de concentration ne sont donc pas faites (obligatoirement) aux mêmes niveaux.
Pour pouvoir moyenner les concentrations de ces 2 vols (ce que tu semble vouloir faire), il faut que les 2 grilles de pression soient les mêmes. C'est le but de l'opération d'interpolation.
Il faut juste choisir une grille de pression commune (et régulière, c'est mieux) pour nos 2 vols.
Avec mon exemple, elle pourrait être quelque chose comme [25000:-5:19000]. Et pour le vol 2 qui n'a pas de mesures entre 25000 et 21600, tu mets à ces niveaux une valeur interdite (indéfinie).

[Jerome Briot]

Je ne suis pas sûr de tout bien comprendre...

Pour tracer la moyenne des concentrations pour chaque valeur de pression :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
[p,c] = textread('noy_brut_200104.txt','%*f%*f%*f%*f%f%*f%f%*f%*f%*f');
 
up = unique(p);
mc = zeros(numel(up),1);
for n = 1:numel(up)
    mc(n) = mean(c(p==up(n)));
end
 
figure
plot(mc,up,'ro')

J'ai tout bon ?

[j_ere_my]

ok, je te remercie