从 Rasterbrick 中的较低层提取数据

Question

所以我从我使用这个问题的方法制作的光栅砖中提取数据：How to extract data from a RasterBrick?

除了从日期给定的图层中获取数据之外，我还想提取几个月前的数据。在我最好的猜测中，我通过执行以下操作来做到这一点：

sapply(1:nrow(pts), function(i){extract(b, cbind(pts$x[i],pts$y[i]), layer=pts$layerindex[i-1], nl=1)})

所以提取时应该查看 layerindex i-1，这应该会给出一个月前的数据。所以 layerindex = 5 的点，应该看第 5-1 = 4 层。
然而，它并没有这样做，并且似乎给出了一些随机数或几个月前的重复。解决这个问题的正确方法是什么？

Answer 1

您的代码从前一点的层获取值，而不是前一层。

想象一下，我们正在查看第 2 行 (i=2) 中的点。您指示该层的代码是pts$layerindex[i-1]，即pts$layerindex[1]。也就是说，第1行的点所在的层。

修复很简单。为了清楚起见，我将单独编写函数：

foo = function(i) extract(b, cbind(pts$x[i],pts$y[i]), layer=pts$layerindex[i]-1, nl=1)
sapply(1:nrow(pts), foo)

我没有测试过，不过应该就这些了。