Implement crs for AbstractGeoTable

JuliaEarth · Nov 7, 2024 · 216cf04 · 216cf04
1 parent 03a1cdb
commit 216cf04
Show file tree

Hide file tree

Showing 3 changed files with 8 additions and 0 deletions.
diff --git a/src/GeoTables.jl b/src/GeoTables.jl
@@ -23,6 +23,7 @@ import DataAPI: nrow, ncol
 import Meshes: partitioninds
 import Meshes: sampleinds
 import Meshes: sortinds
+import Meshes: crs
 import TransformsBase: apply, revert, reapply
 
 # required for VSCode table viewer

diff --git a/src/abstractgeotable.jl b/src/abstractgeotable.jl
@@ -112,6 +112,8 @@ function Base.parentindices(geotable::AbstractGeoTable)
   end
 end
 
+crs(geotable::AbstractGeoTable) = crs(domain(geotable))
+
 # -----------
 # IO METHODS
 # -----------

diff --git a/test/basics.jl b/test/basics.jl
@@ -41,6 +41,11 @@
     @test parentindices(dat) == inds
     @test isnothing(values(pdat))
 
+    # coordinate reference system (crs)
+    dom = CartesianGrid(2, 2)
+    dat = dummy((a=[1, 2, 3, 4], b=[5, 6, 7, 8]), dom)
+    @test crs(dat) == crs(dom)
+
     # Tables interface
     dom = CartesianGrid(2, 2)
     dat = dummy((a=[1, 2, 3, 4], b=[5, 6, 7, 8]), dom)