Intersectar shapefile con raster, ¿luego determinar la proporción de cada valor de raster?

Question

Utilizando el módulo arcpy en Python, quiero encontrar la proporción de cada valor de píxel dentro de los límites de un polígono superpuesto. (El polígono está en un archivo shapefile.)

Por ejemplo, hay 10 valores posibles en el ráster, y al intersectar el polígono el 17% de los píxeles son clase1, el 8% son clase2, etc.

Esto puede implicar el uso de la herramienta Estadísticas Zonales, pero aparentemente estoy teniendo problemas para trabajar la sintaxis de forma programática.

Answer 1

Yo usaría Tabulate Area (Spatial Analyst) en lugar de Zonal Statistics si estás buscando proporciones.

---

import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
TabulateArea("zonedata.shp", "IDStr", "valueraster", "VALUE",
             "C:/sapyexamples/output/areatable.dbf", 2)