PoR_map: Plot data in PoR map

In tectonicr: Analyzing the Orientation of Maximum Horizontal Stress

Plot data in PoR map

Description

Plot data in PoR map

Usage

PoR_map(
  x,
  PoR,
  pb = NULL,
  type = c("none", "in", "out", "right", "left"),
  show.deviation = FALSE,
  ...
)

Arguments

x, pb	sf objects of the data points and the plate boundary geometries in the geographical coordinate system
PoR	Pole of Rotation. "data.frame" or object of class "euler.pole" containing the geographical coordinates of the Pole of Rotation
type	Character. Type of plate boundary (optional). Can be "out", "in", "right", or "left" for outward, inward, right-lateral, or left-lateral moving plate boundaries, respectively. If "none" (the default), only the PoR-equivalent azimuth is returned.
show.deviation	logical. Whether the data should be color-coded according to the deviation from the prediction, or according to the stress regime? Is ignored if type=='none'.
...	optional arguments passed to tectonicr.colors()

Value

plot

See Also

PoR_shmax(), axes(), tectonicr.colors()

Examples

data("nuvel1")
na_pa <- subset(nuvel1, nuvel1$plate.rot == "na")

data("plates")
plate_boundary <- subset(plates, plates$pair == "na-pa")

data("san_andreas")
PoR_map(san_andreas, PoR = na_pa, pb = plate_boundary, type = "right", show.deviation = TRUE)

tectonicr documentation built on Sept. 11, 2024, 6:05 p.m.