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...
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...
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...
Notre organisation envisage de déplacer notre workflow de géotraitement vers PostGIS. Nous utilisons actuellement ArcGIS, avec une pléthore d'outils Python personnalisés utilisés dans ModelBuilder. Nous déplaçons la plupart de nos données dans PostGIS pour qu'elles soient utilisées par une variété...
Je fais un workflow à petite échelle dans un cahier ipython avec des géopandas et bien fait, en tirant un tas de données géospatiales, parfois à partir de fichiers de formes, parfois à partir de postgis (où un traitement plus coûteux est effectué). Maintenant, je tire ensuite les tables postgis en...
J'ai des points en WGS84 lat / long et j'aimerais mesurer de "petites" distances (moins de 5 km) entre eux. Je peux utiliser la formule haversine de http://www.movable-type.co.uk/scripts/latlong.html et cela fonctionne très bien. J'aimerais utiliser des bibliothèques Python Shapely, pour que je...
J'essaie de créer un rectangle tourné en coordonnées lat / lon. Fondamentalement, j'ai les coordonnées du centre et l'angle de rotation autour du centre. Est-il possible de faire en utilisant, galbé ou gdal, ou autre
Fermé. Cette question est hors sujet . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle soit sur le sujet pour l'échange de piles de systèmes d'information géographique. Fermé il y a 3 ans . Je suis assez nouveau sur python et...
Je crée des rasters à virgule flottante dans GDAL avec le pilote GeoTIFF. Lorsque je charge les images résultantes dans QGIS ou Arc, la symbologie par défaut met le min à -3.40282e38 et le max à 3.40282e38, de sorte que le raster semble exactement gris. Existe-t-il un moyen d'écrire directement la...
J'utilise un Ubuntu 14.4 et j'ai installé un environnement virtuel python. Ensuite, j'ai installé plusieurs bibliothèques python (gdal, lxml, psycopg2, etc.) Aussi dans mon système, j'ai le module pour QGIS. Donc, quand je tape python, je peux importer ce module en tant que: import qgis.core Mais...
J'écris une méthode Python qui ajoute deux rasters et génère une seule sortie raster. Pour des raisons indépendantes de ma volonté, les étendues des rasters en entrée sont différentes, mais elles se chevauchent. Est-il possible d'opérer exclusivement sur les pixels du raster en entrée qui se...
Nous avons un grand référentiel (10000 et plus) de données spatiales (vecteurs shp / gdb / tab et rasters ecw / tif et aussi des données las, xyz et DEM) où certains ont des enregistrements de métadonnées en xml en utilisant l'éditeur de métadonnées anzlic. Je veux vérifier si chaque jeu de données...
J'essaie d'apprendre les ficelles du traitement d'image de télédétection en utilisant les liaisons Python GDAL et numpy. Dans un premier temps, je lis un fichier géotiff Landsat8, je fais une manipulation simple et j'écris le résultat dans un nouveau fichier. Le code ci-dessous semble fonctionner...
Il y a une géodatabase fichier sur un site ftp que je voudrais télécharger avec un script Python. En ce moment, je pense qu'une façon de le faire est de copier la géodatabase ftp dans une géodatabase sur mon ordinateur. Voici le script que j'ai commencé. Est-ce que quelqu'un sait comment je peux...
J'implémente un outil qui permet à l'utilisateur de créer le profil d'une entité (dans mon cas parcelles urbaines). Est-il possible de prendre automatiquement une capture d'écran de cette fonctionnalité (affichée dans une fenêtre QGIS) et de l'afficher dans le profil PDF? L'image ne doit montrer...
Étant donné une boîte englobante, comment puis-je compter le nombre de valeurs spécifiques (disons, je suis intéressé par le nombre de value == 1) dans un raster dans 1) arcpy, ou 2) un package raster dans
Avec le code suivant, je peux définir 100% de transparence sur la valeur raster 0: map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0)
Une fois le calque créé, comment le cacher / l'afficher? Je peux activer / désactiver le rendu d'une couche spécifique en cochant la case via QGIS, mais je dois le faire par programmation à partir du code python. Comment puis-je afficher / masquer (pas supprimer) une étiquette du code python? Je...
Quelle est la bonne façon d'insérer un point dans une base de données PostGIS à l'aide de
J'ai trouvé dans l'aide d'ESRI .NET: Python Toutes les applications ArcGIS Desktop incluent un langage de script intégré appelé Python. De nombreux domaines d'ArcGIS, en particulier le géotraitement, sont accessibles via des interfaces de programmation d'application (API) Python simplifiées, ce qui...