Okej, więc mam trochę dość różnych rozwiązań rastrowych w bazie danych, które wypróbowałem przez ostatnie 2 tygodnie, które obecnie nie działają (raster PostGIS) lub nie mogę poradzić sobie z ilością danych (Rasterlite), więc Czytam o formacie VRT. Miejcie na uwadze, chciałbym wasze opinie na ten...
Próbuję zainstalować GDAL 2.X na CentOS7 bez budowania ze źródła, ponieważ będę musiał skonfigurować go na innych serwerach, a instalacja według kodu źródłowego będzie bardzo męczącym i czasochłonnym zadaniem. Znalazłem niektóre repozytoria, takie jak opengeo i ELGIS, ale dotyczą one CentOS7. W...
Krótko: chcę zmienić rozdzielczość rastra i wygładzić szarą rampę, jak pokazano na poniższych zdjęciach. Preferowane jest użycie GDAL, PIL lub Numpy. Opis: Za pomocą High Performance Geostatistic Library dzielę punkty na rastry o rozdzielczości wyjściowej 20 metrów . Nie chcę zmieniać...
Jestem wyjątkowo nowy w GIS. Używam gdaldo czytania na mapie użytkowania terenu / pokrycia terenu i muszę wybrać szerokość / długość niektórych rodzajów pokrycia terenu, aby zindeksować do innego zestawu danych, który jest wyrażony tylko w szerokości / długości. Niestety nie rozumiem formy...
Łączę niektóre pliki kształtów i miałem pewne problemy w QGIS, więc używam bezpośrednio ogr2ogr. Robię to (partią): ogr2ogr -overwrite %destination% %n1% ogr2ogr -update -append %destination% %n2% -nln all_new ogr2ogr -update -append %destination% %n3% -nln all_new ogr2ogr -update -append...
Mam plik obrazu tif i plik świata .tfw (utworzony z danych OSM za pomocą Mapertive). Jak mogę (za pomocą wiersza poleceń systemu Linux) przekonwertować to na plik GeoTIFF? (tzn. chcę mieć jeden plik, a nie 2). Jestem pewien, że jest to możliwe z Gdalem czy coś
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...
Regularnie tworzę własne rastry GeoTIFF z GDAL w Pythonie, np .: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band =...
Właśnie eksportuję tabelę Postgis do shp, korzystając z tych wskazówek, ale nie jestem w stanie zaimportować shp do Postgis przy użyciu tej samej biblioteki (ogr). Dowolny pomysł? Wielkie
Piszę proste narzędzie do przycinania partii wielopasmowych plików rastrowych geotiff do tego samego (mniejszego) obszaru. Korzystając z gdalwarp, mogę z łatwością przyciąć plik za pomocą pliku kształtowego przycinającego jeden wielokąt: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif...
Próbuję przetwarzać niektóre dane rastrowe przy użyciu ogr / gdal i nie mogę uzyskać pełnego wykorzystania wszystkich rdzeni na moim komputerze. Kiedy uruchamiam ten proces tylko na jednym rdzeniu, uzyskuję 100% wykorzystania tego rdzenia. Kiedy próbuję podzielić się na wielordzeniowy (w poniższym...
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...
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...
Zainstalowałem GDAL i sprawdziłem w wierszu poleceń, działa, a kiedy uruchamiam gdal2tiles.py w wierszu poleceń, działa, ale kiedy uruchamiam ten sam plik gdal2tiles.py z python IDLE, otrzymuję komunikat „Brak modułu o nazwie osgeo” błąd „”. Czy możesz
Używam Pythona 2.6.5 (32-bitowego) z Numpy 1.3 i Gdal 1.9.1 zainstalowanym na Windows 7 64-bitowym. Próbuję odczytać zestaw danych rastrowych Imagine (.img) o pojemności 800 MB do tablicy Numpy, aby wykonać algebrę rastrową, ale gdy tylko uruchomię następujący kod, Python.exe ulega awarii. from...
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...
Mam dane rastrowe w tabeli postgres, którą chcę wprowadzić do Pythona jako tablicę numpy. Korzystam z psycopg2, aby połączyć się z bazą danych. Jestem w stanie pobrać dane, ale są one zwracane jako ciąg (prawdopodobnie serializowany plik binarny). Czy ktoś wie, jak pobrać ten ciąg i...
Chcę utworzyć dwa wielokąty. Jeden jest z prostokątnych rozmiarów mojego rastra. Wiem, że mogę uzyskać zasięg użycia rastra, r <- raster("band5.tif") e <- extent(r) plot(e)daje mi to: Jak mogę utworzyć JEDEN wielokąt granicy rastra, jak pokazano poniżej? Tak wygląda mój...
Jestem przeważnie wektorem, ale przy bieżącym projekcie muszę pracować z kilkoma starymi zeskanowanymi mapami papierowymi (mapy szkód bombowych z Londynu, jeśli jesteś zainteresowany!) Przeskanowaliśmy i dokonaliśmy georeferencji map, a teraz chcemy utworzyć kafelkową warstwę kompozytową do...