Road network data from OSM for testing od_disaggregate.
library(sf)
nrow(od_road_network)
#> [1] 667
head(od_road_network)
#> Simple feature collection with 6 features and 11 fields
#> Geometry type: LINESTRING
#> Dimension: XY
#> Bounding box: xmin: -1.576616 ymin: 53.76156 xmax: -1.506878 ymax: 53.81465
#> Geodetic CRS: WGS 84
#> osm_id name highway waterway aerialway barrier man_made maxspeed
#> 5 1709456 Hanover Way tertiary <NA> <NA> <NA> <NA> 30 mph
#> 88 3999930 Wakefield Road primary <NA> <NA> <NA> <NA> 40 mph
#> 89 3999931 <NA> primary <NA> <NA> <NA> <NA> 40 mph
#> 90 3999934 Queen Street tertiary <NA> <NA> <NA> <NA> <NA>
#> 93 4002495 Cardigan Road tertiary <NA> <NA> <NA> <NA> 30 mph
#> 95 4004386 Victoria Road tertiary <NA> <NA> <NA> <NA> <NA>
#> cycleway_left z_order
#> 5 <NA> 4
#> 88 <NA> 7
#> 89 <NA> 7
#> 90 <NA> 4
#> 93 <NA> 4
#> 95 <NA> 4
#> other_tags
#> 5 "lit"=>"yes"
#> 88 "junction"=>"roundabout","lanes"=>"2","lit"=>"yes","ref"=>"A61"
#> 89 "lanes"=>"2","lit"=>"yes","oneway"=>"yes","ref"=>"A61"
#> 90 "lanes"=>"2"
#> 93 <NA>
#> 95 <NA>
#> geometry
#> 5 LINESTRING (-1.55754 53.799...
#> 88 LINESTRING (-1.506881 53.76...
#> 89 LINESTRING (-1.507634 53.76...
#> 90 LINESTRING (-1.516503 53.77...
#> 93 LINESTRING (-1.57529 53.810...
#> 95 LINESTRING (-1.576119 53.81...
plot(od_road_network$geometry)
plot(od_data_zones_min$geometry, lwd = 3, col = NULL, add = TRUE)
#> Warning: 'x' is NULL so the result will be NULL