Skip to contents

Azimuth Conversion From PoR to Geographical Coordinate Reference System

Source: R/model_shmax.R
PoR2Geo_azimuth.Rd

Conversion of PoR azimuths into geographical azimuths

Usage

PoR2Geo_azimuth(x, PoR, axial = TRUE)

Arguments

x

sf object or a data.frame containing the coordinates of the point(s) (lat, lon columns). x must contain the direction of \(\sigma_\text{Hmax}\) as column azi, its standard deviation (column unc) is optional).

PoR

data.frame or object of class euler.pole containing the geographical coordinates of the Euler pole.

axial

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

Value

numeric vector of transformed azimuths (in degrees)

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). doi:10.1038/s41598-023-42433-2 .

See also

PoR_shmax()

Examples

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()
#> [1] 50 54 24 41 30 27