如何裁剪MultiPolygon Shapefile？求助解决R语言st_crop警告问题

核心问题拆解 #

• 经纬度参数错位：你使用的2020人口普查城乡数据坐标系为WGS84（EPSG:4326），其中x对应经度（美国东北部为西经，数值为负），y对应纬度（北纬为正）。你之前的裁剪参数不仅把经纬度的x/y位置搞反，还未给经度加上负号，导致裁剪范围完全偏离目标区域。
• st_crop警告的本质：这个警告是因为st_crop的快速裁剪逻辑默认假设属性字段在裁剪后的几何中保持空间恒定，但MultiPolygon被裁剪后可能拆分出多个子几何，属性会被重复赋值，因此触发提示。改用st_intersection能更规范地处理复杂几何的裁剪，同时规避该警告。

修正后的完整代码 #

# 加载依赖包
library(sf)
library(ggplot2)

# 读取Shapefile
practice <- st_read("tl_rd22_us_uac20.shp")

# 确认坐标系（应为EPSG:4326）
st_crs(practice)

# 定义正确的裁剪范围：西经80°至72°，北纬40°至45°
target_bbox <- st_bbox(
  c(xmin = -80, ymin = 40, xmax = -72, ymax = 45),
  crs = st_crs(practice)
)
# 将边界框转换为多边形，用于交集运算
bbox_polygon <- st_as_sfc(target_bbox)

# 用st_intersection完成裁剪，替代st_crop
practice_cropped <- st_intersection(practice, bbox_polygon)

# 绘制裁剪后的结果
ggplot() + 
  geom_sf(data = practice_cropped, size = 1.5, color = "black", fill = "cyan1") + 
  coord_sf()

补充说明 #

• 若坚持使用st_crop，可修正参数后添加warn = FALSE参数抑制警告，但st_intersection更适配MultiPolygon这类复杂几何，能保证属性与裁剪后的几何正确对应。
• 确保Shapefile文件路径正确，或直接将下载解压后的文件放置在当前R工作目录下。

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