Einzugsgebietsanalyse mit dem Openrouteservice (QGIS3)

Die Einzugsgebietsanalyse ist hilfreich, um die Erreichbarkeit von Standorten auszuwerten. Mit einem gegeben Standort von Feuerwehren, Krankenhäusern, Stationen des öffentlichen Personenverkehrs usw. kann mit Hilfe der Analyse bestimmt werden, welche Gebiete in Abhängigkeit von der aufzuwendenden Zeit oder Strecke bedient werden können. Bis vor kurzem waren solche Analysen mit Hilfe von Open-Source-Tools und freien Daten schwierig. Doch jetzt haben wir durch OpenStreetMap (OSM) Zugriff auf ein globales Straßennetzwerk und freie Webdienste wie Openrouteservice (ORS), der komplexe Routingaufgaben unter Nutzung von OpenStreetMap (OSM) Daten ausführen kann. In diesem Tutorial nutzen wir das ORS Tools Plugin, um eine Einzugsgebietsanalyse mit QGIS durchzuführen.

Überblick über die Aufgabe

Wir nutzen die U-Bahn Stationen der Stadt Kochi (Indien), um das Gebiet zu bestimmen, von dem aus eine Station in 15 Minuten zu Fuß erreicht werden kann.

Weitere Fähigkeiten die wir erlernen

  • Wie man Daten im digitalen Austauschformat General Transit Feed Specification (GTFS) in QGIS lädt.

  • Wie man eine Reihe von Punktdaten in Linien mit Hilfe des Tools Points to Path konvertiert.

Beschaffung der Daten

Kochi Metro Rail Limited (KMRL) stellt frei verfügbare Daten für die U-Bahn in Kochi im Format Global Transit Feed Specification (GTFS-static) bereit. Die Daten können zum download auf der Seite Open Data abgefragt werden.

Der Einfachheit halber können wir eine Kopie des Datensatzes unter folgendem Link herunterladen:

KMRL-Open-Data.zip

Datenquelle [KMRL]

Einrichtung

Openrouteservice API stellt Routingalgorithmen bereit, die mit freien geographischen Daten wie z.B. OpenStreetMap arbeiten. Es handelt sich um einen freien webbasierten Service, auf den mit einem QGIS Plugin zugegriffen werden kann. Obwohl der Service frei ist, ist eine Registrierung zum Erhalt eines API-Schlüssels erforderlich. Der API-Schlüssel dient dazu, Missbrauch zu verhindern und beschränkt die Nutzung.

  1. Gehe zur Seite Openrouteservice Sign Up und erstelle einen Account. Nachdem der Account aktiviert wurde, besuche die Seite Dashboard und fordere ein Token an. Wähle Free als the Token type und gebe ORS Tools QGIS als Token name ein. Klicke auf CREATE TOKEN.

../../_images/setup12.png
  1. Nach der Erstellung kopieren wir die lange Zeichenkette unter Key. Das ist der eindeutige mit dem Account verbundene Identifikator der zur Autorisierung des Services verwendet wird.

../../_images/setup21.png
  1. Öffne QGIS. Gehe zu Erweiterungen ‣ Erweiterungen verwalten und installieren. Suche nach dem Plugin ORS Tools und installiere es. Klicke auf Schließen.

../../_images/setup3.png
  1. Gehe im QGIS Hauptfenster zu Web ‣ ORS Tools ‣ Provider Settings.

../../_images/setup4.png
  1. Erweitere den Bereich openrouteservice und kopiere den Schlüssel (aus Schritt 2) in das Textfeld API Key. Klicke auf OK.

../../_images/setup5.png

Arbeitsablauf

  1. Unzip the downloaded KMRL-Open-Data.zip file to a folder on your computer. You will notice that the unzipped directory contains many text files. Each file contains data for a different aspect of the transit system. The format of the files and their uses are described in GTFS Reference. Out of all the files, 2 files contain geospatial data and are of interest to us. The file shapes.txt contains points that describe a physical path that the vehicle takes, and the file stops.txt contains the location of each transit stop. Both of these are CSV files that can be imported into QGIS. Click the Open Data Source Manager button.

../../_images/1141.png
  1. Im Dialogfenster Data Source Manager gehen wir zum Feld Getrennte Texte. Klicke auf den Knopf neben Dateiname und öffne die Datei shapes.txt. Wähle CSV (comma separated values) als Dateiformat. Das X-Feld und das Y-Feld sollten automatisch gefüllt sein. Klicke auf Hinzufügen.

../../_images/2121.png
  1. Klicke nocheinmal auf den Knopf und wähle die Datei stops.txt. Klicke auf Hinzufügen und danach auf Schließen.

