🐛 correct check on layer inner type (#146)

* 🐛 correct check on layer inner type

* 🟩 prepare for tests

* ⚙️ handle the internal type check gracefully

* ✅ read/write test passing

Project.toml

@@ -1,7 +1,7 @@
 name = "SpeciesDistributionToolkit"
 uuid = "72b53823-5c0b-4575-ad0e-8e97227ad13b"
 authors = ["Timothée Poisot <timothee.poisot@umontreal.ca>"]
-version = "0.0.3"
+version = "0.0.4"
 
 [deps]
 ArchGDAL = "c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3"

src/io/geotiff.jl

@@ -183,9 +183,10 @@ function _prepare_layer_for_burnin(
     layer::SimpleSDMPredictor{T},
     nodata::T,
 ) where {T <: Number}
-    @assert eltype(layer) <: Number
+    K = SimpleSDMLayers._inner_type(layer)
+    @assert K <: Number
     array = replace(layer.grid, nothing => nodata)
-    array = convert(Matrix{eltype(layer)}, array)
+    array = convert(Matrix{K}, array)
     dtype = eltype(array)
     array_t = reverse(permutedims(array, [2, 1]); dims = 2)
     return array_t

test/01_integration_read.jl

@@ -0,0 +1,20 @@
+module TestThatWeCanReadStuff
+
+using SpeciesDistributionToolkit
+using Test
+
+# Get some EarthEnv data
+layer = SimpleSDMPredictor(RasterData(EarthEnv, LandCover); layer = 1)
+D = SimpleSDMLayers._inner_type(layer)
+
+# Write the data
+f = tempname()
+SpeciesDistributionToolkit._write_geotiff(f, [layer]; driver = "GTiff", nodata = typemax(D))
+@test isfile(f)
+
+# Read the data
+layer2 = SpeciesDistributionToolkit._read_geotiff(f, SimpleSDMPredictor)
+@test typeof(layer2) <: SimpleSDMPredictor
+@test all(values(layer) .== values(layer2))
+
+end

test/Project.toml

@@ -0,0 +1,2 @@
+[deps]
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"

test/runtests.jl

@@ -1 +1,26 @@
-nothing
+using SpeciesDistributionToolkit
+using Test
+
+global anyerrors = false
+
+tests = [
+    "read/write layers" => "01_integration_read.jl",
+]
+
+for test in tests
+    try
+        include(test.second)
+        println("\033[1m\033[32m✓\033[0m\t$(test.first)")
+    catch e
+        global anyerrors = true
+        println("\033[1m\033[31m×\033[0m\t$(test.first)")
+        println("\033[1m\033[38m→\033[0m\ttest/$(test.second)")
+        showerror(stdout, e, backtrace())
+        println()
+        break
+    end
+end
+
+if anyerrors
+    throw("Tests failed")
+end