`line_via()` returns geometries that lack a CRS

[Robinlovelace]

library(stplanr)
library(sf)
#> Linking to GEOS 3.9.0, GDAL 3.2.1, PROJ 7.2.1
l <- flowlines_sf[2:4, ]
p <- destinations_sf
lv <- line_via(l, p)
sf::st_crs(l)
#> Coordinate Reference System:
#>   User input: EPSG:4326 
#>   wkt:
#> GEOGCRS["WGS 84",
#>     DATUM["World Geodetic System 1984",
#>         ELLIPSOID["WGS 84",6378137,298.257223563,
#>             LENGTHUNIT["metre",1]]],
#>     PRIMEM["Greenwich",0,
#>         ANGLEUNIT["degree",0.0174532925199433]],
#>     CS[ellipsoidal,2],
#>         AXIS["geodetic latitude (Lat)",north,
#>             ORDER[1],
#>             ANGLEUNIT["degree",0.0174532925199433]],
#>         AXIS["geodetic longitude (Lon)",east,
#>             ORDER[2],
#>             ANGLEUNIT["degree",0.0174532925199433]],
#>     USAGE[
#>         SCOPE["unknown"],
#>         AREA["World"],
#>         BBOX[-90,-180,90,180]],
#>     ID["EPSG",4326]]
sf::st_crs(lv$leg_orig)
#> Coordinate Reference System: NA
sf::st_crs(lv$leg_via)
#> Coordinate Reference System: NA
sf::st_crs(lv$leg_dest)
#> Coordinate Reference System: NA

^{Created on 2021-07-20 by the reprex package (v2.0.0)}