1、多个单波段栅格数据
2、合并成1个多波段
library(raster) library(terra) library(tidyverse) setwd("H:/zxzy") flnames <- list.files(path = 'H:/zxzy/ndvi/aral/', pattern = '.tif$') fl <- paste0("H:/zxzy/ndvi/aral/", flnames) fl yndvi20<- rast(fl) #Terra包读取文件夹下所有TIFF,组成一个多波段数据 names(yndvi20) = seq(1:20) #重命名为1-20 print(yndvi20) plot(yndvi20)