Czy można wczytać surowe dane do pliku geopandas GeoDataFramea la a la pandas DataFrame? Na przykład następujące prace: import pandas as pd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") pd.read_json(io.BytesIO(r.content)) Następujące nie: import geopandas as gpd import requests data = requests.get("https://data.cityofnewyork.us/api/geospatial/arq3-7z49?method=export&format=GeoJSON") gpd.read_file(io.BytesIO(r.content)) Innymi słowy, czy można odczytać dane geoprzestrzenne znajdujące się w …
Zbudowałem wtyczkę do zasilania bazy danych Spatialite za pomocą maski wprowadzania. Istnieją dwa pola reprezentujące minimalny / maksymalny wymiar obiektu. W momencie, gdy wywołuję narzędzie Measure ( iface.actionMeasure().trigger()) za pomocą dwóch przycisków „miara”, kliknij wymiar i wpisz wartość za pomocą klawiatury w polach (QLineEdit). Czy istnieje możliwość, że po wykonaniu …
Podaję kafelki wektorowe za pomocą TileStache , mam wszystko ustawione tak, jak chcę. Moje dane są przechowywane w Postgres i używam dostawcy VecTiles do obsługi kafelków GeoJSON . Chcę buforować wszystkie moje kafelki, aby kafelki działały szybciej. Korzystam z tilestache-seed.py do zapasu mojej pamięci podręcznej. Używam seed-tilestache na kilku maszynach. …
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, …
Zamknięte . To pytanie musi być bardziej skoncentrowane . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby skupiało się tylko na jednym problemie, edytując ten post . Zamknięte 4 lata temu . Szukam dobrego zasobu (MOOC, książka itp.), Aby uczyć Pythona w aplikacjach do zdalnego wykrywania. Znam …
Chciałbym wywołać funkcję wtyczki interpolacji QGIS (metoda TIN) (Raster-> Interpolate) z konsoli python. Nie mogę znaleźć odpowiedniej funkcji w QGIS API lub na liście algorytmów przetwarzania. Znalazłem algorytm triangulacji SAGA, który działa dobrze, ale jest 5-10 x wolniejszy, a szybkość jest ważna w moim przypadku. Masz pomysł, jak to wykonać?
Jestem w tym nowy. „Buduję” wtyczkę Python dla qgis i chcę importować moduły z GRASS i SAGA, ale nie rozumiem, jak to zrobić. Wtyczka, którą chcę zaimportować, to interpolacja Kriging z SAGA lub v.surf.idw z GRASS lub wtyczka interpolacyjna z Raster qgis. Chcę, aby metoda interpolacji przekształciła punkty w siatkę. …
Próbuję zmienić wtyczkę z QGIS 1.8 na 2.0. W wersji 1.8 tworzę tabelę nazw pól takich jak to: layer = vlayer provider = layer.dataProvider() feat = QGSFeatrue() allAttrs = provider.attributeIndexes() provider.select(allAttrs) fieldNames = [] while provider.nextFeature(feat) fields = layer.pendingFields() for (k, field) in fields.iteritems(): fieldNames.append(field.name()) Próbowałem tego dla wersji 2.0 …
Zobacz ten link, aby uzyskać więcej informacji. Problem: Chcę przejść przez ciągły raster (taki, który nie ma tabeli atrybutów), komórka po komórce i uzyskać wartość komórki. Chcę wziąć te wartości i uruchomić na nich warunki warunkowe, emulując kroki algebry mapy wyszczególnione poniżej, bez korzystania z kalkulatora rastrowego. Na prośbę o …
Udało mi się dodać QgsVectorLayerw QGISpośrednictwem wtyczki Python, na przykład: vl = self.iface.addVectorLayer(uri.uri(), layerName, self.dbConn.getProviderName()) Ale jak mogę to usunąć?
Byłem bardzo niezadowolony z obliczania długości pościeli w WGS84 w Milach . Zastanawiałem się, czy istnieje wygodniejszy, Pythoński sposób obliczania długości linii WKT zgodnie z danym SRID. Mam na myśli coś takiego: srid="WGS84" line="LINESTRING(3.0 4.0, 3.1 4.1)" print length(line, srid) Szukam dokładnej odpowiedzi, a nie sin\cosprzybliżeń. Jakieś pomysły?
Skrypt, który piszę, umieszcza dwie ramki danych i ustawia ich zakres. Robiąc to, przerysowuje cały widok Active 4 razy, znacznie spowalniając skrypt. Naciśnięcie klawisza F9 lub kliknięcie przycisku „Wstrzymaj rysowanie” przed uruchomieniem skryptu wyłącza rysunek i powoduje, że skrypt działa znacznie szybciej, ale chcę, aby skrypt zrobił to automatycznie. Chcę …
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?
W dokumentacji gdal_calc podano kalkulator rastrowy wiersza poleceń ze składnią numpy . Później jest kilka przykładów, w których w jednym z nich: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - oznacza ustawione wartości od zera i poniżej na zero Niestety nie ma …
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.