Czytam plik kształtu za pomocą Fiony: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] daje mi to dyktando zawierające listę list współrzędnych i wpisz „MultiPolygon”. Nie jestem jednak pewien, jak użyć danych...
Robię drobny przepływ pracy w notatniku ipython z geopandami i kształtnie, wyciągając garść danych geoprzestrzennych, czasem z plików kształtów, czasem z postgis (gdzie odbywa się droższe przetwarzanie). Teraz przeciągam tabele Postgis do Pythona sqlalchemy, przekształcając geometrię do WKT po...
Jeśli mają dwa punkty, z których chcę utworzyć prosty LineStringobiekt: from shapely.geometry import Point, LineString A = Point(0,0) B = Point(1,1) Podręcznik Shapely dlaLineString stanów: Sekwencja Pointinstancji nie jest poprawnym parametrem konstruktora. A LineStringjest opisany punktami,...
Wykonuję bardzo proste obliczenia długości polilinii za pomocą zgrabnego: 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) Moje współrzędne są w WGS84. Nie mogę znaleźć żadnych informacji o...
Mam ramkę danych geopandas wiązki linii, które mają pewne dane związane z każdym wierzchołkiem / punktem: Point_x = (Lat, Lon, Time, ID, Data1, Data2, Data3) Punkty są konwertowane na linie na podstawie identyfikatora i uporządkowane według czasu. Chcę przełamać pasy w miejscu, w którym...
Próbuję utworzyć obrócony prostokąt we współrzędnych lat / lon. Zasadniczo mam współrzędne środka i kąt obrotu wokół środka. Czy można to zrobić za pomocą, foremnego lub gdal, czy czegoś
Próbuję znaleźć miejsce, w którym wiele linii przecina wielokąt dla dwóch różnych ramek geodanych: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 =...
Próbuję generować wielokąty dla satelitarnych pokosów orbit. Do tej pory mam metodę generowania dwóch linii, które reprezentują krawędź każdego pokosu w [łac, długo]. Niektóre pokosy przekraczają międzynarodową linię danych i dlatego obejmują: Udało mi się to rozwiązać za pomocą ogr2ogr...
Chcę przyciąć zestaw polilinii (czarne linie na obrazku poniżej) do zewnętrznej granicy wielokąta. Wszelkie puste przestrzenie w obrębie wielokąta należy zignorować. Moje idealne wyjście to przerywane żółte linie. Linie początkowe mogą, ale nie muszą być proste. Obraz jest uproszczonym przykładem,...
Mam dwa wielokąty, POLYGON ((- 5 -5, -5 5,5 5,5 -5, -5 -5)) POLYGON ((1 1, 1 7, 7 7, 7 1, 1 1)) i chcę je równomiernie podzielić. Wynik powinien być podzielony wzdłuż czerwonej zaznaczonej linii. Używam qgis / python zgrabnie i chcę to zrobić za pomocą