PoR_distance: Distance to Pole of Rotation

In tectonicr: Analyzing the Orientation of Maximum Horizontal Stress

Distance to Pole of Rotation

Description

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)

# 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)

tectonicr documentation built on April 4, 2025, 12:54 a.m.