Convertir una tabla a shape en R
El objetivo el presente ejercicio es realizar la conversión de una tabla, que contiene coordenadas dadas en longitud y latitud, a un archivo shapefile en R, es un proceso muy común y útil en los SIG; para el presente caso tenemos las coordenadas de incendios en el Valle del Cauca para el periodo 2000 – 2015, la idea es convertir esta tabla a shapefile y guardar este archivo en el disco duro de nuestro pc.
Primero habilitamos las librerías de las cuales haremos uso:
require(raster)
require(rgdal)
Luego abrimos la tabla en R bajo el comando “read.table”, como se muestra a continuación.
Tabla = read.csv(“E:/Blogger/Post_5/_data/incendios_valle.csv”)
Seguido de ello convertimos la tabla a formato shapefile, teniendo en cuenta que las columnas 1 y 2 corresponden a la longitud y latitud respectivamente:
coordinates(tabla) <- ~LONGITUDE+LATITUDE # se escribe de esta manera longitud y latitud pues así se llama el nombre de estas columnas
Podemos ver el tipo de clase a la que queda convertido el objeto “tabla” bajo el siguiente comando:
Class(tabla) #como podemos ver ahora es un “SpatialPointsDataFrame”
plot(tabla) #miramos que en efecto el archivo es una nube de puntos, que representan los incendios en el Valle del Cauca.
Por último, procedemos a realizar el guardado de este archivo.
writeOGR(tabla, dsn = «E:/Blogger/Post_5/_data», layer= «incendios_valle», driver = «ESRI Shapefile»)
De este link se puede descargar el dato que se hace uso en este tutorial.
SIG con R - Tutoriales 
Saludos Fabio, excelente tutorial. Tengo una pregunta como asigno un sistema de coordenadas al shape de salida?Gracias.
Me gustaMe gusta
Problema solucionado ;)proj4string(tabla) = CRS(\»+init=epsg:4326\») ## WGS84
Me gustaMe gusta