require(raster)
## Cargando paquete requerido: raster
## Cargando paquete requerido: sp
require(sp)
require(mapview)
## Cargando paquete requerido: mapview
nombres=list.files("~/R/",full.names = T)
nombres
## [1] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-01.tif"
## [2] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-02.tif"
## [3] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-03.tif"
## [4] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-04.tif"
## [5] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-05.tif"
## [6] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-06.tif"
## [7] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-07.tif"
## [8] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-08.tif"
## [9] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-09.tif"
## [10] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-10.tif"
## [11] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-11.tif"
## [12] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2020-12.tif"
## [13] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-01.tif"
## [14] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-02.tif"
## [15] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-03.tif"
## [16] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-04.tif"
## [17] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-05.tif"
## [18] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-06.tif"
## [19] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-07.tif"
## [20] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-08.tif"
## [21] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-09.tif"
## [22] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-10.tif"
## [23] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-11.tif"
## [24] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2021-12.tif"
## [25] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-01.tif"
## [26] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-02.tif"
## [27] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-03.tif"
## [28] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-04.tif"
## [29] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-05.tif"
## [30] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-06.tif"
## [31] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-07.tif"
## [32] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-08.tif"
## [33] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-09.tif"
## [34] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-10.tif"
## [35] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-11.tif"
## [36] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2022-12.tif"
## [37] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-01.tif"
## [38] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-02.tif"
## [39] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-03.tif"
## [40] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-04.tif"
## [41] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-05.tif"
## [42] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-06.tif"
## [43] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-07.tif"
## [44] "C:/Users/sebas/Documents/R/wc2.1_cruts4.09_2.5m_prec_2023-08.tif"
tmax_2024_9=raster("~/R/wc2.1_cruts4.09_2.5m_prec_2020-01.tif")
tmax_2024_9
## class : RasterLayer
## dimensions : 4320, 8640, 37324800 (nrow, ncol, ncell)
## resolution : 0.04166667, 0.04166667 (x, y)
## extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## source : wc2.1_cruts4.09_2.5m_prec_2020-01.tif
## names : wc2.1_2.5m_prec_01
## values : 0, 1071.8 (min, max)
## z-value : 2020
plot(tmax_2024_9)
tmax_all=stack(nombres)
tmax_all
## class : RasterStack
## dimensions : 4320, 8640, 37324800, 44 (nrow, ncol, ncell, nlayers)
## resolution : 0.04166667, 0.04166667 (x, y)
## extent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## names : wc2.1_2.5m_prec_01.1, wc2.1_2.5m_prec_02.1, wc2.1_2.5m_prec_03.1, wc2.1_2.5m_prec_04.1, wc2.1_2.5m_prec_05.1, wc2.1_2.5m_prec_06.1, wc2.1_2.5m_prec_07.1, wc2.1_2.5m_prec_08.1, wc2.1_2.5m_prec_09.1, wc2.1_2.5m_prec_10.1, wc2.1_2.5m_prec_11.1, wc2.1_2.5m_prec_12.1, wc2.1_2.5m_prec_01.2, wc2.1_2.5m_prec_02.2, wc2.1_2.5m_prec_03.2, ...
## min values : 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ...
## max values : 1071.8, 1299.9, 2209.6, 746.8, 1095.7, 1742.4, 3322.5, 3620.1, 1300.8, 1757.0, 1083.7, 1326.7, 1238.7, 893.8, 1904.9, ...
plot(tmax_all[[1:12]])
loc_puj_cali=cbind(-76.4147,3.8969)
plot(tmax_all[[1]])
points(loc_puj_cali)
extract(tmax_all[[1:24]],loc_puj_cali)
## wc2.1_2.5m_prec_01.1 wc2.1_2.5m_prec_02.1 wc2.1_2.5m_prec_03.1
## [1,] 52.6 103.8 86.4
## wc2.1_2.5m_prec_04.1 wc2.1_2.5m_prec_05.1 wc2.1_2.5m_prec_06.1
## [1,] 116 133.8 152.6
## wc2.1_2.5m_prec_07.1 wc2.1_2.5m_prec_08.1 wc2.1_2.5m_prec_09.1
## [1,] 131.2 103 103.9
## wc2.1_2.5m_prec_10.1 wc2.1_2.5m_prec_11.1 wc2.1_2.5m_prec_12.1
## [1,] 104.6 152.4 154.9
## wc2.1_2.5m_prec_01.2 wc2.1_2.5m_prec_02.2 wc2.1_2.5m_prec_03.2
## [1,] 77.3 128.9 205.5
## wc2.1_2.5m_prec_04.2 wc2.1_2.5m_prec_05.2 wc2.1_2.5m_prec_06.2
## [1,] 181 182.8 173.9
## wc2.1_2.5m_prec_07.2 wc2.1_2.5m_prec_08.2 wc2.1_2.5m_prec_09.2
## [1,] 80.2 157.9 115.7
## wc2.1_2.5m_prec_10.2 wc2.1_2.5m_prec_11.2 wc2.1_2.5m_prec_12.2
## [1,] 214.7 212.5 84.5
##extraer la serian de tiempo tmax en javeriana Cali
#y graficarla
tmax_puj=extract(tmax_all,loc_puj_cali)
tmax_puj
## wc2.1_2.5m_prec_01.1 wc2.1_2.5m_prec_02.1 wc2.1_2.5m_prec_03.1
## [1,] 52.6 103.8 86.4
## wc2.1_2.5m_prec_04.1 wc2.1_2.5m_prec_05.1 wc2.1_2.5m_prec_06.1
## [1,] 116 133.8 152.6
## wc2.1_2.5m_prec_07.1 wc2.1_2.5m_prec_08.1 wc2.1_2.5m_prec_09.1
## [1,] 131.2 103 103.9
## wc2.1_2.5m_prec_10.1 wc2.1_2.5m_prec_11.1 wc2.1_2.5m_prec_12.1
## [1,] 104.6 152.4 154.9
## wc2.1_2.5m_prec_01.2 wc2.1_2.5m_prec_02.2 wc2.1_2.5m_prec_03.2
## [1,] 77.3 128.9 205.5
## wc2.1_2.5m_prec_04.2 wc2.1_2.5m_prec_05.2 wc2.1_2.5m_prec_06.2
## [1,] 181 182.8 173.9
## wc2.1_2.5m_prec_07.2 wc2.1_2.5m_prec_08.2 wc2.1_2.5m_prec_09.2
## [1,] 80.2 157.9 115.7
## wc2.1_2.5m_prec_10.2 wc2.1_2.5m_prec_11.2 wc2.1_2.5m_prec_12.2
## [1,] 214.7 212.5 84.5
## wc2.1_2.5m_prec_01.3 wc2.1_2.5m_prec_02.3 wc2.1_2.5m_prec_03.3
## [1,] 62.7 195.9 174.6
## wc2.1_2.5m_prec_04.3 wc2.1_2.5m_prec_05.3 wc2.1_2.5m_prec_06.3
## [1,] 187.1 150.8 229.5
## wc2.1_2.5m_prec_07.3 wc2.1_2.5m_prec_08.3 wc2.1_2.5m_prec_09.3
## [1,] 99.7 91.9 137.9
## wc2.1_2.5m_prec_10.3 wc2.1_2.5m_prec_11.3 wc2.1_2.5m_prec_12.3
## [1,] 228.3 227.9 88.3
## wc2.1_2.5m_prec_01.4 wc2.1_2.5m_prec_02.4 wc2.1_2.5m_prec_03.4
## [1,] 169.2 69.4 172.3
## wc2.1_2.5m_prec_04.4 wc2.1_2.5m_prec_05.4 wc2.1_2.5m_prec_06.4
## [1,] 155.5 140.4 61
## wc2.1_2.5m_prec_07.4 wc2.1_2.5m_prec_08.4
## [1,] 66.4 88.5
plot(tmax_puj[1,],type="b")
library(ggplot2)
#crear vector de fechas mensuales
fechas_mensuales = seq(from = as.Date("2020-01-01"),
to = as.Date("2023-08-01"),
by = "month")
fecha=fechas_mensuales
tmax=as.numeric(tmax_puj[1,])
df=data.frame(fecha,tmax)
df
g1=ggplot(df,aes(x=fecha,y=tmax))+geom_line()+theme_bw()
g1
