Should Geomorphometry.jl also include the ability for sparse calculations?

[alex-s-gardner]

I have a global high-resolution DEM (cop30) and a ton of profile data (ICESat-2) and I'd like to calculate the slope and curvature at each ICESat-2 measurement point. I could do this by calculating slope and curvature everywhere but it would be much faster (less allocations) to simply calculate slope and curvature at each point. This is a utility that I often need but have not found.

Looking into Stencils.jl it seems like this shouldn't be too difficult... the steps would look like this

1. lazily read in profile data (probably Rasters.jl)
2. eagerly read data cropped to extents of profile data + a buffer for calculation of slope, curvature, ect
3. cast copped data to a SencilArray = A
4. loop though each intersecting index of A and apply functions

I'm guessing that this is functionality that others have built or would want but it's not clear where such functionality should live. Given that Geomorphometry.jl is home to most of the funtions I'm wondering if it should live here.

Thoughts?

[alex-s-gardner]

Closing as this is probably better suited for Slack