zijderveld {StratigrapheR} | R Documentation |
Draws a Zijderveld plot
Description
Draws a Zijderveld plot: it projects 3D points (having declination, inclination and intensity) in 2D, horizontally and vertically.
Usage
zijderveld(
dec,
inc,
int,
xh = "WE",
xv = xh,
centre = F,
xlim = NA,
ylim = NA,
unit = NA,
xlab = "",
ylab = "",
labels = NA,
nlabels = 1,
h = list(pch = 19),
v = list(pch = 21, bg = "white"),
f = list(pch = 21, bg = "white", cex = 1.5),
t = list(pos = 3, offset = 0.5),
l = list(),
anchored = T,
style = "branches",
tcl = 0.2,
orientation = TRUE,
scientific = NA,
decimals = 10,
add = FALSE
)
Arguments
dec |
declination of the data; it is the angle from the north taken on an horizontal plane. It is measured clockwise from North and ranges from 0 to 360° (Tauxe 2010). Values outside this range are corrected by incfix(). |
inc |
inclination of the data; it is the angle from the horizontal, is positive downward, and ranges from +90° for straight down to -90° for straight up (Tauxe, 2010). Values outside this range are corrected by incfix(). |
int |
intensity of the data. |
xh |
orientation of the x axis for the horizontal points: can be 'SN' or 'WE'. |
xv |
orientation of the x axis for the horizontal points: can be 'SN', 'WE' or 'modified' (for the latter the horizontal projection of the vector given by the square root of the addition of the squared horizontal components). |
centre |
logical, whether the [0,0] point should be in the centre of the plot. Is ignored if xlim and/or ylim are defined. |
xlim , ylim |
the x and y minimal limits. The actual limits can change to keep a x/y ratio of 1. |
unit |
the tick interval. |
xlab , ylab |
the titles for the axes. |
labels |
a character vector of labels to add to each point. |
nlabels |
the number of labels to skip (for clarity). |
h , v , f , t , l |
list of graphical parameters to feed the graphical functions: h, v and f are fed to points() for the horizontal, vertical and first points respectively; t is fed to the text() for the labels and l is fed to lines() for the lines joining each horizontal and vertical points. See ?points, ?text and ?lines help page for the possible arguments. See ?merge_list for further information. |
anchored |
logical, whether the lines should be anchored to the [0,0] point. |
style |
the style of the plot: 'branches', 'box0', 'box1', or 'box2'. The boxes are advised when zooming using xlim and/or ylim. |
tcl |
The length of tick marks (see par() help page). |
orientation |
logical, whether to add captions indicating the orientation of the plot. |
scientific |
logical or NA, whether have scientific notation (e.g. -1.0E-06) or not (e.g. 0.00015). If NA, R will be left only judge. |
decimals |
the number of decimals if scientific is T or F. Having not enough decimals can lead to override the unit parameter, but the tick labels will be correctly aligned. |
add |
logical, whether to add the plot to an existing plot. |
Details
By default horizontal projection is made of black points, vertical of white points.
References
Tauxe, L., 2010. Essentials of Paleomagnetism. University of California Press.
See Also
Examples
zd <- zeq_example
ori <- par()$mfrow
par(mfrow = c(1,2))
zijderveld(dec = zd$Dec, inc = zd$Inc, int = zd$Int,
xh = "WE", unit = 10^-5)
zijderveld(dec = zd$Dec, inc = zd$Inc, int = zd$Int,
style = "box1", scientific = FALSE, decimals = 5,
labels = zd$Treat, nlabels = 2)
par(mfrow = ori)