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...
Używam GDAL z pythonem do manipulowania plikami kształtów. Podczas otwierania pliku źródłowego najpierw ustawiam sterownik na „ESRI Shapefile”, zanim wywołam metodę „Open” sterownika z wymaganym plikiem shapefile. Nie miałem z tym żadnych problemów; jestem jednak ciekawy, jaka jest różnica między...
Chciałbym uruchomić skrypt Pythona, qgis --code myscript.pya następnie natychmiast wyjść. Używam iface.actionExit().trigger(), który zabija QGIS, gdy uruchamiam go z konsoli Pythona, ale nie po włożeniu do skryptu, który przekazuję --code. Jaki jest właściwy sposób natychmiastowego wyjścia?...
Mam geodataframe z punktami i niektórymi powiązanymi danymi. Chcę wykreślić go na mapie przy użyciu geopand i mieć rozmiar punktów odpowiadający jednej z kolumn w ramce geodata. Do tej pory mam następujący kod: base = world.plot(color='white', figsize=(20,10)) geo_df.plot(ax=base, marker='.',...
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)...
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...
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...
Chciałbym zapisać zależną od projektu konfigurację mojej wtyczki Python w pliku projektu QGIS. Czy to możliwe i jakie są najlepsze praktyki, aby to
Istnieją 3 prawdopodobne scenariusze, które próbuję uchwycić w bliskiej odległości dla: Stacja metra wymiany, która ma 2 lub więcej stacji sąsiednich. Oznacza to, że dana stacja łączy 2 lub więcej głównych tras i ma 2 lub więcej stacji sąsiednich. Terminalowa stacja metra, która ma tylko 1 stację...
Piszę aplikację internetową intensywnie wykorzystującą dane, która jest dostarczana za pośrednictwem apache. Moje pytanie dotyczy tego, jak najlepiej zorganizować przetwarzanie, biorąc pod uwagę, że istnieje wiele opcji. Mam do dyspozycji OpenLayers / JQuery / JavaScript, PostGIS / Postgresql (z...
Obserwowałem niezwykłe wyniki za pomocą skryptu geoprzetwarzania Pythona. (Załączony) skrypt wykonuje następujące działania: Użyj kursora wyszukiwania, aby wyszukać strefę UTM odpowiadającą cechom wielokąta Utwórz przestrzenny obiekt odniesienia na podstawie wyników kursora wyszukiwania Konwertuj...
Mam do czynienia z różnorodnymi przestrzennymi plikami wielokątów (pliki kształtów, klasy plików / osobiste klasy gdb, klasy cech sde (niektóre z widokami przestrzennymi)) iw zależności od typu pliku otrzymuję wiele nazw pól cech klasy kształtu ( wymagane pole, a nie utworzone) Na przykład nazwy...
Korzystając z konsoli python QGIS 2.4, próbuję wyświetlić tylko podzbiór funkcji warstwy. Udało mi się uzyskać funkcje zarówno z warstwy: featuresIterator = layer.getFeatures(QgsFeatureRequest().setFilterExpression ( u'"field"=4')) i od dataProvider: featuresIterator =...
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 2 lata temu . Buduję skrypt w...
Jestem raczej nowy w świecie GIS, a zwłaszcza PostGIS, więc przepraszam, jeśli odpowiedź wydaje się oczywista ... Chciałbym przeprowadzić analizę szeregu budynków. Interesuje mnie ich powierzchnia elewacji wraz z odpowiednią orientacją. Jak pokazano na poniższym obrazku, chciałbym mieć długość i...