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...
Nowicjusz tutaj zmaga się z GIS. Próbuję wytyczyć totemy dla miasta Milwuakee, używając plików kształtów znalezionych na stronie internetowej hrabstwa . Śledzę ten wątek tutaj z pewnym sukcesem. Mój kod daje: from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire...
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...
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...
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...
Korzystam z poniższego wiersza, aby poinformować użytkownika o statusie: iface.mainWindow().statusBar().showMessage("Status:" + str(i)) Wtyczka działa na moim zestawie danych około 2 minut, ale system Windows wykrywa go jako „nie odpowiada” i przestaje wyświetlać aktualizacje stanu. Dla nowego...
Chcę utworzyć funkcję, która wybiera funkcję i przybliża ją (podobnie w QGIS). Dlatego istnieje następująca funkcja: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) def zoomTo(self): layer = self.vlayer atable = self.ui.table selectList=[] for i in
Próbuję utworzyć plik kształtów wielokątów z listy współrzędnych za pomocą narzędzi Python typu open source. Oto, co mam do tej pory, które zostało zhakowane razem z książki kucharskiej Pythona GDAL / OGR i odpowiedzi GIS SE . Istnieje podobne pytanie Python: Jak utworzyć plik kształtu wielokąta z...
Czy można zmienić kolor podczas dodawania warstwy wektorowej? Mam warstwę wektorową z obszarami wodnymi i chcę załadować ją w stylu koloru niebieskiego. Warstwa jest dodawana z: QgsMapLayerRegistry.instance().addMapLayer(self.vlayer) Czy mogę w jakiś sposób zmienić kolor po lub podczas...
Znalazłem dane o wysokiej precyzji wysokości (rozdzielczość 1/3 - 1/9 sekundy łukowej) z krajowego zestawu danych wysokości dostarczonego przez USGS. Pochodzi w formatach IMG, Arcgrid i GridFloat. Wiem, że ArcGrid odpowiada płatnemu pakietowi oprogramowania, ale staram się trzymać swobodnie...
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...
Czy istnieje sposób na uzyskanie listy zainstalowanych formatów gdal w systemie z pythonem? W wierszu poleceń można to łatwo zrobić za pomocą: gdalinfo --formats ... ale czy można to zrobić w
Chcę stworzyć prostą funkcję if / else w moim Modelu i chciałem wiedzieć, czy jest to możliwe, czy nie. Chcę, aby to zrobić, jeśli wybrano dane wejściowe dla „Segment Layer”, że nie musi ono wykonywać „Segmentation (średnich zmian)” i zamiast tego używa „Segment Layer” jako danych wejściowych...
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...
Chcę georeferencję rastra za pomocą pythoni GDAL. Moje obecne podejście polega na dzwonieniu gdal_translatei gdalwarpkorzystaniu os.systemz brzydkiej listy punktów kontroli naziemnej. Naprawdę chciałbym zrobić to w sposób natywny python. To jest obecny proces, którego używam: import...
Jak dotąd korzystasz z ArcObjects z Pythona? to moje najczęściej czytane i wymieniane pytania i odpowiedzi na temat wymiany stosów GIS. Mimo tego sukcesu jest to prawdopodobnie jeden z moich najsłabszych obszarów, jeśli chodzi o faktyczne wykorzystanie. Duża część tego kiepskiego pokazu wynika z...
Jak mogę uzyskać rzutowane współrzędne, a także rzeczywiste wartości pikseli na tych współrzędnych z pliku GeoTiff, a następnie zapisać je w tablicy liczbowej? Mam plik arsenci020l.tif, a jego współrzędne są w metrach. Poniżej znajduje się skrócona wersja gdalinfo, którą na nim uruchomiłem. ~$...
Zastanawiam się, czy istnieje sposób, aby otworzyć konsolę Python bezpośrednio przy uruchomieniu programu. Istnieje możliwość ustawienia skrótu do konsoli Python, ale nie mogę znaleźć takiej opcji do uruchomienia programu
Mam zestaw linii reprezentujących linie autobusowe. Niektóre linie nachodzą na siebie i biegną tymi samymi drogami. Jestem w stanie wyodrębnić węzły. Jednak interesuje mnie wyodrębnianie tylko takich przejść: W jaki sposób mogę to zrobić? Szukam sposobów na QGIS lub Python. Próbowałem...