# Query without boundary definition
q <- ohsome_query(
"elements/count/groupBy/boundary",
filter = "building=*",
time = "2022-01-01"
)
# Use franconia from the mapview package as bounding polygons
# \donttest{
set_boundary(q, mapview::franconia, digits = 4)
# }
# Use the bounding box of franconia
# \donttest{
set_boundary(q, sf::st_bbox(mapview::franconia))
# }
if (FALSE) {
# Get bounding box of the city of Kigali from OSM
set_boundary(q, osmdata::getbb("Kigali"))
}
# Definition of two named bounding circles
set_boundary(q, c("Circle 1:8.6528,49.3683,1000", "Circle 2:8.7294,49.4376,1000"))
Run the code above in your browser using DataLab