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)...
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...
Parfois, mais trop souvent malheureusement, mes services AGS s'arrêtent sans raison. Je ne le remarque généralement pas jusqu'à ce que les utilisateurs m'appellent et me demandent pourquoi leurs cartes n'ont pas de données. Existe-t-il un moyen (je pensais à Python) d'obtenir une réponse d'AGS...
J'utilise la ligne suivante pour informer l'utilisateur de l'état: iface.mainWindow().statusBar().showMessage("Status:" + str(i)) Le plugin prend environ 2 minutes pour s'exécuter sur mon ensemble de données, mais Windows le détecte comme "ne répondant pas" et arrête d'afficher les mises à jour...
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...
Je veux créer une fonction qui sélectionne une fonction et effectue un zoom sur celle-ci (similaire dans QGIS). Il y a donc la fonction suivante: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) def zoomTo(self): layer = self.vlayer atable = self.ui.table selectList=[] for i in...
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...
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 à...
J'essaie de créer un fichier de formes polygonales à partir d'une liste de coordonnées en utilisant des outils open source python. Ce qui suit est ce que j'ai jusqu'à présent qui a été piraté ensemble à partir du livre de recettes Python GDAL / OGR et ce GIS SE répond . Il y a une question...
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...
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...
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 un ensemble de lignes représentant des lignes de bus. Certaines lignes se chevauchent et empruntent les mêmes routes. Je peux extraire les nœuds. Cependant, je suis intéressé à extraire uniquement les croisements comme celui-ci: Comment puis-je faire ceci? Je cherche des moyens avec QGIS ou...
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 #...
Je me demande s'il existe un moyen d'ouvrir la console Python directement au lancement du programme. Il existe une possibilité de définir un raccourci pour la console Python, mais je ne trouve pas une telle option pour le lancement du programme
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'...