Systemy Informacji Przestrzennej

Adam Inglot
Adiunkt w Katedrze Geodezji
Wydział Inżynierii Lądowej i Środowiska, Politechnika Gdańska
adam.inglot@pg.edu.pl

2022-08-21, aktualizacja 2024-10-04

© ainglot.pl, 2022. Udostępnianie i wykorzystanie zgodnie z licencją Creative Commons Attribution 4.0 International license (CC BY 4.0).
Projekt strony zapożyczony za zgodą autora - Paulo Raposo

CC By 4.0

Wprowadzenie do ArcGIS Pro i QGIS

Celem przygotowanych ćwiczeń jest wprowadzenie użytkownika do prac w oprogramowaniu GIS. Pierwsze zagadnienie przeprowadzi czytelnika przez możliwość pobierania danych z państwowego zasobu geodezji i kartografii. Należy pamiętać, iż Geoportal Krajowy podlega zmianom i poniższe odnośniki mogą różnić się od aktualnego stanu. Zadania do wykonania:

katalog z rozpakowanymi danymi BDOT i PktA

Pobieranie danych z zasobów GUGiK

Wraz z nowelizacją ustawy prawo geodezyjne i kartograficzne, od 31 lipca 2020 część danych z zasobów Głównego Urzędu Geodezji i Kartografii zostały bezpłatnie udostępnione. Dotyczy to ortofotomapy, NMT, NMPT, danych ALS, BDOT10k, BDOO, Podstawowej Osnowy Geodezyjnej, Państwowy Rejestr Granic i model 3D budynków. Udostępnione zbiory danych zawierają ogromną ilość danych przestrzennych którymi będziemy się posługiwać w ramach ćwiczeń. Sposób pobierania danych jest stosunkowo prosty, należy przejść na stronę internetową Geoportal Krajowy a następnie wybrać dane które chcemy pobrać. Poniżej przedstawione zostały, krok po kroku czynności potrzebne do przygotowania środowiska pracy z danymi przestrzennymi.

1. Przejdź do przeglądarki internetowej i uruchom Geoportalu Krajowego.

logo geoportalu

W pierwszej kolejności pobierzmy dane rastrowe, ortofotomape dla terenu Politechniki Gdańskiej i okolicy.

2. Zbliż się do terenu PG za pomocą scrolla lub przycisku "+".

okno geoportalu

3. Przejdź do okna po prawej stronie.

menu boczne

4. Rozwiń Dane do pobrania\Ortofotomapa. Do wyboru są trzy opcje: według aktualności, według rozdzielczości i prawdziwa ortofotomapa.

ortofotomapa do pobrania

5. Rozwiń Ortofotomapa wg rozdzielczości, a następnie zaznacz Piksel do 5 cm - arkusze.

ortofoto wg aktualności

6. Wstaw lewy przycisk myszy (LPM) w okolicach budynku hydro PG.

wskazanire gmachu głównego PG

7. Rozwiń zbiory danych i pobierz dane PPM na wskazany link.

menu ortofotomapy

Dla wskazanej lokalizacji mamy do wyboru dużą paletę danych, według roczników: 2022, 2021, 2020, 2019, 2018, 2016, 2013, 2010, 2004. Co oznacza że na podstawie udostępnionych danych można zaobserwować zmiany przestrzenne w zakresie 17 lat. Dane te są dostępne w postaci zobrazowań multispektralnych: RGB - zakres światła widzialnego, CIR - zakres bliskiej podczerwieni. Rozdzielczość przestrzenna dostępnych danych to: 0.5, 0.25, 0.1, 0.05 metra.

8. Rozwiń pierwszy dostępny zbiór danych ortofotomapy (N-34-50-C-c-4-2 2022 RGB), zobacz jakie są informacje dostępne odnośnie tych danych.

metadane ortofotomapa

W pierwszej kolejności mamy dostępny link do danych, następnie numer godła mapy, data wykonania zobrazowania, rozdzielczość przestrzenna w metrach oraz w jakim układzie przestrzennym zostały zapisane dane. Pozostałem informacje nie wnoszą znaczącej treści w naszych rozważaniach.

