Rose de pollution

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Manon Bouqueniau
Messages : 4
Enregistré le : 14 Mai 2018, 06:51

Rose de pollution

Messagepar Manon Bouqueniau » 29 Mai 2018, 14:19

Bonjour à tous !

J'ai un jeu de données journalières sur des concentrations de polluants dans les pluies. Je voudrais les associer à des données météo pour obtenir des roses de vent qui me permettraient de connaître l'origine géographique des polluants. Le problème est que mes données météo sont elles de fréquence horaires. mais je ne peux pas réaliser des moyennes pour mes données de direction du vent, car ils sont en degrés. Auriez-vous des idées quant à ce que je pourrais faire pour obtenir tout de même ces roses de vents ?

Merci d'avance !

Serge Rapenne
Messages : 1426
Enregistré le : 20 Aoû 2007, 15:17
Contact :

Re: Rose de pollution

Messagepar Serge Rapenne » 29 Mai 2018, 19:29

Bonjour,

Ta question est une question de méthodologie et pas une question sur R, elle est donc hors charte ici.

Toutefois,il n'est pas impossible de calculer des moyennes de vents, il suffit de transformer le vent sous forme d'angle (la direction), normes (la vitesse) en vecteur (u,v).
si tes données sont de la forme suivante
dta un data frame avec les variables suivantes
    date : date au format POSIXct
    ws : vitesse du vent m/s
    wd : direction du vent en degre

Code : Tout sélectionner

dta$u<-dta$ws*sin(2*pi*dta$wd/360)
dta$v<-dta$ws*cos(2*pi*dta$wd/360)

mean.u<-with(dta,tapply(u,format(date,"%Y-%m-%d),mean,na.rm=T))
mean.v<-with(dta,tapply(v,format(date,"%Y-%m-%d),mean,na.rm=T))

mean.vent<-merge(mean.u,mean.v,by="date")

mean.vent$wd<-with(mean.vent,(atan2(mean.u,mean.v) * 360/2/pi)+180)
mean.vent$ws<-with(mean.vent,((mean.u^2+mean.v^2)^.5)


Code non testé car j'ai la flemme de constituer une jeu de données bidon mais l'idée est là .

Serge

Michaël Delorme
Messages : 67
Enregistré le : 04 Avr 2016, 10:21

Re: Rose de pollution

Messagepar Michaël Delorme » 31 Mai 2018, 05:53

Voir les packages CircStats, circular... et par exemple la fonction CircStats::circ.mean()


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 1 invité