Werken met gpkg
Een van de dingen die ik telkens vergeet is hoe je makkelijk kunt omgaan met geopackage-bestanden (.gpkg). Het is het standaarduitwisselformaat van Geonovum dat op de pas-toe-of-leg-uit lijst staat. Zie ook: https://www.geonovum.nl/geo-standaarden/geopackage.
Eigenlijk is het vrij simpel. De sf-library is hiervoor geschikt.
Om te zien welke lagen er in een .gpkg opgeslagen zitten, gebuik je het commando:
library(sf)
layers <- st_layers("museum.gpkg")
Dit levert een list op, met een aantal velden. Eén van die velden is name. Dit geeft de namen weer van de beschikbare lagen in het .gpkg bestand.
layers$name
[1] "museum_afbakening_gebied" "adressen_depont"
[3] "adressen_textielmuseum" "looproutes_depont"
[5] "loopader_depont" "loopader_textielmuseum"
[7] "isochroon_depont" "isochroon_textielmuseum"
[9] "loopaders_gezamelijk" "looproutes_textielmuseum"
[11] "bestemming" "afbakening4326"
[13] "qgis_projects"
Het daadwerkelijk inlezen van een laag in R gebeurt als volgt:
library(sf)
isotextiel <- st_read("museum.gpkg", layer="isochroon_textielmuseum")