Kartendaten mit Python

In Google Maps kann man Dateien einlesen, in denen Orte oder Flächen markiert werden und kann diese verschicken. Dies kann eine Vielzahl an nützlichen Anwendungen haben:

1. Ferienziele: Mache für dich und deine Freunde eine Weltkarte, auf der deine wichtigsten Ferienziele markiert sind.
2. Reise-Route: Plane deine nächste Reise
3. Weltherrschaft: Du willst dein eigenes Land haben? Stecke die Grenzen auf Google Maps ab und teile es mit deinen Freunden.

Typisches Vorgehen:

1. Gehe auf Google My Maps …
2. … und logge dich ein.
3. Klicke auf „Neue Karte erstellen“ und „Importieren“ …
4. … und ziehe dort ein passendes CSV-File hinein.
5. Dann muss man die „Spalten zur Positionierung …“ auswählen: typischerweise „WKT“ …
6. … und dann die „Spalte zur Benennung …“: typischerweise „Name“

WKT steht für Well-Known Text und ist ein Textformat für räumliche Datenmodelle. Es ermöglicht die Beschreibung geometrischer Objekte wie Punkte, Linien, Polygone und mehr in einem Format, welches sowohl für Menschen wie auch für Computer gut lesbar sind.

Einen Punkt mit Koordinaten $(x,y)$ notiert man darin wie folgt:

"POINT (x y)"

Eine Linie, die von $(x1,y1)$ zu $(x2,y2)$ und von da zu $(x3,y3)$ geht, notiert man so:

"LINESTRING (x1 y1, x2 y2, x3 y3)"

Ein geschlossenes Polygon (Vieleckt) notiert man so:

"POLYGON ((x1 y1, x2 y2, x3 y3, x1 y1))"

In diesem Falle erhält man natürlich ein Dreieck. Beachte hier die Verwendung von doppelten Klammern.

Man kann auch gleichzeitig mehrere Punkte, Linien und Polygone festlegen:

"MULTIPOINT ((x1 y1, x2 y2, x3 y3))"
"MULTILINESTRING ((x1 y1, x2 y2), (x3 y3, x4 y4))""
"MULTIPOLYGON (((x1 y1, x2 y2, x3 y3, x1 y1)),((x4 y4, x5 y5, x6 y6, x4 y4)))"

Um Punkte, Linien und Polygone in Google Maps einzuzeichnen, kann man diese Formen in einem CSV-File festlegen. Dazu benötigt man die entsprechenden Koordinaten ($x:$ longitude/Längengrad, $y:$ latitude/Breitengrad).

Wir betrachten die folgenden vier Ortschaften mit zugehörigen Längen- und Breitengrad:

• Frauenfeld (8.89656423219735, 47.558162001326)
• Weinfelden (9.10636576683266, 47.5662241119746)
• Romanshorn (9.37937276486743, 47.5655209998311)
• Kreuzlingen (9.16907715682578, 47.6525502207935)

tg.csv

Name,WKT
Weinfelden,"POINT (9.10636576683266 47.5662241119746)"
Line FF to WF,"LINESTRING (8.89656423219735 47.558162001326, 9.10636576683266 47.5662241119746)"

Wann und wo benötigt man Anführungs- und Schlusszeichen?

• Für die WKT-Befehle werden diese benötigt. Der Grund ist, dass diese Befehle Kommas beinhalten, welche auch als Trennzeichen (Delimiter) des CSV agieren.
• Für die Namen hingegen sind sie optional, ausser die Namen beinhalten Kommas, dann braucht man auch hier Anführungs- und Schlusszeichen.

Beachte, dass das CSV einen passenden Header benötigt: „Name,WKT“. Dank diesem weiss Google Maps, welche Bedeutung die Werte im CSV haben.

Dieses CSV-File kann man herunterladen, in Google Maps importieren und dann sieht man den Punkt und die Linie am richtigen Ort mit der angegebenen Beschriftung.

Arbeite nicht mit Excel / Google Sheets, sondern mit einem reinen Texteditor wie Visual Studio Code.

Erstelle ein leeres CSV-File mit Header „Name,WKT“ und füge Punkte für die vier Gemeinden Frauenfeld, Weinfelden, Romanshorn, Kreuzlingen hinzu. Verwende die Koordinaten aus der Theorie oben:

Importiere das File dann in Google Maps und stelle sicher, dass alles stimmt.

Welches sind deine Top 3 Orte, die du noch besuchen möchtest?

1. Finde deren Koordinaten …
2. … erstelle ein entsprechendes CSV-File …
3. … und importiere sie auf Google Maps.

Spielregel: Es müssen mind. zwei unterschiedliche Kontinenten vorkommen.

Von nun arbeiten wir wieder mit Python und dem folgenden Datensatz, der die Koordinaten aller Bahnhöfe in der Schweiz beinhaltet: sbb_csv_json.zip.

