在R语言中使用多边形裁剪栅格遇mask函数方法匹配错误

问题根源 #

你遇到的错误核心是空间对象格式不兼容：maptools的wrld_simpl返回的是sp包的SpatialPolygonsDataFrame，而geodata获取的栅格是terra包的SpatRaster，两者分属不同的空间工具链，导致mask函数无法识别匹配的方法。

解决方案（推荐使用terra生态） #

步骤1：转换多边形格式 #

用terra包的vect()函数将旧格式的多边形转为terra兼容的SpatVector：

# 加载依赖包
library(maptools)
library(geodata)
library(terra)

# 获取智利多边形并转换格式
data(wrld_simpl)
cl <- wrld_simpl[wrld_simpl$NAME == "Chile", ]
cl_vect <- vect(cl)

# 获取2009年footprint栅格
foot_p <- footprint(year = 2009)

步骤2：裁剪+掩膜贴合边界 #

先通过crop()缩小栅格范围提升效率，再用terra的mask()函数保留多边形边界内的栅格值：

# 裁剪栅格至多边形外接范围
foot_cropped <- crop(foot_p, cl_vect)
# 按多边形边界掩膜
foot_masked <- mask(foot_cropped, cl_vect)

备选方案（使用sp/raster生态） #

若习惯旧工具链，可将栅格转为raster包的RasterLayer格式：

library(raster)
# 转换栅格格式
foot_raster <- raster(foot_p)
# 裁剪+掩膜一步完成
foot_final <- mask(crop(foot_raster, cl), cl)

注意事项 #

• 优先选择terra生态：它是raster/sp包的官方替代，功能更全、处理效率更高，建议后续空间数据处理逐步迁移至该体系。
• CRS自动继承：转换格式后无需重新设置坐标系，原对象的WGS84参数会被保留。

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