J'ai hérité d'un portail Web assez grand / complexe (interne) (exécuté sous python 2.6.6), qui a des capacités de cartographie, en utilisant Mapnik 2.2.0. Le service Web a récemment commencé à générer 500 erreurs lorsqu'il essaie de créer certaines des cartes affichées, avec la seule erreur: $...
J'utilise TileStache pour servir mes propres tuiles vectorielles à l'aide d'une source de données postgis. J'ai d'abord re-projeté les fichiers de formes que j'avais en utilisant ArcMap (reprojeté à Web Mercator (sphère auxiliaire)) et les ai chargés dans Postgres en utilisant shp2pgsql (SRID:...
J'ai une couche avec des polygones (appelons-la "bois") qui est coupée par "n" objets linéaires ("routes"). Chaque route a un attribut roadid unique. Je dois ajouter TOUS les roadid (pas seulement les premières trouvailles de QGis) à la couche de bois-polygones qu'ils croisent pour cette dernière...
J'ai un comportement étrange dans QGIS avec les outils GRASS. J'essaie de créer des bassins dans la région alpine et d'utiliser les données SRTM3. Comme je peux le voir sur le raster de la direction du flux, les données SRTM3 semblent agir correctement: j'utilisais l'algorithme r.watershed avec une...
J'essaie d'écrire un script qui crée une nouvelle classe d'entités surfaciques contenant un seul polygone (carré) avec les coordonnées suivantes: (0, 0), (0, 1 000), (1 000, 0) et (1 000, 1 000) . import arcpy import fileinput import string import os from arcpy import env env.workspace =...
Je veux remplir les valeurs NA en utilisant les valeurs moyennes des voisins les plus proches: r <- raster(matrix(1:16, nrow=8, ncol=8)) r[r==12] <- NA
Tout d'abord, je tiens à dire que je sais qu'une question similaire a été soulevée auparavant, mais qu'elle n'a pas fourni de solution satisfaisante. J'ai besoin d'obtenir le nombre de sommets de chaque objet sur le calque de ligne vectorielle. Sur la base de cet article:
J'essaie d'importer un fichier de formes dans QGIS. Lorsque je l'importe, la limite est affichée dans la liste, mais lorsque j'essaie de la visualiser dans la zone de mappage, elle n'apparaît pas. Je peux le faire afficher les nœuds si je passe en mode édition, donc je sais que les données sont là....
Je développe une application SIG à l'aide d' OpenLayers-3 . Dans cette application, les utilisateurs dessinent des entités et impriment la carte. Supposons que l'utilisateur dessine ses fonctionnalités dans une grande zone. Maintenant, pour imprimer toutes ses fonctionnalités, il effectuerait un...
Je développe un module complémentaire pour QGIS qui doit attribuer un UUID à n'importe quelle fonctionnalité dès qu'il est ajouté à la couche. J'utilise le signal featureAdded pour écrire l'UUID dans le champ correspondant de la fonctionnalité. Voici le code que j'utilise (aussi simplifié que...
Je cherche à automatiser une tâche commune pour nous. Par exemple, Jaromir Jagr est né à Kladno en 1972. Je veux pouvoir prendre cela et produire une chaîne. Né à Kladno, Tchécoslovaquie (actuelle République tchèque). Je connais un peu geonames.org, mais je ne sais pas s'il a la composante...
J'ai un ensemble de données de routes dans une ville stocké dans une base de données compatible PostGIS. Je voudrais identifier et supprimer toute caractéristique qui n'est pas connectée au reste du réseau routier. Ce comportement est courant sur les bords des réseaux, où les données ont été...
J'essaie d'accorder le rôle GeoServer ADMIN aux membres d'un groupe Active Directory. J'ai suivi la documentation ici et ici . Mes utilisateurs Active Directory peuvent se connecter avec succès, mais les membres de mon groupe d'administration Active Directory ne reçoivent pas les privilèges ADMIN...
J'exécute une analyse spécifique où nous utilisons galbée pour créer des tampons autour des points (emplacements de magasin), puis vérifie si d'autres points (emplacements d'utilisateur) sont présents dans cette valeur de tampon. Quand je vérifie la distance en utilisant galbée, elle se révèle être...
Je suis assez nouveau dans tout cela. Je veux avoir un survol / survol ET un clic sur ma fenêtre geojson. Voici le code jusqu'à présent, mais je ne sais pas où inclure le survol sans gâcher la fenêtre contextuelle de clic. Le survol doit afficher le contenu d'un champ de texte et mettre en...
J'essaie de tracer une carte raster en utilisant Pyhton. L'image n'est pas enregistrée, bien qu'elle apparaisse. De plus, les pixels n'apparaissent pas à l'endroit exact. Je reçois l'image comme indiqué ci-dessous (ce qui est faux, cependant, l'étiquetage lat-lon est correct.) Vous trouverez...
J'ai besoin de faire un motif de grille avec des points qui échelonnent pour faire un motif de diamant au lieu du motif régulier qu'une grille en résille fait. Il est parfois fait référence à un "départ décalé". Mes points doivent être à 800 pieds l'un de l'autre. Je crée ma carte dans ArcGIS 10.2....
J'ai une couche polygonale avec de nombreux trous. Je voudrais le généraliser. J'ai utilisé v.generalize d'herbe. Avec lui, j'ai utilisé l'algorithme chaiken. Mais, seules les bordures des polygones ont été généralisées, les trous qu'elles contiennent n'ont pas lissé leurs bords. Que puis-je faire...
J'essaie de mettre en œuvre une vérification très efficace pour voir si deux points sont à moins d'un mile l'un de l'autre, ou non. Mon approche actuelle consiste à calculer la distance Haversine , puis à vérifier si elle est inférieure à un mile. L'efficacité est importante dans ce cas car je dois...
J'ai un DEM raster et plusieurs points. Je veux déterminer l'élévation DEM à chaque point. Hormis l'examen manuel du DEM à chaque point, (cela prend beaucoup de temps avec beaucoup de points), y a-t-il un moyen de le faire sans accéder à l'outil Spatial Analyst Extraire les valeurs en points? J'ai...