Próbuję zmienić wartość pola pliku kształtu. Nie zmienia to jednak wartości. Co jest nie tak z moim kodem? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature = layer.GetNextFeature() dist = 233 while...
W Pythonie, używając GDAL, wyodrębniłem projekcję rastra jako ciąg WKT w następujący sposób: wkt = dataset.GetProjection() # wkt is
Mam plik kształtu z wielokątami. I mam globalny plik rastrowy. Chcę nałożyć wieloboki pliku kształtu na siatkę rastrową i obliczyć średnią wartość rastra dla każdego wielokąta. Jak mogę to zrobić za pomocą GDAL, zapisując wyniki w pliku
Przygotowuję dane do włączenia do węzła geonetwork, w którym mapy są wyświetlane za pośrednictwem geoserver. Problem polega na tym, że komórki, które nie mają danych, są wyświetlane jako geotiff w kolorze czarnym, a ja chciałbym, aby kolor był ustawiony na biały lub przezroczysty. Korzystam z...
Próbuję przetworzyć niektóre bioklimatyczne pliki rastrowe, takie jak te, które można pobrać ze strony http://www.worldclim.org/current (zestaw bioklimatyczny). Wydaje się, że mają ustawione wartości nodata -3.4e+38zgodnie z QGIS (patrząc na wynik gdalinfo, to -3.39999999999999996e+38). Wydaje...
W kodzie Python GDAL dość często widzę, że ludzie zamykają zestawy danych na końcu skryptu. Dlaczego warto zamykać zestaw danych w Pythonie GDAL? Czy są jakieś konsekwencje, jeśli tego nie zrobię? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds =...
Ponownie projektuję rastry w Pythonie za pomocą GDAL. Muszę rzutować kilka tiffów ze współrzędnych geograficznych WGS 84 na WGS 1984 Web Mercator (Sfera Pomocnicza), aby użyć ich później w Openlayers razem z OpenStreetMap i być może mapami Google. Korzystam z Pythona 2.7.5 i GDAL 1.10.1 stąd i...
Przepraszam, jeśli poniższe pytanie jest nieco głupie, ale jestem BARDZO nowy w tej całej sprawie GIS. Próbuję przekonwertować niektóre rzutowane obrazy geoTiff na WGS84 przy użyciu gdal w python. Znalazłem post opisujący proces przekształcania punktów w rzutowanych GeoTiff przy użyciu czegoś...
Czy ktoś wie, jak najłatwiej zainstalować GDAL 2.0 na MacOSX? Zwykła podejrzana dystrybucja binarna kingchaos i homebrew wydaje się utknąć na GDAL 1.11 Potrzebuję wsparcia dla właściwości funkcji TopoJSON, które wydają się być zawarte w 2.0 EDYCJA 2017/01/07: Dystrybucja binarna kingchaos...
Mam obraz .tiff, który próbuję nałożyć na Google Earth. W tym celu potrzebuję najpierw georeferencji. Projekcja obrazu to UTM z NAD83 jako punktem odniesienia. Znam długość / długość czterech rogów obrazu. Robię w następujący sposób: gdalinfo, aby uzyskać cztery współrzędne. gdal_translate, aby...
Obecnie pracuję nad importem danych klimatycznych CANGRID (dostarczanych jako Surfer Grid ascii, pliki „.grd”) do ArcGIS. Siatka ma 95 rzędów na 125 kolumn. Metadane zapewniają lat / lon pochodzenia (lewy dolny róg), rozmiar komórki (50 km), a także odwzorowanie notatek jako stereograficzne...
Właśnie otrzymaliśmy duży zestaw DEM w pracy i chciałbym z nich wygenerować kontury. DEM mają rozdzielczość 1m i rozmiar 1km x 1km. Dane wyjściowe z gdalinfo: Driver: AAIGrid/Arc/Info ASCII Grid Files: 380000_6888000_1k_1m_DEM_ESRI.asc Size is 1000, 1000 Coordinate System is `' Origin =...
Jako przykład wezmę następujący kafelek http://a.tile.openstreetmap.org/3/4/2.png i zapiszę go jako „4_2.png”. Współrzędne WGS84 tej płytki można obliczyć lub czytać tam klikając na odpowiedni żeton: 0 66.51326044311185 45 40.97989806962013 (West North East South) Jak poprawnie dokonać...
Mam wiele źródeł kafelków wygenerowanych za pomocą gdal2tiles, które chciałbym przedstawić na tej samej mapie. Dlatego podczas prezentacji kafelka muszę określić, z którego źródła go podawać, sprawdzając, czy nie ma granic. Czy ktoś wie, jak obliczyć granice długości geograficznej dla pojedynczego...
Próbuję zaimportować gdal dla Python 2.7.3 na Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Python27\lib\site-packages\gdal.py", line 2, in <module> from osgeo.gdal import deprecation_warn File...
Chcę stworzyć portal internetowy, który będzie wyświetlał obrazy dynamicznie tworzone przez QGIS, ale dostarczone przez QGIS obrazy tiff są zbyt duże i nie można ich pokazać na stronie za pomocą przeglądarki zdjęć. Właśnie dlatego zdecydowałem się przekonwertować te pliki tiff na jpg / jpeg, aby...
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ć...
czy możliwe jest podzielenie pliku shapefile na funkcję w Pythonie? (najlepiej byłoby rozwiązaniem, w którym mogę tymczasowo zapisać wynikowe obiekty wektorowe w pamięci zamiast na dysku). Powód: Chcę użyć funkcji gdal rasterizeLayer z kilkoma różnymi podzbiorami pliku kształtu. Ta funkcja wymaga...
Wydaje się, że powinno to być proste, ale nie mogę tego zrozumieć w życiu. Biorąc pod uwagę istniejący GeoTiff zawierający jeden zespół, jak możesz dodać kolejny zespół do tego GeoTiff? Czy jest to jedyny sposób, aby to zrobić przy użyciu GDALCreate (), aby utworzyć nowy plik wyjściowy z...
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...