Question

我有一个关于app()家庭功能的简单问题

我想将一个向量（n 个值）除以一个栅格并得到 n 个结果栅格

s <- rast(system.file("ex/logo.tif", package="terra"))[[1]] 

lapp(s, aa = c(1, 10, 100, 1000), function(x, aa){0.46 * aa/(x + 1)})

同样想知道要使用 SpatRasterDataset 来完成。这可能吗？

这是一个例子：

fun = function(i, j) { 
    v <- c(1, 10, 100, 1000)
    f <- i + j
    v / f
}
srd <- sds(list(i = s, j = s/10))
r <- lapp(srd, fun = fun, recursive = TRUE )

Answer 1

您可以这样做（使用稍微修改的示例，以便更容易地看到发生了什么）：

library(terra)
s <- rast(ncols=10, nrows=10)
s <- init(s, 1:10)
v <- c(1, 10, 100, 1000)    
a <- app(s, function(x, aa){ sapply(aa, function(i) 2 * i/x) }, aa=v)
a
#class       : SpatRaster 
#dimensions  : 10, 10, 3  (nrow, ncol, nlyr)
#resolution  : 36, 18  (x, y)
#extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=longlat +datum=WGS84 +no_defs 
#source      : memory 
#names       : lyr.1, lyr.2, lyr.3 
#min values  :   0.2,   2.0,  20.0 
#max values  :     2,    20,   200

使用 terra 版本 1.3-17（当前为开发版本），您现在还可以：

b <- 2 * v/s

得到相同的结果。当使用比层数长的向量进行算术时，SpatRaster 的层会被回收。

使用 SpatRasterDataSet，您可以执行以下操作。在 fun 之外测试 lapp 以查看返回的值的数量及其顺序是否有意义，这一点很重要。

fun = function(i, j) { 
    v <- c(1, 10, 100, 1000)
    f <- i + j
    f <- rep(f, each=length(v))
    matrix(v / f, ncol=length(v), byrow=TRUE)
}
srd <- sds(list(i = s, j = s/10))
r <- lapp(srd, fun = fun)

使用 terra 包的 SpatRaster 和向量之间的数学运算

1 个答案: