Zadanie : operacje na wektorach - interpolacja TIN

1. Stwórz mapy wektorowe obrazujące gęstość zaludnienia województwa śląskiego w roku 2008. Pokaż (graficznie) miejsca o największym indeksie zmian o podstawie zmiennej (In plus oraz In minus). W jakiej gminie zanotowano największy spadek, a gdzie największy wzrost omawianego wskaźnika.

2. Do bazy danych o ludności gmin podepnij informacje o powierzchni mieszkaniowej na osobę (wyniki przedstaw na mapach rastrowych).

3. W oparciu o moduł TIN wyznacz na podstawie centroidów gmin województwa śląskiego interpolowaną mapę powierzchni mieszkaniowej za rok 2008. Podpowiedź: mapę centroidów należy importować z formatu shp z użyciem systemu referencyjnego UTM34n, dane z odpowiednich pól eksportuj do pliku avl a następnie uzupełnij tabele centroidów.

Przykład:

Interpolacja gęstości zaludnienia

Ad1

Stwórz mapy wektorowe obrazujące gęstość zaludnienia województwa śląskiego. Pokaż (graficznie) miejsca o największym indeksie zmian o podstawie zmiennej (In plus oraz In minus). W jakiej gminie zanotowano największy spadek, a gdzie największy wzrost omawianego wskaźnika.

1. Zmień tzw. working folder na C:\(...)\Z4 aep\tworzenie mapy + tin\,

2. Otwórz DISPLAY Launcher, wybierz zakładkę warstwy wektorów i otwórz plik poligony_gminy_slaska_utm34n. Wybierz opcję Advanced Palette / Symbol Selection. Ponieważ każdy z poligonów identyfikowany jest za pomocą numeru ID nie zawiera informacji ilościowej można wybrać atrybut Qualitative (w zakładce Data relationship) oraz color Variety (black outline). Otworzono mapę wektorową Śląska który jest podzielony na 167 gmin - kliknięcie na poligon za pomocą Cursor Inquiry, ujawnia informację o numerze ID obszaru miejskiego/.

3. Kliknij ikonę Database Workshop (ikonka tabelki na przyborniku) i otwórz bazę danych regionalnych populacji - ludnosc_2008_gminy_slaska.mdb. Otworzy się baza danych zawierająca informację o populacji gmin. Utwórz łącznik (plik *.vlx) [komenda „establish display link], który jest plikiem „łącznikiem” między bazą danych ludnosc_2008_gminy_slaska.mdb a mapą wektorową ludnosc_2008_gminy_slaska.vct) /w terminologii baz danych kolumna nazywana jest polem (fields) zaś wiersz rekordem (record). Przełącz kursor w tryb Cursor Inquiry Mode i klikając w poszczególne poligony zaobserwuj zmiany w bazie danych.

(na łącznik składają się trzy komponenty: nazwa pliku łącznika (.vlx), nazwa warstwy (mapy) wektorowej (.vct) oraz nazwa pola bazy danych które ma niepowtarzalne ID.)

