Shapely ist ein Open-Source-Python-Paket für die Konstruktion und Analyse von planaren geometrischen 2D-Objekten, die in geografischen Informationssystemen üblich sind.
Ich habe mehrere Codebeispiele mit Bibliotheken wie Shapefile, Fiona und Ogr ausprobiert, um zu überprüfen, ob ein Punkt (x, y) innerhalb der Grenzen eines mit ArcMap erstellten Multipolygons (und damit im Shapefile-Format) liegt. Keines der Beispiele funktioniert jedoch gut mit Multipolygonen, obwohl sie mit regulären Einzelpolygon-Shapefiles gut funktionieren. Einige Ausschnitte, …
Ich benutze Shapely in Python und bekomme eine MultiLineStringmit ein paar LinestringObjekten. Ich kann garantieren, dass alle LineStringObjekte einfache Linien mit nur 2 Eckpunkten sind und dass sie alle Teil einer einzelnen Linie sind (keine Verzweigungen). Ich möchte "die Punkte verbinden" und eine einzelne erstellen LineString. Muss ich dafür eine …
Wenn ich versuche, Shapely auf meinem Windows 64-Bit-Computer zu installieren, funktioniert die GEOS-Bibliothek nicht. Bisher habe ich das OSGeo4W-Installationsprogramm ausgeführt, von dem aus ich GDAL installiert habe (ich glaube, die Geos-Bibliothek ist in diesem Paket enthalten). Danach überprüfte ich und ich habe geos_c.dll in meinem C:\OSGeo4W\binVerzeichnis, aber entweder habe ich …
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 …
[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 …
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 mache einen kleinen Workflow in ipython notebook mit Geopandas und Shapely und ziehe eine Reihe von Geodaten, manchmal aus Shapefiles, manchmal aus Postgis (wo teurere Verarbeitung durchgeführt wird). Jetzt ziehe ich dann die Postgis-Tabellen in Python mit sqlalchemy, transformiere die Geometrie auf dem Weg in WKT und erhalte so …
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 durch Punkte beschrieben, besteht jedoch nicht aus Punktinstanzen. Wenn ich also …
Ich mache eine sehr einfache Berechnung der Länge einer Polylinie mit Shapely: from shapely.geometry import LineString ... xy_list = [map(float,e) for e in xy_intm] line = LineString(xy_list) s = '%s,%s,%s' % (fr,to,line.length) Meine Koordinaten sind in WGS84. Ich kann anscheinend keine Informationen über das Längenattribut von Shapely finden. Was ist …
Ich versuche, ein gedrehtes Rechteck in Lat / Lon-Koordinaten zu erstellen. Grundsätzlich habe ich die Koordinaten des Zentrums und des Drehwinkels um das Zentrum. Ist es möglich, formschön oder gdal oder etwas anderes zu verwenden?
Um die äußeren Koordinaten zu erhalten, muss ich ein formschönes MultiPolygon in ein Polygon konvertieren. So mach ich es: if poly.geometry.type == 'Polygon': x, y = poly.geometry.exterior.xy elif poly.geometry.type == 'MultiPolygon': allparts = [p.buffer(0) for p in poly.geometry] poly.geometry = shapely.ops.cascaded_union(allparts) x, y = poly.geometry.exterior.xy # here happens the error …
Ich versuche, Polygone für Satellitenorbitalschwaden zu erzeugen. Bisher habe ich eine Methode, um zwei Linien zu erzeugen, die die Kante jedes Schwads in [lat, long] darstellen. Einige der Schwaden überqueren die internationale Datenlinie und wickeln sich so um: Ich konnte dies lösen mit ogr2ogr -wrapdateline: ogr2ogr -wrapdateline -f "ESRI Shapefile" …
Ich möchte eine Reihe von Polylinien (schwarze Linien im Bild unten) an der Außengrenze eines Polygons befestigen. Alle Hohlräume innerhalb des Polygons sollten ignoriert werden. Meine ideale Ausgabe sind die gestrichelten gelben Linien. Die Anfangslinien können gerade sein oder nicht. Das Bild ist ein vereinfachtes Beispiel, in Wirklichkeit ist das …
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.