使用 R 中的 rasterVis 在堆栈中的每一层上绘制特定点数据

Question

我有 2 个光栅堆叠在一起：

library(rasterVis)
r1 <- raster(system.file("external/test.grd", package="raster")) 
r2 <- r1 / 2
r.stack <- stack(r1, r2)

由于稍后我想为堆栈中的每个特定图层突出显示一些区域，因此我根据栅格值创建了两个点数据集：

pts1 <- rasterToPoints(r1, spatial=T)
idx <- which(as.data.frame(pts1)[, 1] >= 400)
pts1 <- pts1[idx, 1]

pts2 <- rasterToPoints(r2, spatial=T)
idx <- which(as.data.frame(pts2)[, 1] >= 400)
pts2 <- pts2[idx, 1]

现在，我想用 R 中 rasterVis 的 levelplot 绘制光栅堆栈。我还想用 pts1 覆盖 r1，用 pts2 覆盖 r2。

但是，一旦我添加一个点数据集，它就会用于两个层：

levelplot(r.stack) + layer(sp.points(pts1, pch=20, cex=0.1, col="black"))

如何在使用栅格堆栈的同时使用具有特定图层的特定点数据集？

我想通过单独绘制每个图层及其特定点数据集然后使用 print.trellis 来避免制作自己的子图。我试过了，但结果比带栅格堆栈的 levelplot 差。

知道如何实现吗？

Answer 1

使用 panel.number 函数，您可以根据您所在的面板对数据进行子集化：

pts <- list(pts1, pts2)

levelplot(r.stack) +
    layer(sp.points(pts[[panel.number()]],
                    pch=20, cex=0.1, col="black"))

使用 R 中的 rasterVis 在堆栈中的每一层上绘制特定点数据

Plot specific point data on each layer in stack with rasterVis in R

r

gis

lattice

r-raster

levelplot