[R-es] Pregunta sobre manipulación de shapefile

[Francisco J]

Gracias Óscar pruebo la segunda opción mañana y te cuento. El tema es que quiero superponer las coordenadas de los polígonos sobre los mapas de Google vía la librería ggmap, por eso la estructura de datos que decía es ideal a mi propósito

Enviado desde mi smartphone Sony Xperia™

Oscar Perpiñán Lamigueiro <oscar.perpinan en upm.es> wrote:

>
>Hola,
>
>Cuando usas readShapPoly obtienes un objeto de clase
>SpatialPolygonsDataFrame. Si lo que necesitas es añadir nuevas columnas,
>puedes hacerlo como harías con cualquier data.frame. Usando el ejemplo
>que viene en su página de ayuda, añado una nueva variable llamada test:
>
>     library(maptools)
>     xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
>      IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))
>     xx$test <- rnorm(nrow(xx))
>
>También te puede interesar la función spCbind, que te permite unir un
>Spatial*DataFrame con un data.frame.
>
>Finalmente, si necesitas las coordenadas del centroide de cada polígono,
>puedes usar coordinates:
>
>     coordinates(xx)
>
>Pero si necesitas las coordenadas de todos las líneas es más fácil si
>lees el mismo fichero con readShapeLines, y usas de nuevo coordinates. Ahora
>obtendrás una lista con las coordenadas de cada uno de los polígonos.
>
>Saludos.
>
>Oscar.
>
>
>-- 
>Oscar Perpiñán Lamigueiro
>Dpto. Ingeniería Eléctrica (ETSIDI-UPM)
>Grupo de Sistemas Fotovoltaicos (IES-UPM)
>URL: http://oscarperpinan.github.io