Być może to trochę dziwne pytanie, ale pozwólcie, że wyjaśnię to krótko, zanim zacznę moje rzeczywiste pytania: Mikroskopia sił atomowych (AFM) to metoda, która w skrócie (i o mojej ograniczonej wiedzy) pozwala badaczom skanować obszary w mikro- i nanoskali. Działa poprzez „skanowanie” obszaru za...
Mam duży plik rastrowy (129600 na 64800 pikseli) z globalnymi częściami wód (1 bit wartości 0 i 1) i próbuję wydobyć linie brzegowe oceanów i wód śródlądowych. Próbowałem z ArcGIS i QGIS przekonwertować z rastra na polilinię, ale to trwa wieki. Czy ktoś zna lepszy / szybszy sposób (Python lub R)...
Odwiedziłem i ponownie odwiedziłem stronę dotyczącą obsługi geometrii w książce kucharskiej PyQGIS: http://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html, ale wydaje się, że nie mogę dowiedzieć się, jak uzyskać wielokąt do rysowania z konsoli Python. Czy ktoś może...
Aby uzyskać wartość pojedynczego piksela w punkcie rastra za pomocą rasterio, istnieje przykład: https://github.com/mapbox/rasterio/pull/275 Czy jednak istnieje bezpośredni interfejs API w rasterio (a nie cli), którego można użyć do wyodrębnienia wartości w jednym punkcie rastra? --...
Mam liczbę całkowitą rastrową, dla której chciałbym zbudować tabelę atrybutów rastrowych przy użyciu Pythona i GDAL. Mogę utworzyć tabelę atrybutów rastrowych GDAL w Pythonie w następujący sposób: >>> rat = gdal.RasterAttributeTable() Działa to dobrze, jak widzimy: >>>...
Próbuję zoptymalizować mój kod python za pomocą pypy . Chociaż jestem w stanie zaimportować większość bibliotek i uruchomić się pomyślnie, nie jestem w stanie zaimportować bibliotek GDAL. Te same biblioteki po zaimportowaniu do Pythona 2.7 działają poprawnie. Próbowałem dodać lokalizację GDAL do...
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...
Chcę opracować GUI w Pythonie do geoprzetwarzania ArcGIS. Czy mogę używać PyQT do programowania GUI w ArcGIS? Rozważałem także użycie Tkintera i wxPython.
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...
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 =...
mam GeoPanda 0.2.1 py27_0 fiona 1.7.0 np110py27_1 Python 2.7.10 zainstalowany w środowisku Anaconda 2-4.1.1-Windows-x86. Potrafię budować GeoDataFrame, odczytując wejściowe zestawy danych i operując danymi, ale zapisanie wyjściowego zestawu danych nie zachowuje układu współrzędnych. import...
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...
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ść...
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
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...
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...
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ść?...
Było mnóstwo pytań na ten temat, ale żadne z nich nie rozwiązało mojego problemu. Mam maszynę z: Windows 7 x64 Python 3.4.3 x64 gdal zainstalowane: gdalinfo --version --> GDAL 1.11.4, released 2016/01/25; gdal został zainstalowany z kołem GDAL-1.11.4-cp34-none-win_amd64.whlod tutaj GDAL...
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...
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