2017-10-16 13 views
0

Ok, donc j'ai fait maid une carte dans R de l'état du Colorado. Id aiment être en mesure d'ombrer certains comtés. Mes données me donnent les noms des comtés mais je ne sais pas comment l'appliquer à mon ggmap.Comment remplir des comtés spécifiques dans un état des États-Unis R

C'est le nom de tous les comtés données de mes données de polygones

[email protected][["NAME_2"]] 
[1] "Adams"  "Alamosa"  "Arapahoe" "Archuleta" "Baca"  
[6] "Bent"  "Boulder"  "Broomfield" "Chaffee"  "Cheyenne" 
[11] "Clear Creek" "Conejos"  "Costilla" "Crowley"  "Custer"  
[16] "Delta"  "Denver"  "Dolores"  "Douglas"  "Eagle"  
[21] "El Paso"  "Elbert"  "Fremont"  "Garfield" "Gilpin"  
[26] "Grand"  "Gunnison" "Hinsdale" "Huerfano" "Jackson"  
[31] "Jefferson" "Kiowa"  "Kit Carson" "La Plata" "Lake"  
[36] "Larimer"  "Las Animas" "Lincoln"  "Logan"  "Mesa"  
[41] "Mineral"  "Moffat"  "Montezuma" "Montrose" "Morgan"  
[46] "Otero"  "Ouray"  "Park"  "Phillips" "Pitkin"  
[51] "Prowers"  "Pueblo"  "Rio Blanco" "Rio Grande" "Routt"  
[56] "Saguache" "San Juan" "San Miguel" "Sedgwick" "Summit"  
[61] "Teller"  "Washington" "Weld"  "Yuma"  

Je voudrais être en mesure de choisir une des comtés et l'ombre il un peu de couleur.

My map

J'ai aussi les chemins ici pour arriver à mes données. J'ai regardé autour et ne peux pas sembler obtenir ceci ensemble. toute aide serait appréciée

code

library(raster) 
library(ggplot2) 
library(rgdal) 

#calling our state 
states <- c('Colorado') 

#getting our countys and states 
Co <- getData("GADM",country="USA",level=2) 
Colorado <- Co[Co$NAME_1 %in% states,] 

# getting map 
bm <- ggmap(get_map(location = c(-105.56, 39), 
       maptype = "hybrid", zoom = 7)) 

#overlaying our polygon onto ggmap 
gg <- bm + geom_polygon(data = Colorado, aes(y=lat,x=long, group=group), 
alpha = 0, color = "red") 
gg + geom_path() + coord_map() 

Répondre

1

Vous pouvez sous-ensemble de vos données Colorado pour extraire les informations d'un comté spécifique, par exemple "El Paso":

subset(Colorado, Colorado$NAME_2 == "El Paso") 

Utilisez les données de sous-ensembles pour créer un polygone plein.

# Select county and fill color 
county <- "El Paso" 
county_fill_color <- "pink" 

# overlaying our polygon onto ggmap 
gg <- bm + geom_polygon(data = Colorado, aes(y=lat, x=long, group=group), alpha = 0, color = "red") 
gg <- gg + geom_path() + coord_map() 

# overlay selected county polygon with fill and alpha parameters 
gg + geom_polygon(data = subset(Colorado, Colorado$NAME_2 == county), aes(y=lat, x=long, group=group), alpha = 0.75, fill = county_fill_color) 

Vous pouvez ensuite envelopper cela dans une fonction telle county et county_fill_color sont des arguments à votre fonction.

enter image description here