如何使用 leaflet R 绘制基于国家/地区的等值线图

Question

世界边界 geo.json 从此处下载。https://github.com/johan/world.geo.json

我试图突出显示 3 个国家（在世界地图视图中），并根据该国家/地区的项目数量将它们绘制为渐变颜色。

这是我的步骤：

首先下载世界边界geo.json文件并读取作为底图；然后我尝试突出显示数据中的国家/地区多边形。然而事实证明，世界上所有国家都是由这三个国家的信息随机着色和标记的。这是地理数据框子集问题吗？

WorldCountry <-geojsonio::geojson_read("./GeoData/countries.geo.json", what = "sp")

#Dataframe for choropleth map
Country <- c("Bulgaria","Pakistan","Turkey")
Projects <- c(2,1,6)
data <- data.frame(Country,Projects)

#basemap
Map <- leaflet(WorldCountry) %>% addTiles() %>% addPolygons()

#set bin and color for choropleth map
bins <- c(0,1,2,3,4,5,6,7,8,9,10,Inf)
pal <- colorBin("YlOrRd", domain = data$Projects, bins = bins)

#set labels
labels <- sprintf(
  "<strong>%s</strong><br/>%g projects <sup></sup>",
  data$Country, data$Projects) %>% lapply(htmltools::HTML)

#add polygons,labels and mouse over effect
Map %>% addPolygons(
  fillColor = ~pal(data$Projects),
  weight = 2,
  opacity = 1,
  color = 'white',
  dashArray = '3',
  fillOpacity = 0.7,
  highlight = highlightOptions(
     weight = 5,
    color = "#666",
    dashArray = "",
    fillOpacity = 0.7,
    bringToFront = TRUE),
  label = labels,
  labelOptions = labelOptions(
    style = list("font-weight" = "normal", padding = "3px 8px"),
    textsize = "15px",
    direction = "auto")
)

我期待这样的事情：

Answer 1

这样就可以解决问题了！使用以下方法对 WorldCountry 进行子集化：

data_Map <- WorldCountry[WorldCountry$id %in% data$Country, ]
Map <- leaflet(data_Map) %>% addTiles() %>% addPolygons()