Übung 5

Integration von Geodaten

• Bisher haben wir alle Datensätze einzeln betrachtet
• Wenn wir alle Datensätze ins gleiche Bezugssystem bringen, können wir diese integrieren bzw. überlagern
• Überlagern kann heissen:
  • gemeinsam visualisieren
  • Information übertragen

Vorbereitung:

• Starte ein neues Script Uebung_5.R
• Importiere darin alle räumlichen libraries

library("tmap")
library("sf")
library("terra")

# tmap_options(check.and.fix = TRUE)
tmap_mode("plot")

Übung 5.1

Importiere die Datensätze ausserberg.gpkg (aus Übung 2) sowie dhm200_2056.tif (aus Übung 3, data-processed/ausserberg.gpkg bzw. data-processed/dhm200_2056.tif)

Übung 5.2

Überlagere die beiden Datensätze in einem tmap-Plot, indem du diese mit + verkettest.

Übung 5.3

• Mit crop() können wir ein Raster auf den “extent” von einem Vektor Datensatz zuschneiden
• Schneide dhm200 auf den extent von ausserberg zu
• Visualisiere das resultierende Raster mit tmap (wieder gemeinsam mit ausserberg)

Übung 5.4

• die Auflösung des Raster Datensatzes ist zu grob!!
• Lösung: Hoch aufgelöster Datensatz dhm25 (aus Übung 4) und einlesen (zip-File: processed/dhm25_2056.tif)
• wiederhole das Zuschneiden mittels crop sowie das Visualisieren mittels tmap

Input: Rasterwerte extrahieren

• bisher haben wir zwei Datensätze (Raster und Vektor) visuell überlagert
• nächster Schritt: Information von Raster → Punkt Datensatz übertragen
• dazu müssen wir ausserberg von einem sp- in ein SpatVector Objekt konvertieren
• danach können wir das SpatVector Objekt gemeinsam mit extract verwenden

Übung 5.5

• Verwende die Funktion extract mit ausserberg um die Höhenwerte aus dhm25 zu extrahieren
• Speichere den output in einer Variabel und beguteachte diese

Übung 5.6

Spiele die Höheninformation aus extract zurück in ausserberg.

## Übung 5.7

Visualisiere nun ausserberg und Färbe die Punkte nach ihrer Höheninformation ein.

Input: Vektordaten zuschneiden

• nun wollen wir zwei Vektordatensätze miteinander verschneiden
• Ausgangslage:
  • wir verfügen über einen TWW Datensatz der Schweiz
    • URL: https://www.bafu.admin.ch/bafu/de/home/themen/biodiversitaet/zustand/karten/geodaten.html)
    • Shorturl: https://bit.ly/3CqNRKT
  • wir verfügen über den Gemeindelayer der Schweiz
    • URL: https://data.geo.admin.ch/ch.swisstopo.swissboundaries3d/swissboundaries3d_2023-01/swissboundaries3d_2023-01_2056_5728.gpkg.zip)
    • Shorturl: https://bit.ly/47wCOzo)
  • wir wollen alle TWW Flächen innerhalb der Gemeinde Landquart erhalten

Übung 5.8

• Lade diese beiden Datensätze herunter und importiere sie in R (swissboundaries Hoheitsgebiet)
• Transformiere sie in ESPG 2056

Übung 5.9

Erstelle ein neues Objekt landquart, welches nur die Gemeinde Landquart beinhaltet und visualisiere diese.

Übung 5.10

Überlagere die TWW Flächen mit der Gemeindegrenze von Landqart.

Übung 5.11

• Verwende die Funktion st_intersection() um die TWW-Flächen auf die Gemeindegrenze von Landquart zu zuschneiden.
• Visualisiere das Resultat

Input: Vektordaten selektieren

Mit st_intersection haben wir TWW Flächen verschnitten, da st_intersetion die Schnittmenge beider Polygone nimmt

• Alternativ können wir alle TWW Flächen selektieren, die mindestens Teilweise innerhalb des Gemeindegebietes liegen

tww_landquart2 <- tww[landquart, ]

Übung 5.12

• Selektiere die TWW Flächen, welche sich zumindest Teilweise in der Gemeinde Landquart befinden und speichere den Output als tww_landquart2
• Visualisiere das Resultat mit tmap
• Vergleiche tww_landquart2 mit tww_landquart. Wie unterscheiden sich diese?

Übung 5.13