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.

2 respuestas a “Shape a Raster en R

Deja una respuesta

Introduce tus datos o haz clic en un icono para iniciar sesión:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Salir /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Salir /  Cambiar )

Conectando a %s