如何从 terra 栅格中获取一组值？

Question

如何为 terra 栅格重现 getValuesBlock 包的函数 raster？例如，如何将terra包的getValuesBlock文档中给出的例子翻译成raster包？

logo <- brick(system.file("external/rlogo.grd", package="raster"))
getValuesBlock(logo, row=35, nrows=3, col=50, ncols=3, lyrs=2:3)

我可以写：

logoT <- rast(logo*1)
Cells <- rep(50:(50+3-1),3)+rep(((35-1):(35-1+3-1))*ncol(logoT),each=3)
extract(logoT[[2:3]], Cells)

但是公式不是那么容易找到并且容易出现脚本错误。 terra 中是否有更直接的等价物？

Answer 1

我得到了更好的解决方案，这要归功于 cellFromRowColCombine 函数使任务更简单：

Cells <- cellFromRowColCombine(logoT, 35:(35+3-1), 50:(50+3-1))
extract(logoT[[2:3]], Cells)

不再需要 getValuesBlock。

Answer 2

谢谢罗伯特。您的解决方案（当然！）要好得多，因为执行速度要快得多。但是，从帮助页面中不太容易理解 readValues 的行为。经过几次测试，我明白这取决于 SpatRaster 是否在内存中：

logo <- rast(system.file("ex/logo.tif", package="terra"))
# the SpatRaster is in a file, no value returned
readValues(logo,row=35, nrows=3, col=50, ncols=3, mat=TRUE)
#   red green blue
 
logo2 <- logo*1
# the SpatRaster is in memory, readValues works
readValues(logo2,row=35, nrows=3, col=50, ncols=3, mat=TRUE)
#      red green blue
# [1,] 155   168  220
# [2,] 167   176  231
# [3,] 165   175  226
# ...

当源是文件时，SpatRaster 必须先“打开”：

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

readStart(logo)
readValues(logo,row=35, nrows=3, col=50, ncols=3, mat=TRUE)
#      red green blue
# [1,] 155   168  220
# [2,] 167   176  231
# [3,] 165   175  226
# ...
readStop(logo)

是吗？如果是这样，也许应该在帮助页面 SpatRaster 的某处提及文件中和内存中 read and write 处理之间的区别。

此外，在 readValues 帮助页面中提及 extract 作为 see also 函数可能会有所帮助，因为它高度相关（并且更有效，用于提取块）。

最后，以下脚本引发错误：

logo <- rast(system.file("ex/logo.tif", package="terra"))
readValues(logo,row=35, nrows=3, col=50, ncols=3, mat=TRUE)
#   red green blue
 
readStart(logo)
# Error: [readStart] the file is not open for reading

对 readValues 的错误调用似乎对后面的 readStart 有副作用，本不该如此。