J'essaie de traiter certaines données raster à l'aide d'ogr / gdal et je n'arrive pas à obtenir la pleine utilisation de tous les cœurs de ma machine. Lorsque je n'exécute le processus que sur un seul cœur, j'obtiens une utilisation à 100% de ce cœur. Lorsque j'essaie de diviser en multicœurs (dans...
Lors de nos relevés aériens de la faune, nous générons généralement des dizaines de fichiers GPX et convertissons généralement la plupart d'entre eux en fichiers de formes à afficher dans une application SIG. Les fichiers GPX sont un excellent moyen de conserver toutes (ou presque) toutes vos...
J'ai besoin d'aide ... J'ai des problèmes avec QGIS. J'ai des rasters qui semblent avoir des valeurs nulles, même si je les ai reclassés pour obtenir un raster booléen :( Comment puis-je changer les valeurs NULL à
Je suis surtout un vecteur de vecteurs, mais pour un projet en cours, je dois travailler avec un tas de vieilles cartes papier numérisées (cartes de dégâts de bombe ww2 pour Londres, si vous êtes intéressé!) Nous avons numérisé et géoréférencé les cartes, et souhaitons maintenant créer une couche...
J'utilise ArcGIS 10.1 et je souhaite créer un nouveau raster basé sur deux rasters préexistants. Le RasterToNumPyArray a un bon exemple que je veux adapter. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape =...
J'ai des ortophotographies d'une zone dont j'ai besoin pour servir de fond à mes vecteurs. Je l'ai obtenu sous la forme de 2500 fichiers (71,5 Mo chacun) au format TIFF brut avec le fichier de mots TFW correspondant - 180 Go de données. Le système de coordonnées est local et il correspond à mes...
Existe-t-il un moyen de définir tous les pixels de valeur 0 d'un raster sur nodata à l'aide de GDAL? Je voudrais utiliser fillnodata pour remplir certains trous, mais les trous sont définis sur 0, je dois donc les supprimer
J'ai réussi à ajouter QgsVectorLayerà QGIStravers un plugin python, par exemple: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Mais comment puis-je le
Fermé . Cette question doit être plus ciblée . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin qu'elle se concentre sur un problème uniquement en modifiant ce message . Fermé il y a 2 ans . Contexte C'est ma deuxième question...
Je veux projeter cette carte en projection robinson: library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) Et je voudrais changer la projection en "Robinson" (en suivant les conseils de la réponse à ma question précédente: Quelle projection...
J'analyse les chemins probables empruntés par les navetteurs pour se rendre au travail, sur la base des données de flux du recensement, des distributions de densité d'emploi et du réseau de transport. Pour le moment, j'utilise pgRouting et le module complémentaire pgRouting pour QGIS (décrit ici )...
Étant donné un objet spatial dans R, comment puis-je découper tous ses éléments pour qu'ils se trouvent dans un cadre de sélection? Il y a deux choses que j'aimerais faire (idéalement, je saurais faire les deux, mais l'une ou l'autre est une solution acceptable à mon problème actuel - restreindre...
J'ai posté sur le forum arcgis mais je n'ai jamais eu de réponse. Quelqu'un peut-il me dire comment résoudre ce problème? Supposons que vous ayez un dossier de suivi de véhicule du 1er octobre au 31 décembre avec des informations sur la vitesse du vent. Toutes les données sont enregistrées dans la...
J'ai une question concernant le tracé d'une carte de choroplèthes. Que recommande-t-on pour tracer une carte avec des coupures personnalisées? Et rendez-le un peu attrayant. Dans mon script, j'ai utilisé la commande plot and spplot pour visualiser certaines données mais je ne suis pas entièrement...
Existe-t-il une bibliothèque JavaScript open source comme Leaflet ou OpenLayers qui peut être utilisée dans une application Web ou mobile et qui charge des cartes vectorielles au lieu d'images? Je n'ai besoin que des informations routières et de certains POI. La seule bibliothèque que j'ai trouvée...
J'ai utilisé la gIntersectfonction du rgeospackage pour déterminer quels polygones se croisent dans deux SpatialPolygonsDataFrames (spdfs). Le résultat est une matrice logique avec des noms de domaine comprenant les ID de polygone de sfdf # 1 et des noms de cols les ID de polygone de spdf # 2. Je...
J'avais un plugin fonctionnel dans QGIS 1.8 qui lisait les données d'une base de données MSAccess et les ajoutait à une série de couches de mémoire. Il y a un certain traitement impliqué dans l'intérim, donc je ne pense pas que simplement utiliser QGIS pour lire directement à partir de la base de...
Ma question concerne l'extraction des valeurs des rasters à l'emplacement des points. Avec l'extrait de fonction, c'est très facile, et la fonction me donne une trame de données avec les valeurs de toutes les variables dans les points. Je veux avoir dans cette trame de données les coordonnées de...
Fermé . Cette question est basée sur l'opinion . Il n'accepte pas actuellement les réponses. Voulez-vous améliorer cette question? Mettez à jour la question afin d'y répondre avec des faits et des citations en modifiant ce message . Fermé il y a 4 ans . Je suis analyste SIG depuis plus de 10 ans...
J'ai ces deux polygones: library(sp); library(rgeos); library(maptools) coords1 <- matrix(c(-1.841960, -1.823464, -1.838623, -1.841960, 55.663696, 55.659178, 55.650841, 55.663696), ncol=2) coords2 <- matrix(c(-1.822606, -1.816790, -1.832712, -1.822606, 55.657887, 55.646806, 55.650679,...