将多个栅格提取编译到一张表中

Question

我有3个栅格，其中使用shapefile中的折线提取了数据。目前，我将其提取为3个单独的列表。有没有一种方法可以从所有三个栅格中提取数据并将它们编译到一个具有不同列的表中，以用于每个栅格的数据？

这是我正在使用的当前代码

{{1}}

Answer 1

简而言之，您可以stack()提取所有要从中提取数据的栅格，也可以extract从堆栈中提取数据。

这是一个完全可复制的示例，使用两个栅格和一个SpatialLines对象，就像您在问题中遇到的那样。跳到最后一个代码块，直接回答您的问题。

library(sp)
library(raster)

# function to generate random rasters
gen_raster <- function(){
  r <- raster(nrows = 10, ncols = 10, res = 1, 
              xmn = 0, xmx = 10, ymn = 0, ymx = 10, 
              vals = rnorm(100, 5, 1))
  return(r)
}

# generate 2 random rasters
r1 <- gen_raster()
r2 <- gen_raster()

# view
par(mfrow = c(1,2))
  plot(r1, main = "raster 1"); plot(r2, main = "raster 2")
dev.off()

# generate transect (`SpatialLines` object)
m <- as.matrix(data.frame(x = 5.5, y = seq(0, 10, 1)))
l <- list(Lines(Line(m), "m"))
l <- SpatialLines(l)

# view the transect
plot(r1, main = "raster 1 with transect"); lines(l)

在堆积的栅格上运行extract将返回一个列表，其中包含矩阵。您需要做的最后一件事是将其作为data.frame进行提取，这有点棘手。

rs <- stack(r1, r2)            # stack any amount of rasters to extract from
re <- extract(rs, l)           # extract at locations `l`
do.call(rbind.data.frame, re)  # convert to data.frame

    layer.1  layer.2
1  4.586890 5.115136
2  4.780503 5.093281
3  6.877302 3.337345
4  5.913230 3.755099
5  4.907834 4.887160
6  5.576908 5.386136
7  3.572350 5.225392
8  4.778727 5.391765
9  6.600041 4.205841
10 6.946321 5.544172

列的名称是堆栈中栅格图层的名称。您可以使用names(rs)访问这些名称，并使用names(rs) <- c("new_name_1", "new_name_2")对其进行修改。