Add methods to operations like disjoint, crosses etc to accept arbitrary GeoInterface.jl compatible geometries and convert them to ArchGDAL geometries with the shared to_gdal method and geointerface_geomtype.
The returned values are always still ArchGDAL geometries.
We don't use convert(ArchGDAL, geom) here as it has a small overhead to call a pacage method, which is unnecessary to pay locally - especially for points of anything there are a lot of.
rafaqz
commented
1 year ago
Add methods to operations like
disjoint,crossesetc to accept arbitrary GeoInterface.jl compatible geometries and convert them to ArchGDAL geometries with the sharedto_gdalmethod andgeointerface_geomtype.The returned values are always still
ArchGDALgeometries.We don't use
convert(ArchGDAL, geom)here as it has a small overhead to call a pacage method, which is unnecessary to pay locally - especially for points of anything there are a lot of.Needs tests, and https://github.com/JuliaGeo/GeoInterface.jl/pull/78