Appendix A: Working with Austrian data

Table of content for chapter 04

Chapter section list

R Code A.1 : Numbered R Code Title

Listing / Output A.1: Listing title

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/_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 A.2 : Numbered R Code Title

Listing / Output A.2: Listing title

ggplot2::ggplot(data = at_regions) +  
  ggplot2::geom_sf(fill = NA)

R Code A.3 : Numbered R Code Title

Listing / Output A.3: Listing title

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/_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

ggplot2::ggplot(data = at_pol_bezirk) +  
  ggplot2::geom_sf(fill = NA)