From 393a7b78beb337186d2baa6517b964470d6d0f31 Mon Sep 17 00:00:00 2001 From: Quarto GHA Workflow Runner Date: Fri, 7 Feb 2025 14:49:54 +0000 Subject: [PATCH] Built site for gh-pages --- .nojekyll | 2 +- 01-geodata.html | 60 +- 02-geoviz.html | 52 +- 04-geoproc.html | 48 +- 05-transforms.html | 1368 +++++++++++++++++++++---------------------- 07-pipelines.html | 244 ++++---- 10-correlation.html | 52 +- 12-mining.html | 486 +++++++-------- 14-energy.html | 64 +- index.html | 4 +- references.html | 8 +- search.json | 26 +- 12 files changed, 1207 insertions(+), 1207 deletions(-) diff --git a/.nojekyll b/.nojekyll index 3725e04..e55159e 100644 --- a/.nojekyll +++ b/.nojekyll @@ -1 +1 @@ -c9e3b98b \ No newline at end of file +8479b378 \ No newline at end of file diff --git a/01-geodata.html b/01-geodata.html index e6b5d84..f5240f4 100644 --- a/01-geodata.html +++ b/01-geodata.html @@ -1128,63 +1128,63 @@

-0.216003 -0.60917 -0.863141 +0.722209 +0.613786 +0.0994426 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -0.691252 -0.937862 -0.225594 +0.438343 +0.305785 +0.239826 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -0.820834 -0.757224 -0.682132 +0.431154 +0.148168 +0.29044 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -0.697256 -0.45779 -0.66697 +0.371207 +0.472349 +0.0751472 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -0.996161 -0.969452 -0.278064 +0.597996 +0.242541 +0.539035 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -0.860413 -0.209935 -0.709269 +0.205348 +0.68878 +0.0779545 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -0.864347 -0.75301 -0.0745924 +0.965599 +0.858176 +0.0992563 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -0.00386113 -0.704749 -0.387572 +0.297715 +0.513887 +0.307331 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -0.634169 -0.266041 -0.344567 +0.448824 +0.0499551 +0.959651 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -0.286819 -0.0207844 -0.979672 +0.101168 +0.510622 +0.126042 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) diff --git a/02-geoviz.html b/02-geoviz.html index 41e1283..06eb719 100644 --- a/02-geoviz.html +++ b/02-geoviz.html @@ -426,53 +426,53 @@

-0.782612 -0.782563 +0.572878 +0.7025 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -0.439685 -0.630608 +0.573844 +0.660015 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -0.12363 -0.234298 +0.332714 +0.864078 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -0.483283 -0.545023 +0.708591 +0.187223 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -0.161346 -0.0687522 +0.300084 +0.86401 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -0.622278 -0.743368 +0.127223 +0.838724 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -0.228983 -0.177852 +0.621324 +0.992886 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -0.786215 -0.718279 +0.917433 +0.823236 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -0.79631 -0.198231 +0.337756 +0.547082 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -0.238591 -0.106484 +0.509657 +0.792821 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -520,14 +520,14 @@ 

(A=rand(10, 10),)
-
(A = [0.29400635786174834 0.30136313860386343 … 0.850396069835956 0.7564893677059472; 0.5891600875088715 0.4710638478717837 … 0.734748856127191 0.03945956182772836; … ; 0.3117776524678104 0.6773628636686471 … 0.5886532490854754 0.0847992968801613; 0.9116939753552071 0.3382451121210255 … 0.26891034974686256 0.35595459123777595],)
+
(A = [0.5384948124207134 0.23613853906095217 … 0.6948648644828772 0.2569801356859971; 0.30206326166511155 0.8158499910806138 … 0.0002703007746212016 0.5899112203913356; … ; 0.040296573964802684 0.22863706060806555 … 0.15509044334683908 0.011624346969265065; 0.25392327182804386 0.9228296844829125 … 0.2233308962293422 0.3774488045745531],)

or a semicolon before the key/value pair:

(; A=rand(10, 10))
-
(A = [0.10133496725296853 0.9594888116965612 … 0.5887708010722427 0.682060769727528; 0.38405261827877735 0.7897336411274214 … 0.3162251703014688 0.4338471138926461; … ; 0.2642404473597736 0.46550378133843406 … 0.3115977150115633 0.868080255067318; 0.8606586577849369 0.5922089191891257 … 0.04843178973336981 0.04806460399423795],)
+
(A = [0.2740340754330898 0.43711818676827296 … 0.7408610136843572 0.9857691593592899; 0.7167090508413428 0.8990187229633196 … 0.28885272737610834 0.6587513422096297; … ; 0.4946076917399742 0.6139162989809983 … 0.6229704477933816 0.28295025509727023; 0.2588356410933321 0.9309814188120664 … 0.28515666211267054 0.5225756553414658],)
@@ -543,7 +543,7 @@ 

viz(img.geometry, color = img.A)
- +

but any vector with the same length can be passed:

@@ -557,7 +557,7 @@ 

viz(img.geometry, color = img.B, alpha = rand(length(img.B)))
- +

Other aesthetic options are available in the official documentation.

@@ -682,7 +682,7 @@

viewer(geotable)
- +
@@ -704,7 +704,7 @@

geotable |> viewer

- +
diff --git a/04-geoproc.html b/04-geoproc.html index 0373370..9b45e93 100644 --- a/04-geoproc.html +++ b/04-geoproc.html @@ -855,26 +855,26 @@ 

points = rand(Point, 100, crs=Cartesian2D)

100-element Vector{Point{𝔼{2}, Cartesian2D{NoDatum, Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}}}}:
- Point(x: 0.9389356104983692 m, y: 0.15296281346016882 m)
- Point(x: 0.1460818975888738 m, y: 0.43648723644532483 m)
- Point(x: 0.23742231595856478 m, y: 0.6405291044150762 m)
- Point(x: 0.7206559019000279 m, y: 0.4861469744426379 m)
- Point(x: 0.24481125599776943 m, y: 0.7805322875738824 m)
- Point(x: 0.33603743972132105 m, y: 0.7238618945957376 m)
- Point(x: 0.7631476851032564 m, y: 0.2100979048381364 m)
- Point(x: 0.6525703942632162 m, y: 0.10843645938992619 m)
- Point(x: 0.6189305175858374 m, y: 0.35996879268325865 m)
- Point(x: 0.5217370033068044 m, y: 0.14188801197837198 m)
+ Point(x: 0.8756947465880447 m, y: 0.8990537137809087 m)
+ Point(x: 0.1956765101666681 m, y: 0.8336818649334522 m)
+ Point(x: 0.8279897329569817 m, y: 0.11568228336009967 m)
+ Point(x: 0.8288674236423157 m, y: 0.4450355464124124 m)
+ Point(x: 0.06472372246097335 m, y: 0.4379809005496862 m)
+ Point(x: 0.15805967399021315 m, y: 0.16472526591813264 m)
+ Point(x: 0.4314496514072652 m, y: 0.6630688294104694 m)
+ Point(x: 0.04836356906120487 m, y: 0.3314375502895359 m)
+ Point(x: 0.8758800547852885 m, y: 0.008928352674892448 m)
+ Point(x: 0.3636918859743625 m, y: 0.04446603781468206 m)
  ⋮
- Point(x: 0.33840693612728334 m, y: 0.4212174755989695 m)
- Point(x: 0.5495477165134407 m, y: 0.6636507804393046 m)
- Point(x: 0.7138898295553878 m, y: 0.7263374727279495 m)
- Point(x: 0.7114161255008237 m, y: 0.6635461670537387 m)
- Point(x: 0.43454266510211725 m, y: 0.6218504702122563 m)
- Point(x: 0.6279184723857204 m, y: 0.8546149636650487 m)
- Point(x: 0.2953485404647025 m, y: 0.4076587711193762 m)
- Point(x: 0.3919943999447002 m, y: 0.350106902402918 m)
- Point(x: 0.7917676998641788 m, y: 0.08505674384761841 m)
+ Point(x: 0.1285997162616831 m, y: 0.07640123979800262 m) + Point(x: 0.4966182879547494 m, y: 0.6192425829200147 m) + Point(x: 0.05800198309024951 m, y: 0.6370632062270113 m) + Point(x: 0.594366730491882 m, y: 0.6059051267636071 m) + Point(x: 0.8749539071271232 m, y: 0.4254424482751783 m) + Point(x: 0.22377311925764853 m, y: 0.926559318968934 m) + Point(x: 0.011423450893109588 m, y: 0.4393450751380037 m) + Point(x: 0.9981017385417039 m, y: 0.5329531489106377 m) + Point(x: 0.2989786219945655 m, y: 0.9157438575127722 m)
@@ -892,9 +892,9 @@ 

