# Create an initial litholog ----
l <- c(-2,-1,0,1,2)
r <- c(-1,0,1,2,3)
h <- c(4,3,4,3,4)
i <- c("B1","B2","B3","B4","B5")
log <- litholog(l, r, h, i)
# Get a custom bedding joint to specific dimensions using changejoint() ----
liq <- changejoint(oufti99$liquefaction,
yleft = 0, ymax = 0.3,
xmin = 1, xmax = 2)
nlog <- weldlog(log, dt = 0, seg = list(liq = liq), j = c("liq"))
# Plots for visualisation ----
plot.new()
plot.window(xlim = c(0,5), ylim = c(-2,3))
axis(1)
axis(2)
multigons(nlog$i, nlog$xy, nlog$dt)
Run the code above in your browser using DataLab