The goal of isonoi is to provide example data and functions to demonstrate the concept of ‘iso voronoi’ polygons.
It is designed to be a prototype implementation of ideas outlined and demonstrated with example data and reproducible code in https://github.com/isonoi/isonoi
Give it a spin with the reproducible examples shown below.
remotes::install_github("isonoi/isonoi")
# TODO: remove this and keep only library call
devtools::load_all()
#> ℹ Loading isonoi
# library(isonoi)
library(sf)
x = net_oldenburg_raw
points = points_oldenburg
grid = iso_grid(x)
plot(grid, border = "grey")
plot(sf::st_geometry(x), add = TRUE)
plot(points, col = "red", add = TRUE, pch = 18, cex = 2)
#> Warning in plot.sf(points, col = "red", add = TRUE, pch = 18, cex = 2):
#> ignoring all but the first attribute
The 3 figures below show outputs from the packages next to the standard voronoi diagram.
library(tmap)
poly = sf::st_convex_hull(sf::st_union(x))
vnoi = iso_vnoi(points, poly)
#> Warning in st_voronoi.sfc(sf::st_union(points), poly): st_voronoi does not
#> correctly triangulate longitude/latitude data
#> Warning: attribute variables are assumed to be spatially constant throughout
#> all geometries
vgrid = iso_vgrid(points, grid)
#> lon-lat points
net = sfnetworks::as_sfnetwork(x)
vsfn = iso_sfn(points, net, grid)
#> lon-lat points
#> lon-lat points
#> Warning: st_centroid assumes attributes are constant over geometries
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> lon-lat points
#> Warning in shortest_paths(x, from, to, weights = weights, output = "both", : At
#> core/paths/dijkstra.c:446 : Couldn't reach some vertices.
#> Joining with `by = join_by(name)`
m1 = tm_shape(vnoi) + tm_fill("name")
m2 = tm_shape(vgrid) + tm_fill("name")
m3 = tm_shape(vsfn) + tm_fill("name")
tmap_arrange(m1, m2, m3, nrow = 1)
#> Some legend labels were too wide. These labels have been resized to 0.24, 0.62, 0.29. Increase legend.width (argument of tm_layout) to make the legend wider and therefore the labels larger.
#> Some legend labels were too wide. These labels have been resized to 0.24, 0.62, 0.29. Increase legend.width (argument of tm_layout) to make the legend wider and therefore the labels larger.
#> Some legend labels were too wide. These labels have been resized to 0.24, 0.62, 0.29. Increase legend.width (argument of tm_layout) to make the legend wider and therefore the labels larger.