weldjoint: Changes boundaries segments in basic lithologs

In StratigrapheR: Integrated Stratigraphy

Changes boundaries segments in basic lithologs

Description

Adds personalised segments to bed boundaries of lithologs from "litholog()"-like data frames

Usage

weldjoint(
  log,
  dt,
  collection,
  sym,
  yinv = F,
  xinv = F,
  yleft = 0,
  yright = NA,
  ymin = NA,
  ymax = NA,
  xmin = 0,
  xmax = max(log$xy),
  add.dt = 0,
  tolerance = 8
)

Arguments

log	a "litholog()"-like data frame on which the new segments need to be welded.
dt	the position of the n boundaries to change.
collection	a collection object (e.g. oufti)
sym	the name of the symbols in the collection. This should be a symbol that can be considered as a bedding joint; see is.joint.
yinv, xinv	whether to inverse the plotting for x and y values (T or F)
yleft, yright	the depth/height/time value for the extreme point at the right or left of the joint (yleft overruns yright, which overruns ymin and ymax)
ymin, ymax	the extreme values for the y axis (in case of conflict with yleft and/or yright, defaults to the smallest exaggeration)
xmin, xmax	the extreme values for the x axis
add.dt	whether to automatically add the dt value to the dt of the segments (with the add.dt value when it is not zero)
tolerance	the order of tolerance for errors, i.e. the number of decimals considered as being meaningful for matching dt to log

Value

a "litholog()"-like data frame, with new bed boundaries

Examples

# Generate litholog ----

l <- c(0,1,2,3,4)
r <- c(1,2,3,4,5)
h   <- c(4,3,4,3,4)
i <- c("B1","B2","B3","B4","B5")
log  <- litholog(l, r, h, i)

# Modify the boundaries of the litholog ----

nlog <- weldjoint(log, c(1,2,3,4), oufti99,
                  sym = c("1sin", "stylolith", "3sin", "liquefaction"),
                  ymax  = c(NA, NA, NA, 0.2),
                  xmin  = c(0,0,0,1),
                  xmax  = c(4,4,4,1.5))

# Visualise

par(mfrow = c(1,2))

plot.new()
plot.window(xlim = c(0,5), ylim = c(0,5))

axis(1)
axis(2)

multigons(log$i, log$xy, log$dt)

plot.new()
plot.window(xlim = c(0,5), ylim = c(0,5))

axis(1)
axis(2)

multigons(nlog$i, nlog$xy, nlog$dt)

StratigrapheR documentation built on July 9, 2023, 6:02 p.m.