rand(Triangle, 3)

3-element Vector{Triangle{𝔼{3}, Cartesian3D{NoDatum, Quantity{Float64, 𝐋, Unitful.FreeUnits{(m,), 𝐋, nothing}}}}}:
- Triangle((x: 0.926907 m, y: 0.790056 m, z: 0.232759 m), (x: 0.886033 m, y: 0.374097 m, z: 0.911197 m), (x: 0.503991 m, y: 0.937864 m, z: 0.913291 m))
- Triangle((x: 0.523967 m, y: 0.582064 m, z: 0.294668 m), (x: 0.552555 m, y: 0.64402 m, z: 0.939042 m), (x: 0.969351 m, y: 0.167959 m, z: 0.92495 m))
- Triangle((x: 0.0891154 m, y: 0.982618 m, z: 0.271923 m), (x: 0.959133 m, y: 0.628085 m, z: 0.853343 m), (x: 0.700139 m, y: 0.848431 m, z: 0.636408 m))
+ Triangle((x: 0.661352 m, y: 0.815234 m, z: 0.587243 m), (x: 0.763481 m, y: 0.332212 m, z: 0.0666455 m), (x: 0.233197 m, y: 0.724928 m, z: 0.489059 m)) + Triangle((x: 0.0395623 m, y: 0.00905283 m, z: 0.951567 m), (x: 0.0941487 m, y: 0.477524 m, z: 0.693062 m), (x: 0.503129 m, y: 0.963429 m, z: 0.268928 m)) + Triangle((x: 0.524629 m, y: 0.822033 m, z: 0.160219 m), (x: 0.732884 m, y: 0.637768 m, z: 0.559781 m), (x: 0.040715 m, y: 0.602569 m, z: 0.420415 m))
@@ -904,7 +904,7 @@

viz!(points, color = "black") Mke.current_figure()
- +
@@ -912,7 +912,7 @@

viz!(points, color = "black") Mke.current_figure()

- +

Geometries can be discretized into geospatial domains (i.e., collections of geometries). We have seen some of these domains in previous chapters, including the GeometrySet and the CartesianGrid. Here, we will focus on the Mesh subtypes:

diff --git a/05-transforms.html b/05-transforms.html index 0c42715..24d2195 100644 --- a/05-transforms.html +++ b/05-transforms.html @@ -421,73 +421,73 @@

-8.66638 -1.36663 -1.09848 -1.58887 +6.97666 +-5.10709 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 +8.4477 +-2.07371 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 +6.34579 +6.50635 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 +6.03157 +-5.25486 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 +4.81221 +1.04529 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 +3.42317 +-0.369744 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 +7.45059 +-1.6694 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 +8.85509 +0.804656 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 +7.25952 +0.719406 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 +8.70975 +-4.7497 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -505,14 +505,14 @@ 

Mke.hist(gtb.a, color = "gray80")
- +

We can compute the cross-correlation between columns “a” and “b”:

cor(gtb.a, gtb.b)
-
0.0011106632198854395
+
0.008062003117630903

And inspect bivariate distributions of the values of the geotable with PairPlots.jl by Thompson (2023):

@@ -521,7 +521,7 @@

pairplot(values(gtb))
- +

This pattern is useful to answer geoscientific questions via marginal analysis (i.e. entire columns treated as measurements of a single random variable). However, the answers to many questions in geosciences depend on where the measurements were made.

@@ -530,26 +530,26 @@ 

gtb.a
10000-element Vector{Float64}:
-  8.666381809751455
-  7.7010119957059215
-  2.5537153446068386
-  9.048989853817329
-  5.351979872681088
-  5.9383574448864005
-  5.773195011732517
-  3.779320608778172
-  4.772167719345372
-  5.3804908475885975
+  6.976656168217884
+  8.447697767753994
+  6.345789678893408
+  6.031566300246763
+  4.8122133315296125
+  3.423170264528298
+  7.450592926414082
+  8.855092841560099
+  7.259524396222013
+  8.709749432872337
   ⋮
-  0.32149740360756635
- -1.160128995053951
-  1.468191445528531
- -1.2616041401088116
- -0.05702398706854036
- -1.449646415698771
- -0.7673014401364087
-  0.9336966733365696
-  0.43034983850403313
+ -0.26002378155797246 + 2.3798159192086064 + -0.00879557675793819 + -0.5927430966461085 + -0.983600170557722 + -1.4299830648536611 + -0.36674400878799396 + -0.10058903215138397 + 0.8615869982051558
@@ -602,53 +602,53 @@

-8.66638 -1.36663 +6.97666 +-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 +8.4477 +-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 +6.34579 +6.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 +6.03157 +-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 +4.81221 +1.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 +3.42317 +-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 +7.45059 +-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 +8.85509 +0.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 +7.25952 +0.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 +8.70975 +-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -688,63 +688,63 @@

-8.66638 -1.36663 -1.09848 +6.97666 +-5.10709 +-0.647744 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 +8.4477 +-2.07371 +1.02694 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 +6.34579 +6.50635 +-0.291667 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 +6.03157 +-5.25486 +-1.21931 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 +4.81221 +1.04529 +0.7528 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 +3.42317 +-0.369744 +-0.409321 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 +7.45059 +-1.6694 +0.828357 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 +8.85509 +0.804656 +-0.646786 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 +7.25952 +0.719406 +-0.505269 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 +8.70975 +-4.7497 +-0.847882 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -784,63 +784,63 @@

-1.36663 -1.09848 -1.58887 +-5.10709 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -3.37709 -0.257203 -0.353562 +-2.07371 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -4.78739 --0.165323 -0.554613 +6.50635 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) --3.02341 --1.7829 --1.82417 +-5.25486 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) --3.48527 -0.129093 --0.730952 +1.04529 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -1.73545 -0.309309 -1.17752 +-0.369744 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) --0.510192 -0.597076 --0.188414 +-1.6694 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) --7.55738 --0.071525 --0.00570931 +0.804656 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -1.09176 --0.158519 -0.0301378 +0.719406 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -4.46788 --0.62577 --0.815971 +-4.7497 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -878,53 +878,53 @@

-8.66638 -1.36663 +6.97666 +-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 +8.4477 +-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 +6.34579 +6.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 +6.03157 +-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 +4.81221 +1.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 +3.42317 +-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 +7.45059 +-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 +8.85509 +0.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 +7.25952 +0.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 +8.70975 +-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -964,63 +964,63 @@

-8.66638 -1.09848 -1.58887 +6.97666 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -0.257203 -0.353562 +8.4477 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 --0.165323 -0.554613 +6.34579 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --1.7829 --1.82417 +6.03157 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 -0.129093 --0.730952 +4.81221 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -0.309309 -1.17752 +3.42317 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 -0.597076 --0.188414 +7.45059 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --0.071525 --0.00570931 +8.85509 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 --0.158519 -0.0301378 +7.25952 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 --0.62577 --0.815971 +8.70975 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1060,7 +1060,7 @@ 

gtb |> Select("a") |> viewer
- +

@@ -1096,73 +1096,73 @@

-8.66638 -1.36663 -1.09848 -1.58887 +6.97666 +-5.10709 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 +8.4477 +-2.07371 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 +6.34579 +6.50635 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 +6.03157 +-5.25486 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 +4.81221 +1.04529 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 +3.42317 +-0.369744 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 +7.45059 +-1.6694 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 +8.85509 +0.804656 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 +7.25952 +0.719406 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 +8.70975 +-4.7497 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1207,73 +1207,73 @@

-8.66638 -1.36663 -1.09848 -1.58887 +6.97666 +-5.10709 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 +8.4477 +-2.07371 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 +6.34579 +6.50635 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 +6.03157 +-5.25486 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 +4.81221 +1.04529 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 +3.42317 +-0.369744 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 +7.45059 +-1.6694 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 +8.85509 +0.804656 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 +7.25952 +0.719406 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 +8.70975 +-4.7497 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1318,73 +1318,73 @@

