por_transformation function

Conversion between spherical PoR to geographical coordinate system

Transformation from spherical PoR to geographical coordinate system and vice versa

geographical_to_PoR(x, PoR)

PoR_to_geographical(x, PoR)

Arguments

• x: Can be either a "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, a two-column matrix containing the lat and lon coordinates, a sf object, or a raster object.

• PoR: Pole of Rotation. "data.frame" or object of class "euler.pole"

  containing the geographical coordinates of the Euler pole

Returns

object of same type of x with the transformed coordinates. If x

was a data.frame, transformed coordinates are named 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)
head(PoR_to_geographical(san_andreas.por, por))