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 obiektami?
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, chociaż dobrze sobie radzą z normalnymi, pojedynczymi wielokątami …
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 C:\OSGeo4W\binkatalogu, ale albo przegapiłem kilka kroków konfiguracji, albo …
Chciałbym uzyskać przecięcie wielu wielokątów. Korzystając z shapelypakietu Pythona , mogę znaleźć przecięcie dwóch wielokątów za pomocą intersectionfunkcji. Czy istnieje podobna skuteczna funkcja uzyskiwania przecięcia wielu wielokątów? Oto fragment kodu, aby zrozumieć, co mam na myśli: from shapely.geometry import Point coord1 = ( 0,0 ) point1 = Point(coord1) circle1 = …
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 współrzędnych, aby utworzyć Shapely MultiPolygon, ponieważ wymaga a sequence of exterior ring …
[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 „przydatnością” danego wielokąta). Po przeczytaniu pliku kształtu mam listę sekwencji współrzędnych …
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 spełniony jest jakiś warunek. W tej …
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 drodze, otrzymując coś takiego: sql = """ SELECT ST_AsText(ST_Transform(the_geom,4326)) as newgeom,* …
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, ale nie składa się z instancji Point. Więc jeśli mam dwa punkty Ai Bczy …
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 = LineString([(0.9, 0.9), (0.2, 0.6)]) poly_gdf = gpd.GeoDataFrame(geometry=[polygon]) …
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 atrybucie długości foremnej. Jaka jest jednostka atrybutu długości? Czy istnieje …
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ś innego?
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 -wrapdateline: ogr2ogr -wrapdateline -f "ESRI Shapefile" test.shp orbits.shp Co prawdopodobnie …
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, w rzeczywistości wielokąt jest o wiele bardziej złożony i zawiera …
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ą python.
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.