Add image_from_paths/image_from_paths!

Project.toml

@@ -1,7 +1,7 @@
 name = "Thermochron"
 uuid = "b9a8379e-ee1a-4388-b7ca-7852af1cef08"
 authors = ["C. Brenhin Keller and collaborators"]
-version = "0.11.1"
+version = "0.11.2"
 
 [deps]
 BangBang = "198e06fe-97b7-11e9-32a5-e1d131e6ad66"
@@ -30,5 +30,5 @@ LsqFit = "0.15.0"
 NaNStatistics = "0.6"
 ProgressMeter = "1"
 Reexport = "1"
-StatGeochemBase = "0.4.3, 0.5, 0.6"
+StatGeochemBase = "0.6.11"
 julia = "1.11"

src/inversion.jl

@@ -771,3 +771,36 @@
         )
         return ttresult, kineticresult
     end
+
+
+    """
+    ```julia
+    image_from_paths(tT::TTResult; 
+        xresolution::Int=1800, 
+        yresolution::Int=1200, 
+        xrange = nanextrema(tT.tpointdist), 
+        yrange = nanextrema(tT.Tpointdist), 
+    )
+    ```
+    Produce a 2d image (histogram) of path densities given a `TTResult`
+    
+    ```julia
+    julia> imgcounts, xq, yq = image_from_paths(tT; xrange=boundary.agepoints, yrange=boundary.T₀)
+    ```
+    """
+    function StatGeochemBase.image_from_paths(tT::TTResult;
+            xresolution::Int = 1800, 
+            yresolution::Int = 1200, 
+            xrange = nanextrema(tT.tpointdist), 
+            yrange = nanextrema(tT.Tpointdist), 
+        )
+        image_from_paths(tT.tpointdist, tT.Tpointdist; xresolution, yresolution, xrange, yrange)
+    end
+    function StatGeochemBase.image_from_paths!(tT::TTResult;
+            xresolution::Int = 1800, 
+            yresolution::Int = 1200, 
+            xrange = nanextrema(tT.tpointdist), 
+            yrange = nanextrema(tT.Tpointdist), 
+        )
+        image_from_paths!(tT.tpointdist, tT.Tpointdist; xresolution, yresolution, xrange, yrange)
+    end

test/testcomplete.jl

@@ -432,4 +432,22 @@ DN17Eamean = mean(kinetics.zdmdist .|> x-> x.DN17Ea)
 @info "Mean zircon DN17Ea: $DN17Eamean"
 rmr0mean = mean(kinetics.zdmdist .|> x-> x.rmr0)
 @test 0 < rmr0mean < 1
-@info "Mean zircon rmr0: $rmr0mean"
+@info "Mean zircon rmr0: $rmr0mean"
+
+## --- Test conversion of t-T images 
+
+ttimage, xq, yq = image_from_paths(tT, xresolution=200, yresolution=100, yrange=boundary.T₀)
+@test xq == 8.75:17.5:3491.25
+@test yq == 2.0:4.0:398.0
+@test ttimage isa Matrix
+@test axes(ttimage, 1) == 1:100
+@test axes(ttimage, 2) == 1:200
+
+ttimage2, xq2, yq2 = image_from_paths!(tT, xresolution=200, yresolution=100, yrange=boundary.T₀)
+@test xq2 == 8.75:17.5:3491.25
+@test yq2 == 2.0:4.0:398.0
+@test ttimage2 isa Matrix
+@test axes(ttimage2, 1) == 1:100
+@test axes(ttimage2, 2) == 1:200
+
+@test ttimage == ttimage2