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.