returning subset of raster data via API
Hi all: as part of working on OGC API - Coverages support in pygeoapi, I'm writing the equivalent of what WCS GetCoverage would do; that is, clip by bbox and subset by bands if requested.
An initial implementation can be found in this branch.
At the moment, the data is returned as JSON-ified ndarrays. I would like to return the data
in its native format. What would be a general approach here? Would one use a