StratigrapheR (version 1.3.1)

weldjoint: 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
)

Value

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

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

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)