J'ai une liste de quelques centaines de villes avec leur latitude / longitude. Étant donné un autre emplacement (également en lat / long), je dois trouver la ville la plus proche. Comme je n'utilise aucun SIG, l'algorithme évident consiste maintenant à faire une boucle pour toutes les villes, en...
J'essaie de comprendre comment charger un raster dans une base de données PostGIS2.0 (j'ai posé des questions précédentes sur ce sujet ici et ici ). J'essaie d'utiliser le raster2pgsql.exeprogramme fourni avec PostGIS2.0. Après avoir compris que l'invite de commande dans Windows doit être exécutée...
Est-il possible de créer des vues PostgreSQL sur une table contenant des données raster? J'ai réussi à créer la vue mais l'entrée raster_columnsn'est pas correcte (toutes les informations raster sont manquantes). De plus, je n'ai pas pu ouvrir la vue SQL à l'aide de QGIS (probablement à cause des...
J'ai un petit problème avec l'API Python OGR. Ce que j'essaie de faire, c'est d'obtenir toutes les coordonnées de chaque sommet de l'anneau extérieur d'un polygone. Voici ce que j'ai jusqu'à présent: import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file =...
Par coutume, je ne veux pas dire une structure différente du fichier KML, mais dans la <description></description>section, j'ai des données avec des balises HTML de ce type: <?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Document>...
Avec beaucoup d'aide de gis.stackexchange.com, j'ai créé des lignes de repère pour les symboles de point. Les lignes sont basées sur une symbologie de générateur de géométrie: Raccourcir la géométrie des lignes à l'aide de la boîte de dialogue d'expression dans QGIS? Maintenant, je veux implémenter...
J'ai une classe Contours Featureclass, avec un intervalle de 0,2 mètre. Je souhaite sélectionner uniquement les contours qui sont à 1 m d'intervalle? J'ai essayé d'utiliser la syntaxe suivante dans la boîte de dialogue "Sélectionner par attribut" (dans ArcMap), mais elle sélectionne tous les...
J'essaie d'obtenir un nombre en comptant le nombre de points dans un fichier de formes. Et cela fonctionne, sauf que j'ai du mal à utiliser ce numéro ailleurs. Finalement, je vais utiliser ce nombre dans certaines mathématiques (calculatrice de terrain), mais pendant le débogage, je rencontre une...
Ce message est un wiki communautaire . Modifiez les réponses existantes pour améliorer ce post. Il n'accepte pas actuellement de nouvelles réponses. Je développe du code en Python, et je veux utiliser les liaisons python GDAL / OGR. Quel est le minimum que je dois installer pour utiliser les...
Je veux convertir un GeoPDF qui a une entité vecteurs en un format vectoriel et conserver les coordonnées réelles de ces vecteurs. Le GeoPDF a été exporté depuis ArcMap et contient des couches. J'ai essayé GDAL, mais il pixellise le PDF avant de l'exporter. Ce qui n'est pas bon pour moi. De plus,...
J'ai un MNT de terre nue fabriqué à partir de données lidar (A). Il fonctionne bien et a passé tous les tests de précision. Il est utilisé dans les modèles d'ondes de tempête. Le problème que nous rencontrons concerne les passages souterrains. Je dois permettre à l'eau de s'écouler dans les...
J'essaie d'effectuer une union sur un champ commun après avoir fusionné deux fichiers de formes adjacents. Les fichiers de formes se retrouvent avec au moins un mince ruban d'espace entre eux. Lorsque je tente une union, j'obtiens l'erreur de trou orphelin suivante: Erreur dans
J'espère que quelqu'un a déjà compris celui-ci. Je viens d'installer Geoserver 2.9 sur une distribution Ubuntu 16.04 vanilla. La méthode Geoserver 2.8 d'activation de CORS avec la classe shanbe.hezoun ne fonctionne plus avec Jetty 9.2.13. Il est mentionné que la prise en charge de CORS est déjà...
Je voudrais généraliser un SpatialPolygonsDataFrame en R pour rendre le traçage plus efficace. Les polygones dont je dispose se traduisent par d'énormes tracés pdf qui sont lents à charger et difficiles à intégrer dans un document. J'ai essayé d'utiliser la gSimplifyfonction du rgeospackage, mais...
J'ai un raster (DEM USGS en fait) et je dois le diviser en petits morceaux comme le montre l'image ci-dessous. Cela a été accompli dans ArcGIS 10.0 à l'aide de l'outil Fractionner le raster. Je voudrais une méthode FOSS pour ce faire. J'ai regardé GDAL, pensant sûrement qu'il le ferait (en quelque...
Je cherche à calculer des statistiques focales pour chaque cellule d'un raster, dans un voisinage d'un critère spécifié. Contexte - J'ai trois rasters binaires, chacun représentant un seul type de végétation d'intérêt. J'aimerais calculer le pourcentage de couverture de chaque type de végétation...
J'utilise actuellement le script suivant pour ajouter des données d'attribut d'une table à de nombreux fichiers de formes individuels: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in
En regardant cette page d'informations pour un ensemble de données hébergé sur esri.com: http://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ImageCat_NLT/MapServer/layers Essayer de comprendre comment je peux avoir accès aux données brutes (coordonnées brutes lat / lng ou SHP). Il semble...
J'ai un fichier de formes avec des polygones. Et j'ai un fichier raster global. Je veux superposer les polygones du fichier de formes sur la grille raster et calculer la valeur raster moyenne pour chaque polygone. Comment puis-je faire cela en utilisant GDAL, en écrivant les résultats dans le...
J'utilise ArcGIS 10.2 for Desktop, et j'ai parcouru ce site, et je n'arrive pas encore à comprendre ma réponse. Très nouveau pour VBA et Python, mais beaucoup d'années avec ArcGIS. Je sais que je peux le faire lentement avec Select By Attributes, mais cela prend du temps. J'essaie une jonction...