如何合并每个ID的波段，以在R中创建两个多波段TIFF？

Question

我有16个独立的光栅文件，按ID和波段命名，用于8个波段和2个ID。如何合并每个ID的频段以创建两个多频段TIFF？以下是我的数据：

    dput(rastlist)
    c("7398_b10_new.tif", "7398_b2_clip.tif", "7398_b3_clip.tif", 
"7398_b4_clip.tif", "7398_b5_clip.tif", "7398_b6_clip.tif", "7398_b7_clip.tif", 
"7398_pan_new.tif", "9609_b10_new.tif", "9609_b2_clip.tif", "9609_b3_clip.tif", 
"9609_b4_clip.tif", "9609_b5_clip.tif", "9609_b6_clip.tif", "9609_b7_clip.tif", 
"9609_pan_new.tif")

我在列表的前8个元素中做到了这一点。下面是操作步骤：

setwd("mydir")

#first import all files in a single folder as a list 
rlist <- list.files(path = "mydir", pattern='.tif$', all.files=T, full.names=F)

#substract the  first 8 raster of the list
n = tail(rlist,8)

#stack layers
rstack = stack(n)

#substract the first element of the list
one = tail(n, 1)

#get the first 4 letters of the first element of a list (to be used as raster name of the raster stack)
rexport = substr(one, 1, 4)

writeRaster(rstack, filename = rexport, options = "INTERLEAVE=BAND", overwrite = T, format = "GTiff")

基本上，我想要的是创建一个函数，每8个光栅执行上述代码。

Answer 1

循环遍历每个唯一ID：

for(ras in unique(substr(rastlist, 1, 4))){
  rlist <- list.files(path = "mydir", pattern=paste0(ras, '*.tif$'), all.files=T, full.names=F)
  rstack = stack(rlist)
  writeRaster(rstack, filename = paste0(ras, ".tif"), options = "INTERLEAVE=BAND", overwrite = T, format = "GTiff")
}

或者，同样的事情，但作为一个函数，就像你在问题中说的那样：

stacker <- function(ras){
  rlist <- list.files(path = "mydir", pattern=paste0(ras, '*.tif$'), all.files=T, full.names=F)
  rstack = stack(rlist)
  writeRaster(rstack, filename = rexport, options = "INTERLEAVE=BAND", overwrite = T, format = "GTiff")
}

rass <- unique(substring(rastlist, 1, 4))

lapply(rass, stacker)