Systèmes d'information géographique

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

5
Créer un diagramme de voronoi à partir de segments de ligne
Je cherche un moyen de créer un diagramme de voronoi basé sur des segments de ligne au lieu de points. Voir l'exemple ci-dessous (tiré de cette présentation ). Idéalement, j'aimerais quelque chose que je puisse écrire à partir de Python, mais une solution utilisant ArcGIS ou similaire serait également acceptable. …
2
Propriétés et opérateurs SpatialPointsDataFrame dans R
J'ai créé un objet de type en SpatialPointsDataFrameutilisant le sppackage dans R. Cependant, je suis confus au sujet des @, $, . and []opérateurs et quand les utiliser pour accéder aux différentes propriétés de mon objet. Voici mon exemple de code: library(sp) library(rgdal) #creating a SpatialPointsDataFrame with sample points in …
14 r  syntax 
1
PostGIS: Count number of vertices of lines and polygons in tables
I have 23 tables in a PostGIS schema which I need to count the number of vertices of. The tables are a mixture of lines and polygons so realised I need to use ST_NPoints(geom) So I ran the following query SELECT count(ST_NPoints(geom) FROM lines; the result/count equals the number of …
14 postgis 
2
Python, GDAL et construction de tables d'attributs raster
J'ai un raster entier pour lequel je voudrais construire une table d'attributs raster en utilisant Python et GDAL. Je peux créer une table d'attributs raster GDAL en Python comme suit: >>> rat = gdal.RasterAttributeTable() Cela fonctionne bien, comme nous pouvons le voir: >>> rat <osgeo.gdal.RasterAttributeTable; proxy of <Swig Object of …
2
Comment diviser une fonctionnalité non contiguë dans QGIS?
J'ai un fichier de formes avec une fonctionnalité qui n'est pas contiguë. Je veux séparer les zones à gauche de la zone à droite. Cependant, l'outil Diviser les fonctionnalités de QGIS vous oblige à tracer une ligne sur la fonctionnalité. Lorsque je trace une ligne entre les deux zones de …
14 qgis  shapefile 
2
Identifier l'emplacement de l'image?
J'ai une image satellite d'une seule partie d'une ville avec beaucoup de végétation. J'ai besoin de savoir où est prise cette image. Existe-t-il un moyen de savoir où il est pris sur Google Maps, Google Earth ou quelque chose de similaire? C'est une image .jpg donc je ne sais pas …
2
Comment créer une légende d'histogramme dans QGIS print composer?
J'ai trouvé ce billet de blog en cherchant un moyen de construire dynamiquement un graphique (histogramme) dans le compositeur d'impression. La fonctionnalité que je souhaite répliquer à partir du blog consiste à dessiner des rectangles statiques dans le compositeur d'impression qgis, puis à les redimensionner en fonction d'une valeur de …
1
Dépliant: Comment utilisez-vous removeLayer?
Ça me rend fou. Voir le tutoriel Leaflet pour ajouter un contrôle de couche à votre carte: L.Icon.Default.imagePath = "Scripts/images"; var mapUrl = 'http://{s}.tile.openstreetmap.se/hydda/full/{z}/{x}/{y}.png'; var secondMap = 'http://{s}.tile.osm.org/{z}/{x}/{y}.png'; var littleton = L.marker([39.61, -105.02]).bindPopup('This is Littleton, CO.'), denver = L.marker([39.74, -104.99]).bindPopup('This is Denver, CO.'), aurora = L.marker([39.73, -104.8]).bindPopup('This is Aurora, CO.'), …
2
Convertir .tif avec .tfw en GeoTiff?
J'ai un fichier image tif et un fichier mondial .tfw (créé à partir de données OSM avec Mapertive). Comment puis-je (en utilisant la ligne de commande linux) convertir cela en un fichier GeoTIFF? (c'est-à-dire que je veux avoir un fichier, pas 2). Je suis sûr que c'est possible avec gdal …
6
Lecture de données brutes dans des géopandas
Est-il possible de lire des données brutes en a geopandas GeoDataFrame, a la a pandas DataFrame? Par exemple, les travaux suivants: import pandas as pd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") pd.read_json(io.BytesIO(r.content)) Ce qui suit ne: import geopandas as gpd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") gpd.read_file(io.BytesIO(r.content)) En d'autres termes, est-il possible …
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.