Add implementations for Base types.

Author: evetion

docs/src/guides/defaults.md

@@ -45,3 +45,17 @@ npoint(::QuadTrait, geom) = 4
 npoint(::PentagonTrait, geom) = 5
 npoint(::HexagonTrait, geom) = 6
 ```
+
+# Implementations
+GeoInterface is implemented automatically for `NTuple`s, `NamedTuple`s and `Vector`s to behave as Points.
+
+```julia
+a = [1,2,3]
+x(a) == 1
+
+b = (1,2,3)
+x(b) == 2
+
+c = (;X=1, Y=2)
+x(c) == 1
+```

src/GeoInterface.jl

@@ -45,5 +45,6 @@ include("types.jl")
 include("interface.jl")
 include("fallbacks.jl")
 include("utils.jl")
+include("base.jl")
 
 end # module

src/base.jl

@@ -0,0 +1,19 @@
+# Implementation of GeoInterface for Base Types
+
+GeoInterface.isgeometry(::Type{<:AbstractVector{<:Real}}) = true
+GeoInterface.geomtrait(::AbstractVector{<:Real}) = PointTrait()
+GeoInterface.ncoord(::PointTrait, geom::AbstractVector{<:Real}) = Base.length(geom)
+GeoInterface.getcoord(::PointTrait, geom::AbstractVector{<:Real}, i) = getindex(geom, i)
+
+GeoInterface.isgeometry(::Type{<:NTuple{N,<:Real}}) where {N} = true
+GeoInterface.geomtrait(::NTuple{N,<:Real}) where {N} = PointTrait()
+GeoInterface.ncoord(::PointTrait, geom::NTuple{N,<:Real}) where {N} = N
+GeoInterface.getcoord(::PointTrait, geom::NTuple{N,<:Real}, i) where {N} = getindex(geom, i)
+
+for i in 2:4
+    sig = NamedTuple{default_coord_names[1:i],NTuple{i,T}} where {T<:Real}
+    GeoInterface.isgeometry(::Type{<:sig}) = true
+    GeoInterface.geomtrait(::sig) = PointTrait()
+    GeoInterface.ncoord(::PointTrait, geom::sig) = i
+    GeoInterface.getcoord(::PointTrait, geom::sig, i) = getindex(geom, i)
+end

src/utils.jl

@@ -1,7 +1,7 @@
 """Test whether the required interface for your `geom` has been implemented correctly."""
 function testgeometry(geom)
     try
-        @assert isgeometry(geom)
+        @assert isgeometry(geom) "Geom doesn't implement `isgeometry`."
         type = geomtrait(geom)
 
         if type == PointTrait()
@@ -18,7 +18,7 @@ function testgeometry(geom)
                     issub = geomtrait(g2) isa subtype
                     !issub && error("Implemented hierarchy for this geometry type is incorrect. Subgeometry should be a $subtype")
                 end
-                @assert testgeometry(g2)  # recursive testing of subgeometries
+                @assert testgeometry(g2) "Subgeometry implementation is not valid."
             end
         end
     catch e

test/test_primitives.jl

@@ -290,3 +290,37 @@ end
     @test GeoInterface.astext(geom) isa String
     @test GeoInterface.asbinary(geom) isa Vector{UInt8}
 end
+
+@testset "Base Implementations" begin
+
+    @testset "Vector" begin
+        geom = [1, 2]
+        @test testgeometry(geom)
+        @test GeoInterface.x(geom) == 1
+        @test GeoInterface.ncoord(geom) == 2
+        @test collect(GeoInterface.getcoord(geom)) == geom
+    end
+
+    @testset "Tuple" begin
+        geom = (1, 2)
+        @test testgeometry(geom)
+        @test GeoInterface.x(geom) == 1
+        @test GeoInterface.ncoord(geom) == 2
+        @test collect(GeoInterface.getcoord(geom)) == [1, 2]
+    end
+
+    @testset "NamedTuple" begin
+        geom = (; X=1, Y=2)
+        @test testgeometry(geom)
+        geom = (; X=1, Y=2, Z=3)
+        @test testgeometry(geom)
+        geom = (; X=1, Y=2, Z=3, 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]
+
+    end
+end