使用R语言terra包导出多波段栅格时重命名各波段

我通过Terra包对多波段栅格执行裁剪、投影、掩膜、重采样操作后，将每个波段导出为单独的栅格文件。代码能正常运行，但输出的单个栅格文件内部仍沿用原多波段栅格的波段名称（如wc2.1_30s_bioc_GFDL-ESM4_ssp370_2021-2040_1），而非自定义的GFDL_SSP3_2040_BIOCLIM_1这类名称。需要在循环中用terra::names()实现自动重命名。

修改后的完整代码如下：

# 加载包
library(terra)

# 定义多波段栅格路径
multiband_raster_path <- "mypath.tif"

# 定义掩膜栅格路径
mask_raster_path <- "mypath.tif"

# 定义输出文件夹
output_directory <- "mypath_output"

# 加载掩膜栅格
mask_raster <- terra::rast(mask_raster_path)

# 加载多波段栅格
multiband_raster <- terra::rast(multiband_raster_path)

# 裁剪多波段栅格以提升处理速度
multiband_raster <- terra::crop(multiband_raster, mask_raster)

# 投影、掩膜、重采样多波段栅格
multiband_raster <- terra::project(multiband_raster, mask_raster)
masked_raster <- terra::mask(multiband_raster, mask_raster)
resampled_raster <- terra::resample(masked_raster, mask_raster)

# 获取波段数量
num_bands <- terra::nlyr(resampled_raster)

# 循环导出每个波段为单独文件
for (i in 1:num_bands) {
  # 提取第i个波段
  single_band <- resampled_raster[[i]]
  # 自定义波段名称
  names(single_band) <- paste0("GFDL_SSP3_2040_BIOCLIM_", i)
  # 定义输出文件路径
  output_file_path <- file.path(output_directory, paste0("GFDL_SSP3_2040_BIOCLIM_", i, ".tif"))
  # 保存单波段栅格
  writeRaster(single_band, output_file_path, overwrite=TRUE)
}

关键修改说明 #

在循环中提取单波段后，新增names(single_band) <- paste0("GFDL_SSP3_2040_BIOCLIM_", i)这一行，直接为当前单波段栅格设置自定义名称。这样导出的栅格文件内部就会使用你指定的名称，而非原多波段栅格的默认波段名。

内容的提问来源于stack exchange，提问作者O.Cowan