../../_images/364.png
  1. Wir sehen 2 neue Layer stops und shapes, die im Layer Bereich angezeigt werden. Wir werden jetzt den shapes Punktlayer in einen Linienlayer konvertieren. Er repräsentiert den Verlauf der U-Bahnlinie. Gehe zu Verarbeitung ‣ Werkzeugkiste.

../../_images/430.png
  1. Suche das Tool Vektorerzeugung ‣ Punkte zu Weg und klicke doppelt darauf, um es zu starten.

../../_images/528.png
  1. In the Points to Path dialog, select shapes as the Input point layer. As per GTFS specifications, each individual route has a unique shape_id so select that from the drop-down menu as the Path group expression. We can also specify the order of points that will form the line by selecting shape_pt_sequence as the Order expression. Click Run.

../../_images/627.png
  1. Ein neuer Layer Pfade wird im Bereich Layer hinzugefügt. Wir können die Sichtbarkeit des Layers shapes abschalten, um den neu hinzugefügten Layer zu sehen.

../../_images/727.png
  1. Jetzt wo wir die U-Bahnstationen und die Linie hinzugefügt haben, können wir mit der Netzwerkanalyse beginnen. Unter Verarbeitungswerkzeuge suchen wir das Tool ORS Tools ‣ Isochrones ‣ Isochrones From Layer. Klicke doppelt darauf, um es zu starten.

../../_images/826.png
  1. Wir wählen openrouteservice als Provider. Wir werden das Gebiet in 15-minütiger fußläufiger Reichweite zu jeder U-Bahnstation ermitteln. Wähle stops als Input Point Layer und stop_id als Input Layer ID Field. Im Dropdown-Menü Travel mode wählen wir foot-walking. Da wir an der zeitbasierten Fläche interessiert sind, wähle time als Dimension. Zum Schluss geben wir 15 Minuten als ranges vor. Klicke auf Starte.

../../_images/927.png

Bemerkung

Beachte, dass die Openrouteservice API ein Limit von 20 Anfragen je Minute hat. Wenn ein Layer mehr als 20 Punkte hat, erscheinen unter Umständen Fehlermeldungen wegen Überschreitung des Limits. Man kann das Tool weiter laufen lassen, es setzt die Verarbeitung mit 20 Punkten je Minute fort.

  1. Nach Beendigung der Berechnung sehen wir einen neuen Layer Isochrones, der im Layer Bereich geladen wurde. Zu jedem Punkt wurde ein Polygon erzeugt, dass das Gebiet in 15-minütiger fussläufiger Reichweite darstellt. Um das Ergebnis im Kontext mit den zur Erstellung verwendeten Quelldaten zu sehen, fügen wir eine OpenStreetMap Basiskarte hinzu. Wir scrollen im Browser Bereich nach unten und ziehen XYZ Tiles ‣ OpenStreetMap in den Arbeitsbereich.

../../_images/1033.png
  1. Ein neuer Layer OpenStreetMap wird im Layer Bereich angezeigt. Ziehe den Layer im Layerbereich ganz nach unten. Wir können hinein- und hinauszoomen, um zu sehen, ob das Ergebnis zum Straßennetzwerk passt. Wir sehen, dass die Polygone nicht kreisförmig sind, da die Reisezeit entlang der Straßen berechnet wird. In Regionen mit wenigen Straßen wird eine kleinere Fläche ausgewiesen.

../../_images/1142.png
  1. Um das Einzugsgebiet zu berechnen, müssen wir eine letzte Arbeit ausführen. Wir fügen die einzelnen isochron-Polygone zu einem Polygon zusammen, das die erreichbare Fläche repräsentiert. Suche das Tool Vektorgeometrie ‣ Auflösen.

../../_images/1236.png
  1. Wähle Isochrones als Eingabelayer und klicke auf Starte.

../../_images/1334.png
  1. Nach Fertigstellen der Berechnung wird ein neuer Layer Aufgelöst im Layer Bereich angezeigt. Dieses Polygon repräsentiert die Fläche, die mit einem 15-minütigen Fußweg von der U-Bahn aus erreichbar ist.

../../_images/1431.png

Bemerkung

Das ist ein einfaches Beispiel für eine Einzugsgebietsanalyse eines öffentlichen Verkehrsmittels in QGIS. Eine umfassendere Einzugsgebietsanalyse der U-Bahn würde andere Verkehrsmittel einschließen. Wir könnten Zubringerbusse, nahe gelegene Bushaltestellen und Linien mit Anbindung an die Haltestellen berücksichtigen. Wir könnten auch andere Verkehrsmittel wie das Auto oder Taxi berücksichtigen.


If you want to give feedback or share your experience with this tutorial, please comment below. (requires GitHub account)