Rafael Schouten

Ok, it seems there is a bug in `:touches`. Will have a look at fixing it next weekend.

Thanks! That looks a lot closer. I will try tweaking it to see if can be the same, but Im not sure now if the difference is a real problem...

Yes that's interesting. So the problem is occasionally not including a pixel that gdal includes.

Ok I've solved this. The remaining problem is actually in `crop` rather than `rasterize`. That's why the tests against gdal rasterize didn't catch this. The way DimensionalData.jl works selecting areas...

Just noting this isn't actually merged yet so reopening until its fixed and registered

This should be fixed now.

This will be partly handled with lazy `cat` in DiskArrays.jl https://github.com/meggart/DiskArrays.jl/issues/73

Its Neighborhoods.jl, Neighborhood.jl is an unrelated package ;) I need to think of a new name to be able to register it now. Neighborhoods.jl will definately improve the performance. There...

Now registered: https://github.com/rafaqz/Stencils.jl

Hmm I'm tending to agree with you ans @visr now, maybe `GeoInterface.convert` is better.