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吗？
您现在可以使用terra::values
x <- values(logoT, row=35, nrows=3, col=50, ncols=3, mat=TRUE)
head(x)
# red green blue
#[1,] 155 168 220
#[2,] 167 176 231
#[3,] 165 175 226
#[4,] 154 167 219
#[5,] 165 176 230
#[6,] 165 175 226
但是你不能指定你想要的层，不像raster::getValuesBlock
如果您有坐标而不是行号/列号，则使用起来可能更方便values(crop(logoT, ext(xmin, xmax, ymin, ymax)))
谢谢罗伯特。您的解决方案（当然！）要好得多，因为执行速度要快得多。但是，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，但情况并非如此。
cellFromRowColCombine由于使任务更简单的功能，我得到了更好的解决方案：
Cells <- cellFromRowColCombine(logoT, 35:(35+3-1), 50:(50+3-1))
extract(logoT[[2:3]], Cells)
不再需要 getValuesBlock。