Re: [R-es] realizacion de mapas regionales

Mon, 28 Apr 2014 03:50:26 -0700 

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
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es

Responder a