Systèmes d'information géographique

1

Combler les lacunes en utilisant les voisins les plus proches

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

8 raster r

1

Création de bassin GRASS vs. direction d'écoulement

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 taille minimale …

8 qgis grass watershed basin

3

Créer une classe d'entités surfaciques à partir des coordonnées X, Y à l'aide d'ArcPy?

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 = "D:/UW/Winter_2016/501/w5_more_arcpy8_9_10/ex8" env.overwriteOutput …

8 arcpy cursor valueerror

1

Jointure par attribut / spatialite avec SQL / jointure externe gauche dans QGIS

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 utilisation. De préférence, tous …

8 qgis sql spatialite spatial-join

2

Une transformation Helmert est-elle nécessaire pour convertir entre WGS84 et NAD83?

J'ai lu des choses partout sur Internet qui prétendent que NAD83 et WGS84 sont effectivement équivalents (avec des différences ellipsoïdes mineures). Cependant, si je fais des conversions géodésiques (WGS84) -> ECEF -> géodésiques (NAD83), en utilisant uniquement les paramètres ellipsoïdes WGS84 et GRS80 et aucune transformation Helmert, je ne peux …

8 coordinate-system convert wgs84 nad83

3

Existe-t-il un moyen efficace d'envoyer une grande image raster?

Question assez simple: j'ai deux grandes images raster qui sont stockées dans ArcSDE. Je dois envoyer ces images à quelqu'un d'autre et les publier sur leur site ftp. Quelle est la manière la plus efficace d'envoyer ces rasters? J'ai essayé de créer une géodatabase et de copier des images et …

8 arcgis-desktop raster arccatalog aerial-photography ftp

1

Mapnik échoue à travers une webstack

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: $ Premature end of …

8 python mapnik

2

Problème de superposition de tuiles vectorielles TileStache

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: 900913). Je peux générer les tuiles en …

8 coordinate-system overlay tilestache vector-tiles

2

gdalwarp produisant des rasters vides lorsqu'il est appelé à partir d'un script python, mais pas via la ligne de commande

J'essaie d'agréger certaines données raster (toutes .tif) jusqu'à une résolution plus grossière (de 0,05 degré à 0,25 degré) en utilisant gdalwarp en python, mais la commande ne fonctionne pas. Au lieu d'obtenir une sortie avec une large plage de valeurs, toutes les valeurs de la sortie sont 0. La résolution …

8 python gdal gdalwarp

1

ArcObjects avec jointure gauche

J'attends un code ArcObjects équivalent à la jointure gauche SQL. Par exemple: J'ai le code SQL ci-dessous, c'est une simple jointure pour récupérer les enregistrements des deux tables en fonction de la ChecklistIdvaleur. SELECT T1.ChecklistId, T1.ChecklistName, T2.Latitude, T2.Longitude FROM [dbo].[TableOne] T1 INNER JOIN [dbo].[TableTwo] T2 ON T2.ChecklistId = T1.ChecklistId WHERE …

8 arcobjects definition-query

1

Appel des fonctions PyQGIS à partir d'éditeurs externes (Linux)

Est-il possible d'intégrer des éditeurs Python externes (tels que KDevelop) à QGIS, afin qu'il soit possible d'exécuter des fonctions dans qgis.core, qgis.utilsetc. en dehors de la console QGIS Python? En suivant les directives sur le site Web de QGIS ( http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/intro.html ), j'ai essayé cela, mais il renvoie juste 1 …

8 pyqgis linux

4

Comment calculer les coordonnées max et min pour plusieurs polygones?

Je voudrais calculer les coordonnées latidudinales maximales, médianes et minimales pour un grand ensemble de polygones, en utilisant QGIS ou R. Donc, les points les plus hauts et les plus bas atteints par chaque polygone sur une plage latidudinale. Les coordonnées centrales sont faciles à faire, mais la conversion des …

8 qgis r

3

Comment positionner un guide à un endroit précis de Print Composer?

Plutôt que de faire glisser des guides sur un compositeur d'impression, je voudrais saisir un nombre pour les définir à un point spécifique. Est-il possible de faire cela? Par exemple, je veux m'assurer d'aligner les objets sur les guides de marge à 10 mm des bords de la page. Il …

8 qgis print-composer

5

Quels fournisseurs de géocodage permettent le stockage de données géocodées?

Je voudrais fournir une fonctionnalité où l'utilisateur peut rechercher une adresse / un lieu sur une carte Leaflet, et est capable d'ajouter le marqueur résultant à la carte (qui plus tard peut être enregistré sur le serveur et téléchargé en KML / GPX / GeoJSON format). Le fond de carte …

8 openstreetmap geocoding mapbox bing-maps here-maps

1

Obtenir des ID de couche à l'aide d'ArcPy?

Je travaille dans une session ArcMap, j'ai coché l'option pour obtenir un ID pour chacune des couches de mon document, donc tout le monde l'a. Maintenant, ce dont j'ai besoin, c'est d'un moyen d'obtenir ces valeurs dans un script python (pour les copier automatiquement dans un tableau): il n'y a …

8 arcpy arcgis-10.3 layers