R: État et les contours du comté pour choroplethrZip

Question

J'ai un ensemble de données avec des codes postaux que j'ai utilisé choroplethrZip pour représenter les données. Je regarde les données sur l'état et le niveau du comté. Cependant, les codes postaux ne correspondent pas nécessairement aux lignes d'état et de comté. J'ai essayé d'utiliser reference_map = TRUE, mais il n'a pas de lignes de comté, et commence à avoir l'air un peu occupé. Existe-t-il un moyen de remplacer la carte de référence par défaut par une carte comportant des lignes d'état et de comté sans le reste du détail de la carte? En d'autres termes, je ne veux pas de rues et de topographie.

Voici mon code avec les données d'exemple qui sont similaires aux données que j'utilise. Vous pouvez voir le problème avec l'état des frontières du Texas.

#zip.regions metadata file for choroplethrZip 
data(zip.regions) 
head(zip.regions) 

#Test data file:A data.frame containing population estimates 
# for US Zip Code Tabulated Areas (ZCTAs) in 2012. 
data(df_pop_zip) 

#Create a choropleth of US Zip Codes 
zip_choropleth(df_pop_zip, 
       state_zoom="texas", 
       title="2012 Texas State ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) 

#Zoom County 
dd_fips = c(48113, 48121) 
zip_choropleth(df_pop_zip, 
       county_zoom=dd_fips, 
       title="2012 Denton & Dallas ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) 

TexasPlot

DentonDallasPlot

Answer 1

J'ai utilisé ce blog et trouvé comment faire ce travail: http://www.arilamstein.com/blog/2015/07/02/exploring-the-demographics-of-ferguson-missouri/

Voici mon code final:

library(choroplethrZip) 
library(ggplot2) 
library(choroplethr) 


#Pull in zip.regions metadata file for choroplethrZip 
data(zip.regions) 
head(zip.regions) 

#Test data file:A data.frame containing population estimates 
# for US Zip Code Tabulated Areas (ZCTAs) in 2012. 
data(df_pop_zip) 

# highlight a county 
highlight_county = function(county_fips) 
{ 
    library(choroplethrMaps) 
    data(county.map, package="choroplethrMaps", envir=environment()) 
    df = county.map[county.map$region %in% county_fips, ] 
    geom_polygon(data=df, aes(long, lat, group = group), color = "yellow", fill = NA, size = 1) 
} 

#Zoom County 
dd_fips = c(48113, 48121) 
zip_choropleth(df_pop_zip, 
       county_zoom=dd_fips, 
       title="2012 Denton & Dallas ZCTA Population Estimates", 
       legend="Population", 
       reference_map = TRUE) + 
       highlight_county(dd_fips) 

Rplot Counties

Je viens juste d'ajouter la fonction highlight_county et ça a bien fonctionné. Lorsque j'ai testé avec mes données (au lieu des données de population génériques), cela a fonctionné aussi.