La figura muestra el porcentaje de entradas con mediciones sobre las precipitaciones por año. Siendo el valor esperado = 356*24. Se asume que se toma una medicion por hora.
Es decir, asumiendo mediciones cada una hora. Que porcentaje de mediciones se realizaron para cada una de las 24hs (de acerudo alos años de mediciones de cada estacion)
se separan los resultados para las estaciones 13 y 17 por contener anos con valores muy por arriba de los 200mm. (cerca de 30000mm!!!!) en algunos años. Por ejemplo la estacion 17 tiene esos valores en 2010 y 2014.
DAAC_tbl %>% group_by(estacion) %>% summarise(entradas_precip=sum(!is.na(precipitacion)),entradas=n(),percent=entradas_precip/entradas)
ggplot(DAAC_tbl %>% group_by(estacion) %>% summarise(entradas_precip=sum(!is.na(precipitacion)),entradas=n(),percent=entradas_precip/entradas))+
geom_col(aes(x=estacion,y=percent),fill='skyblue')+
theme_bw()
library(lubridate)
DAAC_tbl %>% filter(year(fecha)==2017) %>% select(fecha,precipitacion) %>% group_by(yday(fecha)) %>% summarise(precip_acum=sum(precipitacion))
DAAC_tbl %>% filter(year(fecha)==2017 & is.na(precipitacion)) %>% select(fecha,estacion,precipitacion) %>% group_by(yday(fecha),estacion) %>% summarise(precip_acum=sum(precipitacion,na.rm = T))
precipitacion_2017<-DAAC_tbl %>% filter(year(fecha)==2017 ) %>% select(fecha,estacion,precipitacion) %>% group_by(year(fecha),estacion) %>% summarise(precip_acum=sum(precipitacion,na.rm = T)) %>% select(estacion,precip_acum)
precipitacion_2017<-inner_join(DAAC_estaciones_tbl,precipitacion_2017,by=c('Estacion'='estacion')) %>% select(Estacion,lat,lng,precip_acum)
precipitacion_2017
Se incluye la precipitacion acumulada promedio, salvo para el caso de la estacion 13 (Jocoli) y la 17 por tener casos muy por arriba de los 200mm
#Calculating membership using the seccion information
coord_estaciones=as.matrix(cbind(DAAC_estaciones_tbl$lng,DAAC_estaciones_tbl$lat))
coord_plumerillo=c(-68.7715,-32.8331)
distances=spDistsN1(coord_estaciones,coord_plumerillo,longlat=T)
coord_estaciones<-cbind(coord_estaciones,estacion=DAAC_estaciones_tbl$Estacion,distances)
La estacion numero 34, es la mas cercana. Sin embargo no tiene datos validos. Por lo que se debe usar la segunda mas cercana que es la numero 16. Ubicada en Lavalle.