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