Visualisierung Zeckenstiche

Einführung

In dieser Aufgabe sollen die Studierenden lernen wie man Geodaten mithilfe der Packages geopandas und folium darstellen kann. Dabei werden originale Zeckenstiche und durch eine Monte Carlo simulierte Zeckenstiche dargestellt. Informationen zu den verwendeten Datensätzen finden Sie in Table 1.

Visualisierung einer Zecke. Bild von brgfx auf Freepik.

Methodik

Zuerst werden die Packages geladen und die folgende Datensätze eingelesen.

Table 1: Verwendete Datensätze

Name	Datentyp	Herkunft	Koordinatensystem	Beschreibung
zeckenstiche_original	gpkg	Zecke App	LV95	Diese Zeckenstiche wurden von März 2015 bis Juli 2019 im Raum Zürich über die App Zecke erfasst. Aus diesem Datensatz wurden random 100 Zeckenstiche ausgewählt und für diese Übung um 100 m verschoben.
zeckenstiche_simuliert	gpkg	FG Geoinformatik	LV95	Simulierte Zeckenstiche aus der Monte Carlo Simulation. Dabei wurden pro Zeckenstich 50 Runs mit 3 unterschiedlichen Distanzen (10m, 100m, 1000m) durchgeführt.
wald	gpkg	Swisstopo	LV95	Dieser Datensatz stammt aus dem swissTLM3D und beinhaltet die Waldflächen und nicht-Waldflächen im Raum Zürich.

Mithilfe von Folium wird eine Base Map erstellt, welche als Mittelpunkt die Stadt Zürich hat. Als es nächstes werden die Layer wald, zeckenstiche_original und zeckenstiche_simuliert hinzugeführt.

Base Map erstellen

# Create a map for all radien
## Basemap mit Koordinaten von Zürich
# Zurich latitude longitude: 47.36667, 8.55
map = folium.Map(location=[47.36667, 8.55], zoom_start = 13, tiles = "cartodbpositron")

Wald Layer hinzufügen

# Wald zur Karte hinzufügen
is_wald = wald[wald["Wald_text"] == "ja"]
is_wald.to_crs(4326).explore(
    m = map,           # pass the map object
    color = "green",   # use green color on all polygons
    tooltip = False,   # hide tooltip
    highlight = False, # disable highlighting
    name = "Wald"      # name of the layer in the map
);

Original Zeckenstiche hinzufügen

# Orginal Zeckenstiche zur Karte hinzufügen
zeckenstiche_original_gpd.to_crs(4326).explore(
    m = map,                                     # pass the map object
    color = "red",                               # use red color on all points
    marker_kwds = dict(radius = 1, fill = True), # make marker radius 1px with fill
    tooltip = "ID",                              # show "ID" column in the tooltip
    name = "Original Zeckenstiche"               # name of the layer in the map
);

Simulierte Zeckenstiche hinzufügen

# Visualisieren der Simulierten Zeckenstiche
zeckenstiche_sim_gpd.to_crs(4326).explore(
    m = map,                                     # pass the map object
    column = "Radius",                           # select the column that should be used for coloring
    cmap = "viridis",                            # use color palette viridis for coloring of the column values
    legend = True,                               # show legend
    legend_kwds = dict(colorbar = False),        # remove colorbar for legend
    marker_kwds = dict(radius = 1, fill = True), # make marker radius 1px with fill
    tooltip = ["ID", "Radius", "Run"],           # show ID, Radius and Run in the tooltip
    name = "Simulation Zeckenstiche"             # name of the layer in the map
);

Layer Control hinzufügen

folium.LayerControl().add_to(map); #Hier kannst du ein ; machen, damit die Message von Folium (<folium.map.LayerControl at ..) unterdrückt wird. 

Resultat

Make this Notebook Trusted to load map: File -> Trust Notebook