Michael Sumner
Michael Sumner
the meat of the task is here https://github.com/hypertidy/filearchy/blob/main/R/gdal_tiles.R#L32-L40 I could probably flesh out a basic go at that with osgeo.gdal, if others were interested to help with the broader python...
what's the pythoniest way to generate a table of te and zoom from a bbox and shape? That's the kind of scaffolding I mean, it's a basic reusable component.
Omg, racqet awesome I came here to protest at your ambition, gawd leave Parquet alone it's already perfect! Glad to find racqet though, I see where it encodes the chunks
Nice!! Thanks, I'm having fun with this - appreciate all the detail and functionality here it really helps a (non-native) Python learner.
One thing is that the coordinate values are currently "left"/"top" aligned, not the centre, so here we start at left/top 0,0 and end at 4,4. ```python from xarray.indexes import CoordinateTransformIndex...
On mixed crs collections, suggest looking at the GDAL warper which takes multiple inputs of any kind, and wrappers like odc-geo loader for xarray, and collection formats like GTI, STACIT,...
(Obviously because of the garbage collector this can only be used transiently, there's no point saving a MEM/DATAPOINTER dsn as they go stale pretty quickly)
(of course there's other ways to populate a MEM dataset so maybe this is just a distraction)
Oh yes entirely ~connected~ confected values
The correct values are in this geotiff, fwiw https://github.com/mdsumner/volcano/tree/main/inst/extdata