Learn R Programming

tectonicr (version 0.4.6)

PoR2Geo_azimuth: Azimuth Conversion From PoR to Geographical Coordinate Reference System

Description

Conversion of PoR azimuths into geographical azimuths

Usage

PoR2Geo_azimuth(x, PoR, axial = TRUE)

Value

numeric vector of transformed azimuths (in degrees)

Arguments

x

data.frame containing the PoR equivalent azimuths (azi.PoR), and either the geographical coordinates of the point(s) or the PoR-equivalent coordinates.

PoR

data.frame containing the geographical location of the Euler pole (lat, lon)

axial

logical. Whether the azimuth is axial (0-180) or directional (0-360).

References

Stephan, T., Enkelmann, E., and Kroner, U. "Analyzing the horizontal orientation of the crustal stress adjacent to plate boundaries". Sci Rep 13. 15590 (2023). tools:::Rd_expr_doi("10.1038/s41598-023-42433-2").

See Also

PoR_shmax()

Examples

Run this code
data("nuvel1")
# North America relative to Pacific plate:
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
data("san_andreas")
san_andreas$azi.PoR <- PoR_shmax(san_andreas, PoR)

# convert back to geo CRS
PoR2Geo_azimuth(san_andreas, PoR) |> head()

Run the code above in your browser using DataLab