R语言中作用于栅格图层的paste0()等效方法是什么？

你要实现的跨栅格对应位置像元值拼接，和向量上用paste0()/paste()的逻辑完全一致，目前栅格处理主流用terra包，不需要写逐像元循环，直接支持向量化拼接操作，步骤如下：

• 前置检查：先确认三个栅格的投影、像元大小、空间范围、行列数完全对齐，若存在错位先以其中一个为参考做重采样对齐，避免拼接结果错位。
• 核心实现代码：

# 加载处理包
library(terra)

# 读入三个栅格数据，替换为你本地的文件路径即可
raster1 <- rast("你的raster1文件路径")
raster2 <- rast("你的raster2文件路径")
raster3 <- rast("你的raster3文件路径")

# 直接拼接生成带分隔符的唯一ID栅格，等效于向量上的paste操作
id_raster <- paste(raster1, raster2, raster3, sep = "-")

运行后生成的id_raster就是字符型栅格，每个像元值就是你需要的格式，比如对应位置raster1=10、raster2=3、raster3=5时，像元值就是10-3-5。

如果后续需要扩展更复杂的自定义拼接逻辑，可以用app()函数实现，灵活性更高：

# 自定义拼接逻辑的写法
id_raster <- app(
  c(raster1, raster2, raster3),
  fun = function(pixel_vals) {
    # pixel_vals是对应位置三个栅格的像元值组成的向量
    paste(pixel_vals[1], pixel_vals[2], pixel_vals[3], sep = "-")
  }
)

如果你还在使用已停止维护的旧raster包，可以用overlay()函数实现相同效果：

library(raster)
raster_stack <- stack(raster1, raster2, raster3)
id_raster <- overlay(
  raster_stack,
  fun = function(x) paste(x[1], x[2], x[3], sep = "-")
)

补充提示：如果后续需要做栅格计算，字符型ID的运算效率很低，你可以根据三个栅格的取值范围直接生成数值型唯一ID，比如用num_id <- raster1*100 + raster2*10 + raster3，得到的数值ID和拼接的字符ID一一对应，运算速度会快很多。

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