Thierry Castel

Regional climate, Variability, Climate Change Impacts

UE4: R climate data practices

Comment récupérer et lire des données climatiques issues du site de service climatique Drias de Météo-France. Au préalable il faut ouvrir un compte (gratuit) pour accéder aux données librement. Les détails de l’utilisation seront vus en séance de TP.

Pour rappel :

  • avec R les # permettent de mettre des commentaires dans le code
  • Pour l'assignement dans un objet on peut utiliser <- ou =
  • Parmis les fonctions très utiles pour explorer les données : class, str, dim, head, tail
## file data from drias : 'tasmintasmaxrstr_metro_CNRM_Aladin_histo_QT_REF_19750101-20051231_1810161527930443.KEYu11UB3Ax3u02D2uxu0D0.zip'
## unzip produces the file : 'tasmintasmaxrstr_metro_CNRM_Aladin_histo_QT_REF_19750101-20051231.txt'
## skip value is estimate by counting header lines up to the data
mydrias <- read.csv("tasmintasmaxrstr_metro_CNRM_Aladin_histo_QT_REF_19750101-20051231.txt", header = FALSE, skip=53)
colnames(mydrias) <- c("date","idpt","lat","lon","alti","Tn","Tx","RR") 
## A adapter en fonction de l'ordre que vous avez précisé lors de l'extration des données sur le site Drias
head(mydrias)
A data.frame: 6 × 8
dateidptlatlonaltiTnTxRR
<chr><int><dbl><dbl><int><dbl><dbl><dbl>
101/01/1975937846.22183.94455302 0.2 7.33.58
201/02/1975937846.22183.94455302 6.4 7.50.00
301/03/1975937846.22183.94455302 0.5 8.30.00
401/04/1975937846.22183.94455302 1.6 3.81.00
501/05/1975937846.22183.94455302-0.3 1.50.13
601/06/1975937846.22183.94455302 0.8-0.60.51
hist(mydrias$Tn)

png

Pour la suite nous vous renvoyons sur le notebook