R语言栅格数据处理合集（二）——多个栅格数据合并成一个多波段数据

最新推荐文章于 2024-10-12 17:44:31 发布

原创最新推荐文章于 2024-10-12 17:44:31 发布 · 1.6k 阅读

9 ·

CC 4.0 BY-SA版权

文章标签：

#r语言 #开发语言

本文介绍了如何使用R语言中的raster、terra和tidyverse库，从H:/zxzy/ndvi/aral/文件夹下的多个单波段TIFF文件中读取并合并成一个多波段数据集，展示了yndvi20数据集的创建过程。

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)