Les fonctions st_read() et st_write() permettent d’importer et d’exporter de nombreux types de fichiers.
Import
Les lignes suivantes importent la couche des communes du département du Lot situé dans le fichier geopackage lot46.gpkg.
#> Linking to GEOS 3.9.0, GDAL 3.2.2, PROJ 7.2.1; sf_use_s2() is TRUE
com <- st_read("data/lot46.gpkg", layer = "commune")
#> Reading layer `commune' from data source
#> `/home/tim/Documents/prz/ElementR_Geomatique_2023/tutoriel/data/lot46.gpkg'
#> using driver `GPKG'
#> Simple feature collection with 313 features and 12 fields
#> Geometry type: MULTIPOLYGON
#> Dimension: XY
#> Bounding box: xmin: 539668.5 ymin: 6346290 xmax: 637380.9 ymax: 6439668
#> Projected CRS: RGF93 / Lambert-93
La fonction st_layers() permet d’avoir un aperçu des couches présentes dans un fichier geopackage.
st_layers("data/lot46.gpkg")
#> Driver: GPKG
#> Available layers:
#> layer_name geometry_type features fields crs_name
#> 1 departement Multi Polygon 96 5 RGF93 / Lambert-93
#> 2 route Line String 16096 1 RGF93 / Lambert-93
#> 3 commune Multi Polygon 313 12 RGF93 / Lambert-93
#> 4 restaurant Point 4734 1 RGF93 / Lambert-93
Export
Les lignes suivantes exportent l’objet com dans la couche commune du geopackage com.gpkg dans le dossier data.
st_write(obj = com, dsn = "data/com.gpkg", layer = "commune")
#> Writing layer `commune' to data source `data/com.gpkg' using driver `GPKG'
#> Writing 313 features with 12 fields and geometry type Multi Polygon.