Je voudrais mettre à jour un attribut d'une fonctionnalité. Cependant, je ne sais pas le mettre à jour en utilisant l'objet de fonctionnalité. Je dois utiliser le fournisseur de données pour le mettre à jour. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp =...
Quels sont les pilotes du package python fiona? Quand je vérifie le manuel, il dit [...] and the possible formats are enumerated in the fiona.drivers list. cependant, quand je tape en python from fiona import drivers print drivers > <function drivers at 0x108763050> Comment puis-je...
Je suis intéressé à apprendre des méthodes pour utiliser toute la puissance de traitement multicœur disponible sur un ordinateur de bureau. Arc indique que le géotraitement en arrière-plan permet à l'utilisateur d'utiliser plusieurs cœurs, cependant, les tâches doivent essentiellement attendre en...
Je travaille avec GDAL en Python et j'ai voulu que GDAL lève des exceptions pour que je puisse faire quelque chose d'utile avec elles. Par exemple, si un fichier ne peut pas être ouvert avec gdal.Open (), je voudrais lever une exception au lieu de gdal renvoyant None. Je voudrais juste un moyen de...
Existe-t-il un package Python qui fournit une implémentation d'un algorithme Straight Skeleton? Je suis conscient que le projet open source (C ++) CGAL contient une implémentation mais il semble que cgal-bindings n'inclut pas ce package CGAL . Dans tous les cas, je préférerais une implémentation...
J'essaie de graver un fichier de formes dans un raster en utilisant RasterizeLayer de GDAL. Je pré-crée une raster de zone d'intérêt à partir d'un fichier de formes différent, étant donné une taille de pixel spécifique. Cet AOI sert ensuite de base à toutes les rastérisations suivantes (même nombre...
J'ai installé des bibliothèques python en utilisant pip (pandas, scipy, simplekml) et depuis chaque fois que j'essaie de charger QGIS, j'obtiens cette erreur: Warning: loading of qgis translation failed [/usr/share/qgis/i18n//qgis_en_US] Warning: loading of qt translation failed...
Je rencontre actuellement une erreur «Accès refusé» après avoir utilisé rasterio pour ouvrir les fichiers d'image Sentinel-2 AWS .jp2 (actuellement hébergé par Sinergise et peut être trouvé ici ). Mon code ressemble à ceci: import rasterio access_key = '*****************' secret_access_key =...
Lorsque j'essaie d'installer Shapely sur mon ordinateur Windows 64 bits, je ne parviens pas à faire fonctionner la bibliothèque GEOS. Jusqu'à présent, j'ai exécuté le programme d'installation OSGeo4W à partir duquel j'ai installé GDAL (je pense que la bibliothèque geos est incluse dans ce package)....
J'ai une grille binaire Arc / Info --- spécifiquement, un raster d'accumulation de flux ArcGIS --- et je voudrais identifier toutes les cellules ayant une valeur spécifique (ou dans une plage de valeurs). En fin de compte, je voudrais un fichier de formes de points représentant ces cellules. Je...
Je voudrais obtenir l'intersection de plusieurs polygones. En utilisant le shapelypackage de Python , je peux trouver l'intersection de deux polygones en utilisant la intersectionfonction. Existe-t-il une fonction efficace similaire pour obtenir l'intersection de plusieurs polygones? Voici un...
je cherche une méthode pour calculer la coordonnée X, Y d'un point intérieur (également connu sous le nom labelPoint dans arcmap) d'un polygone (irrégulier). J'essaie de trouver une solution sans utiliser le module ArcGIS la différence entre un "centroïde" et un labelPoint est que le centroïde peut...
Existe-t-il un moyen d'effectuer la même tâche que l'utilitaire gdalbuildvrt à l'aide des liaisons GDAL Python? Jusqu'à présent, je n'ai trouvé aucun autre moyen de le faire que de créer un vrt d'un seul ensemble de données et de modifier manuellement le xml. Je voudrais créer un vrt à partir de...
Quelle est la meilleure méthode pour calculer la distance de la route la plus proche pour une série de points d'adresse dans QGIS, ayant un très grand ensemble de données? (QGIS 2.0.1, Win XP) J'ai lu les réponses aux questions similaires précédentes: Mesurer la distance entre les lignes et les...
J'ai une couche tampon (polygone vert) que je souhaite diviser en deux polygones chaque fois qu'elle franchit une barrière (ligne bleue). J'ai essayé d'utiliser la méthode "splitGeometry", mais je n'arrive pas à la faire fonctionner. Jusqu'à présent, mon code est le suivant: while...
Je lis un fichier de formes avec Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] cela me donne un dict contenant une liste de listes de coordonnées, et tapez "MultiPolygon" Cependant, je ne sais pas comment utiliser les...
J'ai du mal à utiliser la fonction Polygonize en python. L'exemple de livre de cuisine pour cela peut être trouvé ici . La partie pertinente de mon code est: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized"...
Je suis vraiment étonné des énormes mesures prises par le développement de QGIS au cours des derniers mois. Le nouvel ensemble d'icônes, le compositeur d'impression amélioré et bien d'autres changements à voir dans QGIS 2.0 sont vraiment géniaux! Parmi de nombreuses autres améliorations...
Mon manager m'a demandé de proposer une sorte de mesure de compétence pour permettre aux employés de comprendre les scripts / la programmation. Parce que je suis un programmeur autodidacte, je n'ai vraiment aucune compréhension formelle de la façon de mesurer son niveau de compréhension en matière...
Il y a eu des tonnes de questions à ce sujet, mais aucune n'a résolu mon problème. J'ai une machine avec: Windows 7 x64 Python 3.4.3 x64 gdal installé: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal a été installé à partir d'une roue GDAL-1.11.4-cp34-none-win_amd64.whld' ici Le...