Je sais, ma question est similaire à certaines anciennes sur ce site. J'ai beaucoup de fichiers CSV (coordonnées géographiques) à importer dans qgis (puis à les convertir), et la manière habituelle n'est pas la meilleure façon de le faire (trop longue). J'ai près de 500 fichiers CSV (coordonnées...
J'essaie d'apprendre les ficelles du traitement d'image de télédétection en utilisant les liaisons Python GDAL et numpy. Dans un premier temps, je lis un fichier géotiff Landsat8, je fais une manipulation simple et j'écris le résultat dans un nouveau fichier. Le code ci-dessous semble fonctionner...
J'essaie de générer des polygones pour des andains orbitaux satellites. Jusqu'à présent, j'ai une méthode pour générer deux lignes qui représentent le bord de chaque bande en [lat, long]. Certains des andains traversent la ligne de données internationale et s'enroulent donc: J'ai pu résoudre ce...
J'essaie d'utiliser ArcObject en Python. J'ai un problème avec l'attribution d'une interface (ie IZ.InterpolateZsBetween) à mon objet avec des données de géométrie, le code est: from Comtypes.client import CreateObject, GetModule import arcgisscripting gp = arcgisscripting.create(9.3)...
Je souhaite convertir une table PostGIS en un fichier de formes (sans utiliser pgsql2shp). Afin de créer une géométrie dans le fichier de formes, je dois donner les Xmin, Ymin et Xmax, Ymax, et la géométrie que j'ai dans ma table PostGIS est de forme irrégulière (je peux obtenir l'extérieur en...
Est-il possible d'arrêter et de démarrer par programme les services de carte ArcGIS Server à l'aide de
Est-il possible d'intercepter les messages d'erreur de sortie PostgreSQL dans OGR après l'échec d'une requête? Par exemple conn.ExecuteSQL('SELECT * FROM non_existing_table;') Dans Postgres, cela renvoie un message d'erreur: "Error: relation 'non_existing_table' does not exist." Dans l'ogr, nous...
Lors de l'itération sur une couche vectorielle à l'aide du code suivant (résumé de l'exemple dans le livre de recettes pyqgis ), existe-t-il un moyen de vérifier si une fonction est sélectionnée? provider = vlayer.dataProvider() feat = QgsFeature() allAttrs = provider.attributeIndexes()...
J'ai un problème / une question concernant python et QGIS. Est-il possible d'obtenir les données d'une fonctionnalité sélectionnée dans QGIS en utilisant python, et si oui, comment? Dans mon exemple, j'ai une couche qui montre quelques parcelles pour l'urbanisme. Si je sélectionne maintenant l'un...
J'essaie de faire une figure montrant des données azimutales avec une gamme différente d'incertitudes à chaque point. Cette figure oldschool d'un papier de 1991 capture l'idée de "noeud papillon" que je vise: Des suggestions sur la façon dont je pourrais faire une figure similaire? Je suis...
Est-il possible de définir le projet CRS d'un projet dans QGIS via une commande Python? Dans mon cas, l'utilisateur doit ouvrir QGIS, démarrer mon plugin écrit et ce plugin chargera les données spatiales qui ont le code epsg 31467 (Gauss Krueger Zone 3). Malheureusement, le CRS par défaut de QGIS...
J'ai un outil de script pratique qui parcourt un espace de travail et renomme et copie les fichiers de formes dans un jeu de données d'entité. Cependant, s'il existe un fichier de formes corrompu quelque part dans l'espace de travail, le script échoue et arrête le traitement. Comment gérez-vous des...
J'essaie de comprendre comment je peux utiliser le module csv en python pour ouvrir un fichier csv dans le même dossier que le script python, puis créer un fichier de formes à l'aide du module shapefile pyshp. Le fichier csv ressemble à ceci, mais peut avoir quelques milliers de lignes...
J'ai un grand nombre de fichiers de formes représentant des domaines d'intérêt pour une analyse qui sera réalisée à l'aide de différentes sources d'imagerie satellite (IKONOS, RapidEye, etc.). Malheureusement, l'imagerie n'utilise pas de système pathrow comme Landsat par exemple, de sorte que les...
J'essaie d'apprendre à utiliser l'ogr en python en utilisant les ensembles de données sur les pays et les lieux peuplés à partir de http://www.naturalearthdata.com/downloads/50m-cultural-vectors/. J'essaie d'utiliser des filtres et des tampons pour trouver des points (ne_50m_populated_places.shp)...
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...
J'utilise l'algorithme Birch du package scipy-learn Python pour regrouper un ensemble de points dans une petite ville en ensembles de 10. J'utilise le code suivant: no = len(list_of_points)/10 brc = Birch(branching_factor=50, n_clusters=no, threshold=0.05,compute_labels=True) Dans mon idée, je me...
Je voudrais créer un script Python dans QGIS qui prendra certaines données (shp + tif, dans mon exemple), utiliser Tamplate Map Composer (à partir d'un fichier) et exporter la composition créée vers une image png. Avec presque aucune expérience en programmation (moins que la connaissance de base de...