Rafael Schouten

You could do this manually in a loop over each polygon using `coverage`. But yes it would be nice to automate that process. We can just call `coverage` instead of...

We don't really let DiskArrays figure it out for `extract` unfortunately, it's one by one. For the memory backed alg it's optimised so far that allocations are the main cost,...

Master of both or just here?

Absolutely. gdal just proved a lot of options it will be hard to replicate immediately. Especially warping between any projections. We could treat it like GEOS() in GeometryOps to have...

You have to do a nested rebuild because `rebuild` is very low level. `set` will do what you want like `set(dix, somerange)`. It handles the nesting for you.

I think we need to make a plan for set/rebuild etc... I'm hoping one function can do everything. Like we could allow keywords in `set` so it can do what...

May be a bug! I will have a look. But note that hand-rolling Shapefile.jl geometries is not the best test here, they are not made to be manually defined like...

Ok yes there is a serious bug with polygon rasterization in general for this example (its not just the Shapefile and not just `mask`). I'm not sure how this got...

This was fixed

Looks like you having `missing` values and the alg can't handle that? You can maybe replace them with zeros if they are e.g. the sea. `Raster(filename; missingval=0)` should do that....