Szukam sugestii, jak ulepszyć mój kod Python. Zwykle wydajność nie ma dla mnie znaczenia, ale teraz pracuję z plikiem tekstowym lokalizacji w USA z ponad 1,5 miliona punktów. Przy danej konfiguracji uruchomienie operacji w jednym punkcie zajmuje około 5 sekund; Muszę obniżyć tę liczbę. Używam trzech różnych pakietów GIS Pythona, …
Próbuję uruchomić Shapely na nowej maszynie i próbowałem zainstalować go ze źródła, używając instalatora na jego stronie internetowej i używając pip install shapely. Każde odwołanie, które próbuję zrobić z biblioteki, z wyjątkiem from shapely import *błędu WindowsError: [Error 126] The specified module could not be found. Czy ktoś miał podobne …
Czy ktoś może zademonstrować prosty sposób zapisywania struktur danych geometrii z kształtnych do plików kształtów? Szczególnie interesują mnie wielokąty z dziurami i obramowaniami. Byłoby również korzystne, aby trzymać się z daleka od arcpy (więc osgeo, pyshp itp. Byłoby lepiej).
Czy istnieje prosty sposób na przekształcanie obiektów foremnych (a mianowicie wielokątów i wielokątów) z jednego rzutu na drugi bez konieczności ręcznego kopania i wyodrębniania współrzędnych? W rzeczywistości nie dbam nawet o to, czy w tym momencie są one obiektami foremnymi, chcę tylko przekazać funkcje i projekcję i odzyskać ponownie zestaw …
Chcę obliczyć statystyki ogniskowe dla każdej komórki rastra, w granicach określonych kryteriów. Tło - Mam trzy binarne rastry, z których każdy reprezentuje jeden interesujący rodzaj roślinności. Chciałbym obliczyć procent pokrycia każdego rodzaju roślinności w obrębie (np.) 20 km ^ 2 dowolnej komórki w moim obszarze badań (suma / całkowita liczba …
tło Ze znanego punktu wymagam ustalenia najbliższego otaczającego „widocznego obwodu” względem tabeli MultiLineStrings, jak pokazano na schemacie. Przeszukałem tę stronę z wieloma terminami (np. Minimalna krawędź, minimalny obwód, najbliższy sąsiad, klip, zawierający wielokąt, widoczność, przyciąganie, wycinanie węzłów, ray-trace, wypełnienie zalewowe, wewnętrzna granica, routing, wklęsły kadłub), ale nie mogę znaleźć żadnego …
Biorąc pod uwagę geopandas GeoDataFrameserię wielokątów, chciałbym uzyskać powierzchnię w km2 każdej funkcji na mojej liście. Jest to dość powszechny problem, a zwykle proponowane rozwiązanie w przeszłości był w użyciu shapelyi pyprojbezpośrednio (np tutaj i tutaj ). Czy istnieje sposób, aby to zrobić w czystej postaci geopandas?
Próbuję uchwycić wszystkie nie nakładające się wielokąty wskazane poniżej za pomocą Shapely (podane wielokąty A, B i C). Co więcej, mam nadzieję, że zrobię to bez iteracji, testowania przecięcia itp. Przyjęta odpowiedź na to pytanie wyraża metodę PostGIS, ale wydaje się, że „związek” oznacza różne rzeczy dla różnych ludzi.
Próbowałem stworzyć funkcję, która robi w zasadzie to samo, co funkcja „rozpuszczenia” QGIS. Myślałem, że to będzie super łatwe, ale najwyraźniej nie. Więc z tego, co zebrałem, użycie fiona z kształtnym powinno być tutaj najlepszą opcją. Właśnie zacząłem bawić się plikami wektorowymi, więc ten świat jest dla mnie całkiem nowy …
Próbuję przyciągnąć do siebie dwie linie za pomocą Shapely / Geopandas, ale rezultat przyciągania jest bardzo dziwny. Próbowałem : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) res = lines1 lines2_union = …
Szukam rozwiązania python typu open source do konwersji rastra na wielokąt (bez ArcPy). Znałem funkcję GDAL do konwersji rastra na wielokąt, a oto instrukcja: http://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#polygonize-a-raster-band Niemniej jednak oczekuję, że dane wyjściowe mogą być kształtnymi wielokątami lub dowolnym obiektem, tymczasowo w pamięci, a nie zapisanym jako plik. Czy jest jakiś pakiet …
Mam wielokąty ze współrzędnych w (kształtny python), które wyglądają tak POLYGON ((24.8085317 46.8512821, 24.7986952 46.8574619, 24.8088238 46.8664741, 24.8155239 46.8576335, 24.8085317 46.8512821)) Chciałbym obliczyć powierzchnię tego wielokąta w km ^ 2. Jaki byłby najlepszy sposób, aby to zrobić w Pythonie?
Mam geopandy GeoDataFrame zawierające setki zgrabnychPolygon i MultiPolygongeometrycznych kształtów. Wieloboki zachodzą na siebie w wielu miejscach. Chciałbym stworzyć nową geometrię, która zawiera liczbę nakładających się na siebie. Coś takiego: Czy ktoś ma jakieś pomysły, jak do tego podejść? Nie widzę nawet wejścia. W końcu szczególnie chciałbym móc zważyć wielokąty, aby …
Mam problem ze zrozumieniem użycia indeksów przestrzennych w RTree. Przykład: Mam 300 zbuforowanych punktów i muszę znać obszar przecięcia każdego bufora za pomocą wielokąta. Plik kształtu wielokąta ma> 20 000 wielokątów. Zasugerowano użycie wskaźników przestrzennych w celu przyspieszenia procesu. SO ... Jeśli utworzę indeks przestrzenny dla mojego pliku kształtu wielokąta, …
Używam Shapely w pythonie i dostaję MultiLineStringmnóstwo Linestringobiektów. Mogę zagwarantować, że wszystkie LineStringobiekty są prostymi liniami z tylko 2 wierzchołkami i że wszystkie są częścią jednej pojedynczej linii (bez gałęzi). Chcę „połączyć kropki” i stworzyć singiel LineString. Czy muszę do tego napisać metodę spawania rekurencyjnego, czy jest szybszy sposób?
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.