I would like to get the coordinates of a raster cell with a given value. I suppose many cells can take the same value, in which case several coordinates should be returned. In my first problem I'm interested in the location of the maximum value(s), but also seeking other values in planned future analysis.
I will work with very small rasters, which I cut based on polygons before processing, so speed is not so important in this case. Although fast solutions are welcome, also for wider applicability for others.
I do batch processing using python and QGIS, and can of course use other tools such as SAGA, GRASS, numpy, scipy to do this. The problem seems simple, but I have been completely stuck at this for a while now. I'm quite new to using python and QGIS.