Learn R Programming

tectonicr (version 0.4.4)

PoR2Geo_azimuth: Azimuth conversion from PoR to geographical coordinate reference system

Description

Conversion of PoR azimuths into geographical azimuths

Usage

PoR2Geo_azimuth(x, PoR)

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)

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)

Run the code above in your browser using DataLab