使用terra软件包中的rast()函数创建一个由多个图层组成的栅格堆栈,并使用add()函数将前两个图层相加,然后使用rename()函数将结果保存在第三个图层中。
示例代码:
library(terra)
# 创建栅格图层
r1 <- rast(matrix(1:12, nrow = 4, ncol = 3))
r2 <- rast(matrix(13:24, nrow = 4, ncol = 3))
r3 <- rast(matrix(0, nrow = 4, ncol = 3))
# 将两个图层相加
r_sum <- add(r1, r2)
# 将结果保存在第三个图层中
r_sum <- rename(r_sum, "third_layer")