Implementations for custom order NamedTuple.

Author: evetion

src/base.jl

@@ -17,3 +17,10 @@ for i in 2:4
     GeoInterface.ncoord(::PointTrait, geom::sig) = i
     GeoInterface.getcoord(::PointTrait, geom::sig, i) = getindex(geom, i)
 end
+
+# Custom coordinate order/names NamedTuple
+GeoInterface.isgeometry(::Type{<:NamedTuple{Keys,NTuple{N,T}}}) where {Keys,N,T<:Real} = all(in(default_coord_names), Keys)
+GeoInterface.geomtrait(::NamedTuple{Keys,NTuple{N,T}}) where {Keys,N,T<:Real} = PointTrait()
+GeoInterface.ncoord(::PointTrait, geom::NamedTuple{Keys,NTuple{N,T}}) where {Keys,N,T<:Real} = Base.length(geom)
+GeoInterface.getcoord(::PointTrait, geom::NamedTuple{Keys,NTuple{N,T}}, i) where {Keys,N,T<:Real} = getindex(geom, i)
+GeoInterface.coordnames(::PointTrait, geom::NamedTuple{Keys,NTuple{N,T}}) where {Keys,N,T<:Real} = Keys

test/test_primitives.jl

@@ -312,15 +312,20 @@ end
     @testset "NamedTuple" begin
         geom = (; X=1, Y=2)
         @test testgeometry(geom)
+        @test GeoInterface.x(geom) == 1
+        @test collect(GeoInterface.getcoord(geom)) == [1, 2]
+
         geom = (; X=1, Y=2, Z=3)
         @test testgeometry(geom)
         geom = (; X=1, Y=2, Z=3, M=4)
         @test testgeometry(geom)
+        geom = (; Z=3, X=1, Y=2, M=4)
+        @test testgeometry(geom)
 
         @test GeoInterface.x(geom) == 1
         @test GeoInterface.m(geom) == 4
         @test GeoInterface.ncoord(geom) == 4
-        @test collect(GeoInterface.getcoord(geom)) == [1, 2, 3, 4]
+        @test collect(GeoInterface.getcoord(geom)) == [3, 1, 2, 4]
 
     end
 end