PoR_distance function

Distance to Pole of Rotation

Retrieve the (angular) distance to the PoR (Euler pole).

PoR_distance(x, PoR, FUN = orthodrome)

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

• FUN: function to calculate the great-circle distance. orthodrome(), haversine() (the default), or vincenty().

Returns

numeric vector

Examples

data("nuvel1")
por <- subset(nuvel1, nuvel1$plate.rot == "na") # North America relative to Pacific plate
data("san_andreas")

# distance form sf object
PoR_distance(san_andreas, por)

# distance form data.frame
PoR_distance(sf::st_drop_geometry(san_andreas), por)
PoR_distance(sf::st_drop_geometry(san_andreas), por, FUN = orthodrome)
PoR_distance(sf::st_drop_geometry(san_andreas), por, FUN = vincenty)