Ich projiziere Raster in Python mit GDAL neu. Ich muss mehrere Tiffs von geografischen WGS 84-Koordinaten auf WGS 1984 Web Mercator (Auxiliary Sphere) projizieren, um sie später in Openlayers zusammen mit OpenStreetMap und möglicherweise Google Maps zu verwenden. Ich verwende von hier aus Python...
Ich habe ein Shapefile mit Polygonen. Und ich habe eine globale Rasterdatei. Ich möchte die Polygone des Shapefiles auf das Rasterraster legen und den mittleren Rasterwert für jedes Polygon berechnen. Wie kann ich das mit GDAL machen und die Ergebnisse in das Shapefile
Ich bin ein Neuling in Python und möchte wissen, ob es in ArcGIS 10 eine schnelle Methode gibt, um die Zellwerte eines Rasters Pixel für Pixel und die Koordinaten (XY-Koordinate der Pixelmitte zuordnen) mit Python abzurufen. Um dies weiter zu beschreiben, muss ich die Karte X, die Karte Y und den...
Entschuldigung, wenn die folgende Frage etwas dumm ist, aber ich bin nur SEHR neu in dieser ganzen GIS-Sache. Ich versuche, einige projizierte GeoTiff-Bilder mit gdal in Python nach WGS84 zu konvertieren. Ich habe einen Beitrag gefunden, in dem der Prozess zum Transformieren von Punkten in...
Ich möchte lernen, wie NumPy-Arrays zur Optimierung der Geoverarbeitung verwendet werden. Ein Großteil meiner Arbeit umfasst "Big Data", bei dem die Geoverarbeitung oft Tage in Anspruch nimmt, um bestimmte Aufgaben zu erledigen. Ich bin natürlich sehr daran interessiert, diese Routinen zu...
Angesichts einer geopandas GeoDataFrameReihe von Polygonen möchte ich die Fläche der einzelnen Elemente in meiner Liste in km² anzeigen. Dies ist ein ziemlich häufiges Problem, und die in der Vergangenheit übliche Lösung bestand darin, shapelyund pyprojdirekt zu verwenden (z. B. hier und hier )....
Ich arbeite derzeit daran, CANGRID-Klimadaten (bereitgestellt als Surfer Grid-ASCII-Dateien ".grd") in ArcGIS zu importieren. Das Raster ist 95 Zeilen mal 125 Spalten groß. Die Metadaten enthalten die Länge des Ursprungs (linke untere Ecke), die Zellengröße (50 km) sowie die Notenprojektion als...
Ich möchte programmgesteuert eine GeoJSON-Datei mit Polygonen aus einem Shapefile erstellen, aber Attribute aus meiner eigenen Anwendung hinzufügen. Dies ist für ein Shapefile einfach durchzuführen: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes....
Ich habe eine Software (nicht Arc), die jede Nacht auf einem Nicht-Produktions-Computer ausgeführt wird, der Shapefile-Features aus einer externen Datenbank aktualisiert. In regelmäßigen Abständen möchte ich alle Features im Shapefile löschen (nicht die Datei selbst, die verbleiben muss) und die...
Geschlossen . Diese Frage muss gezielter gestellt werden . Derzeit werden keine Antworten akzeptiert. Möchten Sie diese Frage verbessern? Aktualisieren Sie die Frage so, dass sie sich nur auf ein Problem konzentriert, indem Sie diesen Beitrag bearbeiten . Geschlossen vor 8 Monaten . Gibt es mit...
Ich habe angefangen zu lernen, wie man LAS-Daten in Python manipuliert und wollte sehen, wie andere mit LAS-Dateien umgehen. Ich möchte die Punkte lesen (ich verwende ein Numpy-Array) und die Klassen 1 und 2 (nicht klassifiziert und geerdet) in einem separaten Array herausfiltern. Ich habe den...
Kennt jemand eine reine Python-Bibliothek (keine C- oder Java-Abhängigkeiten), um Geometrieoperationen wie Puffer, Vereinigung, Überschneidungen usw. auszuführen? Ich suche nach etwas, das mit Google AppEngine kompatibel
Ich habe versucht, eine Funktion zu erstellen, die im Grunde dasselbe tut, was die QGIS-Funktion "auflöst". Ich dachte es wäre super einfach aber naja anscheinend nicht. Nach allem, was ich gesammelt habe, sollte der Einsatz von Fiona mit Shapely hier die beste Option sein. Ich habe gerade...
Ich arbeite mit einem Add-In von ArcMap in C #. Aus C # -Code habe ich einige Python-Skripte ausgeführt. Um dieses Skript auszuführen, habe ich den Python-Pfad fest programmiert. Dies ist jedoch nicht portabel. Daher möchte ich den Pfad der ausführbaren Python-Datei aus dem Code abrufen und...
Gegebene Datenpunkte mit Längengrad, Breitengrad und einem dritten Eigenschaftswert dieses Punktes. Wie kann ich Punkte basierend auf dem Eigenschaftswert in Gruppen (geografische Unterregionen) gruppieren? Ich suchte bei Google und fand heraus, dass dieses Problem als "räumlich begrenztes...
Ich verwende die GIS-Daten meines Landkreises , um die Straßen in OpenStreetMap zu aktualisieren. Ich habe eine Kopie, die sie letztes Jahr veröffentlicht haben, und eine, die gerade veröffentlicht wurde, und ich möchte alle LineStrings finden, deren Attribute oder Geometrie sich geändert haben....
Ich habe eine Liste von Datenrahmen, alle mit der gleichen Struktur und alle sind Polygone. Wie kann ich sie alle mit Geopandas
Ich benutze GDAL in Python zum Importieren von Shapefile. Ich möchte die Feldnamen für die Datei kennen. Mein aktueller Weg ist: fields = [] for i in range(1, layer.GetFeature(0).GetFieldCount()): field = layer.GetFeature(0).GetDefnRef().GetFieldDefn(i).GetName() fields.append(field) Auf diese...
Ich versuche, alle nicht überlappenden Polygone, die unten angezeigt werden, mit Shapely zu erfassen (gegebene Polygone A, B und C). Darüber hinaus hoffe ich, dies ohne Iteration, Prüfung auf Schnittmenge usw. zu tun. Die akzeptierte Antwort auf diese Frage drückt die PostGIS-Methode aus, aber es...
Versuch, gdal für Python 2.7.3 unter Windows XP zu importieren: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...