-8.66638 -1.36663 -1.09848 -1.58887 +6.97666 +-5.10709 +-0.647744 +-0.778255 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 +8.4477 +-2.07371 +1.02694 +1.2234 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 +6.34579 +6.50635 +-0.291667 +-0.455588 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 +6.03157 +-5.25486 +-1.21931 +-1.98666 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 +4.81221 +1.04529 +0.7528 +2.09002 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 +3.42317 +-0.369744 +-0.409321 +-0.690845 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 +7.45059 +-1.6694 +0.828357 +-0.685142 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 +8.85509 +0.804656 +-0.646786 +-1.27382 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 +7.25952 +0.719406 +-0.505269 +-0.5437 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 +8.70975 +-4.7497 +-0.847882 +-1.49863 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1408,13 +1408,13 @@ 

gtb |> Functional(cos) |> values |> pairplot
- + 
gtb |> Functional("a" => cos, "b" => sin) |> values |> pairplot
- +

The Map transform can be used to create new columns from existing columns in the geotable. It takes a column specification, calls a function on the selected columns row-by-row, and returns the result as a new column:

@@ -1454,93 +1454,93 @@

-8.66638 -1.36663 -1.09848 -1.58887 -0.687758 -0.202752 +6.97666 +-5.10709 +-0.647744 +-0.778255 +0.63921 +0.384536 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 -0.988323 --0.972398 +8.4477 +-2.07371 +1.02694 +1.2234 +0.828867 +-0.481982 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 -0.554596 -0.074927 +6.34579 +6.50635 +-0.291667 +-0.455588 +0.0625635 +0.975202 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 -0.367006 --0.993025 +6.03157 +-5.25486 +-1.21931 +-1.98666 +-0.248972 +0.516253 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 --0.80234 --0.941522 +4.81221 +1.04529 +0.7528 +2.09002 +-0.995022 +0.501654 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 --0.338035 --0.163907 +3.42317 +-0.369744 +-0.409321 +-0.690845 +-0.277871 +0.93242 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 --0.488169 -0.872651 +7.45059 +-1.6694 +0.828357 +-0.685142 +0.919736 +-0.098447 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 --0.595372 -0.292274 +8.85509 +0.804656 +-0.646786 +-1.27382 +0.539367 +0.693359 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 --0.998214 -0.460923 +7.25952 +0.719406 +-0.505269 +-0.5437 +0.828453 +0.752197 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 --0.784999 --0.242078 +8.70975 +-4.7497 +-0.847882 +-1.49863 +0.655639 +0.0372996 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1589,83 +1589,83 @@

-8.66638 -1.36663 -1.09848 -1.58887 -3.83174 +6.97666 +-5.10709 +-0.647744 +-0.778255 +-10.8619 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m)) -7.70101 -3.37709 -0.257203 -0.353562 -7.01139 +8.4477 +-2.07371 +1.02694 +1.2234 +-3.12049 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m)) -2.55372 -4.78739 --0.165323 -0.554613 -9.40945 +6.34579 +6.50635 +-0.291667 +-0.455588 +12.721 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m)) -9.04899 --3.02341 --1.7829 --1.82417 --7.82973 +6.03157 +-5.25486 +-1.21931 +-1.98666 +-11.729 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m)) -5.35198 --3.48527 -0.129093 --0.730952 --6.84145 +4.81221 +1.04529 +0.7528 +2.09002 +2.84337 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m)) -5.93836 -1.73545 -0.309309 -1.17752 -3.7802 +3.42317 +-0.369744 +-0.409321 +-0.690845 +-1.14881 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m)) -5.7732 --0.510192 -0.597076 --0.188414 --0.423308 +7.45059 +-1.6694 +0.828357 +-0.685142 +-2.51045 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m)) -3.77932 --7.55738 --0.071525 --0.00570931 --15.1863 +8.85509 +0.804656 +-0.646786 +-1.27382 +0.962526 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m)) -4.77217 -1.09176 --0.158519 -0.0301378 -2.025 +7.25952 +0.719406 +-0.505269 +-0.5437 +0.933543 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m)) -5.38049 -4.46788 --0.62577 --0.815971 -8.30999 +8.70975 +-4.7497 +-0.847882 +-1.49863 +-10.3473 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m)) @@ -1699,7 +1699,7 @@ 

gtb |> Filter(row -> row.a < 0 && row.b > 0)
- + @@ -1725,74 +1725,74 @@

-

- - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + @@ -1810,7 +1810,7 @@ 

gtb |> Sort("a", "b")

1278×5 GeoTable over 1278 view(::CartesianGrid, [3398, 4210, 5011, 5014, ..., 9988, 9996, 9997, 9998])1221×5 GeoTable over 1221 view(::CartesianGrid, [464, 789, 1662, 4062, ..., 9990, 9992, 9994, 9995])
a-0.5774910.0720301-0.334594-0.399175Quadrangle((x: 97.0 m, y: 33.0 m), ..., (x: 97.0 m, y: 34.0 m))-0.438512.604690.7620360.555292Quadrangle((x: 63.0 m, y: 4.0 m), ..., (x: 63.0 m, y: 5.0 m))
-0.7819191.807281.473552.32929Quadrangle((x: 9.0 m, y: 42.0 m), ..., (x: 9.0 m, y: 43.0 m))-1.900213.87673-1.60526-1.9003Quadrangle((x: 88.0 m, y: 7.0 m), ..., (x: 88.0 m, y: 8.0 m))
-0.1449941.868640.5537470.911871Quadrangle((x: 10.0 m, y: 50.0 m), ..., (x: 10.0 m, y: 51.0 m))-0.433132.56506-0.572316-1.69405Quadrangle((x: 61.0 m, y: 16.0 m), ..., (x: 61.0 m, y: 17.0 m))
-0.3529172.193090.7623130.988773Quadrangle((x: 13.0 m, y: 50.0 m), ..., (x: 13.0 m, y: 51.0 m))-0.3569551.03008-1.59546-2.34302Quadrangle((x: 61.0 m, y: 40.0 m), ..., (x: 61.0 m, y: 41.0 m))
-0.134830.728194-0.473213-0.100176Quadrangle((x: 18.0 m, y: 50.0 m), ..., (x: 18.0 m, y: 51.0 m))-1.318060.03416031.124092.46925Quadrangle((x: 0.0 m, y: 50.0 m), ..., (x: 0.0 m, y: 51.0 m))
-0.4266450.749988-0.140061-1.12454Quadrangle((x: 29.0 m, y: 50.0 m), ..., (x: 29.0 m, y: 51.0 m))-0.8232962.00016-1.1535-0.272036Quadrangle((x: 7.0 m, y: 50.0 m), ..., (x: 7.0 m, y: 51.0 m))
-2.540780.565004-0.471324-0.140717Quadrangle((x: 31.0 m, y: 50.0 m), ..., (x: 31.0 m, y: 51.0 m))-2.426861.35191-0.532955-1.12277Quadrangle((x: 8.0 m, y: 50.0 m), ..., (x: 8.0 m, y: 51.0 m))
-0.06599041.79317-0.045970.247358Quadrangle((x: 36.0 m, y: 50.0 m), ..., (x: 36.0 m, y: 51.0 m))-1.408770.00141980.3710110.474887Quadrangle((x: 10.0 m, y: 50.0 m), ..., (x: 10.0 m, y: 51.0 m))
-0.2848922.16723-0.328940.188241Quadrangle((x: 38.0 m, y: 50.0 m), ..., (x: 38.0 m, y: 51.0 m))-0.2276111.503560.999945-0.174635Quadrangle((x: 16.0 m, y: 50.0 m), ..., (x: 16.0 m, y: 51.0 m))
-0.815310.3791-0.3531520.418528Quadrangle((x: 41.0 m, y: 50.0 m), ..., (x: 41.0 m, y: 51.0 m))-0.1931151.39062-0.206877-0.285259Quadrangle((x: 17.0 m, y: 50.0 m), ..., (x: 17.0 m, y: 51.0 m))
⋮
- + @@ -1836,74 +1836,74 @@

