求R中光栅和矩阵（或两个矩阵）的每一行之间的角度

Question

我有一个有三个带的光栅和一个有三个列的矩阵。

library(raster)
#Make two single band raster with random values
r1 <- raster(ncol=36, nrow=18, xmn=-1000, xmx=1000, ymn=-100, ymx=900)
set.seed(42)
values(r1) <- runif(ncell(r1))
r2 <- raster(ncol=36, nrow=18, xmn=-1000, xmx=1000, ymn=-100, ymx=900)
set.seed(123)
values(r2) <- runif(ncell(r2))
r3 <- raster(ncol=36, nrow=18, xmn=-1000, xmx=1000, ymn=-100, ymx=900)
set.seed(543)
values(r3) <- runif(ncell(r3))
#Stack them into a three-band raster
r <- stack(r1, r2, r3)

#Make a matrix with values    
y <- matrix(ncol=3, nrow=2, seq(1,6,1))

我想找到y的每一行和r的每一个单元格之间的角度。

#Function for angle
    specAngle <- function(a,b){
        acos(sum(a*b)/( sqrt(sum(a*a))*sqrt(sum(b*b))))
    }

我可以通过以下方式将光栅转换为矩阵：

m <- getValues(r)

现在，每个has矩阵都具有相同的形状（每个has矩阵中有3列，每行作为一个像素或来自y的样本）。我已经研究过如何使用嵌套的apply（）函数，但老实说，我不理解语法，因此无法使其正常工作。

我有40个非常大的超光谱光栅（425个波段）和一个50行矩阵，我需要使用它来生成40个新的光栅，每个光栅有50个波段。如何有效地执行此操作？

Answer 1

我想你可以这样做

x <- lapply(1:nrow(y), function(rownr) calc(r, function(i) specAngle(i, y[rownr,])))
s <- stack(x)

或写为循环

x <- list()
for (i in 1:nrow(y)) {
    x[[i]] <- calc(r, function(j) specAngle(j, y[i,]))
}  
s <- stack(x)

您还可以重写specAngle，以便它可以一次计算所有行（并返回一个向量），而不需要循环/重叠