Using the arcpy module in Python, I want to find the proportion of each pixel value within the boundary of an overlaying polygon. (The polygon is in a shapefile.)
For example, there are 10 possible values in the raster, and intersecting the polygon 17% of the pixels are
class1, 8% are
This may involve use of the Zonal Statistics tool, but apparently I'm having trouble working out the syntax programmatically.