Building data from OSM for testing od_disaggregate.
nrow(od_data_buildings)
#> [1] 303
head(od_data_buildings)
#> Simple feature collection with 6 features and 2 fields
#> Geometry type: MULTIPOLYGON
#> Dimension: XY
#> Bounding box: xmin: -1.5565 ymin: 53.79338 xmax: -1.537664 ymax: 53.80027
#> Geodetic CRS: WGS 84
#> osm_way_id building geometry
#> 1 23150505 retail MULTIPOLYGON (((-1.540118 5...
#> 2 23162531 commercial MULTIPOLYGON (((-1.546567 5...
#> 3 27767021 commercial MULTIPOLYGON (((-1.5565 53....
#> 4 28516780 office MULTIPOLYGON (((-1.553886 5...
#> 5 28516878 commercial MULTIPOLYGON (((-1.550857 5...
#> 6 28516879 commercial MULTIPOLYGON (((-1.552283 5...
plot(od_data_buildings$geometry)
plot(od_data_zones_min$geometry, lwd = 3, col = NULL, add = TRUE)
#> Warning: 'x' is NULL so the result will be NULL