- Se deben instalar y cargar los paquetes raster y rgdal.
library(raster): Proporciona funciones para manipular y analizar datos ráster.library(rgdal): Proporciona funciones para manipular y analizar datos vectoriales, especialmente shapefiles.
install.packages("raster")
install.packages("rgdal")
2. Cargar librerías.
library(raster)
library(rgdal)
3. Asignamos un directorio de trabajo.
setwd("D:/Documents")
4. Cargar la imagen ráster. La función raster() carga la imagen ráster desde el archivo especificado en su argumento.
raster_image <- raster("3_raster/wc_alt.tif")
5. Cargar polígono. La función readOGR() carga el archivo shapefile especificado en su argumento.
polygon <- readOGR("2_vector/Jug_neo_prv.shp")
6. Recortar imagen ráster utilizando el polígono. La función crop() recorta la imagen raster utilizando el polígono especificado en su segundo argumento y el resultado se asigna a la variable «cropped_raster».
cropped_raster <- crop(raster_image, polygon)
7. Visualizar imagen recortada. La función plot() se utiliza para visualizar la imagen recortada.
plot(cropped_raster)
8. Guardar imagen ráster
writeRaster(cropped_raster, "3_raster/Jug_neo_cut.tif", format="GTiff")
