Retrieve the PoR equivalent coordinates of an object
Arguments
- x
sf
ordata.frame
containing lat and lon coordinates (lat
,lon
)- PoR
Pole of Rotation.
"data.frame"
or object of class"euler.pole"
containing the geographical coordinates of the Euler pole
Examples
data("nuvel1")
por <- subset(nuvel1, nuvel1$plate.rot == "na") # North America relative to Pacific plate
data("san_andreas")
san_andreas.por_sf <- PoR_coordinates(san_andreas, por)
head(san_andreas.por_sf)
#> lat.PoR lon.PoR
#> 1 56.93219 -82.67182
#> 2 58.01202 -74.48278
#> 3 57.93239 -74.80319
#> 4 57.83707 -75.51138
#> 5 58.17336 -76.62093
#> 6 57.84485 -76.16849
san_andreas.por_df <- PoR_coordinates(sf::st_drop_geometry(san_andreas), por)
head(san_andreas.por_df)
#> lat.PoR lon.PoR
#> 1 56.93219 -82.67182
#> 2 58.01202 -74.48278
#> 3 57.93239 -74.80319
#> 4 57.83707 -75.51138
#> 5 58.17336 -76.62093
#> 6 57.84485 -76.16849