[R-es] Sobreponer dos mapas

[Oscar Perpiñan]

Hola,

Utiliza el argumento "sp.layout" de spplot. Debes pasarle una lista
cuyo primer elemento es el nombre de la función que vas a usar para
dibujar, en segundo lugar el objeto que vas a representar con esta
función, y a continuación (opcional) los argumentos adicionales para
la función.

En tu caso:

library(sp)
[...]
background <- list('sp.polygons', zonas)
spplot(zonas, "industria_electronica", sp.layout=background)

En help('spplot') encontrarás información y ejemplos para utilizarlo.

Saludos.

Oscar
-----------------------------------------------------------------
Oscar Perpiñán Lamigueiro
Dpto. Ingeniería Eléctrica (ETSIDI-UPM)
Grupo de Sistemas Fotovoltaicos (IES-UPM)
URL: http://oscarperpinan.github.io
Twitter: @oscarperpinan


El día 12 de marzo de 2014, 7:00, Rolando Valdez <rvaldezr en gmail.com> escribió:
> Hola,
>
> Estoy tratando de sobreponer un mapa sobre otro, uno que sirva de "background" y otro con datos. Mi primer intento lo hice usando el paquete PBSmapping.
>
> Aquí lo que hice:
>
> library(PBSmapping)
> estados <- importShapefile("Entidades_2013.shp") #Aquí cargamos las entidades federativas (32 entidades)
> zonas <- importShapefile("zm.shp") #Aquí cargamos las zonas metropolitanas del país (345 municipios)
> plotPolys(estados, xlab="Longitud", ylab="Latitud", border="gray")
> addLines(zonas, col="red")
>
> Este es el resultado: https://www.dropbox.com/s/1286j3ppk6nyve6/zonas.tiff
>
> Ahora lo que deseo hacer es agregar datos, por ejemplo el cociente de localización de la industria de la electronica para cada municipio de las zonas. Anteriormente lo he hecho con spplot, pero primero agregando los datos a traves de:
>
> zm <- readOGR(".", "zm")
> datos <- read.csv(ql.csv) #cocientes de localización por sector por municipio (345 obs)
> zonas en data=data.frame(datos)
>
>  y después lo he pintado así:
>
> spplot(zonas, c("industria_electronica"))
>
> Este es el resultado: https://www.dropbox.com/s/wbnp5hm79kvqedq/ql_indelec.tiff
>
> Quisiera saber entonces, cómo puedo hacerle de tal manera que pueda pintar los datos del cociente de localización pero teniendo de "background" la división política de las entidades federativas, ya que utilizando el spplot sólo me pinta los 345 municipios que tengo en la capa.
>
> De antemano gracias.
>
> Saludos
>
> Rolando Valdez
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es