Enable slicing for 3d estimation result about geostats.jl HOT 5 CLOSED

What we can do here is generalize the slice operation and make it work with any spatial data:

1. We can define the slice(d, xrange, yrange, ...) as the filter of d with a new GeometryFilter
2. This GeometryFilter stores a geometry from GeometryBasics.jl (e.g. Rectangle)
3. For the case of RegularGrid domains, we can add a fast method with Julia's multiple dispatch

For regular grids, given a rectangle representing the slice, we just need to do some basic calculations with origin, spacing and size to find out what is inside the rectangle. It should be fun to implement, please let me know if you would like to tackle this one.

from geostats.jl.

Happy to take this after we finish what is still pending.

from geostats.jl.

Thanks for offering help @exepulveda , I can take care of this one as well given that I have some free time in the following weeks. Could you please review JuliaEarth/GeoStatsBase.jl#48 which follows a similar pattern. That way you can get used to the source code organization in the project.

from geostats.jl.

I am almost there with a first implementation of the slice function. The GeometryFilter is in place for general spatial data, and I am now adding a specialization for regular grids + rectangle geometry.

from geostats.jl.

On a second thought I think we can achieve a better API by introducing a set of geometric operations as opposed to a filter. Here is the draft: JuliaEarth/GeoStatsBase.jl#49 I will signal when it is ready for review 👍

from geostats.jl.