1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53
|
#Repertoire
setwd("XXX")
library(ncdf4)
library(stringr)
library(Rcpp)
library(plyr)
library(lubridate)
#fichier
nc <-nc_open("temp_Corse_surf_avril2013_dailymean.nc")
print(nc)
#Variables
ni <- ncvar_get(nc, "ni")
dim(ni)
nj <- ncvar_get(nc, "nj")
dim(nj)
time = ncvar_get(nc, "time")
head(time)
dim(time())
#dataframe
data_temp<-ncvar_get(nc, "TEMP")
df <- data.frame(cbind(time, as.vector(data_temp)))
colnames(df)<-c("time","temperature")
View(df)
#Enlever les NA
df2 = na.omit(df)
head(df2)
View(df2)
class(df2$temperature)
class(df2$time)
df2$time <- as.POSIXct(time,tz="UTC",origin = "1900-01-01")
df3 <- df2[order(df2$time),]
head(df3)
> head(df3)
time temperature
1 2013-04-01 10:30:00 13.30805
31 2013-04-01 10:30:00 13.24063
98 2013-04-01 10:30:00 13.61955
128 2013-04-01 10:30:00 13.38589
158 2013-04-01 10:30:00 13.28095
188 2013-04-01 10:30:00 13.58657 |
Partager