GIS Cvičení 2. Georeferencování Souřadnicové připojení rastrového souboru Georeferencování Transformace obrazového záznamu do libovolné souřadnicové soustavy. MMÍTŠl [•TíítTrauTičTiTiítTOŠTil vyjádřena již v systému mapových souřadnic. Stupeň transformace - vztah mezi souřadnicemi zdrojovými a cílovými je vyjádřen polynomem n-tého stupně. Output image orientation Input image orientation" Provádí se za pomoci tzv. vlícovacích bodů (control points), kdy jsou určitému pixelu rastru určeny konkrétní souřadnice. Každému pixelu georeferencovaneho obrazu je tak přiřazena nová hodnota, vypočtená na základě obrazu \ původního - tzv. pře vzorkování. metoda nejbližšího souseda bilineární interpolace kubická konvoluce Souřadnicové připojení rastrového souboru Pomocí funkce Georeferencing - nutno aktivovat v menu View- Toolbars -Georeferencing, po zatrhnutí se objeví nástrojová lišta. V okénku vrstvy (Layer) vybereme rastr, který chceme souřadnicově připojit, v našem případě 24-32-05.tif klikneme na ikonu B (Add control points). Georeferencing Georeferencing T Layer: 124-32-05.jpg "3 & Nyní levým tlačítkem myši klikneme na bod rastru o známých souřadnicích (tzv. vlícovací bod) a pomocí pravého tlačítka funkcí Input X and Y... tomuto bodu přiřadíme souřadnice X,Y. Po označení potřebného množství bodů dojde k transformaci rastru. V Data View H Layout View Zoom Layout Bookmarks □I Scrollbars B& Rulers t^ Guides 01 Grid Input X and' Cancei ^oint | ^ Status Bar Overflow Annotation [jj1 Data Frame Properties.. | ^ Main Menu 3D Analyst Advanced Editing Annotation ArcPad ArcScan Data Frame Tools Dimensioning Disconnected Editing l| ^ Draw |V Editor Effects I ^ Georeferencing Enter Coordinates X: -585000 Y: -1150000 V případě použité základní mapy ČR je nutné provést transformaci pomocí mapových křížků (pomyslné průsečíky celých hodnot souřadnic Y a X po 1000 metrech), souřadnice křížků najdeme na mapovém okraji (horní a spodní okraj - souřadnice Y, levý a pravý okraj - souřadnice X, uváděno v tisících metrů). Mapové křížky nejblíže rohům zvětšíme na viditelnost jednotlivých pixelů a kliknutím na střed křížku levým tlačítkem označíme vlícovací bod. Pomocí pravého tlačítka funkcí Input X and Y... přiřadíme souřadnice Y a X zjištěné z mapového okraje. Vzhledem k orientaci souřadnicového systému S-JTSK (kladná poloosa Y směřuje na západ, X na jih) je nutné zadávat souřadnice Y namísto souřadnic X a naopak a předřadit znaménko mínus. Po zadání prvního kontrolního bodu dojde k posunu rastru na zadané souřadnice, je proto nutné si jej zobrazit pomocí ikony □ (Full Extent) a poté pokračovat v transformaci. K transformaci rastru dojde již po zadání dvou bodů, pro větší přesnost a pro výpočet střední kvadratické chyby je nutné zadat nejméně čtyři. Případné opravy cílových souřadnic je možné provést v tabulce kliknutím na ikonu 151 (Link Table) případně je možné uložit transformační tabulku ve formě textového souboru *.txt K ukončení transformace rastru----------------------------------------------- klikneme v menu Georeferencing -Update georeferencing, dojde k vytvoření souborů *.tfw (obsahuje informace o umístění v souřadnicovém systému) a *.aux (definuje souřadnicový systém). Při další práci s rastrem je vždy nutné kopírovat všechny 3 soubory (*.tiff, *.tfw, *.aux). -417.100437 -1220.921530 -4073.204273 -3574.965301 -596000.000000 -1142000.0000... -593000.000000 -1143000.0000... -597000.000000 -1145000.0000... -593000.000000 -1145000.0000... P Auto Adjust Transformation: 11 st 0rder Polynomial [Affine] ^ | Total RMS Error: |0.14677 Load... Save... I 0K~ .'■-**n\. X -596fO0Ö PY,,,1 142 0C I *» J '.- ■>i \ ; ^ >' HA. "VA* km i> ]; :■--,— .^:t:> ':. i*M) ;■■■■ ■' i" / ^r. Á, -1145 OOij i i r vy v/ Vr v^gX'r ' !~- ^ ^ii_-.i Souřadnice bodů Levý horní roh Pravý horní roh Levý dolní roh Pravý dolní roh ink Table ÍŤlfxl Link X Source Y Source X Map I Y Map [ 1209.149808 4847.011757 407.604341 5094.139491 -417.188437 ■1220.921530 ■4073.234273 ■3574.965901 ■596000.000000 ■593000.000000 ■597000.000000 ■593000.000000 -1142000.0000.. -1143000.0000.. -1145000.0000.. -1145000.0000.. 0.12662 0.18992 0.09496 0.15827 j Střední kvadratická odchylka P Auto Adjust Transformation: 11 st 0rder Polynomial [Affine) ▼ | Total RMS Error: 10.14677 OK Dalším způsobem georeferencování je připojení na již souřadnicově připojená data např. na jinou rastrovou či vektorovou mapu. Transformaci provedeme vyhledáním identických vlícovacích bodů nejprve v připojovaném rastru a umístěním do již připojeného rastru či vektoru. V případě síťového připojení na web je možné připojit rastr na data stažená pomocí některé z webových mapových služeb: • Web Map Service (WMS) je základní službou pro práci s mapami na Internetu. Tato služba (protokol) byla vyvinuta jako standard pro přenášení map a polohově určených dat přes Internet za účelem vytváření mapových kompozic. Tato služba umožňuje, aby si uživatel ve své aplikaci kombinoval data nejen v různých formátech, ale i data uložená lokálně, s daty poskytovanými z několika internetových mapových serverů. Tyto servery potom dokonce mohou využívat služby jiných serverů. Pro popis standardního protokolu komunikace je použit formát XML. • Web Feature Service (WFS) - rozhraní pro dotazování geografických dat prostřednictvím webu; pro přenos dat využívá GML založený na XML . UHU1 httD://212.158.143.149/index.Dhi CHOPAV Velkoplošná a maloplošná chráněná území Genové základny Semenné porosty Myslivecká zařízení Směr bořivých větrů Pásmo ohrožení imisemi ÚSES Lesní vegetační stupně Cílové hospodářské soubory Typologická mapa PLO Základní vrstva ČR Výškopis Povodí Ortofoto Rastr LHP LČR Rastr SMO CENIA httoV/aeoDortal.cenia.cz/ Podkladová topografická vrstva Barevná ortofotomapa s prostorovým rozlišer CHOPAV Geologická mapa ČR Geomorfologická mapa ČR Biosférické rezervace UNESCO NATURA 2000 USES Zvláště chráněná území - NP, CHKO, NPR, NPP, PR, PP Zonace NP a CHKO Přírodní parky Průměrná rychlost větru v 10 m nad zemí Koridory tahu čápů Podrobný seznam vrstev naleznete na: http://geoportal.cenia.cz/mapmaker/cen ia/oortal/helo/cz/sluzbv.html Otevření dat z WMS serveru Na stránkách www.uhul.cz přejděte na odkaz MAPY, kde v dolní části naleznete charakteristiku a výběr nabízených webových a mapových služeb s celou URL adresou. Otevření pomocí Add data - Look in - GIS Servers - Add WMS Server. Do horního řádku zkopírujeme pelou URL adresu OPRL (Oblastního plánu rozvoje lesa) ze stránek ÚHÚLu a pomocí Get Layers načteme dostupné mapové vrstvy, potvrdíme OK a ve složce Gis serverů se objeví WMS služba ÚHÚLu, kterou otevřeme. Po otevření se objeví varovná hláška z důvodu transformace souřadnicových systémů (data OPRL zpracována ve WGS 84)—► potvrdíme Yes a zobrazí se nám veškeré vrstvy mapového serveru (můžeme si zapnout pouze požadované vrstvy např. ortofoto + rastr SMO - černá). Pokud se při otevírání dat z WMS serveru varovná hláška neobjeví, je nutné znovu definovat souřadnicový systém projektu. Add Data Look in: Name □ igis2 □ Resultsl HUI blan35 _J Rectify) Hlinie.shJ Q polygon O ReadMJ Student ~3 jy ^IMOil üüfm ss] i Database Connections \ Address Locators Coordinate Systems £3 Interoperability Coi S Scalar References ^ Search Results (^ Toolboxes Q^, Tracking Connections Name: Show of type: | Datasets and Layers (Myr) DIKAV - ArcMap - ArcEditor -n x File Edit View Insert Selection Tools Window Help D^yiUi < Display I Source | Selection | m si* o ^ -]. ^ & Q 4 J m Iffi ^^ CD CD ■ŕ- CD __ g if ► Z: k o 0 U M r^iij 1 * * ť* 1 0# U- tf I* 0 ft A ^ X 0 Jí q 1 e n jj lil d Drawing T ^ ( I D v AT í ť? Arial "TJ |10 _^J B 1 u|,At & ^ -^ ^ * ^ -413223.11 -824660.33 Meters Sŕart ö *? & © ® 1Ä1 ± 2L m R ciipbc Po otevření datových vrstev z mapového serveru je již možné transformovat rastr. I Copy X Remove Zobrazení netransformovaného rastru kliknutím na vrstvu pomocí pravého tlačítka a volby Zoom to layer. (HD Open Attribute Table Joins and Relates Save As Lay_er File... Properties... Nyní opět nástrojem Georeferencing pomocí vlícovacích bodů přiřadíme dobře identifikovatelná místa na rastru např. křižovatku silnic, roh budovy apod, k identickým bodům na ortofotu či rastru SMO (opět nejméně 4 body). Tímto způsobem je možné souřadnicově připojit takřka libovolnou mapu v digitální podobě. Po transformaci můžeme přejít k samotné vektorizaci rastru, Export mapy (výřezu) z WMS serveru File D \3 Edit View Insert Selection Tools Windol New... Ctrl+N Open... Ctrl+O El 5j Save Save As... Save A Copy. Ctrl+S * Add Data... a Page and Print Setup. Print Preview... Print... ]fp Document Properties. Import from ArcView project. I C:\Program Files\ArcGIS.. .\jtsk.mxd Exit Alt+F4 LjIJJUJÍ MíljJ Uložit do: ÖDATA 3 Poslední doki\nenty Öbody_slp (Ělcvičení_2007-08 P^ rámeček, blk ^JSLP Název souboru: Uložit jako typ: TIFF r Po nastavení požadovaného výřezu mapy klikneme na File —* Export Map... vybere složku, kam chceme výřez uložit pojmenujeme ho a zvolíme požadovaný formát zatrhneme Word File pro vytvoření georeference uložíme