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...
Jak mogę wybrać funkcje w PyQGIS za pomocą wyrażenia? Próbowałem użyć QgsExpressionmetody, ale wybrana metoda jej nie przyjmuje: exp = QgsExpression("'ogc_fid' = 482") cLayer = canvas.currentLayer() cLayer.select(exp) Czy to możliwe, a jeśli tak, jak to
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
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ść?...
Szukam metody obliczenia współrzędnej X, Y punktu wewnętrznego (znanego również jako labelPoint w arcmap) (nieregularnego) wielokąta. Próbuję znaleźć rozwiązanie bez użycia modułu ArcGIS różnica między „centroidem” a etykietą polega na tym, że centroid może znajdować się poza wielokątem. Używam...
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...
Czy istnieje sposób, w jaki po uruchomieniu niestandardowej wtyczki stworzonej dla QGIS, powłoka Python Shell może pojawiać się po uruchomieniu wtyczki, jak ekran debugowania? Proszę o to, aby wiedzieć, czy istnieje sposób, aby zobaczyć polecenia drukowania za pomocą skryptu Python podczas...
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...
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...
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 shapefile poddziałów krajowych i chciałbym wyodrębnić tablicę ograniczających szerokości i długości geograficznej dla każdej dywizji .. czy można to
Mam raster, z którym chciałbym wykonać interpolację punktową. Oto gdzie jestem: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array =
Chciałbym zbudować plik kształtu z ramki danych Pandas przy użyciu wierszy lon i lat. Mam plik csv i przetwarzam go za pomocą pand, aby utworzyć ramkę danych, która jest łatwiejsza w obsłudze Czy można to zrobić bez tworzenia pętli linia po
Jakiś czas temu przeprowadziłem ankietę dotyczącą orientacji dachu w sąsiedztwie dotyczącą rozmieszczenia paneli słonecznych, wizualnie sprawdzając każdy dach i zaznaczając indywidualnie, czy są one skierowane na południe, czy prawie na południe. Ponieważ jest to bardzo czasochłonne, chciałbym...
Uczyłem się geoprzetwarzania z Pythonem na kursie MSOG GEOG 825 (geoprzetwarzanie pytona). Niestety link już nie działa. Wiem, że to długa szansa, ale czy ktoś ma kopię materiałów z tego kursu? Naprawdę chciałbym to zakończyć. Sprawdziłem już inne internetowe kursy geoprzetwarzania Pythona....
Mam problem z użyciem funkcji Polygonize w Pythonie. Przykład książki kucharskiej na ten temat można znaleźć tutaj . Odpowiednia część mojego kodu to: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver =...
Chciałbym zaktualizować jeden atrybut funkcji. Nie wiem jednak, jak go zaktualizować za pomocą obiektu funkcji. Muszę użyć dostawcy danych, aby go zaktualizować. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it =...
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...
Mam siatkę binarną Arc / Info --- konkretnie raster akumulacji przepływu ArcGIS --- i chciałbym zidentyfikować wszystkie komórki mające określoną wartość (lub zakres wartości). Ostatecznie chciałbym uzyskać plik kształtu reprezentujący te komórki. Mogę użyć QGIS, aby otworzyć plik hdr.adf i...
Czy istnieje sposób wykonania tego samego zadania, co narzędzie gdalbuildvrt przy użyciu powiązań GDAL Python? Do tej pory nie znalazłem innego sposobu niż utworzenie vrt pojedynczego zestawu danych i ręczna edycja pliku XML. Chciałbym utworzyć vrt z wielu rastrów (zasadniczo wykonujących mozaikę)....