Systemy Informacji Geograficznej

Pytania i odpowiedzi dla kartografów, geografów i specjalistów GIS

2
Dopasowane do mapy punkty GPS do sieci drogowej!
Jestem całkiem nowy w Postgres i PostGIS. Próbuję wykonać proste dopasowanie do mapy (lub wnioskowanie o ścieżce, jeśli chcesz) niektórych punktów GPS. Używam postgres (postgis) . Mam bazę danych w mojej bazie danych zawierającą moje punkty GPS: punkty gp ( punkty , długie, czas, prędkość, azymut, geometria, ...) Zaimportowałem również …
4
Zapisać widok wydruku / mapy kompozytora QGIS jako PNG / PDF przy użyciu Pythona (bez zmiany czegokolwiek w widocznym układzie)?
Mam otwarty widok kompozytora / wydruku QGIS ze wszystkimi elementami dostosowanymi zgodnie z życzeniem. Teraz muszę wydrukować / zapisać / wyeksportować to jako plik PNG lub PDF ze konsoli Python / skryptu Python. Nie chcę niczego zmieniać w bieżącym układzie. Większość przykładów, które znalazłem (na przykład: to ), zmienia pozycję …
3
Przetwarzanie obrazu przy użyciu Pythona, GDAL i Scikit-Image
Mam problem z przetwarzaniem i mam nadzieję, że uda mi się tutaj rozwiązać. Pracuję z teledetekcją stosowaną w leśnictwie, szczególnie z danymi LiDAR. Pomysł polega na użyciu obrazu Scikit do wykrywania wierzchołków drzew. Ponieważ jestem nowy w Pythonie, uznałem wielki osobisty triumf, aby wykonać następujące czynności: Zaimportuj CHM (z matplotlib); …
1
Partia ładuje wiele plików kształtów do Postgis
Zgodnie z shp2pgsqlpomocą opcje upuszczania , dołączania , tworzenia i przygotowywania wzajemnie się wykluczają. Tak więc, jeśli chcę utworzyć tabelę z kształtu, a następnie dołączyć wiele innych plików kształtów, robię coś takiego, zachowując licznik wskazujący, czy jesteśmy w trybie tworzenia, czy dołączania. cnt=0 for shp in $(ls *.shp); do if …
1
QGIS przenoszenie warstwy / wszystkich funkcji za pomocą Vector Affine Transformation
Używam QGIS do analizy niektórych śladów lotu, które zostały wyprodukowane przy użyciu innego oprogramowania. Teraz, gdy znajdują się w QGIS, punktem odniesienia wszystkiego jest 0,0, który jest punktem odniesienia lotniska. Ale w prawdziwym świecie ten punkt odniesienia lotniska jest wartością na wschodzie i północy. tj. Easting = 286499.025 Northing = …
2
Robisz dwuwymiarową rampę kolorów za pomocą QGIS?
Wiele najnowszych map, które widziałem, ma te zgrabne dwukierunkowe gradienty kolorów, których można użyć do korelacji kolorów między ciągłymi oszacowaniami. Oto przykład zaczerpnięty z ostatniej publikacji. Christian Levers, Pieter J. Verkerk, Daniel Müller, Peter H. Verburg, Van Butsic, Pedro J. Leitão, Marcus Lindner, Tobias Kuemmerle - Czynniki wzorców intensywności pozyskiwania …
2
Korzystasz z prostej matematyki w pliku SLD na Geoserver?
Czy można używać prostej matematyki w pliku SLD na Geoserver? Załóżmy, że chcę odczytać właściwość rozmiaru tekstu z bazy danych i pomnożyć jej wartość. Jak to zrobić? <sld:TextSymbolizer> ... <sld:CssParameter name="font-size"> <ogc:PropertyName>text_size</ogc:PropertyName> * 2.1 </sld:CssParameter> ... </sld:TextSymbolizer> To proste podejście oczywiście nie działa, więc pomyślałem, że musi być jakiś sposób, …
11 geoserver  sld 
2
Czy istnieje sposób na podzielenie wielokątów przez kilka linii w QGIS?
Mam warstwę wielokąta (.shp) i linii (.shp). Chcę podzielić wielobok na kilka części wzdłuż linii pokrywających się z wielokątem. Następnie chcę edytować każdy wielokąt osobno (tabela atrybutów). Próbowałem obejść: buforowanie linii z niewielkiej odległości, a następnie zadziałało z narzędziami geoprzetwarzania (różnica - „Unterschied”). Następnie „multipart to singlepart”. Działało, ale wolę …
11 qgis  vector  line 
1
Automatyczne wykrywanie śladów
Biorąc pod uwagę obszar z zagłębieniami w otwartej przestrzeni od; powiedzmy wszystkie pojazdy terenowe jadące na miękkiej powierzchni. Pojazdy będą powodować równoległe wgłębienia (ścieżki) o głębokości około 10-20 centymetrów i szerokości około 15-30 cm, przy długościach różniących się w zależności od wytrzymałości powierzchni. Które platformy teledetekcji byłyby odpowiednie do późniejszej …
1
Zmiana hasła administratora geoserver
Jak mogę zmienić domyślne hasło administratora geoserver geoserverna inne? Próbowałem dwóch różnych rzeczy, ale żadna nie działała. Pierwszą rzeczą, którą próbowałem postępować zgodnie z dokumentacją online, jest przejście do katalogu data_dir pod adresem: /usr/local/lib/geoserver-2.5/data_dir/security i mogłem tylko znaleźć <user enabled="true" name="admin" password="digest1:XXXX..."> w usergroup/default/users.xml. Ale nie mogę edytować podsumowania. Drugą …
1
Jak uzyskać lokalizację funkcji w OpenLayers v3
Po użyciu interakcji rysowania lub modyfikacji interakcji mam detektora, który zwraca funkcję, która zmodyfikowała lub dodała. Przykład draw.on('drawend', function (event) { // get the feature var feature = event.element; // ...listen for changes on it logStatus(feature.getId()); }); Wiem, jak uzyskać identyfikator, ale potrzebuję lokalizacji (lon lat) funkcji, ponieważ muszę ją …
3
Dokładniejszy sposób obliczania powierzchni rastrów
W mojej codziennej pracy ciągle jestem proszony o obliczanie obszarów globalnych zestawów danych rastrowych w projekcji geograficznej przy rozdzielczości 30 sekund kątowych. Te zestawy danych są zwykle wynikiem operacji Łączenia (typowym przykładem są klasy roślinności połączone z warstwą kraju). Aby to zrobić, nasza jednostka utworzyła zestaw danych rastrowych z obszarem …
3
1 km okrąża długie punkty w wielu miejscach na świecie
Mam setki lat-długich punktów rozmieszczonych na całym świecie i muszę tworzyć wielokąty kołowe wokół każdego z nich, o promieniu dokładnie 1000 metrów. Rozumiem, że punkty należy najpierw rzutować ze stopni (długości lat) na coś z jednostkami metrowymi, ale jak można to zrobić bez ręcznego wyszukiwania i definiowania stref UTM dla …
2
Czytanie geometrii postgis z kształtową
Robię drobny przepływ pracy w notatniku ipython z geopandami i kształtnie, wyciągając garść danych geoprzestrzennych, czasem z plików kształtów, czasem z postgis (gdzie odbywa się droższe przetwarzanie). Teraz przeciągam tabele Postgis do Pythona sqlalchemy, przekształcając geometrię do WKT po drodze, otrzymując coś takiego: sql = """ SELECT ST_AsText(ST_Transform(the_geom,4326)) as newgeom,* …
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.