Add geometric transforms to extension

ext/GeoTablesTableTransformsExt.jl

@@ -10,12 +10,19 @@ using TableTransforms
 
 import TableTransforms: divide, attach
 import TableTransforms: applymeta, revertmeta
+import TableTransforms: apply, revert, reapply
 
-# table traits
+# -----------------------
+# TABLE TRANSFORM TRAITS
+# -----------------------
 
 divide(geotable::AbstractGeoTable) = values(geotable), domain(geotable)
 attach(table, dom::Domain) = georef(table, dom)
 
+# -------------------
+# FEATURE TRANSFORMS
+# -------------------
+
 # transforms that change the order or number of
 # rows in the table need a special treatment
 
@@ -108,4 +115,24 @@ function revertmeta(::Sample, newdom::Domain, mcache)
   GeometrySet(ugeoms)
 end
 
+# ---------------------
+# GEOMETRIC TRANSFORMS
+# ---------------------
+
+function apply(transform::GeometricTransform, geotable::AbstractGeoTable)
+  newdom, cache = apply(transform, domain(geotable))
+  newdata = georef(values(geotable), newdom)
+  newdata, cache
+end
+
+function revert(transform::GeometricTransform, newdata::AbstractGeoTable, cache)
+  dom = revert(transform, domain(newdata), cache)
+  georef(values(newdata), dom)
+end
+
+function reapply(transform::GeometricTransform, geotable::AbstractGeoTable, cache)
+  newdom = reapply(transform, domain(geotable), cache)
+  georef(values(geotable), newdom)
+end
+
 end

test/geometric.jl

@@ -0,0 +1,10 @@
+@testset "Geometric transforms" begin
+  d = georef((z=rand(100), w=rand(100)))
+  p = StdCoords()
+  n, c = apply(p, d)
+  dom = domain(n)
+  cen = centroid.(dom)
+  xs = first.(coordinates.(cen))
+  @test dom isa CartesianGrid
+  @test all(x -> -0.5 ≤ x ≤ 0.5, xs)
+end

test/runtests.jl

@@ -22,6 +22,7 @@ testfiles = [
   "indices.jl",
   "operations.jl",
   "feature.jl",
+  "geometric.jl",
   "shows.jl"
 ]