Leaflet地图国家名称匹配异常问题排查请求

我帮你排查了代码和数据集，发现核心问题是空间多边形对象和你的统计数据没有建立正确的关联，直接用case$Total给地图上色是按列表顺序来的，不是按国家匹配关系，所以才会对应不上。另外部分地区的地图名称和你数据集里的简称不一致，也会导致匹配失败，下面是具体的解决办法：

为啥会匹配错位？ #

• 你用iso.expand(case$Code)拿到了对应国家/地区的地图，但转成SpatialPolygons后，没把case里的Total这些数据绑定到这个空间对象上。
• Leaflet渲染的时候，是按空间对象的顺序读取颜色值的，不是按国家名称或代码去匹配，这就造成了数据和地图对不上的情况。

修正后的完整代码 #

library(leaflet)
library(maps)
library(maptools)
library(htmltools)

# 读取你的数据集
case <- read.csv("Cases_02072020_v1.csv", stringsAsFactors = FALSE)

# 获取对应国家/地区的地图轮廓
Country_map <- map("world", fill = TRUE, plot = FALSE, regions = iso.expand(case$Code, regex = TRUE))

# 提取地图里的地区名称（去掉冒号后面的子区域）
IDs <- sapply(strsplit(Country_map$names, ":"), function(x) x[1])

# 转成SpatialPolygons空间对象
Country_sp <- map2SpatialPolygons(Country_map, IDs = IDs, proj4string = CRS("+proj=longlat +datum=WGS84"))

# 关键一步：把统计数据和空间对象绑定起来
# 先统一地图名称和数据集里的名称（map里是全称，你用了简称）
case$map_name <- case$Country
case$map_name[case$Country == "S. Korea"] <- "South Korea"
case$map_name[case$Country == "U.A.E."] <- "United Arab Emirates"
case$map_name[case$Country == "U.K."] <- "United Kingdom"
case$map_name[case$Country == "USA"] <- "United States"

# 匹配地图ID和数据集里的名称，确保顺序对应
match_idx <- match(IDs, case$map_name)
case_df <- case[match_idx, ]

# 生成带属性的空间数据框（这是关联的核心）
Country_spdf <- SpatialPolygonsDataFrame(Country_sp, data = case_df)

# 定义颜色渐变方案
pal <- colorNumeric(palette = "Blues", domain = Country_spdf$Total)

# 创建hover标签
Country_spdf$labels <- sprintf("<strong>Country:%s</strong><br/>Total:%g", 
                               Country_spdf$Country, Country_spdf$Total) %>% 
  lapply(HTML)

# 绘制最终的地图
leaflet(Country_spdf) %>% 
  addTiles() %>% 
  addPolygons(fillOpacity = 0.6, 
              smoothFactor = 0.5, 
              stroke = TRUE, 
              weight = 1, 
              color = ~pal(Total),  # 直接调用绑定的属性，不会错位
              label = ~labels)

重点修复说明 #

• 数据绑定：通过SpatialPolygonsDataFrame把空间轮廓和你的统计数据绑在一起，这样每个多边形都对应正确的Total值，渲染时用~pal(Total)直接调用绑定的属性，彻底解决顺序错位问题。
• 名称统一：map("world")里的部分国家是全称（比如你写的"S. Korea"在地图里是"South Korea"），必须手动映射匹配，否则会出现找不到对应地图的情况。
• 港澳台处理：你的ISO码是对的，iso.expand("HKG")、iso.expand("TWN")、iso.expand("MAC")能正确获取到这些地区的地图，不用额外调整。

额外排查小技巧 #

• 如果还有国家不显示，你可以用unique(map("world")$names)查看地图里所有可用的地区名称，对比你数据集里的名称是否一致。
• 确认case里的Code都是正确的ISO 3166-1三位码，编码错了肯定匹配不上。

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