Ich möchte Fokusstatistiken für jede Zelle eines Rasters in der Nähe eines bestimmten Kriteriums berechnen. Hintergrund - Ich habe drei binäre Raster, die jeweils einen bestimmten Vegetationstyp repräsentieren. Ich möchte die prozentuale Abdeckung jedes Vegetationstyps innerhalb von (z. B.) 20 km ^...
Ich möchte einen Karten-Service auf einem ArcGIS-Server veröffentlichen, der auf Google Maps oder Bing Maps eingeblendet wird. Hierzu muss der Kartendienst in EPSG 3857/900913 vorliegen. Bei der Google Maps-Projektion soll es sich um eine sphärische Erdprojektion handeln. Also ging ich zu...
Ich versuche, eine Zahl zu erhalten, indem ich zähle, wie viele Punkte sich in einem Shapefile befinden. Und das funktioniert, außer dass ich Probleme habe, wenn ich diese Nummer woanders benutze. Irgendwann werde ich diese Zählung in etwas Mathematik (Feldrechner) verwenden, aber während des...
Auf dieser Informationsseite finden Sie einen Datensatz, der auf esri.com gehostet wird: http://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ImageCat_NLT/MapServer/layers Ich versuche herauszufinden, wie ich Zugriff auf die Rohdaten bekomme (entweder Rohdaten / Längenkoordinaten oder...
Ich habe eine Liste von ein paar hundert Städten mit ihren Breiten- und Längengraden. Bei einem anderen Ort (auch in Lat / Long) muss ich die nächste Stadt finden. Da ich kein GIS verwende, besteht der naheliegende Algorithmus darin, eine Schleife für alle Städte zu erstellen und den Abstand...
Mit "benutzerdefiniert" meine ich keine andere Struktur der KML-Datei, aber im <description></description>Abschnitt habe ich Daten mit HTML-Tags dieser Art: <?xml version="1.0" encoding="UTF-8"?> <kml xmlns="http://www.opengis.net/kml/2.2"> <Document> <Placemark>...
Ich habe ein kleines Problem mit der Python-OGR-API. Ich versuche, alle Koordinaten jedes Scheitelpunkts des äußeren Rings eines Polygons abzurufen. Das habe ich bisher: import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file = open(out,'w') for filename in...
Ich versuche, eine Vereinigung auf einem gemeinsamen Feld durchzuführen, nachdem zwei benachbarte Shapefiles zusammengeführt wurden. Zwischen den Shapefiles befindet sich mindestens ein dünner Streifen. Wenn ich versuche, eine Gewerkschaft zu gründen, wird der folgende Fehler angezeigt: Fehler in...
Mein Ziel ist es, ein vorhandenes Shapefile durch Zusammenführen bestimmter Polygone zu ändern. Nachdem ich das Shapefile importiert und den Befehl UnionSpatialPolygons verwendet habe, erhalte ich den gewünschten Polygonumriss. Dies ist jedoch jetzt ein SpatialPolygons-Objekt und kein...
Dieser Beitrag ist ein Community-Wiki . Bearbeiten Sie vorhandene Antworten, um diesen Beitrag zu verbessern. Derzeit werden keine neuen Antworten akzeptiert. Ich entwickle Code in Python und möchte die GDAL / OGR-Python-Bindungen verwenden. Was ist das Minimum, das ich installieren muss, um die...
Ich habe ein Bare-Earth-DEM aus LIDAR-Daten (A). Es funktioniert gut und hat alle Genauigkeitstests bestanden. Es wird in Sturmflutmodellen verwendet. Das Problem, das wir haben, ist Unterführungen. Ich muss Wasser durch Unterführungen fließen lassen, wenn eine Straße die andere kreuzt. Mein...
Ich möchte ein GeoPDF mit einer Vektorentität in ein Vektorformat konvertieren und die realen Koordinaten dieser Vektoren beibehalten. Das GeoPDF wurde aus ArcMap exportiert und enthält Layer. Ich habe GDAL ausprobiert, aber es rastert die PDF-Datei, bevor es exportiert wird. Welches ist nicht gut...
What3words baute auf der Erde ein Gitter von 3 Quadratmetern. Weiß jemand, was der Algorithmus ist? Ich habe mit dieser Methode begonnen: Der Längengradbereich ist [-180, 180] und der Breitengradbereich ist [-90,90]. Wenn Sie sie in gleiche Intervalle schneiden, können Sie fast perfekte Quadrate in...
Ich möchte ein SpatialPolygonsDataFrame in R verallgemeinern , um das Zeichnen effizienter zu gestalten. Die Polygone, die ich habe, führen zu riesigen PDF-Plots, die langsam geladen und schwer in ein Dokument zu integrieren sind. Ich habe versucht, die gSimplifyFunktion aus dem rgeosPaket zu...
Ich versuche, den Feldwert eines Shapefiles zu ändern. Der Wert ändert sich jedoch nicht. Was ist los mit meinem Code? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature = layer.GetNextFeature() dist = 233...
Mit viel Hilfe von gis.stackexchange.com habe ich Führungslinien für Punktsymbole erstellt. Die Linien basieren auf einer Geometriegenerator-Symbologie: Liniengeometrie mit Ausdrucksdialog in QGIS kürzen? Jetzt möchte ich Führungslinien für verschobene Beschriftungen implementieren. Mir ist sehr...
Ich versuche herauszufinden, wie ein Raster in eine PostGIS2.0-Datenbank geladen wird (ich habe hier und hier vorhergehende Fragen zu diesem Thema gestellt ). Ich versuche, das raster2pgsql.exemit PostGIS 2.0 bereitgestellte Programm zu verwenden. Nachdem Sie herausgefunden haben, dass die...
Ich möchte den Feldrechner in ArcMap verwenden, um eine vorhandene Spalte auf zwei Dezimalstellen zu runden. Derzeit habe ich eine Spalte mit 6 Dezimalstellen und möchte sie einfach auf 2 Dezimalstellen abrunden. Ich hatte geplant, den Field Calculator (möglicherweise mit Python) zu verwenden, aber...
Ich versuche einfach, eine .shp-Datei in ein geoJSON-Format zu konvertieren: ogr2ogr -f geoJSON output.json input.shp Nach dem Ausführen des Befehls scheint nichts mehr zu stimmen. Hier ist ein Auszug aus der output.json { "type": "FeatureCollection", "features": [ { "type": "Feature",...
Wie können Sie eine Vektorebene mit überlappenden Polygonen in ein Raster umwandeln, in dem jede Zelle die Anzahl der Polygone zählt, die diese Zelle enthalten? Wir untersuchen Polygone von Brandherden zwischen 1976 und 2000 und möchten wissen, wie viele Brände in jeder Rasterzelle enthalten...