J'ai une application django avec un dépliant sur le frontend et j'ai besoin de calculer la bbox d'un objet geojson pour pouvoir transmettre les coordonnées à map.fitBounds() J'ai essayé quelques bibliothèques ainsi que les getBoundsméthodes de certains objets dans Leaflet, comme FeatureGroup, mais...
J'ai visité et revisité la page sur la gestion de la géométrie dans le livre de recettes PyQGIS: http://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html mais je n'arrive pas à comprendre comment obtenir le polygone à dessiner à partir de la console Python. Quelqu'un...
J'ai un gros fichier raster (129600 par 64800 pixels) avec les plans d'eau mondiaux (valeurs 1 bit 0 et 1) et j'essaie d'extraire les rives des océans et des eaux intérieures. J'ai essayé avec ArcGIS et QGIS de convertir de raster en polyligne, mais cela prend des siècles. Quelqu'un connaît-il un...
J'ai un raster entier pour lequel je voudrais construire une table d'attributs raster en utilisant Python et GDAL. Je peux créer une table d'attributs raster GDAL en Python comme suit: >>> rat = gdal.RasterAttributeTable() Cela fonctionne bien, comme nous pouvons le voir: >>> rat...
J'essaie d'optimiser mon code python en utilisant pypy . Bien que je puisse importer la plupart des bibliothèques et l'exécuter avec succès, je ne suis pas en mesure d'importer des bibliothèques GDAL. Les mêmes bibliothèques lors de l'importation en Python 2.7 fonctionnent correctement. J'ai essayé...
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...
Le script que j'écris repositionne deux blocs de données et définit leur étendue. Ce faisant, il redessine l'intégralité d'Active View 4 fois, ce qui ralentit considérablement le script. En appuyant sur F9 ou en cliquant sur le bouton «Suspendre le dessin» avant d'exécuter le script, vous...
Dans la documentation de gdal_calc, il est indiqué Calculatrice raster en ligne de commande avec syntaxe numpy . Plus tard, il y a quelques exemples où dans l'un d'eux: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - signifie que les valeurs définies de...
J'ai rencontré un problème étrange lors de l'exécution de scripts de traitement PyQGIS autonomes. L'ordre des importations dans le script affecte son exécution normale. Vous pouvez reproduire le problème en ouvrant une console Python et en entrant le script suivant (j'utilise GNU / Linux, QGIS...
Je souhaite développer une interface graphique en Python pour le géotraitement ArcGIS. Puis-je utiliser PyQT pour la programmation GUI dans ArcGIS? J'ai également envisagé d'utiliser Tkinter et wxPython.
J'utilise le code ci-dessous pour trouver un pays (et parfois un état) pour des millions de points GPS. Le code prend actuellement environ une seconde par point, ce qui est incroyablement lent. Le fichier de formes fait 6 Mo. J'ai lu que les géopandas utilisent rtrees pour les jointures spatiales,...
J'essaie de changer un plugin de QGIS 1.8 à 2.0. Dans la version 1.8, je crée une table des noms de champs comme ceci: layer = vlayer provider = layer.dataProvider() feat = QGSFeatrue() allAttrs = provider.attributeIndexes() provider.select(allAttrs) fieldNames = [] while provider.nextFeature(feat)...
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...
Je voudrais appeler la fonction du plugin d'interpolation QGIS (méthode TIN) (Raster-> Interpoler) à partir de la console python. Je ne trouve pas la fonction correspondante dans l'API QGIS ou dans la liste des algorithmes de traitement. J'ai trouvé l'algorithme de triangulation SAGA, qui...
Fermé . Cette question doit être plus ciblée . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle se concentre sur un problème uniquement en modifiant ce message . Fermé il y a 4 ans . Je recherche une bonne ressource (MOOC,...
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...
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...
Je suis nouveau dans ce domaine. Je "construis" un plugin python pour qgis et je veux importer des modules depuis GRASS et SAGA mais je ne comprends pas comment le faire. Le plugin que je veux importer est l'interpolation Kriging de SAGA ou v.surf.idw de GRASS ou le plugin d'interpolation de Raster...
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