创建一个循环以从多个栅格中提取数据

Question

我有 8 个栅格，我厌倦了使用一组点堆叠和提取数据。但它没有成功，因为栅格的范围不同。

我真的不擅长创建循环。如何在 R studio 中编写循环函数，我可以循环提取函数并从每个栅格中提取数据并将数据导出到表中。以下是我当前的脚本。

  H1 <- readOGR(dsn = ".", layer = "Shapefile")
  L1= raster('raster.tiff')
  L2= raster('raster.tiff')
  L3= raster('raster.tiff')
  L4= raster('raster.tiff')
  L5= raster('raster.tiff')
  L6= raster('raster.tiff')
  L7= raster('raster.tiff')
  L8= raster('raster.tiff')

   files<- list()
   files[[1]] <- L1
   files[[2]] <- L2
   files[[3]] <- L3
   files[[4]] <- L4
   files[[5]] <- L5
   files[[6]] <- L6
   files[[7]] <- L7
   files[[8]] <- L8

Answer 1

我假设您的意思是函数raster::extract。另外，您的示例不完整，因此我无法提供完整的解决方案。

另外，我假设您有一组要提取的点，定义为称为coord 的矩阵。

如果我的假设是正确的，你可以使用

d <- lapply(files, extract, y = coord)

第三个参数被转发到extract，所以你实际上是在循环

extract(files[[1]], y = coord)
extract(files[[2]], y = coord)
...

您会在files 中获得一个每个元素一个元素的列表。然后你需要将它们合并在一起，例如

library(dplyr)
output <- bind_rows(d, .id = "raster")