Operaciones entre vectores y raster

Aprendiendo a cortar un raster, estadisticas zonales, extraer información de puntos, etc

Vamos a iniciar…😜

Cortar un raster con un vector

library(raster)
library(sf)

# cargar las capas 

col_shp <- read_sf("Colombia_LO.shp")
bio <- raster("DEM.tif")
plot(bio)

# cortando shape con raster

b1 <- crop(bio, col_shp)

Se coloca primero la capa raster luego el shape. Crop lo que hace es que recorta el raster con la silueta de un cuadrado

# cortando shape con raster

b2 <- mask(b1, col_shp)
           # raster, poligono 

# Juntando las funciones para que sea un solo comando

cortado <- bio |> crop(col_shp) |> mask(col_shp)

Tener en cuenta

Se aconseja utilizar primero la función crop para luego ejecutar mask

Para cortar un stack se aplica el mismo proceso anterior.