choroplethr,ggmap或ggplot2

时间:2015-12-08 14:25:42

标签: r google-maps ggplot2 ggmap choroplethr

我使用library(choroplethr)进行了一些市场分析,我对制作county_choropleth并将其覆盖在ggmap()之上或使用{{1}时有疑问在我的代码中。我想要做的就是把我的县城等一些,放置州州际公路/高速公路,并在它上面绘制圆圈/半径。

我现在拥有的是:

reference_map=TRUE

......这给了我这个情节。 enter image description here

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

我使用此代码检索第二张地图:

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())

我在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的主要问题是,当我将county_choropleth置于其上时,它会使标签,道路等变灰。例如,enter image description here

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

这是我的理想&#34; MS Paint解决此问题的方法:enter image description here

提前感谢您提供任何帮助,如果您有任何需要解答的澄清问题,请与我们联系!

1 个答案:

答案 0 :(得分:2)

对于那些后来偶然发现的人。通过将库更改为leaflettigris,我能够实现我希望做的事情。

enter image description here

我计划进行最终调整以供个人使用,但这里是使用的代码:

library(tigris)
library(leaflet)

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("GEOID", "Score")

IA_counties <- counties(state="IA", cb=TRUE, resolution ="20m")
IA_merged <- geo_join(IA_counties,county,"GEOID", "GEOID")
pal <- colorQuantile("Greens",NULL,n=3)
popup <- paste0("Profitability: ", as.character(IA_merged$Score))

yup2 <- data.frame(lon=c(-93.57217,-95.93779),lat=c(41.67269,41.25861),score=c(1,1))

leaflet() %>%
  addProviderTiles("Esri.WorldStreetMap") %>%
      addLegend(pal = pal, 
      values = IA_merged$Score, 
      position = "bottomright", 
      title = "County Profitablity: ") %>%
  addCircles(lng=yup2$lon, lat=yup2$lat,weight=1,fillOpacity=0.05,color="red",
      radius = 96560) %>%
  addCircles(lng=yup2$lon, lat=yup2$lat,weight=1,fillOpacity=0.025,color="blue",
      radius = 193121) %>%
  addPolygons(data = IA_counties, 
      fillColor = ~pal(IA_merged$Score), 
      fillOpacity = 0.15, 
      weight = 0.2, 
      popup = popup)