Geografisches Informationssystem

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

1
Vereinfachen Sie Polygone von sf-Objekten
Wie kann ich ein sfPolygon vereinfachen, ohne Lücken und Splitter einzufügen? Mit einem Shapefile würde ich zum Beispiel Folgendes verwenden rmapshaper::ms_simplify(): library("pryr") library("rgdal") library("rmapshaper") download.file("https://borders.ukdataservice.ac.uk/ukborders/easy_download/prebuilt/shape/England_gor_2011.zip", destfile = "regions.zip") unzip("regions.zip") regions <- readOGR(".", "england_gor_2011") object_size(regions) # ~13MB regions <- ms_simplify(regions) object_size(regions) # < 1MB Ich habe ausprobiert, sf::st_cast()was aus den Manpages …
14 r  simplify  sf 
1
Teilen von Polygonen am Mittelpunkt mit ArcPy?
Ich versuche, ungefähr 4000 Polygone an ihrem Mittelpunkt senkrecht zu ihrer längsten Achse (dh über die Breite am Mittelpunkt) aufzuteilen, wie im folgenden Diagramm gezeigt. Im Idealfall möchte ich dies automatisch tun und vermeiden, jedes Polygon manuell zu teilen. Ich habe den Mittelpunkt von Polygonen durch Konvertieren der längsten Linien …
3
Unterscheidung zwischen pixelbasierter und objektbasierter Klassifikation?
Ich kämpfe darum, die Unterscheidung zwischen pixelbasierter und objektbasierter Klassifizierung im Fernerkundungsbereich klar zu verstehen, und hoffe, dass jemand aus dieser Community Einblicke gewähren kann. Nach den Informationen, die ich bisher habe, geht mein derzeitiges Verständnis in diese Richtung: Pixelbasierte Klassifizierung: Die Klassifizierung erfolgt auf Pixelebene, wobei nur die für …
3
Über widersprüchliche Bezeichnungen iterieren
So wie ich es sehe, gibt es im Wesentlichen drei Möglichkeiten, wenn zwei zu beschriftende Features zu nahe beieinander liegen. Standardmäßig lässt QGIS eine der Beschriftungen weg. Wenn ich die Eigenschaft "Alle Beschriftungen anzeigen" entweder für den betreffenden Layer oder für die globalen Beschriftungseinstellungen aktiviere, werden überlappende Beschriftungen gedruckt. Wenn …
14 qgis  labeling 
2
Wie übersetzt man John Snow-Datensätze in Kartenkoordinaten?
Im HistData-Paket für R ( https://r-forge.r-project.org/R/?group_id=574 ) sind die Datensätze enthalten, die sich auf John Schnees Karte des Cholera-Ausbruchs in London von 1854 beziehen sind maßgebend, sorgfältig digitalisiert unter Aufsicht von Walter Tobler. Einige Details zu diesen Datensätzen werden von John Mackenzie unter http://www1.udel.edu/johnmack/frec480/cholera/cholera2.html beschrieben . Leider verwenden die Koordinaten …
1
Warum funktioniert Shapelys Schnappschuss (GEO-Schnappschuss) nicht wie erwartet?
Ich versuche, mit Shapely / Geopandas zwei Linien miteinander zu verbinden, aber das Ergebnis des Schnappens ist sehr merkwürdig. Ich habe es versucht : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) …
1
Pixelwerte an einem einzelnen Punkt mit Raster abrufen
Um mit rasterio einen einzelnen Pixelwert an einem Punkt in einem Raster zu erhalten, finden Sie hier ein Beispiel: https://github.com/mapbox/rasterio/pull/275 Gibt es jedoch eine direkte API in rasterio (und nicht in der cli), mit der der Wert an einem einzelnen Punkt in einem Raster extrahiert werden kann? - BEARBEITEN with …
14 python  numpy  rasterio 
1
Unterschied zwischen radiometrischer Kalibrierung und radiometrischer Korrektur
Ich bin neu im Fernerkundungsbereich. Ich war manchmal verwirrt, als ich über die Bildvorverarbeitung las. Kann mir bitte jemand den Unterschied und einige Beispielmethoden für die radiometrische Korrektur und radiometrische Kalibrierung erläutern? Wird die Konvertierung von DN-Werten in TOA-Werte (Top Of Atmospheric) von 2 Bildsätzen (z. B. Landsat 7) aus …
1
Automatische Aktualisierung der WMS-Ebene in OpenLayers 3 anhalten?
Ich habe eine Anwendung, die vor dem Abrufen von WMS-Kacheln vom Server eine HTTP-Anforderung ausführen muss (die einige Kartentransformationen verarbeitet). Gibt es eine Möglichkeit, die automatische Aktualisierung von TileWMS nach dem Ziehen der Karte anzuhalten? Das Festlegen von "unsichtbar" oder das Entfernen und Hinzufügen der Ebene funktioniert, aber ich möchte …
14 wms  openlayers 
2
Wie berechnet man den Begrenzungsrahmen eines Geojson-Objekts mit Python oder Javascript?
Ich habe eine Django-App mit Broschüre im Frontend und muss die bbox eines Geojson-Objekts berechnen, damit ich die Koordinaten weitergeben kann map.fitBounds() Ich habe einige Bibliotheken ausprobiert sowie die getBoundsMethoden bestimmter Objekte in Leaflet ausprobiert FeatureGroup, aber es wurde beanstandet, dass die Grenzen nicht definiert wurden. Kann mir jemand eine …
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.