choroplethr,ggmap或ggplot2中的道路和半径圆

我正在使用library(choroplethr)一些市场分析,我有一些关于制作我的问题county_choropleth,或者将其覆盖在我的代码之上ggmap()或reference_map=TRUE在我的代码中使用.我要做的就是把我的县城等一些地方放置州州际公路/高速公路,并在它上面绘制圆圈/半径.

我现在拥有的是:

library(choroplethr)
data <- Data.frame(County.FIPS = c(19153,19163,19153,19153,19153,19153,19153,19113,19007,19169), Score=c(812.6,769.5,757.9,757.2,722.6,712.4,69727,690.2,64539,642.5)
county <-aggregate(data$Score~data$County.FIPS,data=data,sum)
colnames(county) <- c("region", "value")
mp <- county_choropleth(county, state_zoom=c("iowa"), num_colors = 1) +
  theme(legend.position="none")+
  scale_fill_gradient2("Score", 
                       high = "dark green", 
                       low = "red", 
                       na.value = "grey90", 
                       breaks = pretty(county$value, n = 10), 
                       label = scales::dollar_format())

......这给了我这个情节.

从这里开始,我想要做的是将爱荷华州的主要州际公路覆盖在我的地图上,并创建一些半径圆圈以显示距离某些城市的距离.我希望它能从这张地图中获取元素并理想地将它们合并到我的choroplethr地图中,因为在我看来,它看起来比这个例子更清晰:

我用这段代码检索第二张地图:

library(ggmap)
test<-get_map(location = c(lon=-93.57217,lat=41.67269), maptype="roadmap",source="google",zoom=7,scale="auto")
yup <- data.frame(lon=c(-93.57217,-95.87509), lat=c(41.67269,41.23238),score=c(1,1))
ggmap(test) + stat_density2d(aes(x = lon, y = lat, fill = score,alpha=score), 
               size = 2, bins = 2, data = yup, geom = "polygon") +
  theme(legend.position="none")

我reference_map=TRUE在choroplethr库中使用的主要问题是,当我将county_choropleth放在它上面时,它会使标签,道路等变灰.例如,

那么,有没有一个简单的解决方法,包括道路和地图上绘制圆圈或者我需要用放弃choroplethr并移动到ggmap,GGPLOT2或其他什么东西？我也能够在他们的网站上找到爱荷华州DOT形状文件的路径,所以这是一个包含的选项,但我不知道具体如何只在要绘制和读入R时要求它使用主要的州际公路/高速公路.

这是我对这个问题的"理想"MS Paint解决方案:

提前感谢您提供任何帮助,如果您有任何需要解答的澄清问题,请告诉我们以获取帮助!