Do wykonania mapy z zadania potrzebujemy pobrać trzy roczniki ortofotomapy:

9. Chcąc pobrać dane kliknij LPM* na Pobierz plik danych dla tej sekcji, dane zostaną pobrane do folderu domyślnego.

Następnymi danymi będą dane z Bazy Danych Obiektów Topograficznych dla szczegółowości obiektów w skali 1:10 000 - BDOT10k. Baza ta zawiera informacje przestrzenne o obiektach topograficznych w postaci danych wektorowych. Dane mogą być pobrane w dwóch formatach: ESRI Shapefile - SHP oraz w Geography Markup Language - GML, podczas pobierania decydujemy dla jakiego poziomu administracyjnego chcemy pobrać dane: cały kraj, poszczególne województwo, powiat.

10. Następnie, pobierz dane BDOT10k. w Dane do pobrania rozwiń Topografia i zaznacz Baza danych Obiektów Topograficznych.

pobieranie BDOT10k

11. Skieruj kursor na okno mapy w Geoportalu Krajowym LPM, rozwiń i wybierz Powiat Gdańsk - SHP.

pobieranie BDOT10k pobieranie BDOT10k

Przygotowanie miejsca pracy

Zanim jednak przejdziemy do pracy w środowisku GIS, zastanówmy się nad tym jak skutecznie przygotować dane do wykonywać analizy przestrzennych. QGIS i ArcGIS Pro różnią wieloma elementami jak np. interfejsem, narzędziami, obsługą formatów danych. Jednak w obu środowiskach schemat pracy jest podobny - budujemy projekt do którego "podpinamy" dane. Celowo użyłem słowa podpinamy, żeby zrozumieć że dane nie są kopiowane czy importowane do projektu (oczywiście są od tego wyjątki) ale stanowią niejako załącznik. Za pomocą oprogramowania GIS wyświetlamy dane oraz przeprowadzamy analizę tych danych lub w oparciu o nie. Wiedząc o tym należy pamiętać, że projekt nie przechowuje danych a jedynie informację gdzie dane się znajdują i sposób ich wyświetlania (symbolizację, etykiety). W ArcGIS Pro nastąpiła pewna zmian (porównując wersję ArcGIS 10 i starsze), podczas zakładania projektu budowane jest środowisko pracy (w wybranym przez użytkownika miejscu na dysku) w którym docelowo powinny znajdować się wszystkie dane oraz wyniki. Jednak korzystanie z tych samych danych w wielu projektach i powielanie ich wydaje się niepraktyczne.

Przygotowanie miejsca do pracy stanowi pewne wyzwanie przy dużych projektach, gdy mamy wiele źródeł danych i dużą liczbę danych wynikowych. Struktura folderów w których przechowujemy dane zależy od przyzwyczajenia użytkownika. Należy jednak unikać "niekontrolowanego chaosu", starać się porządkować zbiór danych i usuwać niepotrzebne. Dane w folderach można układać pod względem przestrzennym: zgodnie z podziałem administracyjny lub siatką, możemy przyjąć podział ze względu na rodzaj danych: wektorowe i rastrowe, czy też w zależności od produktu: np. ortofotomapa, BDOT10k, BDOT500. Dużo zależy od specyfiki naszych prac, które ciężko przewidzieć na początku naszych analiz ale jeżeli mamy wyznaczone zadanie do wykonania, to starajmy się przygotować foldery oraz dane do zadania.

Zakładając strukturę folderów należy kierować się kilkoma zasadami, które jednak są zaleceniami a nie restrykcjami:

12. Na dysku (jeżeli komputer ma więcej niż jeden dysk np. dysk D, skorzystaj z niego) przygotuj folder "SIP", a w nim "ProjektyArcGIS".

Pierwsze kroki w ArcGIS Pro

