Cómo cargar un shapefile al entorno de R y obtener información relevante

Para incorporar un archivo vectorial al entorno de R necesitamos la librería rgdal. Una vez instalada y cargada en el intérprete, para hacer lo mismo con el vectorial se tiene que utilizar el método readOGR. Para ello, nos movemos primero a nuestro espacio de trabajo, es decir, el que contiene los vectoriales e invocamos, por ejemplo, la siguiente serie de comandos:

setwd('proyectosR')
library(rgdal)
cuenca<-readOGR('.', 'cuenca')
print(cuenca)

La última instrucción despliega toda la metadata correspondiente al vectorial cuenca que se detalla a continuación:

OGR data source with driver: ESRI Shapefile 
Source: ".", layer: "cuenca"
with 1 features and 4 fields
Feature type: wkbPolygon with 2 dimensions
> print(cuenca)
An object of class "SpatialPolygonsDataFrame"
Slot "data":
  id   count        sum     mean
0  1 7354371 1433542798 194.9239

Slot "polygons":
[[1]]
An object of class "Polygons"
Slot "Polygons":
[[1]]
An object of class "Polygon"
Slot "labpt":
[1]  666852.5 1048982.1

Slot "area":
[1] 6891084688

Slot "hole":
[1] FALSE

Slot "ringDir":
[1] 1

Slot "coords":
          [,1]    [,2]
 [1,] 615534.0 1051550
 [2,] 626817.6 1072537
 [3,] 642388.9 1087432
 [4,] 668792.5 1090591
 [5,] 691359.6 1089237
 [6,] 710541.7 1077051
 [7,] 718214.5 1065090
 [8,] 716183.5 1054258
 [9,] 713926.8 1041169
[10,] 709639.0 1019053
[11,] 687974.5 1008898
[12,] 644420.0 1011380
[13,] 619144.8 1021536
[14,] 612374.6 1035527
[15,] 615534.0 1051550



Slot "plotOrder":
[1] 1

Slot "labpt":
[1]  666852.5 1048982.1

Slot "ID":
[1] "0"

Slot "area":
[1] 6891084688



Slot "plotOrder":
[1] 1

Slot "bbox":
        min       max
x  612374.6  718214.5
y 1008898.1 1090591.1

Slot "proj4string":
CRS arguments:
 +proj=utm +zone=19 +datum=WGS84 +units=m +no_defs +ellps=WGS84
+towgs84=0,0,0 

Si se desea obtener un despliegue gráfico del vectorial entonces se podría hacer lo siguiente:

plot(cuenca, col='blue', bg='yellow')
box()
for(i in 1:4) axis(i)

lo cual dará origen a la gráfica que se obtiene en la siguiente imagen:

vector

A continuación, se tiene el mismo vectorial pero visualizado en QGIS conjuntamente con su tabla atributiva.

vector2

Esta entrada fue publicada en Lenguaje R, SIG, Software Libre. Guarda el enlace permanente.

3 respuestas a Cómo cargar un shapefile al entorno de R y obtener información relevante

  1. Pingback: Cómo cargar un shapefile al entorno de R...

  2. Está muy bueno.
    Me gustaría saber como puedo hacer análisis de Nearest Neighbor de un set de puntos que se encuentran en formato shp.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión / Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión / Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión / Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión / Cambiar )

Conectando a %s