tectonicr (version 0.4.8)

PoR_distance: Distance to Pole of Rotation

Description

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

Usage

PoR_distance(x, PoR, FUN = orthodrome)

Value

numeric vector

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().

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) |> head()

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