如何在R Leaflet中绘制分组区域的合并外轮廓？

核心思路 #

利用sf包的空间数据处理能力，将多个目标区域的几何对象合并为单个（或多部分）多边形，再传入Leaflet绘制统一外轮廓。

步骤1：准备并合并空间数据 #

加载依赖包，获取或读取目标区域的空间数据，筛选后完成几何合并：

# 首次使用安装依赖包
install.packages(c("sf", "leaflet", "rnaturalearth"))

# 加载包
library(sf)
library(leaflet)
library(rnaturalearth)

# 示例：合并美国CA、NV、WA三州
# 获取美国各州sf格式数据
us_states <- ne_states(country = "united states of america", returnclass = "sf")

# 筛选目标州
target_states <- us_states[us_states$name %in% c("California", "Nevada", "Washington"), ]

# 合并几何形状（自动融合相邻区域的边界）
merged_states <- st_union(target_states)

# 转换为标准sf对象，适配Leaflet调用
merged_sf <- st_sf(geometry = merged_states)

步骤2：用Leaflet绘制合并后的外轮廓 #

将合并后的sf对象传入Leaflet，自定义样式实现仅显示外轮廓的效果：

leaflet() %>%
  addTiles() %>%  # 添加默认底图
  addPolygons(
    data = merged_sf,
    color = "#ff4500",       # 外轮廓颜色
    weight = 2,              # 轮廓线条粗细
    fillColor = "transparent",# 填充透明，仅展示外轮廓
    fillOpacity = 0
  ) %>%
  setView(lng = -119, lat = 38, zoom = 5)  # 设置初始视图位置

步骤3：扩展到全球国家合并需求 #

逻辑与上述一致，仅需替换数据源和筛选条件：

# 示例：合并G7国家
# 获取全球国家sf格式数据
world_countries <- ne_countries(returnclass = "sf")

# 筛选目标国家
target_countries <- world_countries[world_countries$name %in% c(
  "Canada", "France", "Germany", "Italy", "Japan", "United Kingdom", "United States of America"
), ]

# 合并几何形状
merged_countries <- st_union(target_countries)
merged_countries_sf <- st_sf(geometry = merged_countries)

# 绘制地图
leaflet() %>%
  addTiles() %>%
  addPolygons(
    data = merged_countries_sf,
    color = "#2e8b57",
    weight = 2,
    fillColor = "#90ee90",
    fillOpacity = 0.2
  ) %>%
  setView(lng = 0, lat = 30, zoom = 2)

关键注意事项 #

• 若使用自定义GeoJSON文件，用st_read("your_file.geojson")读取，自动转为sf格式
• st_union()支持合并任意数量区域，相邻区域会自动融合为连续轮廓，不相邻区域生成MultiPolygon，Leaflet可正常渲染
• 确保sf包为最新版本，避免空间处理的兼容性问题

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