return indices as CartesianIndex

[alex-s-gardner]

If indices were returned as CartesianIndex then they can be used to directly index info a geoarray when dealing with Vector of points

[alex-s-gardner]

This is what I'm doing now:

find indices of points

ind = indices.(Ref(ga0), xy);

extract values

ind = [CartesianIndex((ij[1], ij[2])) for ij in ind]
val = ga0[ind]

[evetion]

I think it's a fair point that we should accept and return CartesianIndices, although it would be breaking. Note you can rewrite your code to:

CartesianIndex.(Tuple.(indices.(Ref(ga), xy)))

And you can actually index into the GeoArray directly with floating point numbers: getindex.(Ref(ga), SVector{2}.(xy)).

[evetion]

In GeoArrays 0.8, indices(::GeoArray, x, y) are now returned as CartesianIndex.