Policz wartość rastra w obwiedni (zasięg przestrzenny)

Biorąc pod uwagę obwiednię, jak mogę policzyć liczbę określonych wartości (powiedzmy, jestem zainteresowany liczbą value == 1) w rastrze w 1) arcpy lub 2) pakiecie rastrowym w R?

W Rużyj, cropaby wyodrębnić wartości i ( np. ) tableJe policzyć.

Na przykład stwórzmy siatkę 1 stopnia obejmującą glob:

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

Obwiednia jest konwertowana na extentobiekt w celu użycia crop:

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

Po wykonaniu tej czynności zestawienie tabel jest proste:

table(getValues(y.raster))

W tym wyniku pierwszy wiersz zawiera wartości, a drugi ich odpowiadające liczby:

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

Jako czek możemy wykreślić raster i zakres:

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

Drobny dodatek: możesz również użyć (bezpiecznej pamięci) funkcji „freq”:

Po odpowiedzi Whuber:

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)

Ale teraz wykonaj:

freq(y.raster)

Ma to znaczenie tylko w przypadku bardzo dużych obiektów (raster w pliku). „freq” zwraca macierz dwukolumnową (wartość / liczba), podczas gdy „table” zwraca tabelę.