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

[Beatriz Martínez]

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 <fjroar en hotmail.com>
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: oscar.perpinan en upm.es
> > To: fjroar en hotmail.com
> > CC: r-help-es en r-project.org
> > 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
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>


-- 
---
Beatriz Martínez
@maritrinez <https://twitter.com/maritrinez>

	[[alternative HTML version deleted]]