Débutant aux prises avec le SIG. J'essaie de cartographier les quartiers de la ville de Milwuakee en utilisant des fichiers de formes trouvés sur le site Web du comté . Je suis le fil ici avec un certain succès. Mon code est donne: from pyproj import Proj, transform # wisconsing EPSG:32054 #...
Est-il possible de changer la couleur lors de l'ajout d'un calque vectoriel? J'ai un calque vectoriel avec des zones d'eau et je veux le charger avec un style de couleur bleue. La couche est ajoutée avec: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) Puis-je en quelque sorte changer la...
Existe-t-il un moyen d'obtenir une liste des formats gdal installés sur un système avec python? Sur la ligne de commande, cela peut être fait facilement avec: gdalinfo --formats ... mais cela peut-il être fait en
J'ai trouvé des données d'élévation de haute précision (résolution 1/3 - 1/9 seconde d'arc) à partir du jeu de données d'élévation national fourni par l'USGS. Il est disponible aux formats IMG, Arcgrid et GridFloat. Je sais qu'ArcGrid correspond à un progiciel payant, mais j'essaie de m'en tenir à...
Je viens de commencer à utiliser pycharm avec qgis mais je n'ai pas pu les connecter tous les deux. Pycharm reste toujours dans un état "d'attente de connexion". La plupart des tutoriels disponibles pointent vers les fenêtres, mais j'utilise ubuntu et je n'ai donc pas pu trouver de moyen de...
Je veux faire une simple fonction if / else à l'intérieur de mon modèle et je voulais savoir si c'est possible ou non. Je veux que cette chose fasse que si une entrée pour 'Segment Layer' est sélectionnée, elle n'a pas à exécuter la 'Segmentation (meanshift)' et utilise à la place 'Segment Layer'...
Je souhaite installer QGIS. Je suis sur un système qui utilise par défaut Python 2.5x et j'ai téléchargé python 3.7x avec homebrew. J'ai créé un environnement virtuel avec ce code. brew install pyenv brew install pyenv-virtualenv exec $SHELL eval "$(pyenv init -)" eval "$(pyenv virtualenv-init -)"...
J'essaie de créer un plugin pour charger un compositeur d'impression à partir d'un fichier, générer un atlas et exporter vers l'image. Jusqu'à présent, j'ai réussi à charger le modèle et à l'exporter vers l'image. Je n'ai pas pu ajouter aucune des couches de la légende (qui se trouvent également...
J'ai été chargé de créer un modèle à l'aide d'ArcGIS ModelBuilder, qui exécute initialement quelques tâches de géotraitement, avant de joindre ensuite un tableau Excel au résultat de ces géotraitements. Cependant, le modèle doit réellement ouvrir le fichier Excel en premier, pour que l'utilisateur...
Je voudrais exécuter un script python via qgis --code myscript.pypuis quitter immédiatement. J'utilise iface.actionExit().trigger(), ce qui tue QGIS lorsque je l'exécute à partir de la console python, mais pas lorsqu'il est placé dans le script auquel je passe --code. Quelle est la bonne façon de...
Je recherche un outil du GDAL capable de reclasser un raster contenant plusieurs valeurs discrètes. Jusqu'à présent, j'ai trouvé des rasters Reclassify utilisant GDAL et Python ; cependant, il semble que seules des valeurs uniques puissent être traitées ici. Existe-t-il quelque chose comme un outil...
Existe-t-il un moyen de vérifier si 2 couches raster données ont un contenu identique ? Nous avons un problème sur le volume de stockage partagé de notre entreprise: il est maintenant si important qu'il faut plus de 3 jours pour effectuer une sauvegarde complète. Une enquête préliminaire révèle que...
Je souhaite géoréférencer un raster à l'aide de pythonet GDAL. Mon approche actuelle est d'appeler gdal_translateet d' gdalwarputiliser os.systemet une vilaine liste de points de contrôle au sol. J'aimerais vraiment un moyen de le faire nativement à l'intérieur python. Voici le processus que...
Comment puis-je obtenir les coordonnées projetées ainsi que les valeurs réelles des pixels à ces coordonnées à partir d'un fichier GeoTiff, puis les enregistrer dans un tableau numpy? J'ai le fichier arsenci020l.tif, et ses coordonnées sont en mètres. Voici la sortie abrégée de gdalinfo que j'ai...
De loin, accéder à ArcObjects à partir de Python? est mon Q&A le plus lu et référencé sur GIS Stack Exchange. Malgré ce succès, c'est probablement l'un de mes points faibles en matière d'utilisation réelle. Une grande partie de cette mauvaise performance provient de ma faible capacité à lire et...
Il y a 3 scénarios probables que j'essaie de capturer à proximité des distances: Une station de métro d'échange, qui a 2 stations voisines ou plus. C'est-à-dire que la station en question relie 2 routes principales ou plus et a 2 stations voisines ou plus. Une station de métro terminale, qui n'a...
ArcMap a la possibilité d'afficher une description étendue de chaque élément de la liste des symboles. Comment puis-je remplir cette propriété à partir d'une table? Le but est d'éviter d'avoir à éditer chaque symbole un par un pour ajouter la description .
Quelles sont les règles arithmétiques qui régissent les valeurs de Nodata? Quelle que soit la façon dont ces valeurs nulles sont définies? J'utilise ArcMap Raster Calculator. Je pense spécifiquement aux opérations comme 1 * Nodata (ou 1 * null), 1 / Nodata, 1 + Nodata, 1-Nodata, Nodata-1 etc...
J'ai un fichier de formes avec des énormes multipolygones, avec 100 000 pièces. Quelle serait la façon la plus simple de les diviser en polygones à une seule partie? Je recherche quelque chose comme la fonction QGIS «Multipart to singlepart», mais le fichier est trop volumineux pour que QGIS puisse...
si j'essaie de charger le plugin openlayers dans une application autonome, il échoue toujours avec False, bien que les chemins soient définis et que le plugin soit disponible. (testé sur Ubuntu 14.04, QGIS 2.4, python-2.7) #!/usr/bin/python # -*- coding: utf-8 -*- import sys, os import qgis import...