-

- - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + @@ -1921,7 +1921,7 @@ 

gtb |> Sort("a", "b", rev=true)

10000×5 GeoTable over 10000 view(::CartesianGrid, [7672, 9050, 7313, 7395, ..., 3176, 3253, 525, 4020])10000×5 GeoTable over 10000 view(::CartesianGrid, [9368, 5564, 8448, 5995, ..., 3889, 4185, 4644, 4690])
a-3.46707-0.383294-2.20706-1.9892Quadrangle((x: 71.0 m, y: 76.0 m), ..., (x: 71.0 m, y: 77.0 m))-3.356861.014730.7128341.02863Quadrangle((x: 67.0 m, y: 93.0 m), ..., (x: 67.0 m, y: 94.0 m))
-3.34926-0.06735450.8362590.0953872Quadrangle((x: 49.0 m, y: 90.0 m), ..., (x: 49.0 m, y: 91.0 m))-3.3043-0.3418750.242428-0.670349Quadrangle((x: 63.0 m, y: 55.0 m), ..., (x: 63.0 m, y: 56.0 m))
-3.15256-3.851150.6215720.782224Quadrangle((x: 12.0 m, y: 73.0 m), ..., (x: 12.0 m, y: 74.0 m))-3.20638-0.847291.718112.74252Quadrangle((x: 47.0 m, y: 84.0 m), ..., (x: 47.0 m, y: 85.0 m))
-3.00851-1.73256-0.594382-0.0182037Quadrangle((x: 94.0 m, y: 73.0 m), ..., (x: 94.0 m, y: 74.0 m))-2.900014.01630.2594981.58336Quadrangle((x: 94.0 m, y: 59.0 m), ..., (x: 94.0 m, y: 60.0 m))
-2.98356-1.876390.3298491.07104Quadrangle((x: 66.0 m, y: 63.0 m), ..., (x: 66.0 m, y: 64.0 m))-2.84947-1.23748-1.08622-0.783291Quadrangle((x: 93.0 m, y: 79.0 m), ..., (x: 93.0 m, y: 80.0 m))
-2.982251.22243-1.3959-0.242132Quadrangle((x: 33.0 m, y: 97.0 m), ..., (x: 33.0 m, y: 98.0 m))-2.83591-0.25163-0.159099-0.225349Quadrangle((x: 59.0 m, y: 52.0 m), ..., (x: 59.0 m, y: 53.0 m))
-2.938691.125391.38942.47889Quadrangle((x: 46.0 m, y: 54.0 m), ..., (x: 46.0 m, y: 55.0 m))-2.83259-2.32321-1.1425-1.35694Quadrangle((x: 42.0 m, y: 61.0 m), ..., (x: 42.0 m, y: 62.0 m))
-2.91635-0.001276170.8259510.84968Quadrangle((x: 77.0 m, y: 56.0 m), ..., (x: 77.0 m, y: 57.0 m))-2.802911.71893-1.32188-1.03587Quadrangle((x: 33.0 m, y: 95.0 m), ..., (x: 33.0 m, y: 96.0 m))
-2.91303-1.18189-0.285407-0.190204Quadrangle((x: 75.0 m, y: 99.0 m), ..., (x: 75.0 m, y: 100.0 m))-2.784193.08927-0.725104-1.39562Quadrangle((x: 88.0 m, y: 98.0 m), ..., (x: 88.0 m, y: 99.0 m))
-2.88527-0.2377190.434829-0.35626Quadrangle((x: 48.0 m, y: 57.0 m), ..., (x: 48.0 m, y: 58.0 m))-2.72102-3.58414-1.18211-2.11982Quadrangle((x: 18.0 m, y: 65.0 m), ..., (x: 18.0 m, y: 66.0 m))
⋮
- + @@ -1947,74 +1947,74 @@

-

- - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + @@ -2055,53 +2055,53 @@

-

- + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -2138,53 +2138,53 @@

-

- + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -2221,53 +2221,53 @@

-

- + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -2303,53 +2303,53 @@

-

- + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -2611,7 +2611,7 @@ 

gtb |> Sample(1000, replace=false) |> viewer
- +
@@ -2634,22 +2634,22 @@

Table with 6 columns and 4 rows: variable mean minimum median maximum nmissing ┌─────────────────────────────────────────────────────────────── - 1 │ a 2.99505 -3.46707 1.9536 13.8075 0 - 2 │ b 0.00814059 -10.3081 0.0176344 10.6618 0 - 3 │ c -0.0131867 -3.50267 -0.00274841 3.5377 0 - 4 │ d -0.0142594 -4.56293 -0.0156412 4.63673 0 + 1 │ a 3.01275 -3.35686 2.02636 12.936 0 + 2 │ b -0.0561674 -11.747 -0.0508579 10.4965 0 + 3 │ c -0.0208011 -3.72148 -0.00897106 3.68313 0 + 4 │ d -0.0227169 -4.72084 -0.0257568 4.7739 0 

gtb |> Center("a") |> describe
Table with 6 columns and 4 rows:
-     variable  mean         minimum   median       maximum  nmissing
-   ┌────────────────────────────────────────────────────────────────
- 1 │ a         1.81899e-16  -6.46212  -1.04145     10.8125  0
- 2 │ b         0.00814059   -10.3081  0.0176344    10.6618  0
- 3 │ c         -0.0131867   -3.50267  -0.00274841  3.5377   0
- 4 │ d         -0.0142594   -4.56293  -0.0156412   4.63673  0
+ variable mean minimum median maximum nmissing + ┌───────────────────────────────────────────────────────────────── + 1 │ a -5.45697e-16 -6.36961 -0.986391 9.92326 0 + 2 │ b -0.0561674 -11.747 -0.0508579 10.4965 0 + 3 │ c -0.0208011 -3.72148 -0.00897106 3.68313 0 + 4 │ d -0.0227169 -4.72084 -0.0257568 4.7739 0
@@ -2658,10 +2658,10 @@

Table with 6 columns and 4 rows: variable mean minimum median maximum nmissing ┌───────────────────────────────────────────────────────── - 1 │ a 0.374082 0.0 0.313795 1.0 0 - 2 │ b 0.491955 0.0 0.492408 1.0 0 - 3 │ c 0.495639 0.0 0.497122 1.0 0 - 4 │ d 0.494439 0.0 0.494289 1.0 0 + 1 │ a 0.390945 0.0 0.330403 1.0 0 + 2 │ b 0.525584 0.0 0.525823 1.0 0 + 3 │ c 0.49978 0.0 0.501378 1.0 0 + 4 │ d 0.494813 0.0 0.494493 1.0 0

The ZScore transform is similar to the Scale transform, but it uses the mean and the standard deviation to standardize the range:

@@ -2671,17 +2671,17 @@

Table with 6 columns and 4 rows: variable mean minimum median maximum nmissing ┌───────────────────────────────────────────────────────────────── - 1 │ a 0.0 -1.91599 -0.308784 3.20585 0 - 2 │ b -6.03961e-18 -4.01852 0.00369816 4.14995 0 - 3 │ c -8.88178e-18 -3.48189 0.0104155 3.54315 0 - 4 │ d 2.62901e-17 -3.90469 -0.00118619 3.99252 0 + 1 │ a -9.09495e-17 -1.87385 -0.290183 2.91929 0 + 2 │ b -7.10543e-19 -4.57418 0.00207741 4.12886 0 + 3 │ c -4.9738e-18 -3.71527 0.0118767 3.71855 0 + 4 │ d 2.64677e-17 -3.99774 -0.00258673 4.08155 0

Another important univariate transform is the Quantile transform, which can be used to convert empirical distribution in a column of the geotable to any given distribution from Distributions.jl by Lin et al. (2023). Selected columns are converted to a Normal distribution by default, but more than 60 distributions are available:

gtb |> Quantile() |> values |> pairplot
- +

In data science, scientific traits are used to link data types to adequate statistical algorithms. The most popular scientific traits encountered in geoscientific applications are the Continuous and the Categorical scientific traits. To convert (or coerce) the scientific traits of columns in a geotable, we can use the Coerce transform:

@@ -3080,7 +3080,7 @@ 

gtb |> Rotate(Angle2d(π/4)) |> viewer
- +

