Créer une carte Web en utilisant Python et GDAL

Apprenez à traiter gratuitement les données SIG à l'aide de Python open source et GDAL. Les données utilisées cette fois-ci sont un modèle numérique d'élévation des informations cartographiques de base fourni par le National Land Research Institute.

La carte Web ci-dessus est créée avec une combinaison de Python et GDAL.

Installation

Tout d'abord, utilisez Miniconda pour installer Python et les bibliothèques requises. Une fois le téléchargement terminé, ouvrez la ligne de commande Miniconda et lancez l'environnement virtuel Python (en utilisant la version par défaut de Python 3):

conda create --name myenv
conda activate myenv

Ensuite, installez les bibliothèques utilisées pour les scripts Python:

conda install -c conda-forge gdal numpy beautifulsoup4

Enfin, téléchargez le script Python utilisé pour créer la carte web:

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

Comment utiliser

Tout d'abord, assurez-vous d'avoir un fichier de modèle d'élévation numérique dans le répertoire . / Data / raw. Si vous avez téléchargé d'autres fichiers, placez-les dans `` ``. / Data / raw à ce stade.

Vous devez modifier une partie du code avant d'activer le script Python. Vous devez spécifier le chemin d'accès à la bibliothèque GDAL installée avec conda```. GDAL est un logiciel coopératif, mais l'installation est très compliquée, alors utilisons le chemin absolu cette fois. Pour trouver le binaire GDAL, utilisez le Linux `` find '' ou le champ de recherche Windos. Si vous recherchez le répertoire de l'environnement virtuel myenv créé ci-dessus, il sortira immédiatement. Recherchez gdaldem.exe et gdalbuildvrt.exe dans le répertoire et gdal2tiles.py dans le répertoire. Une fois que vous l'avez trouvé, définissez gdal_bin_path et gdal_tiles_path dans votre script Python.

Lorsque vous êtes prêt, lancez le script et attendez la carte Web.

Comment ça fonctionne

Les scripts Python peuvent être globalement divisés en deux rôles. La première consiste à convertir un fichier de modèle d'élévation numérique brut (fichier XML ZIP) en GeoTiff à l'aide d'une fonction appelée `` convertir ''. Le second rôle est de laisser GDAL traiter les données directement via la ligne de commande. Trois commandes peuvent être envoyées directement à GDAL: Le traitement des données suivant est effectué:

1. Convertissez GeoTiff de chaque DEM en GeoTiff d'ombrage
2. Combinez tous les GeoTiff d'ombrage dans un jeu de données virtuel vrt
3. Convertir l'ensemble de données virtuel en webmap

La dernière commande gdal2tiles crée un nouveau répertoire . / Data / WTMS```. Il y a un fichier appelé Leaflet.html dans le répertoire, alors double-cliquez dessus. Votre propre carte Web est terminée! Cependant, cette carte Web ne peut être consultée que sur votre propre PC. Si vous voulez que d'autres la voient, placez-la sur les pages Github ou AWS S3, ou faites de votre mieux pour la mettre sur le Web subur!