Mam plik ASCII o szerokości i długości geograficznej oraz data_val w następującym formacie. 35-13.643782N, 080-57.190157W, 118.6 ... Mam plik obrazu GeoTiff i mogę go łatwo wyświetlić. Chcę umieścić „pinezkę” (może być kropką / flagą / gwiazdą lub czymkolwiek najłatwiejszym) na obrazie na określonej szerokości / długości geograficznej znalezionej w pliku …
Używam pyqgis w samodzielnej aplikacji. Teraz chciałem utworzyć niestandardowe wyrażenie do użycia jako filterExpression z reguły ruleBasedRendererV2. (Próbowałem wcześniej renderera z istniejącymi funkcjami i wszystko działało dobrze!) Niestandardowa qgisfunction(bitowa funkcja AND - ponieważ takiej funkcji jeszcze nie ma) wygląda następująco: from qgis.utils import qgsfunction @qgsfunction(args="auto", group='Python') def bitwise_and(value1, value2, feature, …
Próbuję wymyślić, jak użyć Pythona do wyodrębnienia wielokątów w jednym wektorze, które nakładają się o> 90% przez inny wektor. Chciałbym wtedy mieć wektor / mapę, która pokaże tylko te wielokąty. Przykładowe zdjęcie pokazuje moje warstwy. Chcę wszystkie szare wielokąty, które są> 90% czerwone. Muszę to wszystko zrobić za pomocą Pythona …
Próbuję spakować skrypt w języku Python, który używa GDAL dla PyPI. Zacząłem od włączenia bezpośredniego odwołania do mojego setup.py: install_requires=['GDAL==1.11.2'], W ten sposób nie udało się zainstalować pakietu w moim testowym środowisku wirtualnym: extensions/gdal_wrap.cpp:2855:22: fatal error: cpl_port.h: No such file or directory #include "cpl_port.h" ^ compilation terminated. error: Setup script …
Ten mały, mały skrypt usuwa wszystkie małe funkcje na wszystkich warstwach. Działa to dobrze, ale rozmiar plików .shp nie jest zmniejszany. Proszę, powiedz mi dlaczego... canvas = iface.mapCanvas() expr = QgsExpression("$length<100") for l in canvas.layers(): l.startEditing() it = l.getFeatures(QgsFeatureRequest(expr)) for ft in it: l.deleteFeature(ft.id()) l.commitChanges()
Co staram się zrobić: pętlę przez punkt shapefile i wybrać każdy punkt, który wchodzi do wielokąta. Poniższy kod został zainspirowany przykładem zapytania przestrzennego znalezionego w książce: mitte_path = r"D:\PythonTesting\SelectByLocation\mitte.shp" punkte_path = r"D:\PythonTesting\SelectByLocation\punkte.shp" polygon = QgsVectorLayer(mitte_path, 'Mitte', 'ogr') points = QgsVectorLayer(punkte_path, 'Berlin Punkte', 'ogr') QgsMapLayerRegistry.instance().addMapLayer(polygon) QgsMapLayerRegistry.instance().addMapLayer(points) polyFeatures = polygon.getFeatures() pointsCount = …
Czy istnieje oficjalna dokumentacja modułu GDAL w języku Python? Istnieją różne książki kucharskie, ale mało informacji na temat każdej funkcji lub argumenty dla najbardziej podstawowych funkcji.
Chcę utworzyć zestaw danych wektora punktowego o wartości 10000 punktów (lub więcej) w binarnym rastrze, w którym punkty powinny być ograniczone do obszarów, w których wartość rastra wynosi 1. Próbowałem następujące kroki. Poligonizuj raster QGIS: Wektor -> Narzędzia badawcze -> Losowe punkty Działa to dobrze do 2000 punktów, ale wszystko …
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 (normalną) orientację wszystkich krawędzi w szeregu wielokątów. W przykładzie wyróżniłem …
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 Pythonie za pomocą OGR / GDAL. Mam zestaw plików kształtów i zestaw …
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 = layer.dataProvider().getFeatures(QgsFeatureRequest().setFilterExpression ( u'"field"=4')) (Nie jestem pewien, czy rozumiem różnicę). Chcę teraz wyświetlić tylko te funkcje. Dokładnie tak samo, jak przy …
czy ktoś z powodzeniem nakłania pandy do pracy z instalacją QGIS w Pythonie w systemie Windows? Użyłem instalatora OSGeo4W do zainstalowania QGIS2.2, który jest dostarczany z Python 2.7 i Numpy 1.5.1. Każda rozsądnie nowa wersja pand wymaga Numpy 1.6.1 lub nowszej. To wydaje się moim głównym problemem. Próbowałem: - Budowanie …
ArcMap ma możliwość wyświetlania rozszerzonego opisu każdego elementu na liście symboli. Jak wypełnić tę właściwość ze stołu? Celem jest uniknięcie konieczności edycji każdego symbolu jeden po drugim, aby dodać opis .
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 pól różnią się od „Shape_Area”, „SHAPE_Area”, …
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 .csv na warstwę elementów, a następnie na punktową klasę elementów Zauważyłem wyraźnie różne czasy przetwarzania w zależności …
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.