Appendix D: Working with Austrian data

Objectives

This chapter tries to find and to apply Austrian data for map plotting.

Watch out! D.1

I can’t remember from where I got the following Austrian data. Month later I tried to find the original source and learned about the following websites:


The Austrian Map is an interesting website with all Austrian geographical and historical details.

R Code D.1 : Load data for Austrian regions

Listing / Output D.1: Load data for Austrian regions
Code
at_regions <- sf::st_read("_archive/austria/austria-with-regions_.geojson")
#> Reading layer `austria-with-regions_' from data source 
#>   `/Users/petzi/Documents/Meine-Repos/GDSWR-notes-new/_archive/austria/austria-with-regions_.geojson' 
#>   using driver `GeoJSON'
#> Simple feature collection with 9 features and 4 fields
#> Geometry type: MULTIPOLYGON
#> Dimension:     XY
#> Bounding box:  xmin: 9.53099 ymin: 46.37265 xmax: 17.16204 ymax: 49.02135
#> Geodetic CRS:  WGS 84

R Code D.2 : Plot Austrian regions

Listing / Output D.2: Austrian regions
Code
ggplot2::ggplot(data = at_regions) +  
  ggplot2::geom_sf(fill = NA)

R Code D.3 : Load Austrian regions with administrative borders

Listing / Output D.3: Austrian regions with administrative borders
Code
at_folder = "_archive/austria"
pol_bez_folder = "/OGDEXT_POLBEZ_1_STATISTIK_AUSTRIA_20240101/"
pol_bez_file = "STATISTIK_AUSTRIA_POLBEZ_20240101.shp"

at_pol_bezirk <- sf::st_read(
  paste0(at_folder, pol_bez_folder, pol_bez_file)
)
#> Reading layer `STATISTIK_AUSTRIA_POLBEZ_20240101' from data source 
#>   `/Users/petzi/Documents/Meine-Repos/GDSWR-notes-new/_archive/austria/OGDEXT_POLBEZ_1_STATISTIK_AUSTRIA_20240101/STATISTIK_AUSTRIA_POLBEZ_20240101.shp' 
#>   using driver `ESRI Shapefile'
#> Simple feature collection with 117 features and 2 fields
#> Geometry type: MULTIPOLYGON
#> Dimension:     XY
#> Bounding box:  xmin: 112518.2 ymin: 275472 xmax: 685444.5 ymax: 570431.1
#> Projected CRS: MGI / Austria Lambert
Code
ggplot2::ggplot(data = at_pol_bezirk) +  
  ggplot2::geom_sf(fill = NA)