In GIS, this new geotable would be called a rotated “raster”. As another example, let’s translate the geotable to the origin of the coordinate system with the Translate transform:

@@ -3089,21 +3089,21 @@

gtb |> Translate(-to(c)...) |> viewer
- +

and scale it with a positive factor for each dimension:

gtb |> Scale(0.1, 0.2) |> viewer
- +

The StdCoords transform combines Translate and Scale to standardize the coordinates of the domain to the interval [-0.5, 0.5]:

gtb |> StdCoords() |> viewer
- +

In GIS, another very important coordinate transform is the Proj transform. We will cover this transform in the next chapter because it depends on the concept of map projection, which deserves more attention.

@@ -3173,14 +3173,14 @@

viewer(t)
- +

We can use the Detrend transform to remove a trend of polynomial degree 2:

t |> Detrend(degree=2) |> viewer
- +

The remaining component can then be modeled with geostatistical models of geospatial correlation, which will be covered in Part IV of the book.

diff --git a/07-pipelines.html b/07-pipelines.html index 0b53f40..cb8ab61 100644 --- a/07-pipelines.html +++ b/07-pipelines.html @@ -401,73 +401,73 @@

-

- - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + @@ -518,63 +518,63 @@

-

 - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + - - - + + + @@ -659,19 +659,19 @@

- + - + - + @@ -683,19 +683,19 @@

- + - + - + @@ -740,53 +740,53 @@

-

 - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -848,80 +848,80 @@

false

- - + + - - + + + - - - + + - - - + + + - - - + + + - - + + - - - + + + - - - + + + + - - - + + - - + + diff --git a/10-correlation.html b/10-correlation.html index 6d664f6..aced92a 100644 --- a/10-correlation.html +++ b/10-correlation.html @@ -442,7 +442,7 @@

10  
sample = img |> Sample(1000, replace=false)

10000×5 GeoTable over 10000 view(::CartesianGrid, [4020, 525, 3253, 3176, ..., 7395, 7313, 9050, 7672])10000×5 GeoTable over 10000 view(::CartesianGrid, [4690, 4644, 4185, 3889, ..., 5995, 8448, 5564, 9368])
a13.8075-0.0944749-0.485679-0.336729Quadrangle((x: 19.0 m, y: 40.0 m), ..., (x: 19.0 m, y: 41.0 m))12.936-2.72305-1.65361-1.57103Quadrangle((x: 89.0 m, y: 46.0 m), ..., (x: 89.0 m, y: 47.0 m))
13.1831-0.470644-1.29313-0.687942Quadrangle((x: 24.0 m, y: 5.0 m), ..., (x: 24.0 m, y: 6.0 m))12.5552.041171.024610.558263Quadrangle((x: 43.0 m, y: 46.0 m), ..., (x: 43.0 m, y: 47.0 m))
13.0283.472322.718242.96301Quadrangle((x: 52.0 m, y: 32.0 m), ..., (x: 52.0 m, y: 33.0 m))12.53524.788660.9352430.695069Quadrangle((x: 84.0 m, y: 41.0 m), ..., (x: 84.0 m, y: 42.0 m))
12.5963-4.3-0.4188280.866281Quadrangle((x: 75.0 m, y: 31.0 m), ..., (x: 75.0 m, y: 32.0 m))12.3475-0.954583-1.07227-1.00695Quadrangle((x: 88.0 m, y: 38.0 m), ..., (x: 88.0 m, y: 39.0 m))
12.41771.01099-0.163681-0.037842Quadrangle((x: 65.0 m, y: 5.0 m), ..., (x: 65.0 m, y: 6.0 m))12.29112.15845-0.9537170.472874Quadrangle((x: 74.0 m, y: 29.0 m), ..., (x: 74.0 m, y: 30.0 m))
12.1269-0.413411-0.335690.580181Quadrangle((x: 32.0 m, y: 48.0 m), ..., (x: 32.0 m, y: 49.0 m))12.0435-0.488439-0.602107-0.210599Quadrangle((x: 55.0 m, y: 39.0 m), ..., (x: 55.0 m, y: 40.0 m))
11.9822-1.62725-1.14082-1.42278Quadrangle((x: 89.0 m, y: 28.0 m), ..., (x: 89.0 m, y: 29.0 m))12.00150.08742110.157607-0.00494087Quadrangle((x: 29.0 m, y: 23.0 m), ..., (x: 29.0 m, y: 24.0 m))
11.9653-1.86105-0.74178-1.38585Quadrangle((x: 60.0 m, y: 3.0 m), ..., (x: 60.0 m, y: 4.0 m))11.8507-0.1369020.6578981.2378Quadrangle((x: 13.0 m, y: 20.0 m), ..., (x: 13.0 m, y: 21.0 m))
11.95773.92175-0.0290192-1.32897Quadrangle((x: 20.0 m, y: 9.0 m), ..., (x: 20.0 m, y: 10.0 m))11.8279-0.172495-0.818378-1.77468Quadrangle((x: 97.0 m, y: 13.0 m), ..., (x: 97.0 m, y: 14.0 m))
11.7395-3.37414-1.37179-0.271044Quadrangle((x: 75.0 m, y: 14.0 m), ..., (x: 75.0 m, y: 15.0 m))11.8136-3.722870.3661660.302904Quadrangle((x: 45.0 m, y: 28.0 m), ..., (x: 45.0 m, y: 29.0 m))
⋮8.666381.366636.97666-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
7.701013.377098.4477-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
2.553724.787396.345796.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
9.04899-3.023416.03157-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
5.35198-3.485274.812211.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
5.938361.735453.42317-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
5.7732-0.5101927.45059-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
3.77932-7.557388.855090.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
4.772171.091767.259520.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
5.380494.467888.70975-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
8.666381.366636.97666-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
7.701013.377098.4477-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
2.553724.787396.345796.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
9.04899-3.023416.03157-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
5.35198-3.485274.812211.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
5.938361.735453.42317-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
5.7732-0.5101927.45059-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
3.77932-7.557388.855090.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
4.772171.091767.259520.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
5.380494.467888.70975-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
8.666381.366636.97666-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
7.701013.377098.4477-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
2.553724.787396.345796.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
9.04899-3.023416.03157-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
5.35198-3.485274.812211.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
5.938361.735453.42317-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
5.7732-0.5101927.45059-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
3.77932-7.557388.855090.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
4.772171.091767.259520.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
5.380494.467888.70975-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
8.666381.366636.97666-5.10709 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
7.701013.377098.4477-2.07371 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
2.553724.787396.345796.50635 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
9.04899-3.023416.03157-5.25486 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
5.35198-3.485274.812211.04529 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
5.938361.735453.42317-0.369744 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
5.7732-0.5101927.45059-1.6694 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
3.77932-7.557388.855090.804656 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
4.772171.091767.259520.719406 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
5.380494.467888.70975-4.7497 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
5.402713.120320.9974841.113267.08711-0.4021990.9375981.12817 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
4.62316-6.20823-0.834121-0.007472872.50585-1.97214-2.17861-2.66912 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
3.97525-0.519399-1.82212-1.892345.88418-5.197150.480288-0.542933 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
6.101130.238265-0.03949090.3811714.608391.317520.2081381.04836 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
6.158182.74625-0.724002-0.5919083.41838-1.36776-1.01328-0.91081 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
4.529112.74331-1.00981-2.084183.18762-1.80649-0.559284-1.44291 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
5.459382.66034-0.5515750.4479254.51168-2.60257-1.84294-1.3519 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
8.9218-0.8966210.24948-0.01380044.77423-1.54421-0.43362-0.517143 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
4.46891-6.56056-1.33678-1.720758.343175.49575-0.440129-0.690431 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
4.052832.947610.252404-0.965873.96962-3.82759-0.1548930.0566007 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
0.5115011.267990.9842351.05287-0.1695040.925704 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
0.309686-2.33775-0.8516650.0361033-0.817575-2.16224 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
0.202405-0.203429-1.791580.636578-1.995390.473299 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
0.7195530.113039-0.06295790.311790.5452230.197802 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
0.7339161.11933-0.7418140.121367-0.582544-1.01522 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
0.2918521.11792-1.018160.0969073-0.757083-0.569525 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
0.5267031.0889-0.5745440.289237-1.05331-1.83302 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
1.79034-0.3581920.2280590.34832-0.65697-0.45404 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
0.281405-2.46243-1.32071.549762.10836-0.460998 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
0.2144451.195730.2311480.192692-1.51174-0.172302 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
falsetruefalse true Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
falsefalsetrue true Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
falsefalsetrue true Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
falsefalsetrue true Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
falsefalsetrue true Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
falsetruefalse true Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
-1.382541.22097-1.33256-0.153478 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
0.642242-2.430793.16026-0.808094 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
2.42858-0.150932-0.0152279-2.04903 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
-0.1941890.102205-0.7629310.518553 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
0.8536791.111521.27572-0.553782 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
1.949831.132031.27496-0.728995 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
0.1092771.062752.12468-1.04387 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
-0.152716-0.3437350.627421-0.618303 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
2.02588-2.543350.760312.14175 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
0.3898261.186320.0759567-1.50862 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
false true-1.382541.22097-1.33256-0.153478 Quadrangle((x: 0.0 m, y: 0.0 m), ..., (x: 0.0 m, y: 1.0 m))
false true true0.642242-2.430793.16026-0.808094 Quadrangle((x: 1.0 m, y: 0.0 m), ..., (x: 1.0 m, y: 1.0 m))
falsefalse truetrue2.42858-0.150932-0.0152279-2.04903 Quadrangle((x: 2.0 m, y: 0.0 m), ..., (x: 2.0 m, y: 1.0 m))
falsefalse true-0.1941890.102205true-0.7629310.518553 Quadrangle((x: 3.0 m, y: 0.0 m), ..., (x: 3.0 m, y: 1.0 m))
falsefalse true0.8536791.11152true1.27572-0.553782 Quadrangle((x: 4.0 m, y: 0.0 m), ..., (x: 4.0 m, y: 1.0 m))
false true true1.949831.132031.27496-0.728995 Quadrangle((x: 5.0 m, y: 0.0 m), ..., (x: 5.0 m, y: 1.0 m))
falsefalse true0.1092771.06275true2.12468-1.04387 Quadrangle((x: 6.0 m, y: 0.0 m), ..., (x: 6.0 m, y: 1.0 m))
falsefalse true-0.152716-0.343735true0.627421-0.618303 Quadrangle((x: 7.0 m, y: 0.0 m), ..., (x: 7.0 m, y: 1.0 m))
falsefalse truetrue2.02588-2.543350.760312.14175 Quadrangle((x: 8.0 m, y: 0.0 m), ..., (x: 8.0 m, y: 1.0 m))
false true true0.3898261.186320.0759567-1.50862 Quadrangle((x: 9.0 m, y: 0.0 m), ..., (x: 9.0 m, y: 1.0 m))
- + @@ -459,44 +459,44 @@

