I'm starting on a problem where I need to find the high point inside a polygon that exists inside the bounds of a DEM tile. I'm teaching myself GDAL and I'll be using Python to access the GDAL library. I've see Zonal Statistics, but isn't quite what I need. I need the elevation, AND longitude, latitude.
My impression is that I need to INTERSECT the polygon with the raster tile (ArcGrid, Geotiff) and extract an array of points (x,y,z) from the DEM. Then analyze the array of points.
In my learning, I need to figure out how to put these words into commands.