R语言调整RasterStack文件像素尺寸并维持宽高比及解决赋值报错

你遇到的报错是因为RasterStack属于raster/terra包定义的S4类对象，不能直接通过nrow()<-、ncol()<-的方式修改维度属性，R没有为这类对象提供对应的赋值方法，所以会提示找不到nrow<-函数。同时你原有代码的维度计算逻辑也存在错误：等比例缩放时新的列数应该是行数乘以宽高比，原代码的运算逻辑写反了。

推荐使用terra包完成需求，它是传统raster包的官方升级版本，栅格处理效率更高，原生支持RasterStack的等比例缩放，处理后可以直接保持原格式。

第一步：安装加载依赖包 #

# 首次使用先安装
install.packages("terra")
# 加载包
library(terra)

第二步：批量处理栅格列表 #

IDlist_processed <- lapply(IDlist, function(cur_rst) {
  # 读取原始维度
  ori_row <- nrow(cur_rst)
  ori_col <- ncol(cur_rst)
  max_ori_dim <- max(ori_row, ori_col)
  # 尺寸符合要求直接返回不处理
  if (max_ori_dim <= 2000) {
    return(cur_rst)
  }
  # 计算等比例缩放参数
  scale <- 2000 / max_ori_dim
  new_row <- round(ori_row * scale)
  new_col <- round(ori_col * scale)
  # 执行缩放，不裁剪内容
  resized_rst <- terra::resize(cur_rst, nrow = new_row, ncol = new_col)
  # 如果你需要兼容老raster包的RasterStack格式，取消注释下行
  # resized_rst <- raster::stack(resized_rst)
  return(resized_rst)
})

方案优势 #

• 用lapply替代for循环代码更简洁，运行效率更高
• 自动跳过尺寸符合要求的栅格，无冗余计算
• 内置的缩放方法自动处理栅格值重采样，不会损失内容
• 输出结果和原格式完全兼容，无需额外转换

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