解决R中缓冲区湿地类型面积计算的group_by报错问题

错误核心原因：你当前的crop2是sp包的SpatialPolygonsDataFrame类型对象，而dplyr的group_by等函数仅支持tidy兼容的数据结构（如普通数据框、sf包空间对象），旧版sp对象不在支持范围内。

直接用sf包重构代码即可解决问题——sf是当前R空间数据处理的标准工具，和tidyverse完全适配，具体步骤如下：

修正后的完整代码 #

library(tidyverse)
library(sf)

# 用sf的st_read读取空间数据，直接生成sf对象
Neosho_wetlands <- st_read("HU8_11070205_watershed", layer = "HU8_11070205_Wetlands")
Spring_wetlands <- st_read("HU8_11070207_watershed", layer = "HU8_11070207_Wetlands")
North_wetlands <- st_read("HU8_10290104_watershed", layer = "HU8_10290104_Wetlands")
South_wetlands <- st_read("HU8_11070206_watershed", layer = "HU8_11070206_Wetlands")

# 合并sf对象
wetlands <- bind_rows(Neosho_wetlands, Spring_wetlands, North_wetlands, South_wetlands)
rm(Neosho_wetlands, Spring_wetlands, North_wetlands, South_wetlands)
plot(st_geometry(wetlands))

# 处理点数据，转换为sf对象
points <- read.csv("SiteCov.csv") %>%
  st_as_sf(coords = c("long", "lat"), crs = 4326)  # 4326对应WGS84坐标系，即原代码的+proj=longlat

# 投影转换为与湿地数据一致的坐标系
points_sf <- st_transform(points, st_crs(wetlands))

# 创建缓冲区并裁剪湿地数据
buffer2 <- st_buffer(points_sf[1:32,], 2000)
crop2 <- st_crop(wetlands, buffer2)
plot(st_geometry(crop2))

# 按湿地类型分组计算总面积
crop2 %>%
  mutate(area = as.numeric(st_area(.))) %>%  # 计算每个多边形面积，转成普通数值方便求和
  group_by(WETLAND_TY) %>%
  summarize(total_area = sum(area)) %>%
  filter(WETLAND_TY %in% c("Freshwater Emergent Wetland",
                           "Freshwater Forested/Shrub Wetland",
                           "Freshwater Pond",
                           "Lake",
                           "Riverine"))  # 筛选你需要统计的湿地类型

关键修改说明 #

1. 替换空间数据处理工具：用st_read替代readOGR，直接生成sf对象，避免后续类型转换；用st_transform、st_buffer、st_crop替代sp包的对应函数，全程保持sf对象格式。
2. 修正面积计算逻辑：新增mutate(area = as.numeric(st_area(.)))——st_area()会计算每个多边形的面积并返回带单位的空间数值，转成普通数值后才能正常求和；原代码的summarize写法错误，现在改为按类型分组求和。
3. 筛选目标湿地类型：用filter精准筛选你需要统计的湿地类型，不需要的话可以直接删除该行。

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