4. Kliknij gdzieś w polu ludn_2008 (populacja w roku 2008). Kliknij Display Current Field as Map Layer (pierwsze polecenie z menu podręcznego dostępnego po naciśnięciu prawego przycisku myszy lub przedostatnia ikonka w menu Database Workshop (tożsama z ikonką Display). /Kiedy mapa wektorowa jest połączona z bazą danych każde pole może być wyświetlone na mapie tematycznej/.

5. Eksport pola do mapy wektorowej. Kliknij gdzieś w polu ludn_2008. Z menu Database Workshop wybierz File/Export/Field/to Vector File. Zauważ, że program sam podpowiada sugestywną nazwę. Jeżeli zgadzasz się na nią kliknij OK. /program sam uzupełni metadane pobrawszy je z wsadowej mapy wektorowej - wyświetl je za pomocą Layer Properities/Properities/View Metadata/.

6. Eksport pola do mapy rastrowej. Kliknij gdzieś w polu pow_gmin_ha a następnie kliknij w ikonę Create IDRISI Raster Image (ikonka szachownicy) /zauważ że obok jest ikonka Create IDRISI Vector File/. Zostanie zaproponowana nazwa rastra. Po wybraniu OK zostanie wyświetlone okno, w którym można dokonać zmian w ustawieniach rastra - bezpieczniej jest jednak skopiować ustawienia z pliku SLASK_PUSTY.

Dodawanie nowych danych do bazy - obliczenie

7. Zamknij wszystkie mapy związane z tabelą (Window List w głownym meny IDRISI polecenie Close All Map Windows) /będziemy dodawać nowe dane do bazy danych, a można to uczynić (ze względów bezpieczeństwa) tylko jeżeli wszystkie powiązane z nią mapy zostaną zamknięte/

8. Wejdź do Edit w menu Database Workshop i wybierz opcję Add Field. Nazwij nowe pole gestosc_2008 i ustaw typ danych na Real, kliknij Ok.

9. Wybierz ikonę Calculate Field Values (+=) z paska Database Workshop. Z listy Set wybierz gestosc_2008 zaś w zakładce „=” wpisz: ([ludn_2008] / [POW_GMIN_HA]) Jeżeli wszystkie spacje zostały uwzględnione program zapyta czy zmodyfikować bazę danych, potwierdź i zobacz wynik/

10. Zapisz bazę danych. Upewnij się, że aktywna komórka znajduje się w polu gestosc_2008 a następnie wyświetl mapę za pomocą Display. Zapisz bazę danych.

AD2

Do bazy danych o ludności gmin podepnij informacje o powierzchni mieszkaniowej na osobę. Przedstaw te dane na mapach rastrowych.

1. Zamknij wszystkie mapy związane z tabelą (Window List w głownym meny IDRISI polecenie Close All Map Windows) /będziemy dodawać nowe dane do bazy danych, a można to uczynić (ze względów bezpieczeństwa) tylko jeżeli wszystkie powiązane z nią mapy zostaną zamknięte/

2. Dokonaj importu danych z arkusza „powierzchnia na osobe.xls”.Polecenie import/table/from external file. Wskaż na arkusz1 i potwierdź, że pierwszy wiersz to nagłówki.

3. Dokonaj połączenia dodanej tabeli z vectorem ludnosc_2008_gminy_slaska.vct

4. Wyświetl dane o powierzchni mieszkań na mapie rastrowej (analogicznie jak we wcześniejszych działaniach.

5. Można dane o powierzchni przenieść do tabeli ludność 2008 gmin śląska eksportując interesujące nas pole do pliku avl z tabeli pierwotnej i zaimportowanie tych danych do tabeli docelowej.

AD 3

W oparciu o moduł TIN wyznacz na podstawie centroidów gmin województwa śląskiego interpolowaną powierzchnią mieszkaniową na osobę za rok 2008. Podpowiedź: mapę centroidów należy importować z formatu shp z użyciem systemu referencyjnego UTM34n, dane z odpowiednich pól eksportuj do pliku avl a następnie uzupełnij tabele centroidów.

Zależy nam przede wszystkim na uzyskaniu mapy z naniesionymi centroidami o wartościach powierzchni mieszkań. Należy wykonać w tym celu szereg działań:

1. Spreparować mapę rastrową gmin śląska obrazującą wartości pola Idr_id

2. Przekształcić utworzoną mapę rastrową w wektorową (dzięki temu otrzymamy obrysy gmin i centroidy). Polecenie rastervector. Konwersja raster to poligon. Nazwy nadać jak na obrazku niżej.

Uwaga: WYKLUCZ POLYGON TŁA!!!

Otrzymujemy trzy mapy-wektory: poligony (obrysy), linie, punkty (locators) zaczepienia poligonów, centroidy

3. Mapę wektorową centroidów eksportujemy do formatu shape (export/software-specyfic formats/ersi formats/shapeidr)

• Ekportujemy mapę wektorową bez bazy danych (zaznaczamy odpowiednioą opcję)

- Przy okazji tworzenie pliku shp tworzona jest baza danych w formacie *.dbf

4. Import centroidów z pliku shape. Polecenie import/software-specyfic formats/ersi formats/shapeidr

Dzięki temu otrzymujemy plik wektorowy + baza danych w access (baza danych o centroidach)

I oto nam chodziło!!!!!

5. Uzupełniamy otrzymaną bazę danych centroidów o dane dotyczące powierzchni mieszkań. (import field from *.avl)

6. Tworzymy mapę wektorową zawierającą punkty (centroidy) o wartościach odpowiadających powierzchni mieszkań. (Polecenie export/field/vector image w database workshop)

Po otrzymaniu mapy użyj modułu TIN do przeprowadzenia interpolacji punktów pomiędzy centroidami.

7. Polecenie GIS Analysis/Interpolation/Tin Interpolation/TIN

(jako wejście punkty centroidów o wartościach powierzchni mieszkań.)

(na wyjściu mapa rastrowa obrazująca interpolację

(odpowiedź - powierzchnia mieszkań - interpolacja)

---