g <- gc(reset = T)
rm(list = ls())
options(scipen = 999, warn = -1)

setwd("C:/Users/Victor/Documents/ELECTIVO_RSTUDIO/CLASE 0")

require(pacman)
pacman::p_load(rgee, mapview, raster, sf, terra, tidyterra, tidyverse,
               patchwork, gifski, geodata, rgdal, gtools, knitr)

shp_chile <- geodata::gadm(country="CHL", level=1,
                   path=tempdir(),
                   version="latest",
                   resolution=1)

head(shp_chile)

shp_chile$NAME_1

Nota: La posición [6] corresponde a Bío-Bío

plot(shp_chile[6])

biobio <- shp_chile[6]

biobio

Nota: Su geometría corresponde a un polígono

st_crs(biobio)

Nota: La proyección es WGS 84

mapview(biobio)

prec <- worldclim_global(var="prec", res=0.5,
                         path="prec", version="2.1")

prec

Nota: Es un SpatRaster de 12 capas (meses)

st_crs(prec)

Nota: La proyección es WGS 84

Transformar SpatVector a sf

bio_bio <- st_as_sf(biobio)

class(bio_bio)

st_write(bio_bio, "BIOBIO.shp", overwrite=TRUE)

prec_stk <- stack(prec)

layers_prec <- raster::crop(prec_stk, bio_bio) %>% 
  raster::mask(bio_bio)

plot(layers_prec[[7]]) 

mapview(layers_prec[[7]])

layers_prec <- unstack(layers_prec)

namesRaster <- paste0('prec_', 1:12)
Map('writeRaster', x = layers_prec,
    filename = paste0('~/ELECTIVO_RSTUDIO/CLASE 0/',
                      namesRaster, '.tif'))