Skip to contents

Distance to Pole of Rotation

Source: R/coordinates.R
PoR_distance.Rd

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

Usage

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

Value

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) |> head()
#> [1] 30.94452 29.49555 30.61975 30.26368 28.31982 26.60286

# distance form data.frame
PoR_distance(sf::st_drop_geometry(san_andreas), por) |> head()
#> [1] 30.94452 29.49555 30.61975 30.26368 28.31982 26.60286
PoR_distance(sf::st_drop_geometry(san_andreas), por, FUN = orthodrome) |> head()
#> [1] 30.94452 29.49555 30.61975 30.26368 28.31982 26.60286
PoR_distance(sf::st_drop_geometry(san_andreas), por, FUN = vincenty) |> head()
#> [1] 30.94452 29.49555 30.61975 30.26368 28.31982 26.60286