J'ai un projet QGIS avec des couches, basé sur des requêtes SQL PostGIS qui utilisent la date actuelle. Étant donné que les données de la base de données changent, ces requêtes renvoient des données différentes chaque jour. Est-il possible de changer dynamiquement le nom du calque dans le panneau...
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...
Je souhaite charger des fichiers .adf dans R. Les données proviennent de cette page: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 J'ai essayé le code suivant que j'ai trouvé après quelques recherches sur Internet. Le problème est que, dans la classe RasterLayer, j'obtiens des valeurs...
Hier, j'ai posé une question sur les problèmes que j'éprouvais avec l'étiquetage (qui me tourmentent toujours), mais il m'a frappé aujourd'hui qu'il serait probablement plus facile pour moi de simplement créer des sections agrandies de ma carte pour pouvoir afficher toutes les étiquettes...
J'essaie de projeter un raster. Dans R, il y a la projectRaster()fonction à cela (ci-dessous un exemple entièrement reproductible): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to...
Je travaille donc sur un projet en développement depuis plus d'un an, mais je ne me suis impliqué que récemment. Ma tâche consiste à convertir environ 2 500 fichiers KML dans un autre format, à importer ces fichiers dans un Geoserver, puis à copier les informations de style associées pour chaque...
Fermé . Cette question a besoin de détails ou de clarté . Il n'accepte pas actuellement les réponses. Vous souhaitez améliorer cette question? Ajoutez des détails et clarifiez le problème en modifiant ce message . Fermé il y a 6 mois . J'essaie de créer un tampon pour plusieurs polygones dans une...
J'ai créé un data.frame dans R. Je veux obtenir SpatialLinesDataFrame en conséquence, donc je comprends que je dois convertir le data.frame en lignes, les lignes en SpatialLines et les SpatialLines en SpatialLinesDataFrame. j'ai fait ça filedata.frame=data.frame(matrix(file),ncol=14, byrow=T))...
Travailler avec QGIS 2.4.0 Chugiak J'ai une question concernant l'écrêtage d'images raster. J'ai une photo aérienne que je souhaite découper en utilisant l'étendue d'un fichier de formes. Jusqu'ici tout va bien, j'ai utilisé l'outil raster / extraction / clipper pour ce faire. Le problème est que...
Je veux créer un jeu de données vectorielles ponctuelles de 10000 points (ou plus) dans un raster binaire, où les points doivent être contraints aux zones où la valeur du raster est 1. J'ai essayé les étapes suivantes. Polygoniser un raster QGIS: Vector -> Outils de recherche -> Points...
J'utilise la méthode GDAL ReadAsArray pour travailler avec des données raster à l'aide de numpy (en particulier la reclassification). Comme mes rasters sont volumineux, je traite les tableaux en blocs, en répétant chaque bloc et en les traitant avec une méthode similaire à l' exemple GeoExamples ....
Les quartiers Queen et Rook sont deux méthodes courantes pour calculer les statistiques d'une cellule focale. Ils sont également connus dans les quartiers Moore et (von) Neumann. Je pense qu'il existe également une "version évêque", c'est-à-dire du coin inférieur gauche au coin supérieur droit....
J'ai besoin de trouver un algorithme ou une méthode qui puisse détecter des latitude longitude points aberrants dans une trajectoire pendant le post-traitement , qui peut ensuite être corrigé (ramené dans le chemin de la trajectoire en fonction de ses voisins). À titre d'exemple du type de points...
J'essaie de développer un plugin pour ajouter des couches postgis aux qgis à partir d'un plugin dans le cadre de l'apprentissage de python. Le code pour ajouter une couche postgis est tiré du livre de recettes et fonctionne correctement si je l'exécute à partir de la console python dans qgis (si le...
J'ai un grand projet qui se compose d'environ 10x fichiers de projet QGIS. Chaque projet contient un compositeur d'impression avec un atlas configuré pour générer jusqu'à 10 cartes. Quelle est la manière la plus simple d'automatiser l'ouverture de chaque projet et de créer un PDF à l'aide de...
J'utilise un serveur Windows Server 2012. PostgreSQL 9.3.5, 64 bits POSTGIS 2.1.3 QGIS Server 2.6.1-2 QGIS Desktop 2.8.3 et 2.12 J'utilise une tablette Microsoft Surface Pro 4 à l'aide de QGIS Desktop 2.12 pour modifier certaines couches stockées dans la base de données Postgresql. Les couches de...
Il s'agit d'une question de géométrie SQL Server apparemment simple qui, selon moi, aurait une solution prête à l'emploi, mais je n'ai pas de chance d'en trouver une. Mon intention est de sélectionner tous les enregistrements d'une table contenant des polygones imbriqués (contenus) dans un polygone...
J'utilise la bibliothèque R maptools pour analyser le fichier de formes dans une liste de polygones. La fonction readShapeSpatial shp <- readShapeSpatial("<path to my shapefile>") me donne un objet SpatailPolygonDataFrame . Dans mon exemple, mon SpatialPolygonsDataFrame a les colonnes...
J'ai un certain nombre de grilles hexagonales de 1 km qui couvrent divers comtés aux États-Unis dans une base de données postgreSQL / postGIS. Chaque grille a le CRS EPSG: 3857, et la couche comtés a EPSG: 3857. Lorsque vous affichez les grilles avec les comtés dans QGIS, tout semble superbe. Mais...
J'essaie de créer une choropleth()carte correcte dans R en utilisant le GISToolspackage. J'utilise choro.legend()pour afficher ma légende. Cependant, je ne peux pas créer de placement de légende correct, ET / OU Je ne comprends pas le rôle du fmtparamètre dans ma choro.legend()fonction. Il me...