Lors de l'itération sur une couche vectorielle à l'aide du code suivant (résumé de l'exemple dans le livre de recettes pyqgis ), existe-t-il un moyen de vérifier si une fonction est sélectionnée? provider = vlayer.dataProvider() feat = QgsFeature() allAttrs = provider.attributeIndexes() provider.select(allAttrs) while provider.nextFeature(feat): geom = feat.geometry() attrs = feat.attributeMap() …
J'essaie de faire une figure montrant des données azimutales avec une gamme différente d'incertitudes à chaque point. Cette figure oldschool d'un papier de 1991 capture l'idée de "noeud papillon" que je vise: Des suggestions sur la façon dont je pourrais faire une figure similaire? Je suis relativement novice dans le …
Est-il possible d'intercepter les messages d'erreur de sortie PostgreSQL dans OGR après l'échec d'une requête? Par exemple conn.ExecuteSQL('SELECT * FROM non_existing_table;') Dans Postgres, cela renvoie un message d'erreur: "Error: relation 'non_existing_table' does not exist." Dans l'ogr, nous obtenons la norme: "Error: current transaction is aborted, commands ignored until end of …
Je me demandais si quelqu'un avait une certaine expérience dans l'obtention de données d'altitude à partir d'un raster sans utiliser ArcGIS , mais plutôt obtenir les informations sous forme de python listou dict? J'obtiens mes données XY sous forme de liste de tuples: xy34 =[perp_obj[j].CalcPnts(float(i.dist), orientation) for j in range …
Je veux faire une image multispectrale de cero pour faire des tests dessus. Quelque chose de vraiment simple comme 5 bandes complètement uniformes avec du bruit de sel et de poivre dessus ou un carré de valeurs différentes au centre. De toute évidence, ce ne serait qu'une pile de matrices, …
J'utilise GeoPy pour géocoder des adresses en lat, lng. Je voudrais également extraire les composants d'adresse détaillés (rue, ville, état, code postal) pour chaque adresse. GeoPy renvoie une chaîne avec l'adresse - mais je ne trouve pas de moyen fiable de séparer chaque composant. Par exemple: {street: '123 Main Street', …
Je suis nouveau dans la programmation géospatiale en python. j'ai exécuté le script suivant et obtenu le message d'erreur correspondant >>> import osgeo >>> import osgeo.ogr >>> shapefile = osgeo.ogr.Open("tl_2009_us_state.shp") >>> numLayers = shapefile.GetLayerCount() Traceback (most recent call last): File "<pyshell#5>", line 1, in <module> numLayers = shapefile.GetLayerCount() AttributeError: 'NoneType' …
J'ai une classe d'entités ponctuelles représentant les bureaux de tous les employeurs d'une certaine industrie. La classe d'entités possède un attribut pour stocker le nombre d'employés travaillant dans chaque bureau. Quelqu'un a demandé à utiliser ces données, jointes spatialement à la plus petite unité géographique possible - blocs de recensement, …
Je souhaite convertir une table PostGIS en un fichier de formes (sans utiliser pgsql2shp). Afin de créer une géométrie dans le fichier de formes, je dois donner les Xmin, Ymin et Xmax, Ymax, et la géométrie que j'ai dans ma table PostGIS est de forme irrégulière (je peux obtenir l'extérieur …
J'essaie d'utiliser ArcObject en Python. J'ai un problème avec l'attribution d'une interface (ie IZ.InterpolateZsBetween) à mon objet avec des données de géométrie, le code est: from Comtypes.client import CreateObject, GetModule import arcgisscripting gp = arcgisscripting.create(9.3) gp.setproduct("ArcInfo") esriSystem = GetModule("C:/Program Files (x86)/ArcGIS/com/esriSystem.olb") esriGeometry = GetModule("C:/Program Files (x86)/ArcGIS/com/esriGeometry.olb") esriDataSourcesGDB = GetModule("C:/Program Files …
J'ai ramassé une copie du livre super méchant «Python Geospatial Development» d'Erik Westra ( lien Amazon ), et je travaille dessus. Actuellement, il m'apprend à charger des données côtières GSHHS à partir d'un fichier de formes dans une base de données PostGIS, en préparation de la construction d'une application Web …
J'écris un script python qui lit plusieurs fichiers XML contenant des coordonnées x et y et les combine tous dans un seul fichier csv. La latitude et la longitude sont des champs obligatoires dans le csv, mais j'ai du mal à convertir les coordonnées x, y dans l'Ohio North State …
Je travaille actuellement sur un projet dans lequel j'ai besoin de construire un réseau topologique à partir des fonctionnalités géométriques que je trouve dans les fichiers de formes. Jusqu'à présent, en utilisant le projet open source de Ben Reilly, j'ai réussi à transformer des chaînes de lignes en bords networkx, …
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.