dylanbeaudette
commented
6 years ago A very simple, perhaps comical example.
library(raster)
library(wesanderson)
# empty template for horizon classes or membership
r <- raster(nrows=10, ncols=10)
# make up some data representing
# mixed-up 1D but otherwise common 2D description
r[1:28] <- 1
r[29:44] <- 2
r[45:71] <- 3
r[71:100] <- 4
r[77:83] <- 5
r[4:8] <- 3
r[22:28] <- 5
r[80:88] <- 2
# check
as(r, 'matrix')
# simple plot of horizon codes
plot(r, col=wes_palette('Zissou', n=5), asp=4, axes=FALSE)
2D vs. 1D horizonation, useful in situation where cryoturbation is a dominant pedogenic driver.
See Sampling Protocols for Permafrost-Affected Soils for some examples.
Recent papers:
Borrowed from Siewert et al., 2016: