Підрахунок растрових значень у обмежувальному полі (просторовий розмір)

З огляду на обмежувальне поле, як я можу порахувати кількість конкретних значень (скажімо, мене цікавить кількість value == 1) у растрі в 1) arcpy або 2) растровому пакеті в R?

В Rвикористовується cropдля вилучення значень і ( наприклад ) , tableщоб порахувати їх.

Як приклад, давайте створимо сітку 1 ступеня, що охоплює земну кулю:

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

Обмежувальне поле перетворюється на extentоб'єкт з метою використання crop:

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

Зробивши це, таблиця проста:

table(getValues(y.raster))

У цьому висновку перший рядок перераховує значення, а другий перераховує їх відповідні підрахунки:

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

В якості перевірки ми можемо побудувати растр і ступінь:

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

Незначне доповнення: ви також можете використовувати (безпечну для пам'яті) функцію "freq":

Після відповіді від шубер:

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)

Але тепер зробіть:

freq(y.raster)

Це має значення лише для дуже великих об'єктів (растр у файлі). 'freq' повертає матрицю з двома стовпцями (значення / кількість), тоді як 'table' повертає таблицю.