48 lines
1.9 KiB
R
48 lines
1.9 KiB
R
library(jsonlite)
|
|
library(httr)
|
|
devtools::load_all()
|
|
|
|
source("secrets")
|
|
# Create a Bearer token header
|
|
headers.default <- add_headers(
|
|
accept = "*/*",
|
|
apikey = token4
|
|
)
|
|
|
|
getAllFromCoord <- function(coord,start_date,end_date,allstations,N=3,headers,base="https://public-api.meteofrance.fr"){
|
|
three_station=getIdFromCoords(coord,allstations,N=N)
|
|
alldata=lapply(three_station$Id_station,function(statid){
|
|
print(paste("recuperer station",statid))
|
|
allstat=tryCatch(getStationData(start_date=start_date,end_date=end_date,station_id=statid,headers=headers,base=base),
|
|
error=function(e){print(e);NULL})
|
|
print(paste("done, sleep 5 sec"))
|
|
print(dim(allstat))
|
|
Sys.sleep(5)
|
|
return(allstat)
|
|
})
|
|
alldata= do.call("rbind.data.frame",alldata)
|
|
ids=three_station$Nom_usuel
|
|
names(ids)=three_station$Id_station
|
|
cbind.data.frame(alldata,Nom_usuel=ids[as.character(alldata[,1])])
|
|
}
|
|
|
|
|
|
allstations=read.csv("allstations.csv") #get all station
|
|
lacouch.coor <- c(45.3722971,5.6387118)
|
|
lamure.coor <- c(44.9167, 5.8000)
|
|
vignass.coor=c(44.8550665,5.8441789)
|
|
foreve=list()
|
|
for(y in 0:3){
|
|
enddate=format(Sys.Date()-y*365, "%Y-%m-%d")
|
|
startdate=format(Sys.Date()-(y*365+364), "%Y-%m-%d")
|
|
foreve=tryCatch(getAllFromCoord(vignass.coor,startdate,enddate,allstations,headers=headers.default,base="https://public-api.meteofrance.fr"),error=function(e)e)
|
|
}
|
|
|
|
test1=getAllFromCoord(vignass.coor,start_date=startdate,end_date=enddate,allstations,headers=headers.default)
|
|
|
|
cols=palette.colors()[1:length(unique(test1$Nom_usuel))]
|
|
names(cols)=unique(test1$Nom_usuel)
|
|
plot(getDate(test1[,2]),test1[,3],pch=20,col=cols[test1$Nom_usuel],cex=2)
|
|
legend("topleft",col=cols,legend=names(cols),pch=20,cex=2)
|
|
|
|
abline(v=as.numeric(as.POSIXlt("2024-08-07",format="%Y-%m-%d")),lwd=3,col="red")
|