Hola: No sé si he entendido bien, pero:
con la función fortify del paquete ggplot, conviertes un SpatialPolygonDataFrame en un data.frame corriente como el que tú has construido (con una columna adicional que agrupa los puntos que corresponden al mismo polígono). A ese data.frame luego le puedes añadir las columnas que quieras y pintarlo encima de un tile de google maps (siempre que la proyección del SpatialPolygon se establezca a lon lat) Aquí tienes un ejemplo https://github.com/maritrinez/R_Maps_workshop/blob/master/maps_with_R.md#2-plot-a-map-on-a-street-map El enlace debería llevarte al segundo punto del markdown, si no va bien busca "2. Plot a map on a street map", que el documento es largo! Saludos para todos El 8 de octubre de 2014, 21:31, Francisco Rodríguez <[email protected]> escribió: > Oscar, al final no he podido esperar hasta mañana > Con: > library(maptools)library(sp) > poligonos <- readShapeLines("C:\\...\\DATOS\\MUNICIPIO.shp") > sink("C:\\ ... \\DATOS\\MUNICIPIO.csv") > coordinates(poligonos) > > Consigo un fichero con la suficiente estructura para tratarlo fácilmente y > de una tacada montar la tabla .csv tal y como he indicado en un correo > anterior > Necesito el .csv aparte para utilizarlo y manipularlo después en otras > librerías. Aunque con tu respuesta soy capaz de crear una función que > permite la generación de la estructura que te he indicado de un paso de > modo totalmente automatizado, ¿Existe alguna función que haga la conversión > directamente? (Simplemente es para que no re-invente yo mismo la rueda) Mi > intención es usar después ggmap con ggplot2 > Un saludo y gracias en todo caso, con tu indicación ya te digo, dar > solución a mi objetivo de crear algo que directamente me lea los ficheros > de modo como yo quiero en un solo paso y todo dentro de R (como debe ser) > > From: [email protected] > > To: [email protected] > > CC: [email protected] > > Subject: Re: [R-es] Pregunta sobre manipulación de shapefile > > Date: Wed, 8 Oct 2014 19:01:29 +0200 > > > > > > 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 > > [[alternative HTML version deleted]] > > > _______________________________________________ > R-help-es mailing list > [email protected] > https://stat.ethz.ch/mailman/listinfo/r-help-es > > -- --- Beatriz Martínez @maritrinez <https://twitter.com/maritrinez> [[alternative HTML version deleted]] _______________________________________________ R-help-es mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-help-es
