Wziąłem kopię super-nikczemnej książki „Python Geospatial Development” Erika Westry ( link do Amazon ) i pracuję nad tym. Obecnie uczy mnie, jak ładować dane wybrzeża GSHHS z pliku kształtu do bazy danych PostGIS, przygotowując się do budowy geoprzestrzennej aplikacji internetowej. Mój problem...
Próbuję użyć ArcObject w Pythonie. Mam problem z przypisaniem interfejsu (tj. IZ.InterpolateZsBetween) do mojego obiektu z danymi geometrycznymi, kod to: from Comtypes.client import CreateObject, GetModule import arcgisscripting gp =
Chcę przekonwertować tabelę PostGIS na plik shapefile (bez użycia pgsql2shp). Aby stworzyć geometrię w pliku kształtu, muszę podać Xmin, Ymin i Xmax, Ymax, a geometria, którą mam w mojej tabeli PostGIS, ma kształt nieregularny (mogę uzyskać zewnętrzną stronę za pomocą ramki granicznej, ale będzie...
Czy można programowo zatrzymać i uruchomić usługi map ArcGIS Server przy użyciu
Chciałbym zainstalować QGIS. Korzystam z systemu domyślnie Python 2.5x i pobrałem Python 3.7x z Homebrew. Stworzyłem wirtualne środowisko z tym kodem. brew install pyenv brew install pyenv-virtualenv exec $SHELL eval "$(pyenv init -)" eval "$(pyenv virtualenv-init -)" pyenv activate...
Czy można złapać komunikaty o błędach wyjściowych PostgreSQL w OGR po nieudanym zapytaniu? Na przykład conn.ExecuteSQL('SELECT * FROM non_existing_table;') W Postgres zwraca komunikat o błędzie: "Error: relation 'non_existing_table' does not exist." W ogrodach otrzymujemy standard: "Error:...
Czy podczas iteracji po warstwie wektorowej przy użyciu następującego kodu (w skrócie z przykładu w książce kucharskiej pyqgis ) istnieje sposób, aby sprawdzić, czy wybrano funkcję? provider = vlayer.dataProvider() feat = QgsFeature() allAttrs =
Mam problem / pytanie dotyczące Pythona i QGIS. Czy można uzyskać dane jednej wybranej funkcji w QGIS przy użyciu Pythona, a jeśli tak, to w jaki sposób? W moim przykładzie mam warstwę, która pokazuje kilka działek do planowania urbanistycznego. Jeśli teraz wybiorę jedną z paczek, chcę uzyskać...
Próbuję zrobić rysunek przedstawiający dane azymutalne z innym zakresem niepewności w każdym punkcie. Ta oldschoolowa postać z gazety z 1991 roku oddaje ideę „fabuły muszki”, do której dążę: Wszelkie sugestie, w jaki sposób mogę zrobić podobną liczbę? Jestem względnym nowicjuszem w GIS, ale mam...
Czy można ustawić CRS projektu w QGIS za pomocą komendy Python? W moim przypadku użytkownik powinien otworzyć QGIS, uruchomić moją pisemną wtyczkę, a ta wtyczka załaduje dane przestrzenne, które mają kod epsg 31467 (Gauss Krueger Zone 3). Niestety domyślnym CRS QGIS jest WGS84, więc mam kłopoty z...
Mam przydatne narzędzie skryptowe, które przechodzi przez obszar roboczy, zmienia nazwy i kopiuje pliki kształtów do zestawu danych funkcji. Jeśli jednak gdzieś w obszarze roboczym znajduje się uszkodzony plik kształtu, skrypt nie działa i przestaje przetwarzać. Jak radzisz sobie z takimi...
Próbuję zrozumieć, w jaki sposób mogę użyć modułu csv w pythonie, aby otworzyć plik csv w tym samym folderze co skrypt Pythona, a następnie utworzyć plik shapefile przy użyciu pyshp modułu shapefile. Plik csv wygląda tak, ale może zawierać kilka tysięcy wierszy
Mam dużą liczbę plików kształtów reprezentujących obszary zainteresowania do analizy, która zostanie przeprowadzona przy użyciu różnych źródeł zdjęć satelitarnych (IKONOS, RapidEye itp.). Niestety, obrazy nie używają na przykład systemu ścieżek takiego jak Landsat, więc zakres jest bardzo...
Próbuję generować wielokąty dla satelitarnych pokosów orbit. Do tej pory mam metodę generowania dwóch linii, które reprezentują krawędź każdego pokosu w [łac, długo]. Niektóre pokosy przekraczają międzynarodową linię danych i dlatego obejmują: Udało mi się to rozwiązać za pomocą ogr2ogr...
Próbuję nauczyć się korzystać z ogrodu w pythonie przy użyciu zestawu danych o kraju i zaludnionych miejscach z http://www.naturalearthdata.com/downloads/50m-cultural-vectors/. Próbuję użyć filtrów i buforów, aby znaleźć punkty (ne_50m_populated_places.shp) w określonym buforze nazwanego kraju...
Próbuję zbudować wtyczkę, aby załadować kompozytora wydruku z pliku, wygenerować atlas i wyeksportować do obrazu. Do tej pory udało mi się załadować szablon i wyeksportować go do obrazu. Nie byłem w stanie dodać żadnej z warstw legendy (które są również w toc) do eksportowanej mapy, co powoduje,...
Zadanie polegało na utworzeniu modelu za pomocą ArcGIS ModelBuilder, który początkowo uruchamia kilka zadań geoprzetwarzania, a następnie dołącza do tabeli Excela w wyniku tych geoprocesów. Jednak model musi najpierw otworzyć plik Excel, aby użytkownik mógł go ręcznie edytować, zanim następnie...
Chciałbym utworzyć skrypt w języku Python w QGIS, który pobierze pewne dane (shp + tif, w moim przykładzie), użyje tamplate Map Composer (z pliku) i wyeksportuje utworzoną kompozycję do obrazu png. Bez prawie żadnego doświadczenia w programowaniu (mniej niż podstawowa znajomość języka Python)...
Używam algorytmu Birch z scipy-learn pakiet Python do grupowania zestawu punktów w jednym małym mieście w zestawy po 10. Używam następującego kodu: no = len(list_of_points)/10 brc = Birch(branching_factor=50, n_clusters=no, threshold=0.05,compute_labels=True) W moim pomyśle zawsze kończyło mi...
Czy istnieje sposób, aby sprawdzić, czy jakieś 2 podane warstwy rastrowe mają identyczną zawartość ? Mamy problem z wolumenem naszej wspólnej pamięci masowej w firmie: jest on teraz tak duży, że wykonanie pełnej kopii zapasowej zajmuje ponad 3 dni. Wstępne dochodzenie ujawnia, że jednym z...