Shape a Raster en R

Convertir un shapefile a raster en R


El objetivo el presente ejercicio es realizar la conversión de un shapefile a archivo raster en R, es un proceso muy común y útil en los SIG, pues permite incorporar un archivo shapefile a un ejercicio de análisis raster, ejemplo algebra de mapas o conteo de área por unidad político-administrativa, entre otros; para el presente caso tenemos el shapefile de departamentos de Colombia y un archivo raster base con resolución espacial aproximada a 1 km2, el archivo resultante tomara la extensión, la alineación y la resolución acorde a este raster base.

Primero habilitamos las librerías de las cuales haremos uso:
require(raster)
require(rgdal)

Luego abrimos el shapefile en R bajo el comando “shapefile”, y el raster base con el comando “raster”, tal como se muestra a continuación.

x = shapefile("E:/Blogger/Post_6/_data/DepartamentosCOL.shp")
y = raster("E:/Blogger/Post_6/_data/ mascara_col.tif")

Para poder visualizar ambos archivos hacemos uso del comando “plot”.

plot(x)
plot(y)

Luego escogemos un campo de los que tiene el shapefile para hacer sobre este la conversión, para este caso el campo sobre el cual haremos la conversión sera el código del departamento según el código del DANE.

head(x) #con este código podemos observar el encabezado con los primeros seis campos de la tabla de atributos.

x2 <- x[,2] #escogemos la columna Nro. 2, que corresponde al ID del departamento.

Seguido de ello hacemos la conversión bajo el comando “rasterize”:

conversion = rasterize(x2, y, FUN= "first")

Por último, guardamos este archivo raster en nuestro pc bajo el comando “writeRaster”.


writeRaster(conversion, "E:/Blogger/Post_6/_data/conversion.tif")
plot(conversion)





De este link se puede descargar el dato que se hace uso en este tutorial.

Comentarios

  1. Quisiera baral el mapa altitudinal de Colombia "Elevation Above Sea Level" desde https://gadm.org/maps/COL.html y poderlo subir como raster en R, como puedo hacerlo?
    saludos,
    Juan Carlos

    ResponderEliminar
  2. raster("path/myraster.tif"), tengo canal de youtube donde ello lo explico, es Fabio Castro-Llanos on youtube

    ResponderEliminar
  3. Hola por favor necesito el comando adecuado para cargar varios archivos ráster de una misma carpeta. Gracias

    ResponderEliminar
  4. Hola. Si lo que quiero es convertir mi shape en raster, obviamente el raster no existe, como creo un raster base, tal como se usa en este ejemplo, de donde salió el raster base que voy a cargar?

    ResponderEliminar

Publicar un comentario

Entradas populares de este blog

Extracción por mascara en R

Convertir una tabla a shape en R

¿Cómo descargar información de GBIF usando R?