J'ai installé GDAL et vérifié sur la ligne de commande, cela fonctionne et lorsque j'exécute gdal2tiles.py sur la ligne de commande, cela fonctionne mais lorsque j'exécute le même fichier gdal2tiles.py avec python IDLE, j'obtiens «Aucun module nommé» osgeo '" Erreur. Pouvez-vous...
J'ai geopandas GeoDataFrame contenant des centaines de galbe Polygonet MultiPolygongéométries. Les polygones se chevauchent à de nombreux endroits. Je voudrais créer une nouvelle géométrie qui contienne le nombre de chevauchements. Quelque chose comme ça: Quelqu'un at-il des idées sur la façon...
J'ai réussi à ajouter QgsVectorLayerà QGIStravers un plugin python, par exemple: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Mais comment puis-je le
cela ressemble à une question idiote, mais je ne trouve pas de réponse claire à cela: quelles unités géopandas / galbées utilisent le calcul de la distance / surface entre les
j'ai GeoPandas 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 installé dans l'environnement Anaconda 2-4.1.1-Windows-x86. Je suis capable de construire GeoDataFrameen lisant des jeux de données d'entrée et d'exploiter les données, mais l'enregistrement d'un jeu de données de sortie ne préserve...
Existe-t-il un moyen d'exécuter un plug-in personnalisé conçu pour QGIS qu'un shell Python puisse apparaître lors de l'exécution du plug-in, comme un écran de débogage? Je demande ceci pour savoir s'il existe un moyen de voir les commandes d'impression via un script Python lors de la programmation...
Je sers des tuiles vectorielles en utilisant TileStache , j'ai tout configuré comme je veux. Mes données sont stockées dans Postgres et j'utilise le fournisseur VecTiles pour servir les tuiles GeoJSON . Je veux mettre en cache toutes mes tuiles pour rendre les tuiles plus rapides. J'utilise...
Voir ce lien pour plus de détails. Le problème: Je veux parcourir un raster continu (un qui n'a pas de table attributaire), cellule par cellule, et obtenir la valeur de la cellule. Je veux prendre ces valeurs et exécuter des conditions dessus, en émulant les étapes d'algèbre de carte détaillées...
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
J'ai un raster avec lequel j'aimerais faire quelques interpolations ponctuelles. Voici où j'en suis: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array =...
J'écris un code python pour convertir X, Y points en excel en shapefile. Dans ce processus, j'ai des points xy à partir de shhet1. Quelle est l'étape que je dois inclure dans mon processus pour lire les données de la feuille1 d'un classeur Excel (97-2003). J'ai écrit un code comme suit .. import...
Il y a quelque temps, j'ai fait une enquête d'orientation du toit du quartier pour le déploiement des panneaux solaires, en vérifiant visuellement chaque toit et en marquant individuellement s'il était orienté au sud ou presque au sud. Comme cela prend beaucoup de temps, j'aimerais trouver un moyen...
J'ai appris le géotraitement avec Python à l'aide du cours GEOG 825 (géotraitement python) de MSU . Malheureusement, le lien est mort maintenant. Je sais que c'est un plan à long terme, mais quelqu'un a-t-il une copie des documents de ce cours? J'aimerais vraiment le terminer. Vous avez déjà...
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"...
J'aimerais construire un fichier de formes à partir d'un cadre de données Pandas en utilisant les lignes lon & lat. J'ai un fichier csv et je le traite avec des pandas pour créer un bloc de données plus facile à manipuler Est-il possible de le faire sans faire de boucle ligne par...
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 =...
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 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...