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_line() |>
profile_length()
Run the code above in your browser using DataLab