Ich habe Probleme bei der Verwendung der Polygonize-Funktion in Python. Das Kochbuchbeispiel hierfür finden Sie hier . Der relevante Teil meines Codes ist: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = driver.CreateDataSource(outShapefile+ ".shp") outLayer = outDatasource.CreateLayer("polygonized", srs=None) …
Ich lese ein Shapefile mit Fiona: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] Dies gibt mir ein Diktat mit einer Liste von Koordinatenlisten und dem Typ 'MultiPolygon'. Ich bin mir jedoch nicht sicher, wie ich die Koordinatendaten zum Erstellen eines Shapely MultiPolygon verwenden …
Gibt es eine Möglichkeit, dieselbe Aufgabe wie das Dienstprogramm gdalbuildvrt mithilfe von GDAL-Python-Bindungen auszuführen? Bisher habe ich keine andere Möglichkeit gefunden, als ein vrt eines einzelnen Datasets zu erstellen und die XML manuell zu bearbeiten. Ich möchte ein vrt aus mehreren Rastern erstellen (im Wesentlichen ein Mosaik ausführen). Ist das …
[EDIT: Die Lösung hierfür bestand einfach darin, OGR zum Lesen von Shapefiles zu verwenden. Siehe das Beispiel von geographika.] In einem ESRI-Shapefile gibt es keinen Unterschied zwischen Polygonen und MultiPolygonen. Darüber hinaus gibt es keine explizite Unterscheidung zwischen Innenlöchern und Außenringen (abgesehen von der "Händigkeit" eines bestimmten Polygons). Nachdem ich …
Ist es möglich, Rohdaten in a geopandas GeoDataFrame, a la a einzulesen pandas DataFrame? Zum Beispiel funktioniert Folgendes: import pandas as pd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") pd.read_json(io.BytesIO(r.content)) Folgendes nicht: import geopandas as gpd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") gpd.read_file(io.BytesIO(r.content)) Mit anderen Worten, ist es möglich, Geodaten zu lesen, die …
Welche Treiber hat das Python-Paket fiona? Wenn ich das Handbuch überprüfe , heißt es [...] and the possible formats are enumerated in the fiona.drivers list. jedoch, wenn ich Python eingebe from fiona import drivers print drivers > <function drivers at 0x108763050> Wie kann ich das "untersuchen"?
Gibt es eine Möglichkeit, durch Ausführen eines für QGIS erstellten benutzerdefinierten Plugins eine Python-Shell anzuzeigen, wenn das Plugin ausgeführt wird, wie z. B. ein Debug-Bildschirm? Ich frage dies, um zu wissen, ob es eine Möglichkeit gibt, Druckbefehle über ein Python-Skript anzuzeigen, während ein Plugin für QGIS programmiert wird.
Gibt es mit Gdal / OGR eine Möglichkeit, eine shp-, kml- oder PostGIS-Leitungsschicht in ein Netzwerk von Knoten und Segmenten umzuwandeln, damit sie in Paketen wie networkx verwendet werden kann?
Ich versuche, ein Shapefile mit Jupyter-Notizbuch zu öffnen, aber ich kann es nicht zum Öffnen bringen. Das Shapefile befindet sich im selben Ordner wie meine Arbeitsmappe, und Python kann es nicht identifizieren, selbst wenn ich den vollständigen Pfad verwende. import shapefile as shp import matplotlib.pyplot as plt sf = shp.Reader(r'C:\Users\Public\Documents\1-11-99n.shp') …
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 qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from …
Ich möchte den Schnittpunkt mehrerer Polygone erhalten. Mit dem Python- shapelyPaket kann ich mithilfe der intersectionFunktion den Schnittpunkt zweier Polygone ermitteln . Gibt es eine ähnlich effiziente Funktion, um den Schnittpunkt mehrerer Polygone zu erhalten? Hier ist ein Code-Snippet, um zu verstehen, was ich meine: from shapely.geometry import Point coord1 …
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 Linestrings an dem Punkt brechen, …
Ich verwende Ubuntu 14.4 und habe eine virtuelle Python-Umgebung installiert. Dann habe ich mehrere Python-Bibliotheken installiert (gdal, lxml, psycopg2 usw.) Auch in meinem System habe ich das Modul für QGIS. Wenn ich also Python eingebe, kann ich dieses Modul importieren als: import qgis.core Aber wenn ich versuche, dasselbe über meine …
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) map.renderer().setRasterTransparency(tr) map.triggerRepaint() Wie kann ich zwei Werte (0 und 8) so definieren, dass sie 100% transparent sind? Wenn ich die Befehle des letzten …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.