10  -

 - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + - - + + @@ -510,32 +510,32 @@

10  
hscatter(sample, :Z, :Z, lag=0.0)
- +

By increasing the value of the lag, we observe that the correlation is no longer equal to one, and that the linear fit through the points approaches the horizontal axis (i.e., zero correlation):

hscatter(sample, :Z, :Z, lag=3.0)
- + 
hscatter(sample, :Z, :Z, lag=5.0)
- + 
hscatter(sample, :Z, :Z, lag=10.0)
- + 
hscatter(sample, :Z, :Z, lag=50.0)
- +

The Pearson correlation coefficient studied as a function of the lag \(h\) is known as the correlogram function. For example, consider the exponential correlogram function given by \(cor(h) = \exp(-h)\):

diff --git a/12-mining.html b/12-mining.html index ec9e568..ec1bbb1 100644 --- a/12-mining.html +++ b/12-mining.html @@ -966,7 +966,7 @@ 

grades |> CLR() |> ProjectionPursuit() |> values |> pairplot
- +

The ProjectionPursuit is an advanced statistical transform that removes non-linear associations between variables using an iterative procedure (Friedman 1987). The result is a set of independent variables that can be interpolated separately.

@@ -1005,73 +1005,73 @@

- - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + @@ -1099,10 +1099,10 @@

γs = [GeoStatsFunctions.fit(Variogram, g, h -> 1 / h^2) for g in gs] 
4-element Vector{Variogram}:
- MaternVariogram(range: 165.448 m, sill: 0.873103, nugget: 0.431565, order: 1.0)
- CubicVariogram(range: 131.38 m, sill: 0.930791, nugget: 0.635372)
- MaternVariogram(range: 156.205 m, sill: 0.894003, nugget: 0.528769, order: 1.0)
- ExponentialVariogram(range: 284.798 m, sill: 0.904187, nugget: 0.440352)
+ MaternVariogram(range: 165.165 m, sill: 0.873085, nugget: 0.431516, order: 1.0) + CubicVariogram(range: 131.49 m, sill: 0.931023, nugget: 0.635639) + MaternVariogram(range: 156.653 m, sill: 0.894464, nugget: 0.529554, order: 1.0) + ExponentialVariogram(range: 284.741 m, sill: 0.904222, nugget: 0.440377)
@@ -1127,25 +1127,25 @@

gammaplot(vs[1], gs[1], γs[1])

- + 
gammaplot(vs[2], gs[2], γs[2])
- + 
gammaplot(vs[3], gs[3], γs[3])
- + 
gammaplot(vs[4], gs[4], γs[4])
- +

Assuming that the variogram models are adequate, we can proceed to interpolation.

@@ -1187,73 +1187,73 @@

-

 - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + @@ -1271,7 +1271,7 @@ 

interp |> Sample(10000) |> values |> pairplot
- +
@@ -1308,73 +1308,73 @@

- - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + - - - - + + + + @@ -1391,7 +1391,7 @@ 

estim |> Select("Cu") |> viewer
- +
@@ -1441,83 +1441,83 @@

- - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + - - - - - + + + + + @@ -1582,93 +1582,93 @@

<

- - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + - - - - - - + + + + + + @@ -1688,14 +1688,14 @@

< 
estim |> Filter(x -> x.value > 0) |> Select("value") |> viewer
- +

Or any criterion of interest such as positive economic value and small fraction of contaminants:

estim |> Filter(x -> x.value > 0 && x.S < 0.25) |> Select("value") |> viewer
- +
diff --git a/14-energy.html b/14-energy.html index 5906664..7088b36 100644 --- a/14-energy.html +++ b/14-energy.html @@ -432,7 +432,7 @@

14 

14.1 Data

-

We will use reservoir simulation results of the Norne benchmark case, a real oil field from the Norwegian Sea. For more information, please check the OPM project. These results were simulated with the JutulDarcy.jl reservoir simulator by Møyner et al. (2025).

+

We will use reservoir simulation results of the Norne benchmark case, a real oil field from the Norwegian Sea. For more information, please check the OPM project. These results were simulated with the JutulDarcy.jl reservoir simulator by Møyner (2024).

In particular, we will consider only two time steps of the simulation, named norne1.vtu and norne2.vtu. The data are stored in the open VTK format with .vtu extension, indicating that it is georeferenced over an unstructured mesh:

norne₁ = GeoIO.load("data/norne1.vtu")
@@ -684,7 +684,7 @@ 

 zones |> viewer
- +
@@ -731,7 +731,7 @@

0.0 kg

- + @@ -739,7 +739,7 @@

2.43976e6 kg

- + @@ -747,7 +747,7 @@

2.24126e6 kg

- + @@ -755,7 +755,7 @@

2.16366e6 kg

- + @@ -763,7 +763,7 @@

2.21737e6 kg

- + @@ -771,7 +771,7 @@

2.25638e6 kg

- + @@ -779,7 +779,7 @@

1.91029e6 kg

- + @@ -787,7 +787,7 @@

1.78157e6 kg

- + @@ -795,7 +795,7 @@

1.70066e6 kg

- + @@ -803,7 +803,7 @@

1.57036e6 kg

- + @@ -847,19 +847,19 @@

