Shapely to pakiet Pythona typu open source do budowy i analizy płaskich obiektów geometrycznych 2D, wspólnych dla systemów informacji geograficznej.
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...
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...
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...
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...
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...
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ść,...
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...
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...
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...
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...
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...
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':...
Próbowałem kilku przykładów kodu wykorzystującego biblioteki takie jak shapefile, fiona i ogr, aby spróbować sprawdzić, czy punkt (x, y) mieści się w granicach wieloboku utworzonego za pomocą ArcMap (a zatem w formacie shapefile). Jednak żaden z przykładów nie działa dobrze z multipoligonami,...
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...
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...
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ść?...
wygląda to na głupie pytanie, ale nie mogę znaleźć na to jednoznacznej odpowiedzi: jakie jednostki geopandy / zgrabnie używają obliczania odległości / powierzchni między
Podczas próby instalacji Shapely na moim 64-bitowym komputerze z systemem Windows nie mogę uruchomić biblioteki GEOS. Do tej pory uruchomiłem instalator OSGeo4W, z którego zainstalowałem GDAL (uważam, że biblioteka geos jest zawarta w tym pakiecie). Potem sprawdziłem i mam geos_c.dll w moim...
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...
[EDYCJA: rozwiązaniem tego było po prostu użycie OGR do odczytu plików kształtów. Zobacz przykład geographika.] W pliku kształtu ESRI nie ma rozróżnienia między wielokątami i wielokątami. Co więcej, nie ma wyraźnego rozróżnienia między wewnętrznymi otworami a zewnętrznymi pierścieniami (poza...