将boot :: inv.logit（）应用于R中的栅格

Question

我有一个Formal Class Raster，我试图将boot :: inv.logit（）函数应用于它的栅格单元格，例如：

r1 <- raster(nrows=25, ncols=25, vals=rtnorm(n = 625, .1, .9))
r2 <- boot::inv.logit(r1)

但是，当我尝试时，它会返回一个错误：

> Error in plogis(x) : Non-numeric argument to mathematical function

如果我将光栅转换为矩阵，然后返回光栅，它就完成了工作，但是我放弃了与开头时的“Formal Class Raster”相关的所有其他信息，这并不理想：

r2 <- boot::inv.logit(as.matrix(r1))
r2 <- as.raster(r2)

是否有一种简单的方法可以恢复之前的Formal Class Raster信息或将inv.logit()应用于没有as.matrix()变换的栅格？先感谢您。

Answer 1

如果要将该函数应用于栅格，请使用calc中的raster方法：

r2 <- calc(r1,boot::inv.logit)

> r2 

# class       : RasterLayer 
# dimensions  : 25, 25, 625  (nrow, ncol, ncell)
# resolution  : 14.4, 7.2  (x, y)
# extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
# data source : in memory
# names       : layer 
# values      : 0.07434905, 0.9498965  (min, max)

或者，你可以制作r1的空副本，只需填写inv.logit发布的值：

r2 <- raster(r1)
r2[] <- boot::inv.logit(as.matrix(r1))