- - - + + + - - - + + + - - - + + +
1000×2 GeoTable over 1000 view(::CartesianGrid, [2723, 7725, 4929, 1842, ..., 7186, 570, 5535, 2838])1000×2 GeoTable over 1000 view(::CartesianGrid, [6461, 9066, 4198, 7343, ..., 8560, 8275, 3102, 8749])
Z0.409393Quadrangle((x: 22.0 m, y: 27.0 m), ..., (x: 22.0 m, y: 28.0 m))1.24246Quadrangle((x: 60.0 m, y: 64.0 m), ..., (x: 60.0 m, y: 65.0 m))
-0.734586Quadrangle((x: 24.0 m, y: 77.0 m), ..., (x: 24.0 m, y: 78.0 m))-0.146836Quadrangle((x: 65.0 m, y: 90.0 m), ..., (x: 65.0 m, y: 91.0 m))
-1.56952Quadrangle((x: 28.0 m, y: 49.0 m), ..., (x: 28.0 m, y: 50.0 m))0.449928Quadrangle((x: 97.0 m, y: 41.0 m), ..., (x: 97.0 m, y: 42.0 m))
0.191135Quadrangle((x: 41.0 m, y: 18.0 m), ..., (x: 41.0 m, y: 19.0 m))-0.444043Quadrangle((x: 42.0 m, y: 73.0 m), ..., (x: 42.0 m, y: 74.0 m))
-0.457468Quadrangle((x: 49.0 m, y: 6.0 m), ..., (x: 49.0 m, y: 7.0 m))0.0575923Quadrangle((x: 44.0 m, y: 64.0 m), ..., (x: 44.0 m, y: 65.0 m))
0.146243Quadrangle((x: 81.0 m, y: 78.0 m), ..., (x: 81.0 m, y: 79.0 m))0.202576Quadrangle((x: 48.0 m, y: 51.0 m), ..., (x: 48.0 m, y: 52.0 m))
-0.0847255Quadrangle((x: 80.0 m, y: 60.0 m), ..., (x: 80.0 m, y: 61.0 m))-0.272168Quadrangle((x: 76.0 m, y: 7.0 m), ..., (x: 76.0 m, y: 8.0 m))
0.143076Quadrangle((x: 95.0 m, y: 5.0 m), ..., (x: 95.0 m, y: 6.0 m))0.000884331Quadrangle((x: 92.0 m, y: 6.0 m), ..., (x: 92.0 m, y: 7.0 m))
-0.154965Quadrangle((x: 29.0 m, y: 42.0 m), ..., (x: 29.0 m, y: 43.0 m))0.429615Quadrangle((x: 73.0 m, y: 60.0 m), ..., (x: 73.0 m, y: 61.0 m))
-0.139083Quadrangle((x: 75.0 m, y: 90.0 m), ..., (x: 75.0 m, y: 91.0 m))0.0292726Quadrangle((x: 29.0 m, y: 88.0 m), ..., (x: 29.0 m, y: 89.0 m))
⋮
-0.7297331.53324-1.239641.23364-0.7239231.52865-1.234191.23527 (x: 559.725 m, y: -513.31 m, z: 252.82 m)
-1.47832-0.312423-1.876111.08543-1.45634-0.308887-1.841161.0797 (x: 558.955 m, y: -515.23 m, z: 246.87 m)
-0.2834420.9332930.2955981.24323-0.2900310.9359330.2925611.24015 (x: 557.225 m, y: -519.61 m, z: 233.37 m)
-1.340351.538320.108343-0.378056-1.321541.534360.115323-0.375419 (x: 555.375 m, y: -524.3 m, z: 218.92 m)
-0.5766831.20655-0.462563-0.412195-0.5747881.2108-0.465672-0.409665 (x: 553.825 m, y: -528.21 m, z: 206.94 m)
-1.037591.85798-0.118710.0580981-1.027511.84606-0.1125960.0600358 (x: 552.075 m, y: -532.56 m, z: 193.91 m)
-0.3548682.01103-1.776741.46575-0.342722.00237-1.768351.46248 (x: 550.305 m, y: -537.08 m, z: 180.8 m)
-1.951080.365021-0.3087270.763003-1.948220.360919-0.3034170.761106 (x: 547.495 m, y: -544.31 m, z: 160.08 m)
-0.135559-1.1933-1.026961.3586-0.132293-1.19916-1.025991.3655 (x: 544.335 m, y: -552.43 m, z: 137.08 m)
-1.34592-0.568661-0.181257-0.194656-1.36001-0.564604-0.182924-0.193876 (x: 541.695 m, y: -559.19 m, z: 117.96 m)
-0.306621-0.2704690.2487360.430026-0.306652-0.2704150.2504890.431586 Hexahedron((x: 724.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -857.9 m, z: -387.42 m))
0.103434-0.1920520.6466180.4747130.104378-0.1932450.6501410.476164 Hexahedron((x: 749.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -857.9 m, z: -387.42 m))
0.103423-0.1920520.6466240.4743670.104367-0.1932450.6501470.475817 Hexahedron((x: 774.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -857.9 m, z: -387.42 m))
0.103416-0.1920520.6466270.4740690.10436-0.1932450.650150.475518 Hexahedron((x: 799.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 799.995 m, y: -857.9 m, z: -387.42 m))
0.0109184-0.1417440.110584-0.1402170.0102892-0.142480.109897-0.139926 Hexahedron((x: 824.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 824.995 m, y: -857.9 m, z: -387.42 m))
-0.385594-0.04425560.3871160.447342-0.385403-0.04540250.3891340.448921 Hexahedron((x: 674.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 674.995 m, y: -832.9 m, z: -387.42 m))
-0.306623-0.2704690.248740.430606-0.306654-0.2704150.2504940.432168 Hexahedron((x: 699.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 699.995 m, y: -832.9 m, z: -387.42 m))
-0.306618-0.2704690.2487620.430236-0.306649-0.2704150.2505170.431797 Hexahedron((x: 724.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -832.9 m, z: -387.42 m))
0.103412-0.1920520.6466310.4748010.104356-0.1932450.6501540.476253 Hexahedron((x: 749.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -832.9 m, z: -387.42 m))
0.103398-0.1920520.6466380.4744050.104342-0.1932450.6501610.475855 Hexahedron((x: 774.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -832.9 m, z: -387.42 m))
0.7416123.10212e-50.0001995190.2581580.7413493.12956e-50.0002001850.258419 Hexahedron((x: 724.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -857.9 m, z: -387.42 m))
0.7490592.54459e-50.0002071720.2507090.7490392.5455e-50.0002076490.250728 Hexahedron((x: 749.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -857.9 m, z: -387.42 m))
0.748882.54014e-50.0002070010.2508880.7490392.54539e-50.0002076510.250728 Hexahedron((x: 774.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -857.9 m, z: -387.42 m))
0.7488932.54247e-50.0002068650.2508740.749042.54535e-50.000207650.250727 Hexahedron((x: 799.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 799.995 m, y: -857.9 m, z: -387.42 m))
0.7521953.36429e-50.0001988630.2475720.7519953.3731e-50.0001988450.247773 Hexahedron((x: 824.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 824.995 m, y: -857.9 m, z: -387.42 m))
0.7405813.16091e-50.0002049530.2591820.7406933.1575e-50.0002049620.259071 Hexahedron((x: 674.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 674.995 m, y: -832.9 m, z: -387.42 m))
0.7416133.10176e-50.0001994780.2581570.7413473.13e-50.000200190.258422 Hexahedron((x: 699.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 699.995 m, y: -832.9 m, z: -387.42 m))
0.7416153.10201e-50.0001994960.2581540.7413483.12982e-50.0002001870.25842 Hexahedron((x: 724.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -832.9 m, z: -387.42 m))
0.7490552.5446e-50.0002071740.2507120.7490382.54572e-50.000207640.250729 Hexahedron((x: 749.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -832.9 m, z: -387.42 m))
0.7489312.54284e-50.0002070480.2508370.7490392.54539e-50.0002076510.250728 Hexahedron((x: 774.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -832.9 m, z: -387.42 m))
0.7416123.10212e-50.0001995190.2581580.9880150.7413493.12956e-50.0002001850.2584190.987914 Hexahedron((x: 724.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -857.9 m, z: -387.42 m))
0.7490592.54459e-50.0002071720.2507090.9913990.7490392.5455e-50.0002076490.2507280.991422 Hexahedron((x: 749.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -857.9 m, z: -387.42 m))
0.748882.54014e-50.0002070010.2508880.991330.7490392.54539e-50.0002076510.2507280.991422 Hexahedron((x: 774.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -857.9 m, z: -387.42 m))
0.7488932.54247e-50.0002068650.2508740.9913350.749042.54535e-50.000207650.2507270.991423 Hexahedron((x: 799.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 799.995 m, y: -857.9 m, z: -387.42 m))
0.7521953.36429e-50.0001988630.2475720.9925230.7519953.3731e-50.0001988450.2477730.992484 Hexahedron((x: 824.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 824.995 m, y: -857.9 m, z: -387.42 m))
0.7405813.16091e-50.0002049530.2591820.9874530.7406933.1575e-50.0002049620.2590710.987555 Hexahedron((x: 674.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 674.995 m, y: -832.9 m, z: -387.42 m))
0.7416133.10176e-50.0001994780.2581570.9880150.7413473.13e-50.000200190.2584220.987913 Hexahedron((x: 699.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 699.995 m, y: -832.9 m, z: -387.42 m))
0.7416153.10201e-50.0001994960.2581540.9880160.7413483.12982e-50.0002001870.258420.987913 Hexahedron((x: 724.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -832.9 m, z: -387.42 m))
0.7490552.5446e-50.0002071740.2507120.9913980.7490382.54572e-50.000207640.2507290.991422 Hexahedron((x: 749.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -832.9 m, z: -387.42 m))
0.7489312.54284e-50.0002070480.2508370.991350.7490392.54539e-50.0002076510.2507280.991422 Hexahedron((x: 774.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -832.9 m, z: -387.42 m))
0.7416123.10212e-50.0001995190.2581580.9880153.28903e50.7413493.12956e-50.0002001850.2584190.9879143.28616e5 Hexahedron((x: 724.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -857.9 m, z: -387.42 m))
0.7490592.54459e-50.0002071720.2507090.9913993.37404e50.7490392.5455e-50.0002076490.2507280.9914223.37403e5 Hexahedron((x: 749.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -857.9 m, z: -387.42 m))
0.748882.54014e-50.0002070010.2508880.991333.37207e50.7490392.54539e-50.0002076510.2507280.9914223.37403e5 Hexahedron((x: 774.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -857.9 m, z: -387.42 m))
0.7488932.54247e-50.0002068650.2508740.9913353.37223e50.749042.54535e-50.000207650.2507270.9914233.37403e5 Hexahedron((x: 799.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 799.995 m, y: -857.9 m, z: -387.42 m))
0.7521953.36429e-50.0001988630.2475720.992523340803.00.7519953.3731e-50.0001988450.2477730.9924843.40607e5 Hexahedron((x: 824.995 m, y: -882.9 m, z: -399.92 m), ..., (x: 824.995 m, y: -857.9 m, z: -387.42 m))
0.7405813.16091e-50.0002049530.2591820.9874533.2767e50.7406933.1575e-50.0002049620.2590710.987555327830.0 Hexahedron((x: 674.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 674.995 m, y: -832.9 m, z: -387.42 m))
0.7416133.10176e-50.0001994780.2581570.9880153.28904e50.7413473.13e-50.000200190.2584220.9879133.28613e5 Hexahedron((x: 699.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 699.995 m, y: -832.9 m, z: -387.42 m))
0.7416153.10201e-50.0001994960.2581540.9880163.28907e50.7413483.12982e-50.0002001870.258420.9879133.28614e5 Hexahedron((x: 724.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 724.995 m, y: -832.9 m, z: -387.42 m))
0.7490552.5446e-50.0002071740.2507120.9913983.37401e50.7490382.54572e-50.000207640.2507290.9914223.37402e5 Hexahedron((x: 749.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 749.995 m, y: -832.9 m, z: -387.42 m))
0.7489312.54284e-50.0002070480.2508370.991353.37264e50.7490392.54539e-50.0002076510.2507280.9914223.37403e5 Hexahedron((x: 774.995 m, y: -857.9 m, z: -399.92 m), ..., (x: 774.995 m, y: -832.9 m, z: -387.42 m))
7.56895e6 kg 0.0 kg21 Hexahedron((x: 455586.0 m, y: 7.32115e6 m, z: 2586.74 m), ..., (x: 455636.0 m, y: 7.32122e6 m, z: 2605.72 m))
7141.61 kg 2.41219e6 kg21 Hexahedron((x: 455608.0 m, y: 7.32107e6 m, z: 2568.85 m), ..., (x: 455659.0 m, y: 7.32114e6 m, z: 2585.64 m))
3835.96 kg 2.21634e6 kg21 Hexahedron((x: 455625.0 m, y: 7.32099e6 m, z: 2563.06 m), ..., (x: 455678.0 m, y: 7.32106e6 m, z: 2573.92 m))
3652.63 kg 2.13943e6 kg21 Hexahedron((x: 455650.0 m, y: 7.32092e6 m, z: 2559.15 m), ..., (x: 455706.0 m, y: 7.32099e6 m, z: 2569.36 m))
3746.6 kg 2.19304e6 kg21 Hexahedron((x: 455677.0 m, y: 7.32085e6 m, z: 2555.98 m), ..., (x: 455731.0 m, y: 7.32091e6 m, z: 2567.17 m))
3821.64 kg 2.23166e6 kg21 Hexahedron((x: 455700.0 m, y: 7.32078e6 m, z: 2553.79 m), ..., (x: 455756.0 m, y: 7.32084e6 m, z: 2562.11 m))
3226.36 kg 1.88936e6 kg21 Hexahedron((x: 455724.0 m, y: 7.3207e6 m, z: 2552.27 m), ..., (x: 455781.0 m, y: 7.32076e6 m, z: 2560.06 m))
2993.03 kg 1.76212e6 kg21 Hexahedron((x: 455751.0 m, y: 7.32063e6 m, z: 2551.07 m), ..., (x: 455804.0 m, y: 7.32068e6 m, z: 2561.31 m))
2868.5 kg 1.68211e6 kg21 Hexahedron((x: 455785.0 m, y: 7.32057e6 m, z: 2549.52 m), ..., (x: 455831.0 m, y: 7.32061e6 m, z: 2559.75 m))
2642.13 kg 1.55327e6 kg21 Hexahedron((x: 455823.0 m, y: 7.3205e6 m, z: 2547.65 m), ..., (x: 455866.0 m, y: 7.32054e6 m, z: 2556.82 m))
22.97817e10 kgMulti(15752×Hexahedron)14.31915e10 kgMulti(19194×Hexahedron)
33.07385e10 kgMulti(11618×Hexahedron)29.35253e9 kgMulti(22111×Hexahedron)
12.49234e9 kgMulti(17061×Hexahedron)31.04685e10 kgMulti(3126×Hexahedron)
@@ -890,19 +890,19 @@

