PoR coordinates

Retrieve the PoR equivalent coordinates of an object

Usage

PoR_coordinates(x, PoR)

Arguments

x

sf or data.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

Value

data.frame with the PoR coordinates (lat.PoR, lon.PoR)

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 59.05548 -85.46364
#> 2 60.50445 -78.16683
#> 3 59.38025 -84.55317
#> 4 59.73632 -85.74431
#> 5 61.68018 -84.31053
#> 6 63.39714 -80.60932
san_andreas.por_df <- PoR_coordinates(sf::st_drop_geometry(san_andreas), por)
head(san_andreas.por_df)
#>    lat.PoR   lon.PoR
#> 1 59.05548 -85.46364
#> 2 60.50445 -78.16683
#> 3 59.38025 -84.55317
#> 4 59.73632 -85.74431
#> 5 61.68018 -84.31053
#> 6 63.39714 -80.60932