如何用R代码从多幅重叠shapefile中提取最大属性值？

必要工具包 #

先确保安装并加载处理空间数据的sf和数据处理的dplyr：

# 首次运行先安装包
install.packages(c("sf", "dplyr"))

# 加载包
library(sf)
library(dplyr)

步骤1：读取并预处理Shapefile #

假设两个Shapefile路径为"path/to/file1.shp"和"path/to/file2.shp"，属性值列名为value（请替换为你实际的列名）：

# 读取Shapefile并添加来源标识
shp1 <- st_read("path/to/file1.shp") %>% mutate(source = "file1")
shp2 <- st_read("path/to/file2.shp") %>% mutate(source = "file2")

# 合并为单个空间对象
combined_shp <- bind_rows(shp1, shp2)

# 修复几何有效性（避免后续操作报错）
combined_shp <- st_make_valid(combined_shp)

步骤2：拆分重叠区域 #

用st_intersection()拆分所有重叠多边形，得到无重叠的最小单元：

# 拆分后会生成`n.overlaps`列，标记当前小多边形被多少个原多边形覆盖
split_polygons <- st_intersection(combined_shp)

步骤3：计算每个区域的最大属性值 #

按拆分后的多边形分组，提取所有覆盖它的原多边形属性值的最大值：

result <- split_polygons %>%
  group_by(geometry) %>%
  summarise(max_value = max(value, na.rm = TRUE),
            .groups = "drop")

# 查看结果
head(result)

可选：导出结果Shapefile #

如果需要保存结果：

st_write(result, "path/to/max_value_result.shp", delete_layer = TRUE)

关键提示 #

• 全程替换value为你实际的属性列名。
• st_make_valid()是必要步骤，可避免因几何无效导致的报错。
• st_intersection()会自动处理所有重叠场景，包括多层重叠区域，确保最终每个区域唯一。

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