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

问题描述

如何为 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 中是否有更直接的等价物？

解决方法

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

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

不再需要 getValuesBlock。

,

谢谢罗伯特。您的解决方案（当然！）要好得多，因为执行速度要快得多。但是，从帮助页面中不太容易理解 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,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,] 165   175  226
# ...
readStop(logo)

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

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

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

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

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