Azimuth conversion from PoR to geographical coordinate reference system
Source:R/model_shmax.R
PoR2Geo_azimuth.Rd
Conversion of PoR azimuths into geographical azimuths
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). doi:10.1038/s41598-023-42433-2 .
Examples
data("nuvel1")
# North America relative to Pacific plate:
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
data("san_andreas")
head(san_andreas$azi)
#> [1] 36 7 169 1 23 160
san_andreas$azi.PoR <- PoR_shmax(san_andreas, PoR)
res.geo <- PoR2Geo_azimuth(san_andreas, PoR)
head(res.geo)
#> [1] 36 7 169 1 23 160