See https://github.com/pedalea/pctSantiago folder for generation code
# u = "https://github.com/pedalea/pctSantiago/releases/download/0.0.1/od_agg_zone_sub.Rds"
# download.file(u, destfile = "od_agg_zone_sub.Rds")
# desire_lines = readRDS("od_agg_zone_sub.Rds")
santiago_zones
#> Simple feature collection with 75 features and 1 field
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: -70.69203 ymin: -33.47766 xmax: -70.58239 ymax: -33.403
#> old-style crs object detected; please recreate object with a recent sf::st_crs()
#> Geodetic CRS: WGS 84
#> First 10 features:
#> old-style crs object detected; please recreate object with a recent sf::st_crs()
#> old-style crs object detected; please recreate object with a recent sf::st_crs()
#> geo_code geometry
#> 259 259 POLYGON ((-70.65174 -33.431...
#> 260 260 POLYGON ((-70.65008 -33.443...
#> 261 261 POLYGON ((-70.64562 -33.437...
#> 262 262 POLYGON ((-70.64643 -33.477...
#> 263 263 POLYGON ((-70.64913 -33.459...
#> 264 264 POLYGON ((-70.64914 -33.459...
#> 265 265 POLYGON ((-70.65058 -33.450...
#> 266 266 POLYGON ((-70.65058 -33.450...
#> 292 292 POLYGON ((-70.60615 -33.419...
#> 293 293 POLYGON ((-70.60673 -33.415...