如何在R的Leaflet中可视化sfc_Multipolygon类型地理数据？

你的代码方向是对的——Leaflet完全支持sf格式的MultiPolygon数据，问题主要出在数据匹配逻辑、样式区分和NA值处理上。下面是调整后的完整可运行代码，附带关键说明：

调整后的代码 #

library(sf)
library(jsonlite)
library(rnaturalearth)
library(dplyr)
library(leaflet)

# 加载国际组织成员数据（请提前下载数据文件到本地，替换为你的文件路径）
df <- fromJSON("本地路径/Countrycodesfull.json")

# 自定义函数：添加指定国际组织的成员标记列
ISO_memcol <- function(df, orgtosearch) {
  ind <- match(orgtosearch, unlist(df[1, "org_id"]))
  or <- lapply(1:nrow(df), function(x) unlist(df[x, "org_member"])[ind])
  or <- data.frame(matrix(unlist(or)), stringsAsFactors = FALSE)
  names(or) <- orgtosearch
  df2 <- cbind(df, or) %>% as.data.frame(stringsAsFactors = FALSE)
  return(df2)
}

# 添加欧盟（EU）和英联邦（C）的成员标记
df_org <- df %>% ISO_memcol("EU") %>% ISO_memcol("C")

# 获取世界国家边界数据（sf格式）
testmap <- ne_countries(scale = 50, type = "countries", returnclass = "sf") %>%
  select(ISO_3166_3 = adm0_a3, NAME.EN = name_en) %>% # 保留匹配用编码和国家名
  left_join(df_org %>% select(ISO_3166_3, EU, C), by = "ISO_3166_3") # 左连接避免丢失边界数据

# 转换为Leaflet默认支持的WGS84投影（EPSG:4326）
testmap_leaf <- st_transform(testmap, 4326)

# 定义颜色映射：区分不同成员身份的国家
pal <- colorFactor(
  palette = c("#1a75ff", "#ff7f0e", "#9933ff", "#cccccc"), # 蓝=欧盟，橙=英联邦，紫=双成员，灰=非成员
  domain = c("EU_member", "C_member", "Both", "Non_member"),
  na.color = "#cccccc"
)

# 生成统一的成员类型列，用于样式映射
testmap_leaf <- testmap_leaf %>%
  mutate(member_type = case_when(
    EU == "Yes" & C == "Yes" ~ "Both",
    EU == "Yes" ~ "EU_member",
    C == "Yes" ~ "C_member",
    TRUE ~ "Non_member"
  ))

# 创建交互地图
m1 <- leaflet(testmap_leaf) %>%
  addTiles() %>% # 添加默认底图
  addPolygons(
    fillColor = ~pal(member_type), # 根据成员类型填充颜色
    fillOpacity = 0.7,
    color = "#333333", # 国家边界颜色
    weight = 1,
    popup = ~paste0( # 弹窗显示国家详情
      "<strong>", NAME.EN, "</strong><br>",
      "欧盟成员: ", ifelse(is.na(EU), "否", EU), "<br>",
      "英联邦成员: ", ifelse(is.na(C), "否", C)
    )
  ) %>%
  addLegend( # 添加图例说明
    position = "bottomright",
    pal = pal,
    values = ~member_type,
    title = "成员身份",
    opacity = 0.7
  )

# 显示地图
m1

关键调整说明 #

• 数据连接优化：用left_join替代full_join，避免因编码不匹配丢失国家边界数据
• 投影确认：确保sf对象转换为EPSG:4326（Leaflet默认的WGS84坐标系），这一步你之前已经做对了
• 样式区分：通过colorFactor定义颜色映射，用case_when生成统一的成员类型列，让不同组织的国家显示差异化颜色
• NA值处理：非成员或匹配失败的国家用灰色填充，弹窗中明确显示"否"
• 交互增强：添加弹窗显示国家名称和成员身份，添加图例提升地图可读性

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