GIS Cwiczenie6, GIS itp
[ Pobierz całość w formacie PDF ]
//-->Ćwiczenie 6Konwersja układów współrzędnych.Analizy przestrzenne: funkcjaSelectiClip1.2.Skopiuj katalog dane zI:\TechnInfZSZsem2\Cwiczenie6\Danedo swojego kataloguroboczego na dysku na dyskuH: lub D:Otwórz program ArcCatalog i przemieść się do katalogu roboczego i podkataloguDaneD:lubH:\Cwiczenia06\NazwiskoImie\Dane.Sprawdź właściwości (Properties)warstwy tematyczne:Zlew_pn_region.shp, Biebrza_drogi.shp, Biebrza_gleby.shp,Biebrza_miasta.shp, Biebrza_wioski.shp.Zauważ, że nie maja zdefiniowanego XYCoordinate SystemNastępne zadanie będzie polegało na konwersji odwzorowań w zbiorach shp. Konwersjataka w ArcMap odbywa się pomiędzy zdefiniowanymi odwzorowaniami, którychparametry zapisane są w plikach .prj, według ogólnych geometrycznych zasadtransformacji. Wczytane do projektu warstwy nie mają przypisanego żadnegoodwzorowania, a współrzędne wyrażone są w jednostkach stopniach geograficznych.Najpierw więc przypiszemy im odwzorowanie geograficzne, w efekcie czego utworzynam się dodatkowy plik z rozszerzeniem prj (np. dla zbioruBiebrza_wioski.shputworzony zostanie dodatkowy plikBiebrza_wioski.prj)przechowujący informację odanym odwzorowaniu. Następnie przekonwertujemy go do nowego zbioru nadając noweodwzorowanie. Będzie to odwzorowanie mapy topograficznej w układzie PUW 1992.3.Definiowanie układu współrzędnych. WyświetlToolsBox. Znajdź funkcjęDataManagement Tools/Projections and Transformations/Define projection.W okienkumenuDefine projectionwybierz jakoInput datasetjeden ze zbiorów shp,przygotowanych do konwersji npBiebrza_wioski.shp.a jakoCoordinate SystemWGS1984 .prj(Select..Geographic Coordinate System/Word/WGS84.prj).Powtórz czynnośćdla wszystkich pozostałych warstw.W oknie Catalog Tree PKM naDanestwórz nową geobazęplikową o nazwie dane.gdbZmianaukładuwspółrzędnych. Znajdź funkcjęProject (Management).Możesz tozrobić innym sposobem niż w pkt.1: Kliknij na ikonę SearchwpiszProject.Wyświetli się spis znalezionych narzędzi – wybierzProject (DataManagment). Jako Input datasetwskażBiebrza_wioski,nadaj nową nazwę plikowiwynikowemuBiebrza_wioski_92,a jakoOutput Coordinate SystemwybierzPUWG1992(jak go nie ma w zakładceFavoritesto dodajz kataloguH:\Cwiczenia06\NazwiskoImie\Dane\PolskieUkladyWspolrzednych).Jako OutputDataset wstaw H:\....\Dane\Dane.gdb\biebrza_wioski_92.Potwierdź przez kliknięcie naklawiszeZastosujiOk.Powtórz opisaną powyżej procedurę definiowania układu współrzędnych geograficznychi konwersji do układu stosując ją kolejno dla tematów:Biebrza_drogi.shp,Biebrza_gleby.shp, Biebrza_miasta.shp, Zlew_pn_region.shp.Zapisuj je jako:Biebrza_drogi_92, Biebrza_gleb_92, Biebrza_miasta_92, Zlew_pn_region_92.4.5.6.7.8.Otwórz program ArcMap i dodaj wszystkie warstwy z geobazy Dane.gdb zapis projektjakoBiebrza.mxdSprawdź czy projekt ArcMap ustawił odwzorowanie na PUWG 1992 w tym celu wewłaściwościach całego projektu PKM na Layers (Okno Table od Contens)9.10.11.12.13.14.15.16.17.18.19.20.Następnie nadaj geobazie Dane.gdb właściwość roboczej PKMMake DefaultGeodatabase,od tej pory wszystkie analizy będą się zapisywały do niejW Searchzastosuj funkcjęSelect(AnalysisTools\Extract\Select)dla tematuZlew_pn_region_92oknieExpression (optional)wpisz równanie "PROFIL" = 'Sztabin',nazwij wyniksztabin_92.W efekcie zostanie zapisana tylko wybrana wcześniej zlewniajako oddzielny plik.Tworzenie nowych warstw. Analizy przestrzenne – Clip Zastosuj funkcjęClip(SearchlubToolsBoxAnalysis Tools\Extract\Clip)dla tematówBiebrza_drogi_92,Biebrza_gleby_92, Biebrza_miasta_92, Biebrza_wioski _92ograniczając ich obszar dogranic (wycinając) zbiorusztabin_92.Wyniki zapisuj każdorazowo jakoSztabin_drog_92, Sztabin_gleby_92, Sztabin_miasta_92, Sztabin_wioski_92.W wynikudziałania funkcjiClipz tematów zostają usunięte elementy znajdujące się poza obszaremtematu-nakładki:sztabin_92.Z aktualnego projektu Biebrza.mxd usuń wszystkie tematy o nazwieBiebrza…_92,pozostawiając w widoku tylko mapy w granicach zlewni Sztabin i zapisać jakoSztabin.mxd.Skoryguj wizualizację zbiorów tworząc mapy tematyczne: dla dróg (typ legendy:Uniquevalue,kolumna:klasa_drog)i dla gleb (typ:Unique value,kolumna:przydatność).Zapisz projektSztabin.mxd.Obliczenia powierzchni poligonów. W zbiorzeSztabin_gleby_92stwórz nową kolumnę onazwie:Area2,typie:Float.Oblicz powierzchnię obiektów znajdujących się w tym zbiorze. Porównaj z wynikamizgromadzonymi w kolumnieArea.Jak myślisz, co jest przyczyną tak dużych różnic wpowierzchni w obu kolumnach? Odpowiedź zapisz pliku wynik.txt w swoim kataloguroboczym.Otwórz czysty projekt w ArcMap i zapisz go w katalogu roboczym jako Rejestracja.mxdRejestracja mapy rastrowej. Do projektu dodaj skan mapy Mapa_do_rejestracji92.jpg zkataloguDaneSprawdź układ współrzędnych, w którym jest zapisany w treści skanu mapy (tj. znajdźzapis na skanie mapy), żeby ją zobaczyć PKM Zoom to LayerUstaw Odwzorowanie (Coordinate System) projektu na PUWG 1992 (PKM na Layers woknie Table do Contents)21.WArcMapieCustomize/Toolbars/Georeferencingdodajpaseknarzędzi22.Stań na nazwie warstwy rastrowej i PKM Zoom to layer23.Powiększprawy górny róg mapy i narzędziemwskaż (raz kliknij)skrzyżowaniesiatki kilometrowej tak jak na poniższym rysunku . Następnie nie przesuwając myszki zklikniętego miejsca naciśnij PKM i wybierz. Input X and Y i wpisz X=605000 iY=33900024.Stań na nazwie warstwy rastrowej i PKM Zoom to layer25.Znajdź jeszcze dwa przecięcia siatki i powtórz dla nich czynności z punktu 23.26.Po zarejestrowaniu minimum trzech punktów (nie mogą być w jednej linii) należywyeksportować mapę korzystając z opcji Rectify (umieszczonej w PaskuGeoreferencing)27.Następnie skorzystaj z opcji Reset transformation oraz usuńMapa_do_rejestracji.png28.Dodaj w/w zapisaną (w geobazie Dane.gdb) mapę do projektu ArcMap29.Dodaj map podkładową z Internetu (Add Basemap) o nazieTopographic, żeby sprawdzić czy poprawie dokonałej rejestracjimapy. Porównaj je wzrokowo czy nachodzą na siebie.zwidoku
[ Pobierz całość w formacie PDF ]