使用terra::mask用SpatVector掩膜SpatRaster时遇‘无法创建数据集’错误

我尝试用矢量.shp文件通过terra::mask对SpatRaster做掩膜操作，执行代码时出现以下错误：

Error: [mask] cannot create dataset

执行的代码如下：

LCC84 <- rast("C:/Users_forest_VLCE2_1984.tif") 
vec <- vect("C:/Users/Land_Management_Units.shp")
vec_proj <- project(vec, LCC84)
LCC84_masked <- terra::mask(LCC84, vec_proj)

相关数据详情 #

vec
#class       : SpatVector 
#geometry    : polygons 
#dimensions  : 1, 8  (geometries, attributes)
#extent      : -117.3165, -115.1691, 50.70613, 52.27127  (xmin, xmax, ymin, ymax)
#coord. ref. : lon/lat NAD83 (EPSG:4269) 

LCC84 
#class      : SpatRaster 
#dimensions : 128340, 193936, 1 (nrow, ncol, nlyr) 
#resolution : 30, 30 (x, y) 
#extent     : -2660911, 3157169, -851351.9, 2998848 (xmin, xmax, ymin, ymax) 
#coord. ref.: Lambert_Conformal_Conic_2SP 
#source     : CA_forest_VLCE2_1984.tif 
#name       : CA_forest_VLCE2_1984


crs(LCC84, proj=TRUE)
[1] "+proj=lcc +lat_0=49 +lon_0=-95 +lat_1=49 +lat_2=77 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs"

解决方法 #

• 指定输出文件路径，避免内存溢出
  你的栅格尺寸极大（128340×193936），直接在内存生成结果会超出内存限制。terra::mask默认尝试在内存创建数据集，内存不足就会触发报错。解决方式是显式指定输出文件路径，让结果写入磁盘：

  LCC84_masked <- terra::mask(LCC84, vec_proj, filename = "C:/Your/Valid/Output/Path/masked_result.tif", overwrite=TRUE)
  

  替换为你有写入权限的路径，overwrite=TRUE参数允许覆盖已存在的同名文件。

• 验证投影后矢量与栅格的范围重叠
  虽然执行了投影转换，但需确认转换后的矢量是否与栅格有空间重叠，无重叠的掩膜操作可能触发错误：

  # 查看转换后矢量的范围
  print(ext(vec_proj))
  # 检查与栅格范围的交集
  print(intersect(ext(LCC84), ext(vec_proj)))
  

  如果交集为空，需确认矢量文件是否对应正确区域，或投影转换参数是否有误。

• 简化矢量几何以降低计算负载
  复杂的多边形边界会增加计算压力，可尝试简化矢量几何：

  # tolerance单位与投影一致（此处为米），可根据需求调整
  vec_simplified <- simplify(vec_proj, tolerance=10)
  LCC84_masked <- terra::mask(LCC84, vec_simplified, filename="C:/Output/masked_simplified.tif", overwrite=TRUE)
  

• 检查文件与路径权限
  确保R对输入的栅格、矢量文件有读取权限，对输出路径有写入权限。路径尽量使用纯英文，避免中文或特殊字符导致的异常。

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