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

Francisco Rodríguez fjroar en hotmail.com
Mie Oct 8 21:31:18 CEST 2014


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]]



Más información sobre la lista de distribución R-help-es