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 =...
Przejście z listy mailingowej gdal-dev: W poniedziałek, 2 września 2013 o 19:09, David Shean napisał: Cześć lista, próbuję spakować rastry GTiff o identycznej projekcji / zasięgu / rozdzielczości jak pojedynczy plik NetCDF do dystrybucji. Spędziłem ostatnią godzinę konsultując się z internetowym...
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...
Mam warstwę buforową (zielony wielokąt), którą chcę podzielić na dwa wielokąty, gdy tylko przekroczy barierę (niebieska linia). Próbowałem użyć metody „splitGeometry”, ale po prostu nie mogę jej uruchomić. Mój kod do tej pory jest następujący: while ldbuffprovider.nextFeature(feat): while...
Jaka jest najlepsza metoda obliczania odległości od najbliższej drogi dla szeregu punktów adresowych w QGIS, posiadających bardzo duży zestaw danych? (QGIS 2.0.1, Win XP) Czytam odpowiedzi na poprzednie podobne pytania: Mierzenie odległości między liniami i punktami w QGIS Rysowanie linii...
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...
Mam punkty w WGS84 lat / long i chciałbym zmierzyć „małe” (mniej niż powiedzmy 5 km) odległości między nimi. Mogę użyć formuły haversine z http://www.movable-type.co.uk/scripts/latlong.html i działa bardzo dobrze. Chciałbym jednak użyć bibliotek Python Shapely, aby móc wykonywać więcej operacji...
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ś
Zamknięte. To pytanie jest nie na temat . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby było na temat wymiany stosów systemów informacji geograficznej. Zamknięte 3 lata temu . Jestem całkiem nowy w Pythonie i...
Tworzę niektóre zmiennoprzecinkowe rastry w GDAL ze sterownikiem GeoTIFF. Kiedy ładuję powstałe obrazy do QGIS lub Arc, domyślna symbolika ustawia min na -3,40282e38 i max na 3,40282e38, więc raster wygląda dokładnie na szary. Czy istnieje sposób, w jaki mogę zapisać rzeczywisty zakres bezpośrednio...
Piszę metodę Python, która dodaje dwa rastry i generuje jeden wynik rastrowy. Z przyczyn niezależnych ode mnie zakresy rastrów wejściowych są różne, ale się pokrywają. Czy możliwe jest działanie wyłącznie na wejściowych pikselach rastrowych, które nakładają się na 2 zachodzące na siebie regiony,...
Mamy duże (ponad 10 000) repozytorium danych przestrzennych (wektory shp / gdb / tab i rastry ecw / tif, a także dane las, xyz i DEM), w których niektóre mają rekordy metadanych w formacie XML przy użyciu anzlic edytora metadanych. Chcę sprawdzić, czy do każdego zestawu danych przestrzennych jest...
Próbuję nauczyć się obsługi przetwarzania obrazu za pomocą teledetekcji za pomocą powiązań i numpy w Pythonie GDAL. Jako pierwszą próbę czytam plik geotiff Landsat8, wykonuję prostą manipulację i zapisuję wynik w nowym pliku. Poniższy kod wydaje się działać dobrze, z tym wyjątkiem, że oryginalny...
Na stronie ftp znajduje się geobaza plików, którą chciałbym pobrać za pomocą skryptu Python. Obecnie myślę, że jednym ze sposobów jest skopiowanie geobazy ftp do geobazy na moim komputerze. Poniżej znajduje się skrypt, który rozpocząłem. Czy ktoś wie, jak mogę zmienić ten skrypt, aby uzyskać ftp...
Wdrażam narzędzie, które pozwala użytkownikowi utworzyć profil funkcji (w moim przypadku działki miejskie). Czy możliwe jest automatyczne wykonanie zrzutu ekranu tej funkcji (wyświetlanego w oknie QGIS) i wyświetlenie go w profilu PDF? Zdjęcie powinno pokazywać jedynie zasięg funkcji i być może...
Biorąc pod uwagę obwiednię, jak mogę policzyć liczbę określonych wartości (powiedzmy, jestem zainteresowany liczbą value == 1) w rastrze w 1) arcpy lub 2) pakiecie rastrowym w
Za pomocą następującego kodu mogę ustawić 100% przezroczystości na wartość rastrową 0: map=None for lyr in QgsMapLayerRegistry.instance().mapLayers().values(): if lyr.name() == "some_layer": map=
Po utworzeniu warstwy, jak mogę ją ukryć / pokazać? Mogę włączyć / wyłączyć renderowanie konkretnej warstwy, zaznaczając pole wyboru za pomocą QGIS, ale muszę to zrobić programowo z kodu python. Jak mogę pokazać / ukryć (nie usuwać) etykietę z kodu python? Szukam czegoś takiego: aLayer =...
Jaki jest właściwy sposób wstawienia punktu do bazy danych PostGIS przy użyciu
Chcę przestudiować matematykę, aby zobaczyć, jak daną szerokość i długość geograficzną w stopniach dziesiętnych można przekonwertować na EPSG: 3857. Czy ktoś może wskazać dobre odniesienie lub możliwy kod źródłowy biblioteki Ppython / C / C ++ / Java do sprawdzania...