See https://github.com/pedalea/pctSantiago folder for generation code

Examples

# 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...