Übung 4

Komplexe Rasterdaten

Satelliten und Drohnen nehmen meist verschiedene Spektren von Elektromagnetischen Wellen auf
diese Spektren werden in unterschiedlichen Datensätzen abgespeichert
diese Datensätze müssen wieder zusammengefügt werden um ein Gesamtbild zu erhalten
Beispiel: Rot, Grün und Blau werte fügen sich zu einem Farbluftbild zusammen

Übung 4.1

Ladet euch einen Ausschnitt aus dem Datensatz swissimage 10 von Swisstopo herunter: https://www.swisstopo.admin.ch/de/geodata/images/ortho/swissimage10.html
- Shortlink: https://bit.ly/40Fy0Wj
Enzipped den Inhalt in euer RStudio Projekt und schaut den Inhalt an
Öffnet das File in QGIS und schaut euch die Metadaten an
Was ist das Koordinatenbezugssystem? Wie hoch ist die räumliche Auflösung?

Übung 4.2

Erstelle ein neues R Script mit dem Namen Uebung_4.R
Lade die libraries sf, tmap und terra.
Importiere den Swissimage Datensatz
Weise dem importierten Datensatz das Korrekte Koordinatenbezugssystem zu
Schau dir den Datensatz in der Konsole sowie mit plot() an

plot(swissimage)

swissimage

class       : SpatRaster 
dimensions  : 9480, 14000, 3  (nrow, ncol, nlyr)
resolution  : 25, 25  (x, y)
extent      : 2484375, 2834375, 1062000, 1299000  (xmin, xmax, ymin, ymax)
coord. ref. : CH1903+ / LV95 (EPSG:2056) 
source      : SI25-2012-2013-2014.tif 
colors RGB  : 1, 2, 3 
names       : SI25-2012-2013-2014_1, SI25-2012-2013-2014_2, SI25-2012-2013-2014_3

Input: RGB Plots mit `tmap`

Um ein rgb Datensatz mit tmap zu plotten, verwenden wir nicht mehr tm_raster() sondern tm_rgb

tmap_mode("plot")

tmap mode set to plotting

swissimage_10 <- aggregate(swissimage, fact = 10)

tm_shape(swissimage_10) + 
  tm_rgb()

stars object downsampled to 1215 by 823 cells. See tm_shape manual (argument raster.downsample)

Input

Heute haben wir das Höhenmodell dhm200 importiert
Höhenmodell mit 200m Auflösung (→ grob! )
swisstopo stellt zusätzlich das mit 25m Auflösung zur Verfügung (https://www.swisstopo.admin.ch/de/geodata/height/dhm25.html)
- Shortlink: https://www.swisstopo.admin.ch/de/geodata/height/dhm25.html
durch die höhere Auflösung dauert das transformieren in ein neues Koordinatensystem etwas länger

Übung 4.4 (Optional)

Ladet euch das dhm25 mit 25m Auflösung herunter (https://www.swisstopo.admin.ch/de/geodata/height/dhm25.html)
importiert es in R
setzt das korrekte CRS
transformiert es in EPSG 2056 und verwendet dabei folgende Optionen:
- mit filename = den Output direkt in ein File speichern
- mit progress = TRUE den Fortschritt anzeigen lassen
visualisiert es mit tmap

Übung 4.5 (Optional und Open End)

Suche dir auf den gängigen Portalen (s.u.) einen spannenden Datensatz und visualisiere diesen