Sınırlayıcı bir kutudaki raster değerini say (uzamsal boyut)

Bir sınırlayıcı kutu value == 1verildiğinde, 1'deki bir rasterdeki belirli değerlerin (sayıyla ilgileniyorum ) sayısını nasıl sayabilirim?

İçinde R, cropdeğerleri ayıklamak ve ( örneğin ) tablesaymak için kullanın.

Örnek olarak, dünyayı kapsayan 1 derecelik bir ızgara oluşturalım:

library(raster)
x.raster <- raster(outer(179:0, 0:359, `+`), xmn=-180, xmx=180, ymn=-90, ymx=90)

Sınırlayıcı kutu extentkullanmak için bir nesneye dönüştürülür crop:

y.extent <- extent(cbind(c(-125,20), c(-60,50)))
y.raster <- crop(x.raster, y.extent)

Bunu yaptıktan sonra, tablolama basittir:

table(getValues(y.raster))

Bu çıktıda ilk satır değerleri listeler ve ikinci satır karşılık gelen sayılarını listeler:

165 166 167 ... 257 258
  1   2   3 ...   2   1

Bir kontrol olarak raster ve kapsamı çizebiliriz:

plot(x.raster)
plot(y.extent, add=T)

Küçük bir ek: ayrıca (bellek güvenli) fonksiyonu "freq" kullanabilirsiniz:

Whuber'ın cevabını takiben:

library(raster)
x.raster <- raster(outer(179:0, 0:359, '+'), xmn=-180, xmx=180, ymn=-90, ymx=90)
y.extent <- extent(cbind(c(-125,20), c(-60,50)))
y.raster <- crop(x.raster, y.extent)

Ama şimdi yap:

freq(y.raster)

Yalnızca çok büyük nesneler için önemlidir (dosyadaki tarama). 'freq' iki sütunlu bir matris (değer / sayı) döndürürken, 'table' bir tablo döndürür.