Zählen Sie den Rasterwert in einem Begrenzungsrahmen (räumliche Ausdehnung).

Wie kann ich bei einem Begrenzungsrahmen die Anzahl bestimmter Werte (z. B. interessiert mich die Anzahl value == 1) in einem Raster in 1) arcpy oder 2) Rasterpaket in R zählen?

In RVerwenden Sie cropdie Werte und (zu extrahieren zB ) , tablesie zu zählen.

Erstellen wir als Beispiel ein 1-Grad-Raster, das den Globus abdeckt:

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

Der Begrenzungsrahmen wird in ein extentObjekt konvertiert, um Folgendes zu verwenden crop:

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

Danach ist die Tabellierung einfach:

table(getValues(y.raster))

In dieser Ausgabe listet die erste Zeile die Werte und die zweite die entsprechenden Zählungen auf:

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

Zur Kontrolle können wir das Raster und den Umfang zeichnen:

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

Ein kleiner Zusatz: Sie können auch die (speichersichere) Funktion "freq" verwenden:

Nach der Antwort von 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)

Aber jetzt mach:

freq(y.raster)

Dies ist nur für sehr große Objekte von Bedeutung (Raster in Datei). 'freq' gibt eine zweispaltige Matrix (Wert / Anzahl) zurück, während 'table' eine Tabelle zurückgibt.