Oprogramowanie ArcGIS Pro łączy kilka osobnych platform z wersji ArcGIS for Desktop 10 - ArcCatalog, ArcMap, ArcScene, ArcGlobe. Przed instalacją oprogramowania należy porównać parametry komputera na którym chcemy zainstalować ArcGIS do minimalnych wymagań zapewniających poprawne funkcjonowanie programu. W przypadku problemów technicznych z sprawnym funkcjonowaniem środowiska ArcGIS Pro należy sięgnąć po rozwiązania serwerowe dla studentów ePIONIER - najlepiej wybrać Windows10-ANSYS i tam zainstalować ArcGIS.

13. Po instalacji uruchomi się okno logowania do ArcGIS, użyj loginu i hasła do konta ArcGIS Online. Jeżeli nie posiadasz danych do logowania przeczytaj informacje zawarte na stronie internetowej WILiŚ.

logowanie do ArcGIS

14. Utwórz nowy projekt w którym wyświetlisz pobrane dane.

panel projektów

15. Umieść projekt w folderze "ProjektyArcGIS" oraz wprowadź nazwę "GDA".

nowy projekt

Okno katalogu Catalog w ArcGIS można uruchomić przechodząc w menu głównym (lewy górny róg okna projektu) do View a następnie uruchamiając Catalog Pane, (okno katalogu może być już uruchomione):

uruchom okno katalogu

Okno Catalog:

uruchom okno katalogu

16. Przyglądnij zawartość folderu. Plik GDA.aprx to plik projektu, GDA.gdb to domyślna geobaza przygotowana wraz z budową nowego projektu (co wykonane zostało kilka punktów wcześniej).

pusty projekt

17. Przygotuj nowy folder na dane pobrane dla części obszaru Politechniki Gdańskiej. Można wykonać to z poziomu katalogów Windows lub w ArcGIS. Kliknij PPM na folder GDA w Catalog, wybierz New, a następnie Folder. Podaj nazwę "danePG".

nowy katalog

18. Jeżeli chcesz uruchomić katalog z poziomu Windows wskaż (naprowadź kursor) na katalog, pokaże się okno z ścieżką do katalogu. Kliknij LPM na link.

otwórz katalog

19. Przenieś dane do katalogu danePG.

katalog w windows

20. Rozwiń katalog danePG w ArcGIS, nie wszystkie dane są widoczne. Dane ZIP należy rozpakować.

widok katalogu

21. Pobierz kolejne dane, tym razem punkty adresowe. Przejdź ponownie do Geoportalu Krajowego i w Dane do pobrania rozwiń Państwowy Rejestr Granic, a następnie Adresy i ulice wg powiatów.

katalog z rozpakowanymi danymi BDOT i PktA katalog z rozpakowanymi danymi BDOT i PktA katalog z rozpakowanymi danymi BDOT i PktA

22. Zastosuj rady dotyczące nazw folderów i zmień nazwę katalogu SHP2261 na PktA_2261, przenieś dane do katalogu projektu ArcGIS. Sprawdź co się zmieniło w ArcGIS. Jeżeli nic się nie zmieniło, kliknij PPM na katalog i wybierz Refresh.

23. Dodaj ortofotomape z roku 2022 do projektu, przeciągnij z katalogu na okno mapy i upuść.

katalog z rozpakowanymi danymi BDOT i PktA

24. Następnie rozwiń katalog z danymi BDOT10k (2261_SHP) i dodaj do projektu dane PL.PZGiK.336.2261__OT_BUBD_A.shp, jest to warstwa budynków, zbliż się do budynku hydro.

katalog z rozpakowanymi danymi BDOT i PktA

25. Mając wskazaną warstwę PL.PZGiK.336.2261__OT_BUBD_A.shp w tabeli zawartości, w menu głównym LPM na Feature Layer, następnie uruchom okno Symbology i przejdź do zmiany symbolu.

katalog z rozpakowanymi danymi BDOT i PktA

26. Zmień symbol na drugi - Black Outline (2 pts), jest to symbol samego konturu bez wypełnienia.

katalog z rozpakowanymi danymi BDOT i PktA

Dodamy teraz trzy okna mapy, tak aby można było w oknie wydruku jednocześnie zaprezentować dane z różnych okresów czasowych.

27. Przejdź do zakładki Insert w głównym menu programu, następnie rozwiń New Map i wybierz New Map.

