Anshul Singhvi

That's the plan! It's almost there but there are some indexing order issues.

ah fixed it, was just overcomplicating it in my head! ```julia @time zonal(precip_cellarea_stack; of = all_countries, bylayer = true, spatialslices = true, skipmissing = false, threaded = false) do stack...

Yeah, I'm actually thinking maybe we should have that be a preferences.jl thing you can configure, but that's a later thing

But there's a lot of recompilation time here (50%) and it's probably pretty slow in general, this could use a pass by you @rafaqz to clean up the code and...

here's a test script I wrote btw, some of the things still don't work but this is the idea ```julia using Rasters, RasterDataSources, ArchGDAL using NaturalEarth import Proj precip =...

Added a new test / example script to the top comment. But I'm worried runtimes will have dropped significantly with all this complexity.

Some interesting benchmarks and timer runs: ```julia @be zonal($sum, $(precip); of = $all_countries, spatialslices = $false, progress = $false, threaded = $false) seconds=5 ────────────────────────────────────────────────────────────────────────────────────────────── Time Allocations ─────────────────────── ──────────────────────── Tot /...

Same timing, I checked

Closing since this is a pain to actually do, we can just release this in breaking

Here's what the namedtuple warning looks like: ```julia julia> zonal(mean, st; of = [polygon], bylayer=false) ┌ Warning: The function you passed to `zonal` cannot handle a RasterStack's values directly, │...