Leaflet Choropleth地图绘制异常及分箱调整问题咨询

步骤1：数据预处理与地理数据清洗 #

处理空几何集合（GEOMETRYCOLLECTION EMPTY） #

mapSpain获取的省份数据可能包含空几何对象，直接过滤即可：

library(mapSpain)
library(dplyr)
library(sf)
library(leaflet)

# 获取西班牙省份地理数据并过滤空几何
spain_provs <- esp_get_prov() %>%
  filter(!st_is_empty(.)) %>%
  # 统一省份名称格式，匹配房价数据的命名规则
  mutate(prov_name = tolower(ine.prov.name))

坐标系统一 #

mapSpain默认采用ETRS89坐标系（EPSG:4258），而Leaflet仅支持WGS84（EPSG:4326），需转换坐标系：

spain_provs_wgs84 <- st_transform(spain_provs, crs = 4326)

模糊连接匹配省份名称 #

用模糊匹配关联地理数据与房价数据，控制匹配精度避免错误关联：

library(fuzzyjoin)

# 假设你的房价数据框为house_prices，含province（省份名）和avg_price（均价）列
house_prices_clean <- house_prices %>%
  mutate(province = tolower(province))

# 模糊连接，设置匹配阈值（值越小匹配越严格）
merged_data <- stringdist_join(spain_provs_wgs84, house_prices_clean,
                               by = c("prov_name" = "province"),
                               max_dist = 0.2,
                               method = "jw",
                               mode = "inner") %>%
  # 去重，保留每个省份的唯一匹配结果
  group_by(cod.prov) %>%
  slice(1) %>%
  ungroup()

步骤2：解决分箱唯一值问题 #

设置10箱出现breaks非唯一错误，是因为房价数据分布集中导致自动分箱生成重复断点，可通过以下两种方案处理：

方案1：手动指定断点 #

先查看数据分布，再手动设置断点确保唯一性：

# 查看房价数据的极值与四分位数
summary(merged_data$avg_price)

# 手动设置10个断点（示例值，需根据实际数据调整）
custom_breaks <- c(0, 1000, 1500, 2000, 2500, 3000, 3500, 4000, 4500, 5000, Inf)

方案2：使用k-means聚类分箱 #

用聚类方法生成唯一断点，若去重后数量不足，适当调整分箱数：

library(classInt)

# 用k-means生成10个断点
k_breaks <- classIntervals(merged_data$avg_price, n = 10, style = "kmeans")$brks
# 去重确保断点唯一
k_breaks_unique <- unique(k_breaks)
# 若去重后少于10个，减少分箱数重新生成
if(length(k_breaks_unique) < 10) {
  k_breaks <- classIntervals(merged_data$avg_price, n = 8, style = "kmeans")$brks
}

步骤3：绘制Leaflet Choropleth地图 #

用处理好的数据和断点生成地图：

# 选择颜色方案（如viridis）
pal <- colorBin(palette = "viridis", domain = merged_data$avg_price, bins = k_breaks)

# 绘制交互地图
leaflet(merged_data) %>%
  addTiles() %>%
  addPolygons(
    fillColor = ~pal(avg_price),
    weight = 2,
    opacity = 1,
    color = "white",
    dashArray = "3",
    fillOpacity = 0.7,
    popup = ~paste0(prov_name, "<br>均价: ", avg_price, " €/m²")
  ) %>%
  addLegend(pal = pal, values = ~avg_price, opacity = 0.7, title = "平均房价(€/m²)",
            position = "bottomright")

关键问题复盘 #

• 空几何集合：mapSpain的部分海外飞地可能存在空几何，用st_is_empty()过滤即可，不影响本土省份绘制。
• 坐标系不一致：Leaflet仅支持WGS84坐标系，必须通过st_transform()转换地理数据的坐标系。
• 分箱唯一值错误：数据分布集中时自动分箱会生成重复断点，手动指定或聚类分箱可解决；若数据过于集中，建议适当减少分箱数保证断点唯一性。

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