使用R Shiny绘制比利时choropleth：解决geojson_read混合几何报错问题

遇到这种混合wkbPolygon和wkbGeometryCollection的情况，最简单的方法是用sf包来处理几何数据，它能轻松拆解几何集合并保留属性，让Leaflet可以直接绘制同一图层的所有多边形。

具体步骤： #

1. 用sf读取GeoJSON数据
   放弃geojsonio::geojson_read()，改用sf::st_read()，它对复杂几何类型的支持更友好：

   library(sf)
   # 读取本地或在线的GeoJSON文件
   belgium_mun <- st_read("belgium_municipalities.json")
   

2. 拆解GeometryCollection为单一Polygon
   使用st_collection_extract()函数，把所有几何集合中的多边形提取出来，同时完整保留每个市镇的属性数据（比如名称、ID等）：

   # 提取所有多边形类型，统一为可被Leaflet识别的格式
   belgium_mun_clean <- st_collection_extract(belgium_mun, "POLYGON")
   

   这一步会把包含飞地的市镇（比如Baarle Hertog）的几何集合拆分成多个Polygon条目，但每个条目都带着原市镇的属性——也就是说，同一个市镇的所有飞地会共享同一套属性值，绘制时会自动用相同颜色渲染。

3. 按常规流程绘制Choropleth地图
   现在你可以按照RStudio Leaflet教程的步骤继续，比如合并你的数据（比如人口、统计值），然后调用leaflet()绘制：

   library(leaflet)
   
   # 假设你有一个包含市镇ID和统计值的数据框mun_data
   merged_data <- merge(belgium_mun_clean, mun_data, by = "municipality_id")
   
   # 绘制choropleth
   leaflet(merged_data) %>%
     addTiles() %>%
     addPolygons(
       fillColor = ~colorQuantile("YlOrRd", your_value_column)(your_value_column),
       weight = 2,
       opacity = 1,
       color = "white",
       dashArray = "3",
       fillOpacity = 0.7,
       popup = ~paste0(municipality_name, ": ", your_value_column)
     )
   

为什么这方法有效？ #

st_collection_extract()会自动处理几何集合的拆分，同时复制对应的属性信息，确保飞地和主市镇的属性一致。Leaflet可以直接接受处理后的sf对象，不需要额外创建多个图层，完美解决混合几何类型的问题。

内容的提问来源于stack exchange，提问作者Robin VdE