enumerate(GI.getpoint(geom)) is not type inferred on polygons

[asinghvi17]

julia> eltype(enumerate(GI.getpoint(p56)))
Tuple{Int64, Any}

julia> const p5 = p56
GeoInterface.Wrappers.Polygon{false, false, Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, Nothing, Nothing}(GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}[GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}([(0.0, 0.0), (0.0, 10.0), (10.0, 10.0), (10.0, 0.0)], nothing, nothing)], nothing, nothing)

julia> eltype(enumerate(GI.getpoint(p5)))
Tuple{Int64, Any}

julia> typeof(enumerate(GI.getpoint(p5)))
Base.Iterators.Enumerate{Base.Iterators.Flatten{Base.Generator{Vector{GeoInterface.Wrappers.LinearRing{false, false, Vector{Tuple{Float64, Float64}}, Nothing, Nothing}}, GeoInterface.var"#6#7"}}}

Not sure how bad this is for performance - maybe we need a point enumerator wrapper type thing which is known to be type stable (and maybe also known to emit tuple/svector type points)?

[rafaqz]

What about without enumerate

[asinghvi17]

That's a vector so type stability is fine I think

Ah, there's the problem. I guess getpoint for polygons returns a Iterators.Flatten which has eltype of Any.

getpoint(GI.getexterior(p56)) is a vector so no worries about eltype there.

[rafaqz]

Yeah it's double nested for Polygons. We could make a better iterator in GeoInterface.