设置相同CRS后用Shapefile裁剪SpatRaster仍提示范围不重叠

出现这个错误的核心原因是你强行给栅格/矢量设置CRS的操作导致两者坐标范围不匹配，而非数据本身不重叠。以下是具体排查和解决步骤：

1. 先检查原始CRS，不要手动赋值 #

你当前代码中直接用crs(r) <- "epsg:4326"和crs(shp) <- "epsg:4326"的操作是错误的——如果原始数据已有正确CRS，强行覆盖会导致坐标错位。先读取数据后查看原始CRS：

# 读取比利时边界，保留原始CRS
ISO <- countrycode("Belgium", origin = 'country.name', destination = 'iso3c')
shp <- gadm(country=ISO, level=0, path = tempdir(), version="latest")
print("Shapefile原始CRS:")
print(crs(shp))

# 读取ForestClim栅格，保留原始CRS
r <- rast("~/Downloads/ForestClim_01.tif")
print("栅格原始CRS:")
print(crs(r))

2. 统一CRS后再裁剪 #

如果两者CRS不一致，将矢量投影到栅格的CRS（优先选择，避免栅格重采样损失精度）：

# 将比利时边界转换为栅格的CRS
shp_proj <- project(shp, crs(r))

# 执行裁剪+掩膜
cr <- crop(r, shp_proj, mask=TRUE)

如果需要将栅格转换为矢量的CRS（适合栅格分辨率较低的场景）：

# 将栅格转换为WGS84（EPSG:4326）
r_proj <- project(r, "epsg:4326")
cr <- crop(r_proj, shp, mask=TRUE)

为什么处理CHELSA数据时正常？ #

CHELSA的栅格数据原始CRS通常就是WGS84（EPSG:4326），你强行赋值的操作刚好匹配原始CRS，所以坐标范围没有错位；但ForestClim的栅格原始CRS可能不是4326，强行赋值后导致栅格的坐标范围和比利时边界完全不重叠，触发错误。

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