README.md

swissgd

This R package is an interface to parts of the Geo-Information Platform of the Swiss Confederation and work-in-progress. It provides functions to search and download geodata from data.geo.admin.ch and wraps around the Spatial Temporal Asset Catalog (STAC) API.

The acquisition and use of data or services is free of charge, subject to the provisions on fair use. For more information, please see the Terms of Use.

Installation

Install from GitHub for a regularly updated version (latest: 0.1.0):

install.packages("devtools")
devtools::install_github("zumbov2/swissgd")

Functions

get_available_geodata

Retrieves the names of all available datasets (currently 545) and returns the function (download_geodata() or get_stac_assets()) to obtain them.

swissgd::get_available_geodata()

#> # A tibble: 545 x 2
#>    name                                     retrieval_function         
#>    <chr>                                    <chr>                      
#>  1 ch.are.agglomerationsverkehr             swissgd::download_geodata()
#>  2 ch.are.alpenkonvention                   swissgd::download_geodata()
#>  3 ch.are.belastung-personenverkehr-bahn    swissgd::download_geodata()
#>  4 ch.are.belastung-personenverkehr-strasse swissgd::download_geodata()
#>  5 ch.are.erreichbarkeit-miv                swissgd::download_geodata()
#>  6 ch.are.erreichbarkeit-oev                swissgd::download_geodata()
#>  7 ch.are.gemeindetypen                     swissgd::download_geodata()
#>  8 ch.are.gueteklassen_oev                  swissgd::download_geodata()
#>  9 ch.are.landschaftstypen                  swissgd::download_geodata()
#> 10 ch.are.reisezeit-agglomerationen-miv     swissgd::download_geodata()
#> # ... with 535 more rows

search_geodata

Searches for matches to argument pattern within the names of available datasets on the geo-information platform.

swissgd::search_geodata(pattern = "ÖV")

#> # A tibble: 6 x 2
#>   name                                     retrieval_function         
#>   <chr>                                    <chr>                      
#> 1 ch.are.erreichbarkeit-oev                swissgd::download_geodata()
#> 2 ch.are.gueteklassen_oev                  swissgd::download_geodata()
#> 3 ch.are.reisezeit-agglomerationen-oev     swissgd::download_geodata()
#> 4 ch.are.reisezeit-oev                     swissgd::download_geodata()
#> 5 ch.bav.haltestellen-oev                  swissgd::download_geodata()
#> 6 ch.bav.kataster-belasteter-standorte-oev swissgd::download_geodata()

swissgd::search_geodata("Alti")

#> # A tibble: 1 x 2
#>   name                     retrieval_function        
#>   <chr>                    <chr>                     
#> 1 ch.swisstopo.swissalti3d swissgd::get_stac_assets()

show_metadata and show_preview

These functions can be used to visit the entry of a dataset in the Swiss Geometadata Catalogue and preview the data on mapping platform of the Swiss Confederation.

swissgd::show_metadata("ch.are.erreichbarkeit-oev")
swissgd::show_preview("ch.are.erreichbarkeit-oev")

download_geodata

Downloads datasets available directly from the geo-information platform. Other data can be obtained via the STAC API. See the functions below.

swissgd::download_geodata("ch.swisstopo.swissboundaries3d-gemeinde-flaeche.fill")

get_stac_collections

Displays a description of all data provided via the Spatial Temporal Asset Catalog (STAC) REST Interface on the geo-information platform (currently 10).

swissgd::get_stac_collections() %>% 
  dplyr::select(1:3)

#> # A tibble: 11 x 3
#>    title                   id                  description                      
#>    <chr>                   <chr>               <chr>                            
#>  1 swissALTI3D             ch.swisstopo.swiss~ swissALTI3D is an extremely prec~
#>  2 swissTLM3D              ch.swisstopo.swiss~ swissTLM3D is the large-scale to~
#>  3 swissBUILDINGS3D 2.0    ch.swisstopo.swiss~ swissBUILDINGS3D 2.0 is a vector~
#>  4 SWISSIMAGE 10 cm, digi~ ch.swisstopo.swiss~ The orthophoto mosaic SWISSIMAGE~
#>  5 National Map 1:10'000 ~ ch.swisstopo.lande~ The 1:10,000 national map is swi~
#>  6 National Map 1:50'000   ch.swisstopo.pixel~ The National Map 1:50,000 is a t~
#>  7 National Map 1:25'000   ch.swisstopo.pixel~ The National Map 1:25,000 is a t~
#>  8 National Map 1:100'000  ch.swisstopo.pixel~ The National Map 1:100,000 is a ~
#>  9 National Map 1:200'000  ch.swisstopo.pixel~ The National Map 1:200,000 is a ~
#> 10 swissSURFACE3D          ch.swisstopo.swiss~ swissSURFACE3D models all natura~
#> 11 swissSURFACE3D Raster   ch.swisstopo.swiss~ swissSURFACE3D Raster is a digit~

get_stac_assets

Calls the STAC API on the geo-information platform and returns the download links to geo-specific assets. Here the aerial photo with a ground resolution of 10 cm of a part of my hometown Aarau is queried. WGS84, LV03 and LV95 coordinates are possible.

swissgd::get_stac_assets(
    collection_id = "ch.swisstopo.swissimage-dop10", 
    lon = 645685, 
    lat = 249287
    )

#> # A tibble: 2 x 7
#>   type      href      created  updated  `eo:gsd` `proj:epsg` `checksum:multihas~
#>   <chr>     <chr>     <chr>    <chr>       <dbl>       <int> <chr>              
#> 1 image/ti~ https://~ 2021-02~ 2021-02~      2          2056 122053DCAE0197F524~
#> 2 image/ti~ https://~ 2021-02~ 2021-02~      0.1        2056 12208825F4C065E9FA~

download_stac_assets

Downloads assets from the STAC API that were previously queried with a get_stac_assets() call.

res <- swissgd::get_stac_assets(
  collection_id = "ch.swisstopo.swissimage-dop10",
  lon = 645685,
  lat = 249287
  )

swissgd::download_stac_assets(res)

Examples

Spatial distribution of place name suffixes

Setup

Idea: Examine and visualise the spatial distribution of place name suffixes using spatial kernel density estimation. Datasets: ch.swisstopo.swissnames3d, ch.swisstopo.swissboundaries3d-land-flaeche.fill
Packages: sf, raster, btb, ggplot2 and friends
Script: ex1_swissnames.R

Some Results