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 <[email protected]> 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 _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
