計算區域中的R

Question

在光柵不同的土地覆蓋類我怎樣才能給出下面的例子中我的土地覆蓋類別1,2,3,4,5的總面積：

library("raster") 
r <- raster(nrow=10, ncol=10) 
r<-setValues(r,c(rep(1,20),rep(2,20),rep(3,20),rep(4,20),rep(5,20))) 

一方法可能是基於這些值對柵格進行子集化處理，但是像a <- area(r[getValues(r)==1])這樣的一些例子也沒有工作，也沒有a <- area(r[r==1])。

這個website提供了一個解決方案，但它需要爲每個要分析的值創建新的柵格圖層。我寧願不這樣做，因爲我的原始柵格包含許多不同的值，並且非常大。提出了類似的方法here，但它只適用於小區域。

Answer 1

您可以在由r中的值分組的區域上使用基本函數aggregate並對它們進行求和。

aggregate(getValues(area(r, weights=FALSE)), by=list(getValues(r)), sum) 
    Group.1   x 
1  1 48166136 
2  2 126933351 
3  3 320336528