Most objects defined in GeoStats.jl can be plotted directly using the plot command from Plots.jl. For visualization of 3D objects, however, we recommend the experimental MeshViz.jl package. Additional plots are listed below that can be useful for geostatistical analysis.


A hscatter plot between two variables var1 and var2 (possibly with var2 = var1) is a simple scatter plot in which the dots represent all ordered pairs of values of var1 and var2 at a given lag h.

using GeoStats
using Plots

๐’Ÿ = georef((Z=[10sin(i/10) + j for i in 1:100, j in 1:200],))

๐’ฎ = sample(๐’Ÿ, 500)

p1 = hscatter(๐’ฎ, :Z, lag=0)
p2 = hscatter(๐’ฎ, :Z, lag=20)
p3 = hscatter(๐’ฎ, :Z, lag=40)
p4 = hscatter(๐’ฎ, :Z, lag=60)

plot(p1, p2, p3, p4)


The PairPlots.jl package provides the corner plot that can be used with any table, including tables of attributes obtained with the values function.


The StatsPlots.jl package provides various statistical plots such as boxplot, dotplot, violin and other plots commonly used in statistical workflows.