在R语言中裁剪全球Shapefile以仅保留欧洲区域

错误提示st_crs(x) == st_crs(y) is not TRUE的核心原因是冰盖图层和国家图层的坐标参考系（CRS）不一致，空间操作要求参与运算的图层必须使用相同的CRS。

修复步骤： #

• 先查看两个图层的CRS，确认差异
• 将冰盖图层转换为与国家图层一致的CRS
• 再执行裁剪/交集操作

修改后的完整代码： #

library(sf)
library(rnaturalearth)

# 读取冰盖Shapefile
ice_extent <- st_read(shapefile_path)

# 获取世界国家数据（默认CRS是WGS84，EPSG:4326）
world <- ne_countries(scale = "medium", returnclass = "sf")

# 查看两个图层的CRS
print("冰盖图层CRS:")
print(st_crs(ice_extent))
print("国家图层CRS:")
print(st_crs(world))

# 将冰盖图层转换为国家图层的CRS
ice_extent_wgs84 <- st_transform(ice_extent, st_crs(world))

# 定义欧洲边界框并裁剪国家图层
europe_bbox <- st_bbox(c(xmin = -30, xmax = 70, ymin = 30, ymax = 90))
europe_and_surrounding <- st_crop(world, europe_bbox)

# 计算冰盖与欧洲区域的交集
ice_extent_europe <- st_intersection(ice_extent_wgs84, europe_and_surrounding)

# 绘图验证
plot(st_geometry(ice_extent_europe))

额外简化方案： #

如果只需要用边界框直接裁剪冰盖，无需先裁剪国家图层，可直接将边界框转换为对应CRS后操作：

# 将边界框转换为WGS84坐标系的空间对象
europe_bbox_sf <- st_as_sfc(europe_bbox, crs = st_crs(world))
# 直接裁剪转换后的冰盖图层
ice_extent_europe <- st_crop(ice_extent_wgs84, europe_bbox_sf)

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