Je reprojecte des rasters en python en utilisant GDAL. J'ai besoin de projeter plusieurs tiffs à partir des coordonnées géographiques WGS 84 vers WGS 1984 Web Mercator (Sphère auxiliaire), afin de les utiliser plus tard dans Openlayers avec OpenStreetMap et peut-être Google Maps. J'utilise Python...
Je vois assez souvent dans le code Python GDAL que les gens ferment des ensembles de données à la fin de leur script. Pourquoi est-il judicieux de fermer un ensemble de données dans Python GDAL? Y a-t-il des conséquences si je ne le fais pas? import gdal # open dataset ds = gdal.Open('test.tif') #...
Toutes mes excuses si la question suivante est quelque peu stupide, mais je ne suis que TRÈS novice dans ce domaine SIG. J'essaie de convertir des images geoTiff projetées en WGS84 en utilisant gdal en python. J'ai trouvé un article qui décrit le processus de transformation des points dans les...
Étant donné une geopandas GeoDataFramecontenant une série de polygones, je voudrais obtenir la superficie en km2 de chaque entité de ma liste. C'est un problème assez courant, et la solution habituelle suggérée dans le passé a été d'utiliser shapelyet pyprojdirectement (par exemple ici et ici )....
Fermé . Cette question doit être plus ciblée . Il n'accepte pas actuellement de réponses. Vous souhaitez améliorer cette question? Mettez à jour la question pour qu'elle se concentre sur un seul problème en modifiant ce post . Fermé il y a 9 mois . À l'exception d'ArcPy, existe-t-il une...
Le package de géopopie de Python propose deux techniques de mesure de distance: les formules Great Circle et Vincenty . >>> from geopy.distance import great_circle >>> from geopy.distance import vincenty >>> p1 = (31.8300167,35.0662833) # (lat, lon) -
Je travaille actuellement pour importer des données climatiques CANGRID (fournies sous forme de fichiers Surfer Grid ascii, ".grd") dans ArcGIS. La grille est de 95 lignes par 125 colonnes. Les métadonnées fournissent le lat / lon d'origine (coin inférieur gauche), la taille des cellules (50 km)...
Je travaille avec un complément d'ArcMap en C #. A partir du code C #, j'ai exécuté quelques scripts Python. Maintenant, pour exécuter ces scripts, j'ai un chemin python codé en dur. Mais ce n'est pas portable. Donc, je veux obtenir le chemin de l'exécutable Python à partir du code et l'utiliser....
Je suis vraiment nouveau sur Python et j'aimerais savoir s'il existe une méthode rapide pour obtenir les valeurs de cellule d'un raster pixel par pixel et les coordonnées (mapper les coordonnées XY du centre de chaque pixel) en utilisant Python dans ArcGIS 10? Pour décrire cela plus en détail, je...
Je veux créer par programme un fichier GeoJSON en utilisant des polygones à partir d'un fichier de formes mais en ajoutant des attributs de ma propre application. Cela se fait facilement pour un fichier de formes: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data...
Je souhaite apprendre à utiliser les tableaux NumPy pour optimiser le géotraitement. Une grande partie de mon travail implique des «mégadonnées», où le géotraitement prend souvent des jours pour accomplir certaines tâches. Inutile de dire que je suis très intéressé par l'optimisation de ces...
J'ai un logiciel (pas Arc) qui fonctionne tous les soirs sur une machine hors production qui met à jour les fonctionnalités du fichier de formes à partir d'une base de données externe. Je voudrais périodiquement supprimer toutes les fonctionnalités du fichier de formes (pas le fichier lui-même, qui...
est-il possible de diviser un fichier de formes par fonctionnalité en python? (Le mieux serait une solution où je peux enregistrer temporairement les objets vectoriels résultants dans la mémoire à la place sur le disque). La raison: je veux utiliser la fonction gdal rasterizeLayer avec plusieurs...
J'essaie de capturer tous les polygones non superposés indiqués ci-dessous en utilisant Shapely (les polygones A, B et C donnés). De plus, j'espère pouvoir le faire sans itération, sans test d'intersection, etc. La réponse acceptée à cette question exprime la méthode PostGIS mais il semblerait que...
J'ai commencé à apprendre à manipuler les données LAS en python et je voulais voir comment les autres gèrent les fichiers LAS. Je voudrais lire les points (j'utilise un tableau numpy) et filtrer les classes 1 et 2 (non classifiées et au sol) dans un tableau séparé. J'ai le code suivant mais je...
J'ai une liste de dataframes, tous de la même structure et tous sont des polygones. Comment puis-je les ajouter tous ensemble à l'aide de
Quelqu'un connaît-il une bibliothèque Python pure (pas de dépendances C ou Java) pour effectuer des opérations de géométrie telles que tampon, union, intersections, etc.? Je recherche quelque chose qui est compatible avec Google
Étant donné les points de données avec la longitude, la latitude et une troisième valeur de propriété de ce point. Comment puis-je regrouper des points en groupes (sous-régions géographiques) en fonction de la valeur de la propriété? J'ai recherché par google et compris que ce problème semble être...
Essayer d'importer gdal pour Python 2.7.3 sous Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...
Je recherche une solution différente et plus élégante à un problème de statistique spatiale. Les données brutes consistent en une coordonnée xy pour chaque arbre individuel (c'est-à-dire converties en un fichier point .shp). Bien qu'il ne soit pas utilisé dans cet exemple, chaque arbre possède...