if (FALSE) {
# load packages
library(duckspatial)
library(sf)
# create a duckdb database in memory (with spatial extension)
conn <- ddbs_create_conn(dbdir = "memory")
# read data
countries_sf <- st_read(system.file("spatial/countries.geojson", package = "duckspatial"))
# store in duckdb
ddbs_write_vector(conn, countries_sf, "countries")
# combine all geometries into one
ddbs_combine(conn = conn, "countries")
# combine without using a connection
ddbs_combine(countries_sf)
# store result in a new table
ddbs_combine(conn = conn, "countries", name = "countries_combined")
}
Run the code above in your browser using DataLab