Rafael Schouten

Yeah JLD is a pretty bad data format. We are allowed to break your files even on minor version bumps, and its julia only. Kind of worst of all worlds....

`zonal` should be doing a broadcast to read. But we don't test performance on DiskArrays. We just need to go over all the methods and disallow scalar indexing and test...

@asinghvi17 Don't you want a view? That disk array will be materialised. But yeah, RasterSeries should just do that automatically if you pass an empty Ti() dim

Good idea we could totally detect singletons

We have `GI.geometrycolumns` to clarify these in many cases. GeometryOps should use it everywhere, we don't use `:geom` ? except in one tutorial because GADM.jl did. `:geometry` is the standard...

We should as much as possible standardise table columns, it's worth the long term simplicity

Yeah... good demonstratiion that `first(GI.geometrycolumns(geomA))` is waay too long to actually use

Where is it still :geom? This is just GDAL? There is not much more we can do here, GeoInterface has both a standard name (:geometry) and a function to get...

Right yes that's annoying. Do you want a zero? Maybe an `init` keyword would help. Assumed to be some kind of identity like 0 for `sum` or 1 for `prod`,...

Yes so in this case using `init` in `maximum` fixes it. Should we allow passing that? It pretty simple to do. If set to `nothing` we won't pass it in.