Geodatenabfrage aus Google Maps mittels Python

Während klassische Geodatenbankabfragen recht pingelig sind, wenn es darum geht, einen Ort festzulegen, zeigt Google Maps, wie selbst unscharfe Beschreibungen und Tippfehler zur richtigen Fundstelle führen können.

Netterweise gibt es eine Programmierschnittstelle (API), die die Verwendung dieses Komforts in eigenen Programmen erlaubt. Wer diese Schnittstelle intensiv nutzt, wird lediglich gebeten, einen API-Key zu beantragen, um gegebenenfalls die CPU-Zeit in Rechnung gestellt zu bekommen. Für gelegentliche Nutzung kommt man aber auch ganz gut ohne einen API-Key aus. Man muss nur etwas gegen hin und wieder aufkommende Fehlermeldungen tun, die eine zu intensive Nutzung bemängeln. Eine Sekunde zu warten, reicht oft aus.

Das Python-Modul „geopy“ stellt eine ganze Palette an Funktionen zur Verfügung, um neben der Google-Maps-API auch die Programmierschnittstellen vieler anderer Geodatenanbieter abzufragen. Es gehört leider nicht zum Standardinstallationsumfang von Python 3.6, kann aber schnell mittels „pip install geopy“ oder über die Paketverwaltung besserer Betriebssysteme installiert werden.

Die Klasse geopy.GoogleV3 besitzt beispielsweise die Methode geocode, der man nicht mehr als einen ungefähren Ortsnamen („Banohf Bilefelt“) oder eine Beschreibung („höchster Berg in NRW“) übergeben muss, um die exakte Adresse und die Geokoordinaten zu erhalten. Die Klasse geopy.distance.distance berechnet die Entfernung zweier Punkte und berücksichtigt dabei sogar die Erdkrümmung nach verschiedenen Modellen.

Hier ein kleines Python-3.6-Beispielprogramm, das zu einem Ortsnamen oder einer Beschreibung die genaue Adresse heraussucht und gleich noch die nächste Tankstelle, den nächsten Bioladen und den nächsten Aldi-Markt heraussucht :-)

#!/usr/bin/env python3
 
# Info: https://geopy.readthedocs.io/en/stable/
 
# Das geopy-Modul gehört nicht zum Standardumfang und
# muss nachinstalliert werden, z. B. mittels
# pip install geopy
 
import geopy, geopy.distance
import time
 
# Ohne Schlüssel ist die Zahl der Zugriffe pro Zeitintervall begrenzt.
# https://developers.google.com/maps/documentation/geolocation/get-api-key
apikey = None
 
# Aufruf der Google-API. Bei Fehlern wegen API-Quota-Überschreitung
# wird der Versuch nach einer kleinen Wartezeit wiederholt.
def GoSlow(s):
    antwort = ""
    while antwort == "":
        try:
            antwort = geopy.GoogleV3(apikey).geocode(s)
        except geopy.exc.GeocoderQuotaExceeded:
            time.sleep(1)
    return antwort
 
o1 = None
while not o1:
    o1 = GoSlow(input("Gib einen Ortsnamen oder eine Beschreibung ein: "))
    if not o1:
        print("\nDiesen Ort kennt Google nicht.\n")
 
print("Die Adresse lautet:",o1.address,
      f"\n ({o1.latitude:.6f}, {o1.longitude:.6f})")
 
o2 = GoSlow("Aldi nahe "+o1.address)
d = geopy.distance.distance(o1.point, o2.point)
 
print("\nDer nächste Aldi-Markt ist",round(d.m),
      "Meter entfernt:\n",o2.address,
      f"\n ({o2.latitude:.6f}, {o2.longitude:.6f})")
 
o2 = GoSlow("Bioladen nahe "+o1.address)
d = geopy.distance.distance(o1.point, o2.point)
 
print("\nDer nächste Bioladen ist",round(d.m),
      "Meter entfernt:\n",o2.address,
      f"\n ({o2.latitude:.6f}, {o2.longitude:.6f})")
 
o2 = GoSlow("Tankstelle nahe "+o1.address)
d = geopy.distance.distance(o1.point, o2.point)
 
print("\nDie nächste Tankstelle ist",round(d.m),
      "Meter entfernt:\n",o2.address,
      f"\n ({o2.latitude:.6f}, {o2.longitude:.6f})")