return tuple points from GeoInterface methods

Author: rafaqz

src/geo_interface.jl

@@ -41,11 +41,40 @@ GeoInterface.getgeom(::MultiLineStringTrait, geom::MultiLineString, i) =
     getGeometry(geom, i)::LineString
 GeoInterface.getgeom(::MultiPolygonTrait, geom::MultiPolygon, i) =
     getGeometry(geom, i)::Polygon
-GeoInterface.getgeom(
-    ::Union{LineStringTrait,LinearRingTrait},
-    geom::Union{LineString,LinearRing},
-    i,
-) = Point(getPoint(geom, i))
+function GeoInterface.getgeom(::AbstractGeometryTrait, geom::Union{LineString,LinearRing}, i)
+    ref = Ref{Float64}()
+    seq = getCoordSeq(geom::Union{LineString, LinearRing})
+    _get_tuple_point(geom, seq, ref, i)
+end
+function GeoInterface.getgeom(::AbstractGeometryTrait, geom::Union{LineString,LinearRing})
+    n = GeoInterface.ngeom(geom)
+    seq = getCoordSeq(geom::Union{LineString, LinearRing})
+    ref = Ref{Float64}()
+    return (_get_tuple_point(geom, seq, ref, i) for i in 1:n)
+ end
+
+function _get_tuple_point(geom, seq, ref, i)
+    context = get_global_context()
+    if GeoInterface.is3d(geom)
+        GC.@preserve ref begin
+            GEOSCoordSeq_getX_r(context, seq, i - 1, ref)
+            x = ref[]
+            GEOSCoordSeq_getY_r(context, seq, i - 1, ref)
+            y = ref[]
+            GEOSCoordSeq_getZ_r(context, seq, i - 1, ref)
+            z = ref[]
+        end
+        return (x, y, z)
+    else
+        GC.@preserve ref begin
+            GEOSCoordSeq_getX_r(context, seq, i - 1, ref)
+            x = ref[]
+            GEOSCoordSeq_getY_r(context, seq, i - 1, ref)
+            y = ref[]
+        end
+        return (x, y)
+    end
+end
 GeoInterface.getgeom(t::AbstractPointTrait, geom::PreparedGeometry) = nothing
 function GeoInterface.getgeom(::PolygonTrait, geom::Polygon, i::Int)
     if i == 1

test/test_geo_interface.jl

@@ -65,9 +65,10 @@ const LG = LibGEOS
     @test GeoInterface.geomtrait(ls) == LineStringTrait()
     p = GeoInterface.ngeom(ls) == 4
     p = GeoInterface.getgeom(ls, 3)
-    @test p isa LibGEOS.Point
+    @test p isa Tuple{Float64,Float64}
     @test GeoInterface.coordinates(p) == [9, 2]
     @test GeoInterface.testgeometry(ls)
+    @test collect(GeoInterface.getpoint(ls)) == [(8, 1), (9, 1), (9, 2), (8, 2)]
     plot(ls)
 
     @inferred GeoInterface.ncoord(ls)
@@ -98,9 +99,10 @@ const LG = LibGEOS
     @test GeoInterface.geomtrait(lr) == LinearRingTrait()
     @test GeoInterface.ngeom(lr) == 5
     p = GeoInterface.getgeom(lr, 3)
-    @test p isa LibGEOS.Point
+    @test p isa Tuple{Float64,Float64}
     @test GeoInterface.coordinates(p) == [9, 2]
     @test GeoInterface.testgeometry(lr)
+    @test collect(GeoInterface.getpoint(lr)) == [(8, 1), (9, 1), (9, 2), (8, 2), (8, 1)]
     # Cannot convert LinearRingTrait to series data for plotting
     # plot(lr)