Geografisches Informationssystem

Fragen und Antworten für Kartographen, Geographen und GIS-Experten

3
Geopandas räumliche Verbindung extrem langsam
Ich verwende den folgenden Code, um ein Land (und manchmal einen Bundesstaat) für Millionen von GPS-Punkten zu finden. Der Code benötigt derzeit ungefähr eine Sekunde pro Punkt, was unglaublich langsam ist. Das Shapefile ist 6 MB groß. Ich habe gelesen, dass Geopandas rtrees für räumliche Verknüpfungen verwenden, was sie unglaublich …
3
QGIS Print Composer ändert sich nicht von Querformat zu Hochformat?
Mit 2.10.1 öffne ich einen neuen Print Composer und ändere sofort die Ausrichtung der Seiteneinrichtung von Querformat in Hochformat. Die Ausrichtung des Bildschirms ändert sich jedoch nicht. Es zeigt weiterhin Landschaft. Das Ändern des Papierformats wirkt sich ebenfalls nicht auf die Größe / Ausrichtung des Print Composer-Fensters aus. Wenn ich …
3
Erstellen eines Voronoi-Diagramms in PostGIS
Ich versuche von hier aus ein Voronoi-Diagramm aus einem Gitter von Punkten mit modifiziertem Code zu konstruieren . Dies ist eine SQL-Abfrage nach meinen Änderungen: 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"), -- Build edges …
5
Methoden zum Konvertieren von .SVG-Dateien in GeoJSON
Frage Wie konvertiere ich eine SVG-Datei in das GeoJSON-Format? Gewünschtes Ergebnis: Ich möchte eine Choropleth-Map mit etwas wie D3.js oder einer der anderen Bibliotheken da draußen erstellen. Herausforderung: Ich habe eine SVG-Datei, die ich für eine bestimmte Provinz eines bestimmten Landes aus dem öffentlichen Bereich der Verwaltungsregionen erhalten konnte. Ich …
12 geojson  svg  d3 
1
Neuprojektion des Rasters von Lat / Lon nach UTM in R?
Ich muss es in eine UTM umwandeln, um den Puffer funktionsfähig zu machen. wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) Während der Puffererstellung, bei der es sich um die letzte Codezeile handelt, wird folgende Warnung ausgegeben: Warning message: In couldBeLonLat(x) : raster …
5
Fischaugenartige Projektion
Ich benötige eine spezielle Projektion, um sie auf einem Plakat wirkungsvoll präsentieren zu können. Ich möchte das Untersuchungsgebiet (Tschechische Republik) in einem europaweiten Kontext so darstellen, dass die Tschechische Republik die meiste Fläche (wie 70%) in der Mitte einnimmt und das übrige Europa in der Nähe ist. aber lächerlich klein. …
2
writeOGR mit einem durch gSimplify vereinfachten räumlichen Polygon
Ich benutze gSimplify(RGEOS-Paket), um die Geometrien eines Shapefile zu vereinfachen. Die Funktion funktioniert gut, aber jetzt kann ich die Ausgabe nicht in ein neues Shapefile schreiben. Ich habe einige Möglichkeiten ausprobiert: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') Ich bekomme obj muss ein SpatialPointsDataFrame, SpatialLinesDataFrame oder SpatialPolygonsDataFrame sein und mit: writePolyShape(simplyshape, file) …
12 shapefile  r  simplify  rgdal 
3
Geologisches 3D-Modell aus Querschnitten, vorzugsweise mit Open Source-Software
Ich beschäftige mich mit Ingenieurgeologie, und mein Workflow umfasst das Erstellen vieler Querschnitte eines bestimmten Standorts mit Daten aus Karten, Feldarbeiten (Aufschlüssen) und Bohrlöchern. Was ich suche, ist eine Möglichkeit, viele Querschnitte zu nehmen und ein 3D-Modell mit Open-Source-Software zu erstellen, selbst mit viel (CAD-) Zeichnungsarbeit, aber vorzugsweise mit wenigen …
2
Können Beschriftungen für überlappende Punkte zu einer Beschriftung zusammengefasst werden?
Ich habe Punkte, die Beispielorte darstellen. Oft werden mehrere Proben an demselben Ort genommen: mehrere Punkte mit demselben Ort, aber unterschiedlichen Proben-IDs und anderen Attributen. Ich möchte alle Punkte, die sich zusammen mit einem einzelnen Etikett befinden, mit gestapeltem Text kennzeichnen, in dem alle Beispiel-IDs aller Punkte an diesem Punkt …
1
Warum enthält der Slot @ data @ values ​​für RasterLayer nur logische (0) und nicht die tatsächlichen Werte?
Ich versuche herauszufinden, warum der Slot @ data @ values ​​beim Einlesen eines NDVI-Rasters erst dann die tatsächlichen Werte enthält, wenn ich sie manuell einstelle. Beispielsweise: NDVI <- raster("./filename.tif", crs="+proj=longlat +datum=WGS84") NDVI@data@values ## returns: logical(0) Dies ist bei anderen Rastern, die ich mit der gleichen Methode geladen habe, nicht geschehen, …
12 raster  r  memory 
2
Wie funktioniert das räumliche Polygon% über% Polygon, wenn Werte in r aggregiert werden?
Ich arbeite an einem Umweltepidemiologieprojekt, bei dem ich Punktexpositionen habe (~ 2.000 industrielle Schweineoperationen - IHOs). Diese IHOs sprühen auf Felder in der Nähe, aber die Fäkalien, Wassertropfen und Gerüche können kilometerweit wandern. Diese Punktbelichtungen erhalten also 3-mi-Puffer, und ich möchte die Anzahl der IHO-Belichtungen (verschiedener Arten - Summe der …
2
Punktüberzeichnung auf Karten in R (zB mit ggplot2) - wie kann man Punkte zur Seite schieben, ursprüngliche Positionen markieren, Nahpunkte kombinieren ...?
Ich habe eine Datenbank mit Sprachen, deren Längen- und Breitengraden und einem Merkmalswert (entweder Kategorie 1, Kategorie 2 oder beides - im Plot sind diese rot, blau und grün markiert). Es kann bis zu drei Punkte pro Sprache geben und natürlich können zwei Sprachpunkte sehr nahe beieinander liegen. name longitude …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.