Systèmes d'information géographique

Questions-réponses pour les cartographes, les géographes et les professionnels des SIG

1
Géoréférencer des fichiers LAS en coordonnées géographiques (lat / lon) ou cartésiennes (projetées)?
Tout d'abord, je suis un débutant dans le SIG. J'ai un nuage de points lidar avec des données GPS et de télémétrie associées. Ensuite, je transforme les coordonnées de chaque point en un système de coordonnées ECEF et enfin je transforme les coordonnées angulaires en un système de coordonnées cartésiennes, …
2
Démarrage de QGIS avec virtualenv Python?
existe-t-il un moyen de dire à QGIS d'utiliser une installation Python différente de celle par défaut à l'échelle du système? Par exemple, utiliser du Python virtualenv /myvirtualenv/bin/pythonau lieu de /usr/bin/python? Contexte: je développe un plugin personnalisé qui utilise certains modules qui ne sont installés que dans un certain virtualenv, mais …
1
Importer plusieurs images raster empilées dans R?
Je voudrais importer tous les rasters de mosaïques multiples (chaque mosaïque est composée de 7 couches) du dossier indiqué dans mon R. Ensuite, accédez-y en tant que rasters multicouches uniques. J'ai donc fait: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- …
2
tampon et union dans R sans réduire à un multi-polygone
Je voudrais effectuer une superposition spatiale pour identifier le polygone dans lequel se situe un ensemble de points. Cependant, je veux d'abord tamponner et dissoudre les polygones, de sorte que les points situés n'importe où dans les polygones fusionnés (mais pas dans les trous) seront étiquetés de la même manière …
9 r  sp 
2
Filtrage de l'API Overpass par pays
Est-ce que quelqu'un sait comment effectuer une demande d'overpass-api comme celle-ci? "Montrez-moi toutes les banques en Allemagne"? Je pourrais simplement montrer tous les 'amenity = bank' dans une boîte de délimitation spécifique. Mais l'Allemagne (ou la plupart des pays) ne s'intègre pas très bien dans une telle bbox. Y a-t-il …
3
Comment mesurer la similitude des objets SpatialLines
J'ai créé deux SpatialLinesobjets R: . Ces objets ont été créés de cette façon: library(sp) xy <- cbind(x,y) xy.sp = sp::SpatialPoints(xy) spl1 <- sp::SpatialLines(list(Lines(Line(xy.sp), ID="a"))) Maintenant, je veux en quelque sorte conclure qu'il s'agit de la même ligne tournée et inversée, et que la différence entre eux est égale à …
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.