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 extrait de code pour comprendre ce que je veux dire: from shapely.geometry …
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 GDAL est installé …
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 Python pure que je pourrais …
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 = layer.dataProvider() it = dp.getFeatures() for i in …
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" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = …
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 données de coordonnées pour créer …
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 être en dehors …
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 …
[EDIT: la solution à cela était simplement d'utiliser OGR pour lire les fichiers de formes. Voir l'exemple de geographika.] Dans un fichier de formes ESRI, il n'y a pas de distinction entre Polygones et MultiPolygons. De plus, il n'y a pas de distinction explicite entre les trous intérieurs et les …
j'ai un shapefile de subdivisions de pays et je voudrais extraire un tableau de délimitation lat et lon pour chaque division .. est-il possible de faire?
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 examiner "cela"?
en utilisant Gdal / OGR, existe-t-il un moyen de convertir une couche de ligne shp, kml ou PostGIS en un réseau de nœuds et de segments afin qu'elle puisse être utilisée dans des packages tels que networkx?
Je voudrais savoir comment obtenir toutes les valeurs raster dans un polygone en utilisant gdal ou pygeoprocessing, sans lire la grille entière comme un tableau. pygeoprocessing et gdal peuvent faire des statistiques zonales mais seuls les min, max, moyenne, stdev ou count sont disponibles à partir d'une telle fonction. Étant …
J'essaie d'ouvrir un fichier de formes avec le bloc-notes Jupyter, mais je n'arrive pas à l'ouvrir. Le fichier de formes se trouve dans le même dossier que mon classeur et Python ne peut pas l'identifier même si j'utilise le chemin complet. import shapefile as shp import matplotlib.pyplot as plt sf …
J'ai un projet QGIS qui contient un compositeur. Bien sûr, si je vais dans l'interface graphique, je peux l'exporter sous forme de fichier pdf, etc. De la lecture en ligne, voici ma tentative: #!/usr/bin/env python import sys from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from PyQt4.QtCore import QFileInfo from PyQt4.QtXml import QDomDocument …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.