katalog z rozpakowanymi danymi BDOT i PktA

28. Przejdź do nowego okna mapy Map1 a następnie w oknie zawartości 2LPM na Map1 - zaznaczone na czerwono na poniższym obrazku.

katalog z rozpakowanymi danymi BDOT i PktA

29. W zakładce General, zmień nazwę na 2022, zatwierdź OK.

katalog z rozpakowanymi danymi BDOT i PktA

30. Wykonaj analogicznie dla pozostałych roczników.

katalog z rozpakowanymi danymi BDOT i PktA

31. Skopiuj warstwę budynków do pozostałych okien.

katalog z rozpakowanymi danymi BDOT i PktA

32. W każdym oknie umieść odpowiednią ortofotomapę - z danego rocznika.

katalog z rozpakowanymi danymi BDOT i PktA

Przejdziemy teraz do okna w którym będziemy mogli komponować grafikę czy wydruk mapy - Layout. W tym oknie możemy dodawać elementy które nie są dostępne w oknie mapy jak: skala, podziałka liniowa czy strzałka północy. Elementy umieszczane w tym oknie połączone są z wybranym oknem mapy lub z danymi w projekcie np. podziałka liniowa dynamicznie zmienia się względem zmian skali w oknie mapy powiązanym z tą podziałką. Do okna wydruku możemy dodać też elementy nie powiązane z danymi jak tekst, elementy graficzne np. ramka okrąg, również możemy dodać obrazek, zdjęcie czy tabelkę. Po skomponowaniu grafiki czy wydruku mapy, możemy wyeksportować widoczną grafikę w obrębie arkusza wydruku.

33. Ponownie przejdź do menu głównego, zakładka Insert, rozwiń New Layout i wybierz arkusz A4 ISO - Landscape.

katalog z rozpakowanymi danymi BDOT i PktA

34. Otrzymamy nowe okno z widocznym obszarem.

katalog z rozpakowanymi danymi BDOT i PktA

35. W menu głównym, przejdź do Insert a następnie rozwiń Map Frame, wybierz jedną z map.

katalog z rozpakowanymi danymi BDOT i PktA

36. Przeciągając LKM w oknie wydruku mapy, wstaw ramkę mapy.

katalog z rozpakowanymi danymi BDOT i PktA

37. Analogicznie wykonaj dla pozostałych okien mapy.

katalog z rozpakowanymi danymi BDOT i PktA

Jeżeli skala w Twoich opracowaniach nie jest podobna do powyższego rysunku, nie przejmuj się, skale będziemy wprowadzać na późniejszym etapie. Ustalmy teraz rozmieszczenie okiem mapy w oknie wydruku.

38. Wstaw LPM na pierwszą mapę (rok 2018) i przejdź co Properties.

katalog z rozpakowanymi danymi BDOT i PktA

39. Przejdź do ustawień rozmieszczenia obiektu Placement.

katalog z rozpakowanymi danymi BDOT i PktA

40. Wprowadź wartości widoczne na rysunku poniżej.

katalog z rozpakowanymi danymi BDOT i PktA

katalog z rozpakowanymi danymi BDOT i PktA

Widzimy iż dłuższa ściana budynku Centrum Ekoinnowacji PG jest skierowana pod kątem względem kierunku północy. Możemy wykonać obrót okna mapy, dostosowując do okna wydruku.

41. Przejdź do właściwości okna mapy 2LPM na nazwę okna mapy np. 2022 w tabeli zawartości - Contents. Wykonaj analogicznie dla pozostałych map.

katalog z rozpakowanymi danymi BDOT i PktA

Wyszukajmy teraz obiektu reprezentującego już nieistniejącą halę przy budynku hydro. Posłużymy się tym obiektem aby wyśrodkować wszystkie rysunki do tego obiektu.

42. Przejdź do tabeli atrybutowej warstwy budynków jednego z okna mapy np. 2022, warstwa PL.PZGiK.336.2261__OT_BUBD_A. Użyj PPM i przejdź do Attribute Table.

katalog z rozpakowanymi danymi BDOT i PktA

43. W tabeli atrybutowej wyszukaj obiektu o numerze FID = 9988.

