Erstellen Sie eine Webmap mit Python und GDAL

Erfahren Sie, wie Sie GIS-Daten mit Open Source Python und GDAL kostenlos verarbeiten. Bei den diesmal verwendeten Daten handelt es sich um ein numerisches Höhenmodell von grundlegenden Karteninformationen , das vom National Land Research Institute bereitgestellt wird.

Die obige Webmap wird mit einer Kombination aus Python und GDAL erstellt.

Installation

Verwenden Sie zunächst Miniconda , um Python und die erforderlichen Bibliotheken zu installieren. Öffnen Sie nach Abschluss des Downloads die Miniconda-Befehlszeile und starten Sie die virtuelle Python-Umgebung (mit der Standardversion von Python 3):

conda create --name myenv
conda activate myenv

Installieren Sie dann die für Python-Skripte verwendeten Bibliotheken:

conda install -c conda-forge gdal numpy beautifulsoup4

Laden Sie abschließend das Python-Skript herunter, mit dem die Webmap erstellt wurde:

git clone https://github.com/danielhoshizaki/hillshade

Wie benutzt man

Stellen Sie zunächst sicher, dass Sie eine numerische Höhenmodelldatei im Verzeichnis `. / Data / raw``` haben. Wenn Sie andere Dateien heruntergeladen haben, legen Sie sie zu diesem Zeitpunkt in `. / Data / raw``` ab.

Sie müssen einen Teil des Codes ändern, bevor Sie das Python-Skript drehen. Sie müssen den Pfad zur Bibliothek GDAL angeben, die mit `` `conda``` installiert wurde. GDAL ist eine kooperative Software, aber die Installation ist sehr kompliziert. Verwenden wir diesmal den absoluten Pfad. Um die GDAL-Binärdatei zu finden, verwenden Sie das Linux-Suchfeld oder das Windos-Suchfeld. Wenn Sie nach dem oben erstellten Verzeichnis der virtuellen Umgebung von myenv suchen, wird es sofort angezeigt. Suchen Sie im Verzeichnis nach gdaldem.exe und gdalbuildvrt.exe und im Verzeichnis nach gdal2tiles.py. Sobald Sie es gefunden haben, legen Sie den gdal_bin_path und den gdal_tiles_path in Ihrem Python-Skript fest.

Wenn Sie fertig sind, starten Sie das Skript und warten Sie auf die Webmap.

Wie es funktioniert

Python-Skripte können grob in zwei Rollen unterteilt werden. Eine besteht darin, eine rohe numerische Höhenmodelldatei (ZIPed XML-Datei) mit einer Funktion namens `` `convert``` in GeoTiff zu konvertieren. Die zweite Rolle besteht darin, dass GDAL die Daten direkt über die Befehlszeile verarbeitet. Es gibt drei Befehle, die direkt an GDAL gesendet werden können: Die folgende Datenverarbeitung wird ausgeführt:

1. Konvertieren Sie GeoTiff jedes DEM in Shading GeoTiff
2. Kombinieren Sie alle schattierenden GeoTiffs zu einem virtuellen Datensatz vrt
3. Konvertieren Sie das virtuelle Dataset in eine Webmap

Der letzte Befehl gdal2tiles erstellt ein neues Verzeichnis `` `. / Data / WTMS```. Im Verzeichnis befindet sich eine Datei mit dem Namen Leaflet.html. Doppelklicken Sie also darauf. Ihre eigene Webkarte ist fertig! Diese Webkarte kann jedoch nur auf Ihrem eigenen PC angezeigt werden. Wenn Sie möchten, dass andere sie sehen, stellen Sie sie auf Github Pages oder AWS S3 oder geben Sie Ihr Bestes, um sie im Web-Vorort zu veröffentlichen!