Merge pull request #5 from JuliaEarth/natural-earth

Add NaturalEarth backend

Project.toml

@@ -10,6 +10,7 @@ GeoStatsImages = "7cd16168-b42c-5e7d-a585-4f59d326662d"
 GeoTables = "e502b557-6362-48c1-8219-d30d308dcdb0"
 INMET = "1fb021e9-9ab7-46c1-8a23-23f3c2783cf7"
 Meshes = "eacbb407-ea5a-433e-ab97-5258b1ca43fa"
+NaturalEarth = "436b0209-26ab-4e65-94a9-6526d86fea76"
 Unitful = "1986cc42-f94f-5a68-af5c-568840ba703d"
 
 [compat]
@@ -19,5 +20,6 @@ GeoStatsImages = "1.1"
 GeoTables = "1.19"
 INMET = "1.0"
 Meshes = "0.42"
+NaturalEarth = "0.1"
 Unitful = "1.17"
 julia = "1.9"

src/GeoArtifacts.jl

@@ -7,6 +7,7 @@ module GeoArtifacts
 include("image.jl")
 include("gadm.jl")
 include("inmet.jl")
+include("naturalearth.jl")
 
 export GADM, INMET
 

src/naturalearth.jl

@@ -0,0 +1,28 @@
+# -----------------------------------------------------------------
+# Licensed under the MIT License. See LICENSE in the project root.
+# -----------------------------------------------------------------
+
+using GeoIO
+
+import NaturalEarth as NE
+
+"""
+    GeoArtifacts.naturalearth(name::String; fix=true, version = v"5.1.2")
+    GeoArtifacts.naturalearth(name::String, scale::Int; fix=true, version = v"5.1.2")
+
+Load a NaturalEarth dataset as a `GeoTable`.
+
+The `name` should not include the `ne_` prefix, and if providing a `scale`
+should also not include a scale. No suffix should be added.
+
+The option `fix` can be used to fix orientation and degeneracy
+issues with polygons.
+
+# Examples
+
+```julia
+GeoArtifacts.naturalearth("admin_0_countries", 110)
+GeoArtifacts.naturalearth("110m_admin_0_countries")
+```
+"""
+naturalearth(args...; fix=true, kwargs...) = GeoIO.asgeotable(NE.naturalearth(args...; kwargs...), fix)

test/runtests.jl

@@ -32,4 +32,14 @@ using Test
     @test gtb.geometry isa PointSet
     @test embeddim(gtb.geometry) == 3
   end
+
+  @testset "NaturalEarth.jl" begin
+    gtb = GeoArtifacts.naturalearth("admin_0_countries", 110)
+    @test gtb.geometry isa GeometrySet
+    @test embeddim(gtb.geometry) == 2
+
+    gtb = GeoArtifacts.naturalearth("110m_admin_0_countries")
+    @test gtb.geometry isa GeometrySet
+    @test embeddim(gtb.geometry) == 2
+  end
 end