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

