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 = units::set_units(10, "km"),
crs = sf::st_crs("EPSG:26915")
)
Run the code above in your browser using DataLab