katalog z rozpakowanymi danymi BDOT i PktA

44. Nie zapomnij wyczyścić selekcję.

katalog z rozpakowanymi danymi BDOT i PktA

Teraz dla wszystkich okien mapy ustawimy taką samą skale - 1:1000. Wykonamy to dla każdego okna, ważne aby te okna były scentrowane na budynek hali przy hydro.

45. Wskaż LPM jedno z okien mapy w oknie wydruku, przejdź do skali w lewym dolnym rogu okna wydruku - wprowadź skalę 1:1000.

katalog z rozpakowanymi danymi BDOT i PktA

Na koniec wstawmy elementy mapy takie jak: strzałka północy, podziałkę liniową, tytuł wydruku i autora pracy.

46. W menu głównym przejdź do wstawiania tekstu i wprowadź odpowiedni tytuł grafiki.

katalog z rozpakowanymi danymi BDOT i PktA

47. Następnie rozwiń North Arrow i wstaw jedną z wybranych.

katalog z rozpakowanymi danymi BDOT i PktA

48. Wstaw podziałkę liniową.

katalog z rozpakowanymi danymi BDOT i PktA

Dla każdego wstawianego obiektu można wprowadzić zmiany takie jak: wprowadzić kolor tła, dodać ramkę, zaokrąglić naroża. Dla obiektów dynamicznych takich jak podziałka liniowa możemy zmienić rodzaj miar - metry, kilometry, mile itp., zmienić liczbę podziału podziałki.

49. Przejdź do właściwości podziałki liniowej i dodaj białe tło. Kliknij PPM na podziałkę liniową i wybierz Properties.

katalog z rozpakowanymi danymi BDOT i PktA

50. Przejdź do właściwości podziałki liniowej ale w menu głównym i zmień jednostki na metry.

katalog z rozpakowanymi danymi BDOT i PktA

51. Wstaw jako tekst imię i nazwisko - dodaj efekt Hallo tak żeby tekst był widoczny.

katalog z rozpakowanymi danymi BDOT i PktA

Na koniec wprowadzimy roczniki dla wprowadzonych ortofotomap i wyeksportujemy grafikę do obrazu lub do pliku PDF.

52. Korzystając już ze zdobytej wiedzy wprowadź roczniki dla map, możesz pogrubić efekt Hallo.

katalog z rozpakowanymi danymi BDOT i PktA

53. Przejdź do menu głównego do zakładki Share, rozwiń Export Layout i wybierz np. Web JPEG.

katalog z rozpakowanymi danymi BDOT i PktA

54. W zakładce Export Layout można zmienić miejsce i nazwę pliku eksportowanego oraz zmienić rozdzielczość (zostawmy wartości domyślne), żeby zapisać wybierz w prawym dolnym rogu okna Export.

katalog z rozpakowanymi danymi BDOT i PktA

Wynik naszej pracy jest wyeksportowany plik, wygenerowana grafika może być wstawiana w dokumentach lub przesyłana jako pojedynczy plik. Dzięki oknie Layout można generować grafikę do artykułów naukowych, posterów czy map które możemy wydrukować na ploterze.

katalog z rozpakowanymi danymi BDOT i PktA

Pierwsze kroki w QGIS

Program QGIS jest darmowym oprogramowaniem do pracy z danymi przestrzennymi, stanowi poważny zamiennik dla platformy ArcGIS Pro i w niektórych elementach sprawdza się lepiej, np. w prędkości odświeżania wyświetlania dużych zbiorów danych. QGIS również wspiera wolne formaty danych takie jak GeoJSON czy GML.

Zadanie do wykonania: przygotuj wydruk mapy z gmachem głównym Politechniki Gdańskiej w środowisku QGIS, na wydruku mapy umieść dwa okna mapy w jednym ortofotomapę a na drugim Numeryczny Model Pokrycia Terenu (NMPT) pobrany z Geoportalu Krajowego.

Pobierz dane:

55. W przeglądarce internetowej przejdź do Geoportalu Krajowego, w zakładce.

usługi przeglądania WMS i WMTS

