gpt4 book ai didi

r - 在 R 中插入光栅砖的时间序列的有效方法

转载 作者:行者123 更新时间:2023-12-04 01:59:58 25 4
gpt4 key购买 nike

我有一些代表时间序列的大栅格砖(每一层是一个时间点)。我希望通过在层之间进行插值来提高时间分辨率(在下面的示例中,插值使用 zoo::na.spline,但如果有帮助,我愿意接受其他插值方法)。我不能直接在积木上使用 na.spline - 它不是为此而设计的,只是在 as.array.default(X) 中给出 Error : 'dimnames' applied to non-array

目前我正在使用一系列循环来执行此操作(第一个循环生成具有更多时间层的砖 block ,第二个循环将原始砖 block 插入正确的层中,第三个循环使用 na .spline 在每个单元格上)。然而,这在大块上非常慢,而且似乎是一种相当低效的方法。

一个最小的可重现示例。

首先让我们制作一个代表低时间分辨率数据的初始栅格砖。请注意,左上角的单元格始终为 NA,因为任何解决方案都必须对仅包含 NA 的单元格具有鲁棒性:

library(raster); library(rasterVis); library(zoo)
r1 = raster(matrix(c(NA, rep(1,8)), 3, 3))
r2 = raster(matrix(c(NA, rep(2,8)), 3, 3))
r3 = raster(matrix(c(NA, rep(3,8)), 3, 3))
b1 = brick(list(r1,r2,r3))
levelplot(b1)

enter image description here

现在,让我们创建一个空的光栅 block ,在其中插入值:

b2 = brick(lapply(1:9, function(x) raster(matrix(NA, 3, 3))))

接下来我们将 b1 的层插入到 b2 的适当层中。请注意,我什至在这里使用了一个循环,因为 subassignment of selected layers into raster bricks does not work .

old.layers = c(1,4,7)
for (i in 1:nlayers(b1)) {
b2[[old.layers[i]]] = b1[[i]]
}
levelplot(b2, layout=c(9,1))

enter image description here

最后我们遍历每个单元格以进行时间序列插值并将结果插入回 b2:

for (cell in 1:ncell(b2)) {
if (all(is.na(as.vector(b2[cell])))) {
# na.spline wil fail if all are NA
new.values = as.numeric(rep(NA, dim(b2)[3]))
} else {
new.values = na.spline(as.vector(b2[cell]), na.rm = F)
}
b2[cell][] = new.values
}
levelplot(b2, layout=c(9,1))

enter image description here

这可行,但似乎效率极低且不够优雅。有没有更快(最好也更优雅)的方法来做到这一点?

最佳答案

您可以将 raster::calc 与如下函数一起使用:

yy <- rep(NA, 9)
fun <- function(y) {
if (all(is.na(y))) yy else na.spline((yy[c(1,4,7)] <- y), xout=1:9)
}
b2 <- calc(b1, fun)

作为引用,raster::approxNA 接近您想要的(但它不添加层,并使用线性插值)。

关于r - 在 R 中插入光栅砖的时间序列的有效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/47949175/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com