单位为米和度的栅格文件裁剪时范围不重叠报错如何解决

问题根源 #

extents do not overlap报错的核心原因是两个栅格坐标系不匹配：

• 参考裁剪栅格Covs为地理坐标系，坐标单位为十进制度
• 待裁剪栅格Mgpred为投影坐标系，坐标单位为米
  不同坐标系下的范围数值无直接可比性，直接用度单位的范围裁剪米单位的栅格，自然会识别为范围无重叠。

可用解决方案 #

根据你后续数据使用需求，可任选以下两种方案操作：

方案1：重投影待裁剪栅格到参考坐标系后裁剪 #

输出结果和参考栅格Covs坐标系、范围完全匹配，适合后续需要和Covs联合计算的场景：

# 提取参考栅格坐标系
cov_crs <- crs(Covs)
# 重投影待裁剪栅格，res参数可指定为Covs的分辨率保证对齐
Mgpred_proj <- projectRaster(Mgpred, crs = cov_crs, res = res(Covs))
# 执行裁剪
cr <- crop(Mgpred_proj, extent(Covs))
# 如需完全匹配参考栅格的行列数，可追加重采样步骤
cr_aligned <- resample(cr, Covs)

方案2：转换参考范围坐标系后裁剪 #

输出结果保留原待裁剪栅格Mgpred的投影和分辨率，精度损失更小，适合后续需要和其他同投影栅格联合使用的场景：

# 将参考范围转为空间多边形对象
cov_ext <- as(extent(Covs), "SpatialPolygons")
# 给范围指定和Covs一致的坐标系
crs(cov_ext) <- crs(Covs)
# 将范围重投影到Mgpred的坐标系
cov_ext_proj <- spTransform(cov_ext, crs(Mgpred))
# 用转换后的范围裁剪原待裁剪栅格
cr <- crop(Mgpred, cov_ext_proj)

校验建议 #

操作完成后可执行extent(cr)、crs(cr)核对输出结果的范围和坐标系是否符合预期。

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