Systèmes d'information géographique

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

1
PyQGIS enregistre le raster en tant qu'image rendue, puis utilise les outils GDAL dessus
J'ai essayé d'utiliser le code de QGis Save Raster as Rendered Image . iest une couche raster: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs = pipelayer.crs().toWkt() pipe = QgsRasterPipe() pipe.set(pipeprovider.clone()) pipe.set(piperenderer.clone()) pipedFile = os.path.join(tempfile.gettempdir(), safeLayerName + '_pipe.tif') print pipedFile …
3
Construire un diagramme de Voronoi dans PostGIS
J'essaie de construire un diagramme voronoi à partir d'une grille de points en utilisant du code modifié d' ici . Il s'agit d'une requête SQL après mes modifications: DROP TABLE IF EXISTS example.voronoi; WITH -- Sample set of points to work with Sample AS (SELECT ST_SetSRID(ST_Union(geom), 0) geom FROM example."MeshPoints2d"), …
5
Méthodes pour convertir un fichier .SVG en GeoJSON
Question Comment convertir un fichier SVG au format GeoJSON? Résultat souhaité: Je veux créer une carte choroplèth en utilisant quelque chose comme D3.js, ou n'importe quelle autre bibliothèque. Défi: J'ai un fichier SVG que j'ai pu obtenir du domaine public des régions administratives pour une province particulière d'un pays particulier. …
12 geojson  svg  d3 
5
Projection fisheye
J'ai besoin d'une projection spéciale pour une présentation efficace sur une affiche. Je voudrais présenter la zone d'étude (République tchèque) dans un contexte de toute l'Europe de telle manière que la République tchèque occupe la plus grande surface (comme 70%) au centre, et le reste de l'Europe est autour d'elle, …
2
Les étiquettes des points qui se chevauchent peuvent-elles être combinées / fusionnées en une seule étiquette?
J'ai des points représentant des exemples d'emplacements. Souvent, plusieurs échantillons seront prélevés au même emplacement: plusieurs points avec le même emplacement mais différents ID d'échantillons et d'autres attributs. Je voudrais étiqueter tous les points qui sont co-localisés avec une seule étiquette, avec du texte empilé répertoriant tous les exemples d'ID …
3
Modèle géologique 3D à partir de coupes transversales, de préférence à l'aide d'un logiciel Open Source
Je travaille avec la géologie de l'ingénieur et mon flux de travail comprend la réalisation de nombreuses coupes transversales d'un site donné, avec des données de cartes, de travaux sur le terrain (affleurements) et de forages. Ce que je recherche, c'est un moyen de prendre de nombreuses coupes transversales et …
2
writeOGR avec un polygone spatial simplifié par gSimplify
J'utilise gSimplify(package rgeos) pour simplifier les géométries d'un fichier de formes. La fonction fonctionne bien, mais maintenant je ne peux pas écrire la sortie dans un nouveau fichier de formes. J'ai essayé plusieurs façons: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') Je reçois obj doit être un SpatialPointsDataFrame, SpatialLinesDataFrame ou SpatialPolygonsDataFrame et …
12 shapefile  r  simplify  rgdal 
1
Pourquoi l'emplacement @ data @ values ​​pour RasterLayer contient-il uniquement des valeurs logiques (0) et non les valeurs réelles?
En essayant de comprendre pourquoi, lorsque je lis dans un raster de NDVI, l'emplacement @ data @ values ​​ne contient pas les valeurs réelles jusqu'à ce que je les définisse manuellement. Par exemple: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns: logical(0) Cela ne s'est pas produit avec d'autres rasters …
12 raster  r  memory 
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.