Convertir una tabla a shape en R

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.

Comentarios

  1. Saludos Fabio, excelente tutorial. Tengo una pregunta como asigno un sistema de coordenadas al shape de salida?
    Gracias.

    ResponderEliminar
    Respuestas
    1. Problema solucionado ;)

      proj4string(tabla) = CRS("+init=epsg:4326") ## WGS84

      Eliminar

Publicar un comentario

Entradas populares de este blog

Extracción por mascara en R

Creando un mapa en R