Conversion between spherical PoR to geographical coordinate system

Transformation from spherical PoR to geographical coordinate system and vice versa

Usage

geographical_to_PoR(x, PoR)

PoR_to_geographical(x, PoR)

Arguments

x

"data.frame" containing lat and lon coordinates of a point in the geographical CRS or the lat.PoR, lon.PoR) of the point in the PoR CRS.

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 transformed coordinates (lat.PoR and lon.PoR for PoR CRS, or lat and lon for geographical CRS).

Examples

data("nuvel1")
por <- subset(nuvel1, nuvel1$plate.rot == "na") # North America relative to Pacific plate
data("san_andreas")
san_andreas.por <- geographical_to_PoR(san_andreas, por)
head(san_andreas.por)
#>    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
head(PoR_to_geographical(san_andreas.por, por))
#>     lat     lon
#> 1 38.14 -118.84
#> 2 35.97 -114.71
#> 3 37.93 -118.17
#> 4 38.63 -118.21
#> 5 39.08 -115.62
#> 6 38.58 -112.58