Azimuth Between Profile Points
Details
If only two points are given, the azimuth is calculated using triangulation
from the tectonicr
package.
If more than two points are given, the azimuth is calculated using linear
interpolation in the coordinate reference frame given by profile
.
Examples
p1 <- data.frame(lon = -90.8, lat = 48.6) |>
sf::st_as_sf(coords = c("lon", "lat"), crs = "WGS84")
profile_points(p1,
profile.azimuth = 135, profile.length = 10000,
crs = sf::st_crs("EPSG:26915")
) |>
profile_azimuth()
#> Warning: Unit of profile.length not specified. Assuming unit is in meters.
#> 136.7341 [°]