I am trying to perform zonal statistics on a input DEM and a polygon shapefile using ArcPy in Python 2.7. I have done this in Python 3, but want to be able to do it using ArcPy tools in Python 2.7. So far, I successfully got a zonal statistics table, but would like to add another feature.
I would like to create a list of all the raster points that are inside each of the polygons. So for example, one of the polygons has 1631 raster points inside its bounds. I would like to append a column on the table (or some other place, doesn't have to be in the zonal statistics table) that is a list of all these points with the value being the elevation from the DEM.
Figure shows what I would like to get, the list of elevation values for all the raster points in each polygon.
Has anyone done this before or know of a way to do it using zonal statistics in ArcPy for Python 2.7 and ArcGIS 10.6.1?