Auftrag:

• Generiere mithilfe von Python ein CSV-File für Google Maps, welches die Koordinaten einiger Ortschaften St. Gallen, Bern und Schaffhausen beinhaltet.
• Verwende den SBB-Datensatz, um die Koordinaten zu ermitteln.
• Schreibe Code so, dass er auch noch funktioniert, wenn man z.B. St. Gallen durch Winterthur ersetzt.

Tipps

1. Wähle einen Anfangsbuchstaben, z.B. „W“, und speichere diesen in einer Variable.
2. Ermittle nun alle Bahnhöfe, die mit diesem Buchstaben beginnen …
3. und generiere das entsprechende CSV-File.
4. Importiere es in Google Maps, um diese anzuzeigen.

Identifiziere (möglichst) alle Städte mit $10'000$ Einwohnern oder mehr und zeige sie auf einer Google Map an.

Das Problem ist, dass der SBB-Datensatz die Anzahl Einwohner nicht beinhaltet. Deshalb greifen wir auf den Gemeinden-Datensatz vom Thema Datenverarbeitung mit Python zu. Das Vorgehen ist in etwa so:

1. Gehe jede Gemeinde im Gemeinden-Datensatz durch und finde heraus, ob sie mehr als $10'000$ Einwohner hat.
2. Falls ja versuchst du im SBB-Datensatz die Koordinaten dieser Gemeinde zu ermitteln.
3. Achtung: Das geht nicht immer! Warum? Verwende deshalb try-except, damit der Code nicht crashed.
4. Falls du einen passenden Eintrag gefunden hast im SBB-Datensatz, schreibst du einen entsprechenden Eintrag ins CSV.
5. Optional kann man als „Name“ im CSV-File nicht nur den <Gemeindenamen> sondern z.B. <Gemeindenamen (Einwohner)> anzeigen. Dann hat man eine Karte, bei der man auf die Gemeinden klicken kann und gleich die Anzahl Einwohner sieht.
   

Erweitere die letzte Aufgabe: Es fehlen noch einige Ortschaften, auch wichtige wie Zürich und Basel. Warum? Versuche den Code so zu erweitern, dass möglichst viele dieser 'problematischer' Ortschaften auch integriert werden.

Zu Abstimmungen kann man jeweils die detaillierten Daten im JSON-Format herunterladen. Hier der Link zu den Daten der eidgenössischen Abstimmungen

Hier direkt zur Abstimmung vom 18.06.2023

1. Lade den Datensatz zur aktuellsten eidgenössischen Abstimmung herunter.
2. Verschaffe dir einen Überblick: Wie ist der Datensatz aufgebaut? Tipps:
   1. Mit .keys() kannst du sehen, welches die Schlüssel eines Dicts sind.
   2. Mit type(...) kannst du ermitteln, um was für einen Typ es sich handelt.
   3. Achtung: Der Datensatz beinhaltet ein Dict, welches Listen enthält, welches wiederum Dicts enthält u.s.w.
   4. Falls du überfordert bist damit, findest du unten die Antwort
3. Finde in diesem Datensatz die Anzahl Ja-Stimmen zu allen Vorlagen in deinem Wohnort.
4. Finde die extremsten Gemeinden zu den jeweiligen Vorgaben: Höchste und tiefste Anzahl Ja- resp. Nein-Stimmen
5. Wähle eine Vorlage aus. Ermittle von jedem Kanton die beiden Gemeinden mit dem höchsten Ja- resp- Nein Anteil. Erstelle mithilfe von WKT eine Karte, die dies darstellt. Falls im Datensatz sbb.json die jeweilige Gemeinde nicht eingetragen ist, musst du diese halt nachschlagen.

Tipps:Aufbau Datensatz inkl. Beispiel

# lade Datensatz (JSON) in ein dict mit Name data
 
vorlagen = data['schweiz']['vorlagen'] # vorlagen ist Liste, Länge entspricht Anzahl Vorlagen
vorlage = vorlagen[0] # vorlage beinhaltet alle Infos zur ersten Vorlage (erstes Element der Liste vorlagen
print(vorlage['vorlagenTitel'][0]['text']) #  Name der Vorlage (0 für Deutsch)
 
kantone = vorlage['kantone'] # Liste mit Infos zu allen Kantonen (zur gewählten Vorlage)
kanton = kantone[3] # Infos zu einem Kanton (Zahl 0 - 25, weil 26 Kantone)
print(kanton['geoLevelname']) # Name dieses Kantons
 
gemeinden = kanton['gemeinden'] # Liste mit Infos zu allen Gemeinden in diesem Kanton (zur gewählten Vorlage)
gemeinde = gemeinden[7] # Alle Infos zu einer Gemeinde
 
resultat = gemeinde['resultat'] # Resultat
jaStimmenInProzent = resultat['jaStimmenInProzent']
"""