Should subsetByValue work on domain axes or CRS dimensions?

[letmaik]

The Coverage.subsetByValue operation is currently defined on domain axes having an ordering, which essentially means those that correspond directly to a CRS axis. This would make it impossible to subset things like MultiPoint/Polygon coverages based on coordinates. WCS defines value based subsetting strictly based on CRS dimensions, which may also be beneficial for tiling purposes. On the other hand, it is way more complicated to implement since it requires specific algorithms for taking apart (and combining it later again) a polygon for example.

I created this issue just so that we are aware of those differences, and in case we want to change the behaviour.