Cuaderno 6

library(AOI)
library(climateR)
library(sf)
library(raster)
library(rasterVis)
library(dplyr)

tc_prcp <- getTerraClim(valle, param = "prcp", startDate = "2019-02-01")

tc_prcp = getTerraClim(valle, param = "prcp", startDate = "2019-02-01")

tc_tmp <- tc_prcp[[1]]

tc_tmp

class      : RasterStack 
dimensions : 48, 46, 2208, 1  (nrow, ncol, ncell, nlayers)
resolution : 0.04166667, 0.04166667  (x, y)
extent     : -77.58333, -75.66667, 3.083333, 5.083333  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
names      : X2019.02 
min values :     40.7 
max values :    583.2 

library(leaflet)
library(RColorBrewer)
library(rgdal)
pal <- colorNumeric(c("RdYlBu"), values(tc_tmp$X2019.02), 
                    na.color = "transparent")

leaflet() %>% addTiles() %>%
  addRasterImage(tc_tmp$X2019.02 , colors = pal, opacity = 0.8) %>%
  addLegend(pal = pal, values = values(tc_tmp$X2019.02),
    title = "Precip-Feb.2019 [mm]")

Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definitionDiscarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definition

head(param_meta$terraclim)

tc_palmer = getTerraClim(valle, param = "palmer", startDate = "2019-02-01")

tc_tmp <- tc_palmer[[1]]

tc_tmp

class      : RasterStack 
dimensions : 48, 46, 2208, 1  (nrow, ncol, ncell, nlayers)
resolution : 0.04166667, 0.04166667  (x, y)
extent     : -77.58333, -75.66667, 3.083333, 5.083333  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
names      : X2019.02 
min values :     -2.2 
max values :      2.1 

pal <- colorNumeric(c("PuOr"), values(tc_tmp$X2019.02), 
                    na.color = "transparent")

leaflet() %>% addTiles() %>%
  addRasterImage(tc_tmp$X2019.02, colors = pal, opacity = 0.8) %>%
  addLegend(pal = pal, values = values(tc_tmp$X2019.02),
    title = "PDSI-Feb.2019")

Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definitionDiscarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definition

wat_def = getTerraClimNormals(valle, param = "water_deficit", period = "19812010", month=2)

wat_def

$terraclim_19812010_water_deficit
class      : RasterStack 
dimensions : 48, 46, 2208, 1  (nrow, ncol, ncell, nlayers)
resolution : 0.04166667, 0.04166667  (x, y)
extent     : -77.58333, -75.66667, 3.083333, 5.083333  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
names      :  X02 
min values :    0 
max values : 32.6 

tc_tmp <- wat_def[[1]]

pal <- colorNumeric(c("BrBG"), values(tc_tmp$X02), 
                    na.color = "transparent")
leaflet() %>% addTiles() %>%
  addRasterImage(tc_tmp$X02, colors = pal, opacity = 0.8) %>%
  addLegend(pal = pal, values = values(tc_tmp$X02),
    title = "Déficit agua-February")

Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definitionDiscarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null +wktext +no_defs +type=crsDiscarded datum World Geodetic System 1984 in Proj4 definitionSome values were outside the color scale and will be treated as NA