56. Pobierz dane 70934_850784_N-34-50-C-c-4-2.asc.

usługi przeglądania WMS i WMTS

Przejdziemy teraz do programu QGIS, wersje programu mogą się różnic na prezentowanych grafikach, jednak narzędzia z których korzystamy na zajęciach stanowią podstawowe funkcje programu często stabilne i niezmienne pomiędzy wersjami programu. Nie mniej, należy mieć na uwadze że czasami elementy w programie mogą się różnić z tymi na grafice.

57. Uruchom program QGIS. Wybierz Nowy pusty projekt.

usługi przeglądania WMS i WMTS

Dostęp do danych w QGIS jest możliwy w oknie Przeglądarka. Sposób przeglądania danych jest intuicyjny, rozwijamy kolejne foldery w drzewie katalogów. Dane które możemy przeglądać i dodawać są to nie tylko dane przestrzenne. QGIS posiada większą paletę dostępnych formatów pliku do odczytu niż ArcGIS. Jednak przy dodawaniu danych do projektu należy zwrócić uwagę na rozszerzenia plików.

58. Dodaj dane do projektu.

usługi przeglądania WMS i WMTS

Uruchamianie dwóch widoków mapy

59. Wyłącz z wyświetlania NMPT a następnie Dodaj motyw ..

usługi przeglądania WMS i WMTS

60. Wprowadź nazwę orto.

usługi przeglądania WMS i WMTS

61. Analogicznie wykonaj dla NMP.

usługi przeglądania WMS i WMTS

62. Dodajmy teraz drugie okno do projektu. Przejdź do menu głównego Widok a następnie wybierz Nowy widok mapy.

usługi przeglądania WMS i WMTS

63. W nowym widoku mapy wybierz motyw NMPT.

usługi przeglądania WMS i WMTS

usługi przeglądania WMS i WMTS

Zmieńmy teraz sposób wyświetlania warstwy rastrowej - NMPT.

64. Użyj 2LPM na gradient biało-czarny przy warstwie NMPT aby przejść do właściwości warstwy.

usługi przeglądania WMS i WMTS

65. W oknie symbolizacji zmień Sposób wyświetlania na Jednokanałowy pseudokolor.

usługi przeglądania WMS i WMTS

66. Następnie rozwiń paletę kolorów i wybierz Utwórz nową paletę kolorów ..

usługi przeglądania WMS i WMTS

67. Wskaż katalog: cpt-city, zatwierdź OK.

usługi przeglądania WMS i WMTS

68. W oknie po lewej stronie wskaż Topography a następnie wybierz i zatwierdź wkp/schwarzwald/wiki-schwarzwald-d050.

usługi przeglądania WMS i WMTS

usługi przeglądania WMS i WMTS

69. Przejdź do wydruku mapy w menu głównym Projekt a następnie Nowy wydruk..

usługi przeglądania WMS i WMTS

70. Wprowadź nazwę PG.

usługi przeglądania WMS i WMTS

Będziemy teraz pracować w oknie wydruku, tak jak w ArcGIS Pro w Layout tak i tu w QGIS możemy dodawać różne elementy istotne dla grafiki czy wydruku mapy.

71. W menu głównym okna wydruku przejdź do Dodaj obiekt wstaw ``.

usługi przeglądania WMS i WMTS

72. Przeciągnij okno mapy LPM w oknie wydruku.

73. Analogicznie wstaw drugie okno mapy a następnie przydziel dla okna po lewej motyw ortoa dla okna po prawej NMPT zgodnie z grafiką poniżej.

usługi przeglądania WMS i WMTS

74. Dodaj i edytuj elementy mapy: strzałka północy, podziałka liniowa oraz etykiety jako tytuł i jako informacja o autorze.

usługi przeglądania WMS i WMTS

75. Mapę wynikową wyeksportuj do pliku PNG i zapisz domyślne ustawienia.

usługi przeglądania WMS i WMTS

usługi przeglądania WMS i WMTS

76. Zakończ pracę w QGIS nie zapisując projektu.

Używane skróty:

Autor: Adam Inglot, mail: adam.inglot@pg.edu.pl