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

如何解决使用 terra 包的 SpatRaster 和向量之间的数学运算

我有一个关于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,1000)
    f <- i + j
    v / f
}
srd <- sds(list(i = s,j = s/10))
r <- lapp(srd,fun = fun,recursive = TRUE )

解决方法

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

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,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,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)