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 = None
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 przekształcam współrzędne za …
Mam prostą ramkę danych GeoPandas: Chciałbym załadować tę GeoDataframe do tabeli PostGIS. Mam już konfigurację bazy danych z rozszerzeniem PostGIS, ale wydaje się, że nie mogę dodać tej ramki danych jako tabeli. Próbowałem następujące: engine = <> meta = MetaData(engine) eld_test = Table('eld_test', meta, Column('id', Integer, primary_key=True), Column('key_comb_drvr', Text), Column('geometry', …
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 biegunowe z południkiem środkowym (110 stopni W) i …
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 9 miesięcy temu . Wyłączając ArcPy, czy istnieje biblioteka Pythona, która może wykonywać geoprzetwarzanie, takie jak bufor / przecinanie, …
Chciałbym nauczyć się korzystać z tablic NumPy w celu optymalizacji geoprzetwarzania. Duża część mojej pracy dotyczy „dużych zbiorów danych”, w których geoprzetwarzanie często zajmuje kilka dni, aby zrealizować określone zadania. Nie trzeba dodawać, że jestem bardzo zainteresowany optymalizacją tych procedur. ArcGIS 10.1 ma wiele funkcji NumPy, do których można uzyskać …
Chcę programowo utworzyć plik GeoJSON przy użyciu wielokątów z pliku kształtu, ale dodając atrybuty z własnej aplikacji. Można to łatwo zrobić dla pliku kształtu: def create_data_dayer(self,varlist, data): """ Creates a new shape to contain data about nodes. varlist is the list of fields names associated with the nodes. data is …
Mam oprogramowanie (nie Arc), które działa co noc na nieprodukcyjnej maszynie, która aktualizuje funkcje pliku kształtu z zewnętrznej bazy danych. Okresowo chciałbym usunąć wszystkie funkcje z pliku shapefile (nie sam plik, który musi pozostać) i pozwolić programowi „odbudować” plik shapefile od zera. Chciałbym zautomatyzować ten proces. Nie mam obecnie zainstalowanego …
Pracuję z dodatkiem ArcMap w C #. Z kodu C # wykonałem kilka skryptów Python. Teraz, aby uruchomić ten skrypt, mam zakodowaną ścieżkę do Pythona. Ale to nie jest przenośne. Tak więc chcę uzyskać ścieżkę do pliku wykonywalnego Python z kodu i użyć go. Pytanie: Jak mogę uzyskać ścieżkę do …
Jestem naprawdę nowy w Pythonie i chciałbym wiedzieć, czy istnieje szybka metoda uzyskiwania wartości komórek piksela rastrowego piksel po pikselu i współrzędnych (współrzędna mapy XY środka każdego piksela) za pomocą Pythona w ArcGIS 10? Aby opisać to dokładniej, muszę uzyskać mapę X, mapę Y i wartość komórki pierwszego piksela, przypisać …
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ś podobnego do następującego: from osgeo import osr, gdal # …
Biorąc pod uwagę geopandas GeoDataFrameserię wielokątów, chciałbym uzyskać powierzchnię w km2 każdej funkcji na mojej liście. Jest to dość powszechny problem, a zwykle proponowane rozwiązanie w przeszłości był w użyciu shapelyi pyprojbezpośrednio (np tutaj i tutaj ). Czy istnieje sposób, aby to zrobić w czystej postaci geopandas?
Korzystam z danych GIS mojego hrabstwa do aktualizacji dróg w OpenStreetMap. Mam kopię, którą opublikowali w zeszłym roku, i jedną, która właśnie została opublikowana, i chciałbym znaleźć wszystkie LineStrings, które miały zmienione atrybuty lub geometrię. Pomoże mi to upewnić się, że wszystkie nowe i zmodyfikowane ulice są aktualizowane w OpenStreetMap. …
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ć rozdzielczości wyjściowej, ponieważ czas interpolacji rośnie wykładniczo. Przy tej …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.