You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
I am currently testing GeometryOps package to move away from R and the sf package for GIS work. I am facing a minor issue (misusage or bug ?) while using the spatial join function.
The following example with GO.within is working as expected :
using FlexiJoins, DataFrames
import GeoInterface as GI
import GeometryOps as GO
pl = GI.Polygon([GI.LinearRing([(0, 0), (1, 0), (1, 1), (0, 0)])])
pu = GI.Polygon([GI.LinearRing([(0, 0), (0, 1), (1, 1), (0, 0)])])
poly_df = DataFrame(geometry = [pl, pu], color = [:red, :blue])
points = tuple.(rand(1000), rand(1000))
points_df = DataFrame(geometry = points)
joined_df_within = FlexiJoins.innerjoin(
(points_df, poly_df),
by_pred(:geometry, GO.within, :geometry)
)
I am currently testing
GeometryOps
package to move away fromR
and thesf
package for GIS work. I am facing a minor issue (misusage or bug ?) while using the spatial join function.The following example with
GO.within
is working as expected :But the
GO.intersects
predicate fails :Giving the following error:
While the predicate is otherwise working as expected
The versions used are the following :
Thanks
The text was updated successfully, but these errors were encountered: