Ich habe einen Geopandas-Datenrahmen in EPSG: 4326 und würde einen neuen Datenrahmen erstellen, der aus allen Zeilen besteht, die in einen bestimmten Begrenzungsrahmen fallen. Zuerst erhalte ich den Begrenzungsrahmen, der mir wichtig ist (der eigentlich der Begrenzungsrahmen eines anderen...
Ich möchte wissen, wie man alle Rasterwerte innerhalb eines Polygons mit gdal oder pygeoprocessing erhält, ohne das gesamte Raster als Array zu lesen. pygeoprocessing und gdal können zonale Statistiken erstellen, aber von einer solchen Funktion sind nur die Werte min, max, mean, stdev oder count...
Ich habe ein QGIS-Projekt, das einen Komponisten enthält. Wenn ich in die GUI gehe, kann ich sie natürlich als PDF usw. exportieren. Für mein Make-Skript möchte ich dies jedoch über ein Python-Skript erledigen lassen. Vom Online-Lesen hier ist mein Versuch: #!/usr/bin/env python import sys from...
Wie kann ich bei einem Begrenzungsrahmen die Anzahl bestimmter Werte (z. B. interessiert mich die Anzahl value == 1) in einem Raster in 1) arcpy oder 2) Rasterpaket in R
Ich muss Daten aus einer NetCDF mit einem bestimmten Shapefile unterteilen. Die Daten sind die Meeresoberflächentemperatur und die Ozeanfarbe bei einer Auflösung von 1/4 Grad. Ich habe 4 Polygone, die die USA beschreiben. Großes marines Ökosystem im nordöstlichen Festlandsockel und seine...
Ich versuche, die Seile der Fernerkundungs-Bildverarbeitung mit Python GDAL-Bindungen und Numpy zu lernen. Als ersten Versuch lese ich eine Landsat8-Geotiff-Datei, führe eine einfache Manipulation durch und schreibe das Ergebnis in eine neue Datei. Der folgende Code scheint einwandfrei zu...
Wie kann ich eine Ebene nach dem Erstellen ausblenden / anzeigen? Ich kann das Rendern einer bestimmten Ebene aktivieren / deaktivieren, indem ich das Kontrollkästchen über QGIS aktiviere. Dies muss jedoch programmgesteuert über Python-Code erfolgen. Wie kann ich ein Label im Python-Code ein- /...
Ich suche nach einer Methode, um die X, Y-Koordinate eines inneren Punktes (in der Bogenkarte auch als labelPoint bezeichnet) eines (unregelmäßigen) Polygons zu berechnen. Ich versuche, eine Lösung ohne die Verwendung des ArcGIS-Moduls zu finden Der Unterschied zwischen einem "Schwerpunkt" und...
Was ist der richtige Weg, um mit Python einen Punkt in eine PostGIS-Datenbank
Wenn Sie zwei Punkte haben, aus denen ich ein gerades LineStringObjekt erstellen möchte : from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) Das Shapely-Handbuch fürLineString Staaten: Eine Folge von PointInstanzen ist kein gültiger Konstruktorparameter. A LineStringwird...
Ich habe einige Stapelverarbeitungswerkzeuge als Python-Plugins für QGIS 1.8 entwickelt. Ich habe festgestellt, dass die GUI während der Ausführung meiner Tools nicht mehr reagiert. Die allgemeine Weisheit ist, dass die Arbeit an einem Arbeitsthread ausgeführt werden sollte, wobei die Status- /...
Ich bin Pythonist in einem Büro voller ENVI / IDL-Benutzer. Ich habe einige Tools erstellt, um diese Lücke zu schließen, aber ein großer Stolperstein sind ENVI Classic .roi-Dateien, die meine Kollegen häufig verwenden. Ich kann ENVI 5-ROIs problemlos in Python importieren, da das textbasierte...
Ich weiß, meine Frage ähnelt einigen alten auf dieser Seite. Ich habe viele CSV-Dateien (Geokoordinaten) zum Importieren in qgis (und zum anschließenden Konvertieren), und der übliche Weg ist nicht der beste (zu lange). Ich habe fast 500 CSV-Dateien (wgs84-Koordinaten) und dies ist, was ich tun...
Nachdem ich viele Skripte in der Python-Konsole von QGIS ausgeführt habe, habe ich mehr als 10 Wörterbücher und Listen, die ich für die weitere Analyse benötige. Gibt es eine Möglichkeit, eine Python-Sitzung mit all diesen Variablen zu speichern? Im Allgemeinen ist meine Frage, wie ich meine...
Ich arbeite an einem Tool, das Rasterebenen für die Verarbeitung mit NumPy in Arrays konvertiert. Idealerweise möchte ich mit Rastern arbeiten können, die in einer GDB enthalten sind, ohne sie alle zu exportieren (insbesondere, wenn ArcGIS oder ArcPy aktiviert werden müssen) ). Ist dies mit dem...
Mit dem folgenden Code kann ich 100% Transparenz auf den Rasterwert 0 setzen: map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map= lyr tr=None tr=QgsRasterTransparency() tr.initializeTransparentPixelList(0)
Ich habe in der ESRI .NET-Hilfe Folgendes gefunden: Python Alle ArcGIS Desktop-Anwendungen enthalten eine eingebettete Skriptsprache namens Python. Auf viele Bereiche von ArcGIS - insbesondere die Geoverarbeitung - kann über vereinfachte Python-Anwendungsprogrammierschnittstellen (APIs) zugegriffen...
Ich habe eine Datei im ASCII-Rasterformat. Beispielsweise: ncols 480 nrows 450 xllcorner 378923 yllcorner 4072345 cellsize 30 nodata_value -32768 43 2 45 7 3 56 2 5 23 65 34 6 32 54 57 34 2 2 54 6 35 45 65 34 2 6 78 4 2 6 89 3 2 7 45 23 5 8 4 1 62 ... Wie kann ich es mit Python in TIFF oder ein...
Ich habe einen Geopandas-Datenrahmen aus einer Reihe von Linestrings, denen einige Daten mit jedem Scheitelpunkt / Punkt zugeordnet sind: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Die Punkte werden basierend auf der ID in Linestrings konvertiert und nach Zeit sortiert. Ich möchte die...
Ich versuche, mit Python einen Histogrammabgleich durchzuführen, um den Mosaikprozess mehrerer überlappender Raster zu verbessern. Ich stütze meinen Code auf den folgenden: http://www.idlcoyote.com/ip_tips/histomatch.html Bisher habe ich es geschafft, den überlappenden Bereich zweier benachbarter...