有条件地从缓冲区内的栅格中提取

Question

我有一个栅格，存储的值在1到10之间。

如何为每个值提取和求和缓冲区中的像素数量？这样，在指定大小的缓冲区内，我的像素总数为值1、2等。

python decryptRSA.py Hh6+rJdFA0SPWvbLU8gxbrZXTnYXv3M/XlSU2IHgfGvIMXckrJk/3w7OSjadhNeyIHqzfXNXRexn721lmCh7QZbGXB/cKzuEDr9pAZU6kbrc1BWDLkTuOC5e+vAcV21sebuYQUyWjGGkuMrTtXw9nlT0+h9/GAzFS7wVTFE859w= private_key.der

uB▒#▒▒▒{4elem:1234567890:0987654321&i▒_+▒▒,▒I%▒▒▒▒▒▒(d>

Answer 1

在提问时，请提供一些示例数据。

library(raster)
r <- raster(ncol=36, nrow=18, vals=sample(5, 36*18, replace=TRUE))
xy <- cbind(-50, seq(-80, 80, by=20))[1:3,]

您可以将extract与tabulate一起使用

maxvalue <- maxValue(r)  
e <- extract(r, xy[1:3,], buffer=1000000, fun=function(x) tabulate(x, maxvalue))
matrix(e, ncol=maxvalue, byrow=T)

或

maxvalue <- maxValue(r)  
e <- extract(r, xy[1:3,], buffer=1000000)
t(sapply(e, function(x) tabulate(x, maxvalue)))

#     [,1] [,2] [,3] [,4] [,5]
#[1,]    1    3    5    5    6
#[2,]    4    2    0    0    0
#[3,]    1    2    0    1    0