Hola a todos,
Tenía intenciones de haber contestado antes pero ando un poco liado...
Aprovechando la contestación de Javier, yo también utilizaría el
paquete sp y me descargaría los polígonos desde http://www.gadm.org/ (el
objeto de datos es de la clase SpatialPolygons de sp, realmente un
SpatialPolygonsDataFrame porque tiene además información sobre las
regiones: ID_0, ISO, NAME_0, ID_1, NAME_1,...). Este objeto solo
contiene las fronteras con un determinado nivel de "resolución
regional"
(niveles 0, 1 o 2) del tipo: país, autonomía, provincia. Puedes
emplearlo para representar valores regionales (coloreando dependiendo de
población p.e.) o simplemente para pintar por encima de un mapa de una
variable continua (o hacer "overlay" utilizando como máscara; me
pareció
que estabas interesado en la altura...) .
Como comenta Ronaldo, otros paquetes proporcionan funcionalidad
adicional para importar/exportar, transformar, representar datos
espaciales,...: maptools, rgdal, plotGoogleMaps,...
La referencia para el paquete sp sería: Bivand, R.S., Pebesma, E.J.
y Gómez-Rubio, V. (2008). Applied Spatial Data Analysis with R. Use R!
Series. Springer. http://www.asdar-book.org/ y como también comenta
Javier el "CRAN Task View: Analysis of Spatial Data"
http://cran.r-project.org/web/views/Spatial.html.
Un saludo,
Rubén FC
El 28/04/2014 12:49, Javier Villacampa González
escribió:> Se que lo has solucionado pero para otra vez siempre puedes utilizar
> comandos de este estilo. Son rudimentarios, pero para gente menos experta
> le pueden valer. Yo los hice solo por el ocio de ver que podÃa. Asà que
> seguro que existen maneras mejores y más eficientes de tratar mapas.
>
>
> if( !require(maptools)){
> install.packages("maptools")
> }
> require(maptools)
>
> if( !require(raster)){
> install.packages("raster")
> }
> require(raster)
>
>
> # PAra conseguir las coordenadas de las ciudades
> if( !require(maps)){
> install.packages("maps")
> }
> require(maps)
> # Cargamos todas la ciudades del mundo
> data(world.cities)
>
> par(las= 1)
> x11()
> # http://www.gadm.org/
> # http://cran.r-project.org/web/views/Spatial.html ( Descargar el pais ir
> cargar el ficher rData)
>
>
>
####################################################################################################
> # Esto hay que cargarlo antes
> # load("C:/Users/usuario/AppData/Local/Temp/ESP_adm4.RData") #
Poner el
> sitio donde este el fichero
>
####################################################################################################
>
> x11()
> adm <- getData("GADM", country= "Spain", level = 2)
> mar <-
> adm[adm$NAME_0 == "Spain" & adm$NAME_1 !="Islas
Canarias",]
>
> mar$coso <- rep(1, length(mar$NAME_2))
> mar$coso[ mar$NAME_2 =="Huesca"] <- 2
>
>
> plot(mar, bg = "dodgerblue", axes= T)
> plot(mar, lwd = 10, border= "skyblue", add= T)
> plot(mar, col= c("green4","red")[mar$coso], add=T)
> grid()
> box()
>
>
> Ciudades <- world.cities[ world.cities$country.etc == "Spain"
&
> (world.cities$name == "Sabinanigo" | world.cities$pop >
200000),]
> head(Ciudades)
>
> # Dibujo ciudades por población
> points(x= Ciudades$long, y= Ciudades$lat, cex = 5 *
> (Ciudades$pop/max(Ciudades$pop)), pch = 19, col = "indianred2" )
>
> # Selecciono algunas ciudades para poner el nombre
> Ciudades2 <- Ciudades[ Ciudades$name %in% c( "A Coruna",
"Alacant",
> "Barcelona", "Bilbao", "Cartagena",
"Cordoba", "Elx" , "Gijon", "Granada",
> "Madrid" ,
"Malaga","Murcia",
> "Oviedo", "Palma" , "Sabinanigo",
"Sevilla", "Valencia",
> "Valladolid",
"Vigo",
> "Vitoria" , "Zaragoza") , ]
> text(Ciudades2$long, Ciudades2$lat + 0.6*
> (Ciudades2$pop/max(Ciudades2$pop)), Ciudades2$name, cex = 1.2)
>
>
>
>
>
> # CANARIAS
> x11()
> mar <-
> adm[adm$NAME_0 == "Spain" & adm$NAME_1 =="Islas
Canarias",]
>
> get
> plot(mar, bg = "dodgerblue", axes= T)
> plot(mar, lwd = 10, border= "skyblue", add= T)
> plot(mar, col= "green4", add=T)
> grid()
> box()
>
>
> Ciudades <-
> world.cities[ world.cities$country.etc == "Canary Islands"
&
> (world.cities$pop > 50000),]
>
> head(Ciudades)
> summary(Ciudades)
>
> # Dibujo ciudades por población
> points(x= Ciudades$long, y= Ciudades$lat, cex = 5 *
> (Ciudades$pop/max(Ciudades$pop)), pch = 19, col = "indianred2" )
>
> Ciudades2 <- Ciudades[ Ciudades$name %in% c("Santa Cruz de
Tenerife", "Las
> Palmas", "Santa Lucia", "Arrecife",
"Arona"), ]
> # Selecciono algunas ciudades para poner el nombre
> text(Ciudades2$long, Ciudades2$lat + 0.2*
> (Ciudades2$pop/max(Ciudades2$pop)), Ciudades2$name, cex = 1.2)
> --
>
> [[alternative HTML version deleted]]
>
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
[[alternative HTML version deleted]]