R 删除geo_连接中未使用的大地水准面

我正试图在明尼苏达州、爱荷华州和内布拉斯加州的几个县的顶部规划盈利能力。使用

传单

和

底格里斯

，我已经能够绘制出所有的县，无论我是否有相关数据。这就给我留下了几个有颜色的县，其余的都标记为NA。我有没有办法从我的

geo_join

数据中删除所有NA，这样它就不会在威斯康星州未使用的地区使用？我尝试过使用

fortify

，但当我将TIGER边界线与我的县FIPS文件合并以删除它们时，我不知道如何确定要查看的县边界

以下是我的

传单

目前的样子：

我获取地图的代码如下：

library(tigris)
library(leaflet)

pal <- colorNumeric(c("yellow","dark red"),county$Construction.Cost,na.color="white")
IA_counties <- counties(state="IA", cb=TRUE, resolution ="20m")
MN_counties <- counties(state="MN",cb=TRUE,resolution="20m")
NE_counties <- counties(state="NE",cb=TRUE,resolution="20m")
IA_merged <- geo_join(IA_counties,county,"GEOID", "GEOID")
MN_merged <- geo_join(MN_counties,county,"GEOID","GEOID")
NE_merged <- geo_join(NE_counties,county,"GEOID","GEOID")
popupIA <- paste0("County Projects: ", as.character(paste('$',formatC(format(round(IA_merged$Construction.Cost, 0), big.mark=',', format = 'f')))))
popupMN <- paste0("County Projects: ", as.character(paste('$',formatC(format(round(MN_merged$Construction.Cost, 0), big.mark=',', format = 'f')))))
popupNE <- paste0("County Projects: ", as.character(paste('$',formatC(format(round(NE_merged$Construction.Cost, 0), big.mark=',', format = 'f')))))


leaflet() %>%
      addProviderTiles("MapQuestOpen.OSM") %>%
      addLegend(pal = pal, 
                values = IA_merged$Construction.Cost, 
                position = "bottomright", 
                title = "County Projects",
                labFormat=labelFormat(prefix="$")) %>%
      addCircles(lng=yup2$lon, lat=yup2$lat,weight=.75,fillOpacity=0.01,color="red",
                 radius = 96560) %>%
      addCircles(lng=yup2$lon, lat=yup2$lat,weight=.75,fillOpacity=0.01,color="blue",
                 radius = 193121) %>%
      addPolygons(data = IA_counties, 
                  fillColor = ~pal(IA_merged$Construction.Cost), 
                  layerId=1,
                  fillOpacity = .25, 
                  weight = 0.05, 
                  popup = popupIA)%>%
      addPolygons(data=MN_counties,
                  fillColor=~pal(MN_merged$Construction.Cost),
                  fillOpacity=0.25,
                  weight=0.05, 
                  popup = popupMN) %>%
      addPolygons(data=NE_counties,
                  fillColor=~pal(NE_merged$Construction.Cost),
                  fillOpacity=0.25,
                  weight=0.05, 
                  popup = popupNE) 

库（底格里斯）
图书馆（单张）
所以，令人尴尬的是，这个问题的答案和我希望的一样简单。我调整了下面的
na.color
code，它完全符合我的要求

pal <- colorNumeric(c("yellow","dark red"),county$Construction.Cost,na.color="transparent")

pal我是
tigris
软件包的创建者。非常感谢您使用它！在GitHub（）上的
tigris
开发版本中，我在
geo_join
中添加了一个选项，以适应内部连接，这将从生成的空间数据帧中完全删除不匹配的数据（如果这是您正在寻找的）。如果需要，还可以将公共合并列名作为命名参数提供给新的
by
参数。例如：

IA_merged <- geo_join(IA_counties, county, by = "GEOID", how = "inner")

IA\u谢谢你的提醒！我只想说，我真的很感谢你的
tigris
软件包以及它在
传单
中的实用性。继续努力！