I'm performing a fairly routine task, which is to add height information (from a 50m DEM) to survey polygons to work out mean height, as well as max and min.
I've noticed that zonal statistics (Spatial Analyst) is very slow, and depending on the raster analysis settings, can miss out polygons due to raster sampling size. For example, in the image below, the red polygon is missed out of zonal statistics because it's too thin to adequately intersect enough (or indeed any) raster cells.
Interpolate shape -> Add surface information, by contrast, seems much much faster. But I don't understand what's being interpolated here for a polygon feature. Is it just the vertices, or is the interior of the polygon also being compared with my underlying DEM? At what resolution?
Which method would provide more accurate results? And, are there any other solutions. I often find that terrain operations work equally well with rasters containing non-height information, and vice versa.