R：州和縣概述choroplethrZip

Question

我有一個帶有郵政編碼的數據集，我用choroplethrZip繪製數據。我正在查看州和縣一級的數據。但是，郵政編碼並不一定與州和郡線相符。我試圖使用reference_map = TRUE，但它沒有縣線，並開始看起來有點忙。有沒有辦法將默認參考地圖更改爲具有州和縣行而沒有地圖細節其餘部分的地圖？換句話說，我不想要街道和地形。

下面是我的代碼，其中的示例數據與我正在使用的數據類似。你可以看到德克薩斯州邊界的問題。

#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

我用這個博客，發現瞭如何使這項工作：http://www.arilamstein.com/blog/2015/07/02/exploring-the-demographics-of-ferguson-missouri/

這是我的最終代碼：

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

我只是不得不在highlight_county函數中添加它，它工作得很好。當我用我的數據（而不是通用人口數據）進行測試時，它也起作用。