-2 -2.97817e7 Mg -Multi(15752×Hexahedron) +1 +4.31915e7 Mg +Multi(19194×Hexahedron) -3 -3.07385e7 Mg -Multi(11618×Hexahedron) +2 +9.35253e6 Mg +Multi(22111×Hexahedron) -1 -2.49234e6 Mg -Multi(17061×Hexahedron) +3 +1.04685e7 Mg +Multi(3126×Hexahedron) @@ -924,8 +924,8 @@

Fouedjio, Francky. 2016. “A Hierarchical Clustering Method for Multivariate Geostatistical Data.” Spatial Statistics 18: 333–51. https://doi.org/https://doi.org/10.1016/j.spasta.2016.07.003.

-
-Møyner, Olav, Jakob Torben, Øystein Klemetsdal, Bruno M. Pacheco, Andrés Riedemann, Grant Bruer, andreas-brostrom, et al. 2025. “Sintefmath/JutulDarcy.jl: V0.2.40.” Zenodo. https://doi.org/10.5281/zenodo.14671781. +
+Møyner, O. 2024. “JutulDarcy.jl - a Fully Differentiable High-Performance Reservoir Simulator Based on Automatic Differentiation” 2024 (1): 1–9. https://doi.org/https://doi.org/10.3997/2214-4609.202437111.
diff --git a/index.html b/index.html index fca7b7b..689b9c1 100644 --- a/index.html +++ b/index.html @@ -7,7 +7,7 @@ - + Geospatial Data Science with Julia