r - 如何根据另一个二进制文件中的类计算一个二进制文件中变量的平均值？

Question

我有两个具有相同尺寸的二进制文件(栅格):第一个表示 2 个数据之间的相关性，第二个表示具有 10 个类别的土地覆盖图。我想根据土地覆盖类别取相关性的平均值。所以最终我们将得到一个与土地覆盖图相同的 map ，但具有相关性平均值而不是类别数。

这是两个文件:

  1- to read the first file correlation map:

   conne <- file("C:\\corr.bin","rb")
  corr<- readBin(conne, numeric(), size=4,  n=1440*720, signed=TRUE)
  y<-t(matrix((data=corr), ncol=720, nrow=1440))
    r = raster(y)

2-读取第二个文件土地覆盖图:

  conne <- file("C:\\land cover.bin","rb")
  over<- readBin(conne, integer(), size=1,  n=1440*720, signed=F)
  y1<-t(matrix((data=over), ncol=720, nrow=1440))
  r1 = raster(y1)

3-写入结果:

     to.write = file(paste("/orcomplete.bin",sep=""),"wb")
     writeBin(as.double(results), to.write, size = 4)

Answer 1

这不就是这样吗:

 tapply(y, y1, mean, na.rm=TRUE)

如果您希望类均值与输入矩阵的排列相同，请执行以下操作:

outmat <- matrix( ave( y, y1, FUN=mean, na.rm=TRUE), nrow(y), ncol(y) )