ČESKÉ VYSOKÉ UČENÍ TECHNICKÉ V PRAZE FAKULTA STAVEBNÍ
DIPLOMOVÁ PRÁCE
2015
Petra DÍFKOVÁ
ČESKÉ VYSOKÉ UČENÍ TECHNICKÉ V PRAZE FAKULTA STAVEBNÍ OBOR GEODÉZIE A KARTOGRAFIE
DIPLOMOVÁ PRÁCE Zaměření hradu Helfenburk u Úštěka a vytvoření prostorového modelu hradeb
Vedoucí práce: Ing. Tomáš Křemen, Ph.D. Katedra speciální geodézie
červen 2015
Bc. Petra DÍFKOVÁ
ABSTRAKT Tato diplomová práce se zabývá zaměřením hradu Helfenburk u Úštěka a následným zpracováním naměřených dat. Hrad byl podrobně zaměřen metodou laserového skenování. Zaměření bylo provedeno skenovacím systémem Trimble TX5. Práce se zabývá postupem při zaměření objektu, volbou stanovisek a identických bodů. Dále je zde popsán postup při zpracování pořízených mračen bodů v programech Geomagic Studio a Cyclone. Výsledkem zpracování je 3D model hradu Helfenburk ve formě mračna bodů, které je v souřadnicovém systému S-JTSK a výškovém systému Bpv. Výsledný model hradu byl dále použit pro tvorbu 3D CAD modelu pásu vnějších hradeb.
KLÍČOVÁ SLOVA Hrad Helfenburk u Úštěka, laserové skenování, mračno bodů, registrace, CAD model
ABSTRACT The master thesis describes measurement of the castle Helfenburk near Úštěk and consequent processing of the measured data. Technique of laser scanning was used for measuring of the castle. The laser scanner Trimble TX5 was used. The first part of the thesis deals with measurement of the object, choice of standpoints and choice of control points. The second part of the thesis deals with processing of measured data. Point clouds were processed in the program Geomagic Studio and Cyclone. The result of processing is 3D model of the Helfenburk castle in the form of the point cloud. The final point cloud is in the coordinate system S-JTSK and height system Bpv. The resulting model of the castle was used for creating 3D CAD model of the strip outer walls.
KEYWORDS Helfenburk castle near Úštěk, laser scanning, point cloud, registration, CAD model
PROHLÁŠENÍ Prohlašuji, že tuto diplomovou práci jsem vypracovala samostatně na základě odborných konzultací s vedoucím práce. Použitou literaturu a podklady uvádím v seznamu zdrojů.
V Praze dne . . . . . . . . . . . .
...................... (podpis autora)
PODĚKOVÁNÍ Touto cestou bych chtěla poděkovat Ing. Tomáši Křemenovi, Ph.D. za odborné konzultace a vedení této diplomové práce. Dále patří poděkování studentům, kteří spolupracovali hlavně na měřické části této práce. Jedná se o Bc. Alžbětu Prokopovou, Bc. Janu Poesovou, Bc. Martina Touška a Bc. Lukáše Vosyku. Poděkování patří také panu Ing. Tomáši Hončovi z firmy Geotronics Praha, s.r.o za půjčení skeneru Trimble TX5, bez kterého by zaměření hradu nebylo možno provést. V neposlední řadě bych chtěla poděkovat rodině a přátelům za podporu a trpělivost během celého mého studia.
Obsah Úvod ...................................................................................................................... 10 1
Hrad Helfenburk ........................................................................................... 11 1.1 Popis hradu ............................................................................................... 11 1.2 Historie hradu ........................................................................................... 13 1.3 Předchozí zaměření hradu ........................................................................ 13
2
Zaměření hradu ............................................................................................ 15 2.1 Přípravné práce ......................................................................................... 15 2.1.1
Rekognoskace .................................................................................... 15
2.1.2
Budování bodového pole................................................................... 15
2.2 Zvolená metoda měření ............................................................................ 17 2.2.1
Dělení skenovacích systémů .............................................................. 17
2.2.2
Princip metody laserového skenování .............................................. 18
2.2.3
Výhody a nevýhody laserového skenování ....................................... 19
2.3 Použitý přístroj – Trimble TX5................................................................... 19 2.4 Skenování objektu..................................................................................... 20 2.4.1
Průběh skenování .............................................................................. 20
2.4.1.1 Použité terče ............................................................................... 21 2.4.1.2 Rozmístění stanovisek ................................................................ 23 2.4.2 3
Nastavení skeneru při měření na stanovisku .................................... 23
Kontrolní měření........................................................................................... 25 3.1 Kontrolní měření v průběhu skenování .................................................... 25 3.2 Doplňující kontrolní měření ...................................................................... 25 3.2.1
Výpočet souřadnic kontrolních bodů ................................................ 26
3.3 Výsledný soubor kontrolního měření ....................................................... 27 4
Zpracování dat z laserového skenování ....................................................... 28 4.1 Export pořízených dat ............................................................................... 28 4.2 Redukce dat .............................................................................................. 29 4.2.1
Popis redukce dat .............................................................................. 29
4.3 Spojování skenů ........................................................................................ 31 4.3.1
Import dat do softwaru Cyclone........................................................ 31
4.3.2
Modelování identických bodů ........................................................... 33
4.3.3
Registrace .......................................................................................... 35
4.3.3.1 Popis registrace .......................................................................... 36 4.3.3.2 Postup registrace ........................................................................ 37 4.3.3.3 Tvorba podmínek na základě překrytu ....................................... 40 4.3.3.4 Hrubé očištění skenů .................................................................. 41 4.3.3.5 Kontrola registrovaných mračen ................................................ 42 4.4 Zhodnocení modelu na základě kontrolního měření................................ 42
5
4.4.1
Modelové souřadnice kontrolních bodů ........................................... 43
4.4.2
Zhodnocení modelu ........................................................................... 44
Transformace modelu do S-JTSK a Bpv ........................................................ 46 5.1 Sjednocení do jednoho mračna ................................................................ 47 5.2 Výsledný model hradu .............................................................................. 48
6
Drátový model hradeb.................................................................................. 50 6.1 Tvorba CAD modelu .................................................................................. 51
Závěr ...................................................................................................................... 53 Použité zdroje ........................................................................................................ 55 Použité zkratky ...................................................................................................... 56
Seznam obrázků .................................................................................................... 57 Seznam tabulek ..................................................................................................... 59 Seznam příloh ........................................................................................................ 60
Úvod Helfenburk je zřícenina hradu, která se nachází v Severočeském kraji přibližně 2,5 km od města Úštěk. Ze středověkého hradního komplexu jsou zachovány vnější hradby, věž a zdivo vnitřního paláce. O provoz a údržbu hradu se starají členové dobrovolnické skupiny Hrádek. Na základě domluvy zástupců ČVUT a členy skupiny Hrádek vznikl požadavek na tvorbu aktuální dokumentace. Dokumentace v podobě 3D modelu hradu bude sloužit jako podklad pro dobové studie a tvorbu různých plánů nebo 2D výkresů. Dokumentace je vytvářena pro celý rozsáhlý hradní komplex a je rozdělena na zaměření a zpracování naměřených dat. Zaměření hradu provedla skupina pěti studentů, zpracování dat je rozděleno mezi členy této skupiny. Hlavními tématy pro tuto diplomovou práci je zaměření hradu, zpracování dat a tvorba prostorového modelu vnějších hradeb. Pro připojení výsledného modelu hradu do souřadnicového systému S-JTSK a výškového systému Bpv bylo nutno vytvořit nové bodové pole. Jako nejvhodnější metoda pro podrobné zaměření hradu bylo zvoleno laserové skenování. Práce je členěna do několika částí. V první části je uvedeno seznámení s hradem, jeho historií a předchozím zaměřením. Dále je popsán průběh podrobného zaměření, rozmístění stanovisek, volba identických bodů, práce se skenerem a provedené kontrolní měření. V rámci přípravných prací je zde zmíněna i tvorba bodového pole. Zaměření a výpočet bodového pole je naplní jiné diplomové práce. Druhá část práce popisuje zpracování naměřených dat a tvorbu 3D modelu hradu. Zpracování naměřených dat bylo provedeno v programech Geomagic Studio a Cyclone. Popsána je příprava dat, registrace mračen bodů a transformace výsledného 3D modelu do požadovaného souřadnicového systému. Na závěr zpracování je popsán postup ověření přesnosti vytvořeného modelu. V poslední části je popsána tvorba CAD modelu pásu vnějších hradeb z výsledného 3D modelu. Zpracování CAD modelu proběhlo v programech Cyclone a Microstation. Výstupem této práce je spojené 3D mračno bodů zobrazující celý hrad a CAD model vnějších hradeb. Tyto výstupy jsou jednou z částí dokumentace skutečného stavu hradu. 10
1 Hrad Helfenburk Hradní komplex se nachází v Severočeském kraji, okres Litoměřice, katastrální území Rašovice u Kalovic. Rozkládá se na pozemkové parcele č. 776/3 a jeho věž na stavební parcele č. 81. Vlastníkem objektu je město Úštěk, ale o opravy a běžný provoz hradu se stará občanské sdružení Hrádek, které je zde v nájmu.
1.1 Popis hradu Hrad Helfenburk (Obr. 1), jinak v průběhu svých dějin nazývaný též Hrádek, Hradec nebo Hradišťko, se nachází 2,5 km jihovýchodně od města Úštěk. Ze středověkého hradu se zde zachovaly rozsáhlé zříceniny, které se rozkládají na skalnatém pískovcovém hřebeni nad Rašovickým potokem. Hrad patřil v severních Čechách mezi největší, vznikl nejspíše v druhé polovině 14. století. Je tvořen ze dvou časově i stavebně odlišných částí. Kterými jsou vnější opevnění s věží a vnitřní hrad.
Obr. 1: Letecký pohled na Helfenburk [2]
Vnější opevnění bylo postaveno později než vnitřní část hradu a je v současné době spolu s hradní věží nejzachovalejší částí hradu. Vnitřní hrad byl vystavěn na vrcholu pískovcového hřebene, který se rozkládá na třech skalních blocích. Dnes je zachováno pouze zdivo na středním skalním bloku. Celá vnitřní stavba je chráněna skalami a také samostatným pásem hradeb. 11
Vnější opevnění hradu je tvořeno mohutnými hradbami, které jsou zakončeny cimbuřím se střílnami. Hradby jsou místy vysoké až dvanáct metrů a po jejich vnitřní straně probíhá ochoz. Do hradu bylo možno vstoupit branou, ke které vedl padací most přes příkop. Napravo od této vstupní brány byla branka pouze pro pěší. Vstup do vnitřního hradu vedl odsud dále přes druhý padací most. Vstoupit do prostor hradního komplexu bylo možno také dolní bránou, která nyní umožňuje vstup do prostoru mezi vnějším opevněním a vnitřním hradem. Tyto vstupní brány jsou umístěny ve východní straně hradeb spolu s věží (Obr. 2). Hradní pískovcová věž stojící na samostatné skále je součástí opevnění. Věž má čtvercový půdorys a prošla rekonstrukcí v 19. století. Západní strana hradeb je přerušena vysokou otesanou skálou, která je součástí opevnění. U této skály je umístěna v hradbách malá vstupní (výpadní) branka vedoucí taktéž do prostoru mezi vnějším opevněním a vnitřním hradem [1].
Obr. 2: Vlevo: Bývalá hlavní vstupní brána s brankou pro pěší. Vpravo: Dolní brána, dnes sloužící jako vstupní.
12
1.2 Historie hradu Ve druhé polovině 14. století Jan z Klinštějna (z rodu Ronovců) postavil malý hrad, který pojmenoval Helfenburk. V roce 1375 Pražský arcibiskup Jan Očko z Vlašimi tento hrad odkoupil a učinil jej středem nového panství pro všechna arcibiskupská panství na pravém břehu řeky Labe. K první velké přestavbě hradu došlo v letech 1375 – 1379. Zejména byly vystavěny obvodové hradby vnitřního hradu. V letech 1390 – 1395 pokračoval s přestavbou nástupce Jana Očka Jan z Jenštejna. Nechal především postavit hlavní věž a nové vnější opevnění. V pozdějších obdobích byl hrad upravován, ale jeho základní podoba se již podstatně nezměnila. Další významnou osobou v historii Helfenburku byl poslední zde žijící arcibiskup Konrád z Vechty. Jeho působení na Hrádku spadá do období husitských válek. Na pevný Hrádek byli převezeny k uchování některé klenoty, relikvie a knihy ze Svatovítského chrámu. Roku 1421 se Konrád přiklonil k husitům, hrad přešel do rukou šlechty a vystřídal řadu majitelů. V roce 1592 koupil Hrádek Jan Sezima ze Sezimova Ústí na Úštěku, kterému bylo později zkonfiskováno jmění za účast na protihabsburském povstání. Hrad poté přešel do vlastnictví jezuitů. Počátkem třicetileté války byl hrad zcela opuštěn, během války pustošen a brzy bylo zapomenuto i jeho jméno. Počátkem 18. století zde byla zřízena myslivna a časem se zříceniny staly cílem poutníků. Roku 1839 koupil celé panství Ferdinand Lobkovic, jehož roku 1871 vystřídal textilní průmyslník Josef Schroll. V roce 1872 zde byla postavena nová hájenka, která poskytovala také pohoštění pro návštěvníky. V dalších letech byl hrad dále opravován, ale později se začal opět rozpadat [1].
1.3 Předchozí zaměření hradu První zaměření hradu a okolí v roce 1983 provedli M. Záveský a J. Krupka. Nejdříve v místním souřadnicovém systému zaměřili polygonový pořad uvnitř a kolem hradu. Měření bylo prováděno teodolitem s dálkoměrem, body byly stabilizovány pouze dočasně hřeby a dřevěnými kolíky. Poté proběhlo tachymetrické měření, jehož výsledkem byl plán vykreslený v měřítku 1:200 (Obr. 3). 13
Obr. 3: Výřez z výsledného tachymetrického plánu [3]
V druhé polovině 80. let pokračoval v zaměřování Ing. Vladimír Kotrejch. Jeho představou bylo zaměření nejen hradu, ale také stavebních úprav, jednotlivých kapes ve skále a podrobného tvaru zdiva. Nejprve bylo potřeba napojit se na původní zaměření. Část bodů polygonových pořadů již byla zničena. Na hradě se průběžně prováděly terénní úpravy, opravy zdiva a tím se body polygonu postupně ničily. Ing. Kotrejch tedy vytyčil nový polygonový pořad a začal protínáním vpřed zaměřovat terén a kapsy u studny. Jelikož nebylo dostupné dostatečné vybavení jak toto měření vhodně zobrazit, další zaměřování nepokračovalo. V roce 1988 bylo Ing. Pavlem Hlavenkou provedeno fotogrammetrické snímkování hradu. Vyhodnocování proběhlo v letech 1989 – 1990. Výsledný půdorysný plán hradu byl po domluvě připojen do místního souřadnicového systému pomocí vlícovacích bodů, které byly zaměřeny Ing. Kotrejchem [3].
14
2 Zaměření hradu Předmětem měření byl celý hrad Helfenburk. Konkrétně se jedná o vnější hradby, hradní věž, pozůstatky vnitřního paláce a bezprostřední okolí hradu. Přehled provedených prací na hradě je uveden v Tab. 1. Pro podrobné zaměření hradu byla již od začátku zvolena metoda laserového skenování. Zaměření hradu bylo rozděleno na dvě etapy. První etapa proběhla v období 22. a 23. března 2014, druhá etapa proběhla 3. a 4. května 2014. Skenování objektu provedla skupina studentů ve složení Bc. Alžběta Prokopová, Bc. Jana Poesová, Bc. Martin Toušek a Bc. Petra Dífková. Současně se skenováním hradu probíhalo zaměření bodového pole. Bodové pole má sloužit hlavně pro připojení výsledného modelu do souřadnicového systému S-JTSK a výškového systému Bpv. Zpracováním a zaměřením nového bodového pole se zabýval ve své diplomové práci Bc. Lukáš Vosyka. datum 11. 2. 2014 28. 2. 2014 22 - 23. 3. 2014 3. - 4. 5. 2014 22. 2. 2015
provedené práce rekognoskace objektu stabilizace nového bodového pole skenování objektu, zaměření bodového pole kontrolní měření Tab. 1: Rozpis provedených prací
2.1 Přípravné práce 2.1.1 Rekognoskace Rekognoskace objektu byla provedena před začátkem vlastního zaměření celého hradu. Dne 11. února 2014 proběhla prohlídka celého komplexu a také seznámení s objektem. Byl rozvrhnut postup měření a zváženy možnosti budování chybějícího bodového pole.
2.1.2 Budování bodového pole Nové bodové pole na hradě bylo budováno pro potřeby připojení výsledného modelu hradu do souřadnicového systému S-JTSK a výškového systému Bpv. Připojení modelu do těchto souřadnicových systémů je zamýšleno prostřednictvím identických 15
bodů, které budou zaměřeny z nově vybudovaného bodového pole. Nejprve bylo provedeno rozvržení umístění nových měřických bodů. Dále bylo potřeba nové měřické body stabilizovat. Tyto práce proběhli 28. února 2014. Většina bodů byla stabilizována roxory o délce 400 mm a průměru 10 mm. Roxor byl obložen kameny a zabetonován v předem vyhloubené díře. Tam, kde nebylo možno umístit roxory, byly body stabilizovány mosaznými hřeby o délce 80 mm a průměru 4 mm (Obr. 4). Tyto hřeby byly osazeny do vyvrtaného otvoru ve zdivu nebo skále a zalepeny speciálním lepidlem.
Obr. 4: Ukázka použité stabilizace
Celé bodové pole je připojeno GNSS měřením na dvou nejvyšších bodech sítě (č. 601 a 602). Jeden z těchto bodů se nachází na hradní věži a druhý na samostatném skalním bloku, který je situován na západní straně hradu. Z těchto připojovacích bodů následovalo proměření vnitřní a vnější měřické sítě. Obě měřické sítě jsou také vzájemně propojeny měřením. Více informací o hradním bodovém poli lze nalézt v diplomové práci Bc. Lukáše Vosyky [4].
16
2.2 Zvolená metoda měření Pro velký rozsah měřických prací byla zvolena metoda laserového skenování. Tato metoda je v tomto konkrétním případě výhodná pro svoji podrobnost, přesnost a rychlost zaměření celého komplexu hradu. Laserové skenování umožňuje bezkontaktní určování prostorových souřadnic. Prostřednictvím získaných souřadnic lze provádět 3D modelování a vizualizaci složitých staveb, konstrukcí, interiérů, historických a památkových objektů, atd. Nasnímaný objekt je možno pomocí softwaru zobrazit ve formě mračna bodů, ze kterého může být dále vytvořen model objektu, který je možno převést do CAD systému [5].
2.2.1 Dělení skenovacích systémů Skenovací systémy lze dělit pomocí několika kritérií. Dle umístění skeneru při měření jsou rozděleny na pozemní a letecké skenovací systémy. Pozemní skenovací systémy dělíme na statické a kinematické. Dále jsou uvedena pouze kritéria pro rozdělení pozemních statických skenovacích systémů podle: 1) Měřícího principu a. Polární skenery měření tranzitního času měření fázového rozdílu b. Skenery se základnou s jednou CCD kamerou se dvěma CCD kamerami 2) Dosahu a. Systémy s velmi krátkým dosahem D1 (0,1 – 2,0 m) b. Systémy s krátkým dosahem D2 (2 – 10 m) c. Systémy se středním dosahem D3 (10 – 100 m) d. Systémy s dlouhým dosahem D4 (100 – stovky m)
17
3) Přesnosti (přesnost je závislá na vzdálenosti, kterou je skenováno) a. Skenery s malou základnou P1 (0,01 – 1 mm) b. Skenery s velkou základnou a polární skenery s krátkým dosahem skenování P2 (0,5 – 2 mm) c. Polární skenery se středním dosahem P3 (2 – 6 mm) d. Polární skenery s dlouhým dosahem P4 (10 – 100 mm) 4) Rychlosti skenování a. Systémy s velmi vysokou rychlostí R1 (více jak 1 000 000 bodů za sekundu) b. Systémy s vysokou rychlostí R2 (100 000 – 1 000 000 bodů za sekundu) c. Systémy se střední rychlostí R3 (10 000 – 100 000 bodů za sekundu) d. Systémy s nízkou rychlostí R4 (do 10 000 bodů za sekundu) Podrobnější informace jsou uvedeny v [5].
2.2.2 Princip metody laserového skenování V případě polárního skenovacího systému jsou souřadnice bodů určovány prostorovou polární metodou, na základě měřených délek a úhlů. Měření vzdáleností je založeno na elektronickém způsobu měření. Pro určení vzdáleností skenery využívají dálkoměry s různým principem měření (impulsní, fázové). Sběr dat probíhá neselektivně, bezkontaktně a automatizovaně. Prostorové souřadnice měřeného objektu jsou v místním souřadnicovém systému skeneru a tvoří tzv. mračno bodů. Následné zpracování spočívá v očištění jednotlivých mračen bodů od nepotřebných bodů (vegetace, stojící lidé, atd.) V případě, že je objekt zachycen z více než jednoho postavení skeneru, je nutno tato mračna propojit. Propojení do jednoho mračna bodů probíhá tzv. procesem registrace. Jedná se o transformaci dat z jednotlivých stanovisek skeneru do jednoho souřadnicového systému. Registrace je realizována pomocí identických bodů nebo překrytů. Dále je možno provést transformaci mračen bodů do požadovaného souřadnicového systému [5].
18
2.2.3 Výhody a nevýhody laserového skenování Mezi hlavní výhody této metody patří:
Přesné zaměření stávajícího stavu s výrazně vyšší produktivitou práce a také finanční úsporou
Podrobné zaměření s velkou mírou zachycení detailů
Zkrácení práce v terénu s vyšší bezpečností práce
Měření může probíhat za plného provozu, případně s minimálním omezením provozu
Rychlé zpracování digitálního modelu terénu
Zpracování komplexních modelů složitých objektů [5].
Nevýhody této metody:
Časová náročnost na zpracování pořízených dat
Požadavek na počasí – přístroj není vodotěsný, při dešti může být naskenována padající voda
Značné finanční náklady na pořízení přístroje a speciálního softwaru
2.3 Použitý přístroj – Trimble TX5 Pro skenování hradu byl použit statický pozemní laserový skenovací systém Trimble TX5 (Obr. 5).
Obr. 5: Skener Trimble TX5
19
Jedná se o polární skener, který určuje měřenou vzdálenost na základě fázového rozdílu vyslaného a přijatého signálu. Vysílaný infračervený laserový svazek patří do třídy 3R. Přístroj se řadí mezi systémy s vysokou rychlostí skenování R2 a středním dosahem D3 (může sejmout až 976 000 bodů za sekundu v rozsahu až 120 m). Skener je vybaven dvouosým kompenzátorem a rozsahem zorného pole 360° x 300°. Přesnost přístroje je definována chybou určení vzdálenosti mezi přístrojem a kolmo umístěným rovinným terčem. Pro 90% odrazivost je uváděna chyba 0,95 mm na vzdálenost 25 m [6], [7]. Skener je ovládán pomocí dotykové obrazovky. Ovládací software umožňuje nastavení parametrů skenování a také prohlížení pořízených mračen bodů. Naměřená data jsou ukládána na SD kartu, kterou je do skeneru nutno vložit. Detailnější popis ovládání viz kapitola 2.4.2.
2.4 Skenování objektu Skenování hradu bylo provedeno ve dvou etapách. První část skenování proběhla 22. a 23. března 2014, druhá část proběhla 3. a 4. května 2014. Skenování pro pokrytí celého hradu bylo provedeno na 91 stanoviscích ve venkovních prostorách a 36 stanoviscích uvnitř věže.
2.4.1 Průběh skenování Zaměření celého hradu bylo rozděleno do několika bloků. Rozdílný postup byl při zaměření interiéru hradní věže a členitého exteriéru celého komplexu. V průběhu postupu souběžného zaměřování nového bodového pole byly totální stanicí zaměřeny také vybrané identické body pro potřeby tvorby modelu hradu. Přehled provedeného skenování a klimatické podmínky jsou uvedeny v Tab. 2. skenovaná oblast hradní palác (interiér, travnatá plocha) interiér věže nádvoří (vnitřní hrad) prostor mezi hradbami a vnitřním hradem, vnější část hradby
datum 22. 3. 2014 23. 3. 2014 3. 5. 2014 3. 5. 2014 4. 5. 2014
Tab. 2: Rozpis provedeného skenování
20
počasí oblačno zataženo, přeháňky oblačno slunečno slunečno
Stanoviska pro skenování bylo nutno volit s ohledem na další zpracování pořízených dat (registraci jednotlivých mračen). Mezi sousedními stanovisky bylo třeba dbát na oblast společného překrytu mračen a viditelnost několika společných terčů. Seznam pomůcek použitých při skenování: skener Trimble TX5 stativ s adaptérem pro připevnění skeneru kulové a šachovnicové terče
2.4.1.1 Použité terče Pro potřeby následné registrace jednotlivých mračen a transformaci výsledného modelu do konečných souřadnicových systémů byly v terénu umístěny kulové a šachovnicové terče. Kulové terče: Byly použity bílé koule o průměru 200 mm. Tyto terče byly umisťovány tak, aby byla zajištěna jejich viditelnost na více skenech1. K provedení registrace mračna by
měly
být
dostačující
dvě
identické
koule,
protože
pořízené
skeny
jsou již horizontovány (přistroj je vybaven dvouosým kompenzátorem). Při měření jsme měli k dispozici pouze šest těchto koulí, proto bylo nutno pečlivě rozmýšlet jejich umístění a přesouvání. V průběhu měření bylo kvůli nutnosti přemisťování koulí obtížné kontrolovat, zda sousední mračna mají alespoň dvě společné koule. Proto byla snaha mít na jednotlivých skenech zachyceny alespoň tři koule. Na vhodné rozmístění kulových terčů bylo třeba dbát také z hlediska následného propojení jednotlivých zaměřovaných bloků. Například pro propojení vnitřní a vnější části hradeb byly terče umístěny mezi jednotlivá cimbuří hradeb. V pořízených skenech bylo zachyceno 71 identických kulových terčů.
1
Termínem „sken“ je označováno mračno bodů získané z jednoho měření skeneru na stanovisku.
21
Šachovnicové terče: Využity byly šachovnicové terče určené pro použitý skenovací systém. Umístěné šachovnicové
terče
měly
sloužit
především
pro
připojení
celého
modelu
do souřadnicového systému S-JTSK a výškového systému Bpv. Tyto terče byly také využity pro tvorbu identických bodů při registraci jednotlivých mračen. Rozmístěny byly na hradbách, skalních stěnách a uvnitř paláce. Jejich umístění bylo problematické vzhledem ke stavebnímu materiálu hradu (pískovec, hrubě opracované kameny). Bylo tedy velmi obtížné terče kamkoliv připevnit na delší dobu. Zamýšlené připevnění lepicí páskou bylo nefunkční, proto byly terče umisťovány do štěrbin nebo případně zaklíněny mezi kameny. Umístění těchto terčů nebylo tak četné jako v případě umístění terčů kulových. Z vybudovaného bodového pole byly tyto šachovnicové terče zaměřeny totální stanicí pro určení souřadnic jejich středů. Zaměřeno bylo 20 umístěných terčů.
Obr. 6: Ukázka umístění použitých terčů
Přehledka znázorňující rozmístění stanovisek a terčů v průběhu měření je zařazena mezi přílohami této práce (příloha č. 1). Podkladem pro tuto přehledku je půdorysný plán z fotogrammetrického vyhodnocení (viz kapitola 1). Stanoviska skeneru jsou znázorněna černou barvou a křížkem, rozmístění kulových terčů je znázorněno zelenou barvou a tečkou. Rozmístění šachovnicových terčů je uvedeno na přehledce s kontrolními body.
22
2.4.1.2 Rozmístění stanovisek Rozmístění stanovisek pro skenování exteriéru hradu bylo voleno tak, aby na sousedních skenech bylo vždy pokud možno vidět několik společných kulových terčů a také byla zajištěna oblast dostatečného překrytu pořizovaných skenů. Vhodný překryt skenů je důležitý při následné registraci mračen na základě překrytu. Pro pokrytí celé venkovní oblasti bylo skenováno z 91 stanovisek. Hrad je výškově členitý a některé jeho části jsou obtížně dostupné. Pro bezpečný přístup měřičů a dopravu vybavení na některá stanoviska bylo potřeba použít žebřík nebo lano. V interiéru věže bylo třeba volit stanoviska tak, aby byly zajištěny dostatečné překryty jednotlivých skenů. Z důvodu stísněných podmínek zde nebylo možno umístit kulové terče. V interiéru věže bylo pořízeno 36 skenů. Registrace těchto mračen byla prováděna pouze na základě překrytů. Zpracováním modelu hradní věže se zabýval ve své diplomové práci Bc. Martin Toušek. Přehled pořízených skenů v jednotlivých částech hradu je uveden v Tab. 3. skenovaná oblast hradní palác nádvoří (vnitřní hrad) prostor mezi hradbami a vnitřním hradem vnější část hradby interiér věže celkem skenů
čísla pořízených skenů 1 - 17, 211, 201, 202 220 - 226, 203 - 210, 212 - 219
počet pořízených skenů 20 23
227, 256 - 265, 239 - 247, 251 - 255
25
228 - 238, 248 - 250, 266 - 274 101 - 136
23 36 127
Tab. 3: Přehled pořízených skenů
2.4.2 Nastavení skeneru při měření na stanovisku Práce se skenerem probíhala na každém stanovisku obdobně. Na zvoleném stanovisku byl nejprve postaven stativ s trojnožkou, do které byl skener připevněn. Dále byl skener horizontován pomocí krabicové libely, která se nachází na trojnožce. Urovnaný skener byl zapnut tlačítkem On/Off. Před spuštěním vlastního skenování na stanovisku byly nastaveny požadované parametry měření, které ovlivňují výsledný sken. Nastavované parametry byly rozdílné pro skenování venkovní části hradu a pro skenování vnitřních prostor věže. Parametry byly nastaveny v nabídce 23
Parameters – Resolution and Quality. Použité parametry při skenování jsou uvedeny v Tab. 4, náhled menu nastavovaných parametrů skeneru je zobrazen na Obr. 7. skenovaný prostor exteriér hradu interiér místností interiér schodiště
nastaveno kvalita rozlišení 1/4 3x 1/5 3x 1/8 3x
vzdálenost bodů [mm] 6,136 7,670 12,272
Tab. 4: Použité parametry skenování
Nastavení hodnot pro rozlišení a kvalitu je zde realizováno pomocí stupnic. Vhodné nastavení těchto parametrů ovlivňuje výslednou kvalitu skenování. Čím vyšší kvalita je nastavena, tím více je potlačen šum ve skenovaných datech. Mezi stupnicemi pro nastavení kvality a rozlišení jsou vypsány další údaje. Jsou zde uvedeny informace o délce skenování, velikosti pořizovaného skenu a vzájemné vzdálenosti bodů. Vzdálenost bodů udává rozestup dvou naskenovaných bodů, které se nacházejí ve vzdálenosti 10 metrů od skeneru. Velikost skenu je udávána vertikálním a horizontálním počtem bodů. Velikost rozlišení je udávána jednotkami MPts (Mega Points). Tyto jednotky udávají, kolik miliónů bodů přibližně bude naměřený sken obsahovat [6].
Obr. 7: Menu nastavení parametrů měření [6]
24
3 Kontrolní měření Kontrolní měření probíhalo již v průběhu skenování. V průběhu zpracování dat bylo ještě doplněno dalším kontrolním měřením. Celé kontrolní měření bylo provedeno hlavně za účelem zhodnocení přesnosti výsledného modelu hradu. Zhodnocení přesnosti bylo provedeno porovnáním prostorových délek získaných měřením a odečtených z výsledného modelu hradu. Zhodnocení přesnosti modelu je podrobně uvedeno v kapitole 4.4. Dalším účelem měření bylo určit identické body pro transformaci modelu do globálního souřadnicového systému.
3.1 Kontrolní měření v průběhu skenování Některé vybrané kontrolní body byly zaměřeny dříve v průběhu skenování hradu. Jejich zaměření proběhlo v rámci proměřování bodového pole a zaměření středů šachovnicových terčů, které byly rovněž zařazeny mezi kontrolní body. Zaměřeno bylo 10 vybraných bodů a 20 středů terčů. Terče představují jediné kontrolní body z vnější části hradeb. Souřadnice těchto kontrolních bodů a středů terčů jsou jedním z výstupů diplomové práce Bc. Lukáše Vosyky [4].
3.2 Doplňující kontrolní měření Dne 22. 2. 2015 proběhlo kontrolní měření za účelem určení prostorových délek napříč celým hradem. Prostorové délky budou určeny ze souřadnic kontrolních bodů. Kontrolní měření provedly studentky Bc. Alžběta Prokopová, Bc. Jana Poesová a Bc. Petra Dífková. Z hradního bodového pole byly polární metodou zaměřeny vybrané kontrolní body. Kontrolní body byly voleny na nepohyblivých a jasně identifikovatelných místech. Voleny byly např. v rozích oken, klenbě brány, na hranách cimbuří, na výrazných spárách zdiva a výčnělcích kamenů zdiva. Pro pozdější usnadnění identifikace ve zpracovávaném digitálním modelu byla při měření pořizována fotodokumentace zaměřovaných bodů. Na pořizovaných fotografiích byl zaměřovaný bod signalizován pomocí laserové stopy totální stanice (Obr. 8). 25
Zaměřeno bylo celkem 72 kontrolních bodů, které byly rozloženy tak, aby z nich bylo možno realizovat prostorové délky, které budou napříč celým hradem. Tyto délky mají za úkol ověřit, zda vytvářený model hradu odpovídá realitě a není zdeformovaný.
Obr. 8: Ukázka pořízené fotodokumentace kontrolních bodů
Seznam pomůcek použitých při kontrolním měření: totální stanice Trimble M3 (DR5'') v. č. 652352 dřevěný stativ 2x hranol Topcon s výtyčkou 2x hranol mini Leica se stojánkem pásmo teploměr, barometr
3.2.1 Výpočet souřadnic kontrolních bodů Souřadnice bodů bodového pole jsou převzaty z diplomové práce Bc. Lukáše Vosyky [4]. Výpočet souřadnic z námi provedeného kontrolního měření proběhl v programu Groma 8. Body bodového pole byly v rámci kontrolního měření číslovány odlišně. Souřadnice použitých bodů bodového pole a jejich číslování je uvedeno v tabulce, která je zařazena v přílohách této práce (příloha č. 3). Do výpočtu souřadnic kontrolních bodů nebyla zavedena oprava z nadmořské výšky ani kartografického zkreslení. Použitý měřítkový koeficient je roven jedné. Vypočtené souřadnice budou použity pro výpočet prostorových délek, které budou porovnávány s reálným modelem, tudíž je zavádění oprav nežádoucí. Seznam souřadnic 26
vypočtených bodů je zařazen mezi přílohy této práce (příloha č. 4). Protokol o výpočtu je elektronickou přílohou této práce (složka č. 1: Kontrolní měření).
3.3 Výsledný soubor kontrolního měření Z provedených kontrolních měření máme k dispozici 79 kontrolních bodů a 20 zaměřených středů šachovnicových terčů. Pro snazší identifikaci kontrolních bodů byla vyhotovena přehledka s jejich umístěním, tabulka s jejich popisem a
fotodokumentace.
Podkladem
pro
přehledku
je
půdorysný
plán
hradu
z fotogrammetrického vyhodnocení (viz kapitola 1). V tabulce je podrobně popsáno jejich umístění v terénu, a zda je k dispozici fotografie zvoleného kontrolního bodu. V přehledce jsou zeleně vyznačeny body z průběžného kontrolního měření a červeně kontrolní body z měření doplňujícího. Přehledka je zařazena mezi přílohy této práce (příloha č. 2) Tabulka s popisem kontrolních bodů a fotodokumentace je mezi elektronickými přílohami této práce (složka č. 1: Kontrolní měření). V Tab. 5 je ukázka části přehledové tabulky s popisem kontrolních bodů. Číslo bodu
popis
umístění
foto
901
lom oblouku brány
bývalá vstupní brána, východní strana hradeb
ano
902
levý roh kamenného rámu brány
bývalá vstupní brána, východní strana hradeb
ano
903
pravý roh kamenného rámu brány
bývalá vstupní brána, východní strana hradeb
ano
904
střed horního křížení dřevěného zábradlí
bývalá vstupní brána, východní strana hradeb
ano
905
roh prvního pravého kamene od lomu oblouku brány
vstupní brána, východní strana hradeb
ano
Tab. 5: Ukázka z přehledové tabulky
27
4 Zpracování dat z laserového skenování Zpracování pořízených dat bylo rozděleno mezi studenty, kteří se podíleli na měření celého objektu. Bc. Martin Toušek provedl zpracování hradní věže nezávisle na finálním modelu celého hradu. Studentky Bc. Alžběta Prokopová, Bc. Jana Poesová a Bc. Petra Dífková provedly zpracování modelu celého hradu. Společný model je dále využit k dalšímu rozdílnému zpracování. Zpracování pořízených dat vzhledem k jejich objemu proběhlo ve speciálních programech pro laserové skenování. Tyto programy jsou vytvářeny pro zpracování velkých objemů dat, které běžné grafické programy nejsou schopny zpracovat a zobrazit. Konkrétně byly využity programy Geomagic Studio a Leica Cyclone. Prováděné operace v těchto programech jsou náročné na výkon počítače. Proto zpracování probíhalo na výkonném počítači, který nám poskytla Katedra speciální geodézie (Tab. 6). operační systém operační paměť procesor grafická karta
Windows 7 Enterprise 64 GB RAM Intel Core i7-3820 CPU @ 3,60 GHz AMD Radeon HD 7900 Series Tab. 6: Konfigurace školního PC
Následující text se zabývá postupem zpracování dat při tvorbě prostorového modelu hradu. Nejprve je popsána příprava pořízených dat před započetím vytváření modelu hradu. Jedná se o export dat ze skeneru a jejich import do softwaru, kde probíhá zpracování modelu. Dále je popisována registrace jednotlivých mračen bodů do jednoho celku. Nakonec je popsána transformace výsledného modelu do souřadnicového systému S-JTSK a výškového systému Bpv.
4.1 Export pořízených dat Pořízená data přístrojem Trimble TX5 jsou ukládána na vloženou SD kartu. Stažená data z SD karty lze otevřít pouze v programu Scene, který je ke skeneru dodáván. V programu Scene byl proveden export dat do formátu PTX. Tento formát ukládá sedm údajů o každém bodě, který byl naskenován. Ukládány byly souřadnice X, Y, Z [m], intenzita [0,1] a hodnoty RGB [0, 255]. 28
4.2 Redukce dat Skenování dat proběhlo s hustotou odpovídající rozestupu bodů 6,1 mm na vzdálenost 10 m. Hustota skenovaných bodů je závislá na vzdálenosti skenovaných předmětů od skeneru. Blízké předměty byly skenovány s nadbytečnou hustotou. Proto bylo možno před vlastním zpracováním data redukovat. Touto redukcí došlo ke snížení objemu dat a tím i k usnadnění dalšího zpracování. Za vhodnou hustotu, na kterou byla mračna redukována, byl zvolen rozestup mezi body 10 mm. Procesem redukce se celkový objem dat zmenšil ze 157 GB na 25 GB. Redukce dat proběhla v programu Geomagic Studio, který umožňuje dávkové zpracování mračen bodů.
4.2.1 Popis redukce dat Pro usnadnění práce byl proces redukce dat automatizován pomocí dávkového zpracování (Batch Processing). Redukováno bylo automaticky všech 127 pořízených mračen bodů. Pro tento proces bylo nutno vytvořit tzv. Macro, pomocí kterého byl nahrán celý postup ředění bodů. Nastavení makra proběhlo:
spuštění nahrávaní makra: Tools – Macro – Record (zapíší se všechny provedené operace)
ukončení nahrávání makra: Tools – Macro – Stop.
V rámci nahrávání makra bylo použito:
Funkce pro redukci počtu bodů mračna: Points – Sample – Uniform (Obr. 9) rozestup bodů byl nastaven na 10 mm nastavení maximální priority křivosti (Curvature Priority - Max) pro zachování bodů tam, kde je povrch zakřiven a je potřeba pro jeho popsání více bodů než v případě rovné plochy.
29
Obr. 9: Nastavení funkce Uniform
Nahrané makro bylo spuštěno funkcí Tools – Advanced – Batch processing. Pro automatický proces byly nastaveny tyto parametry (Obr. 10):
adresář kde jsou uloženy skeny, které je třeba zpracovat
název makra, podle kterého bude proces pobíhat: Macro Name
adresář, kam budou ukládána zpracovaná mračna
formát, ve kterém budou zpracovaná mračna ukládána (*.vtx) tento formát ukládá souřadnice bodu X, Y, Z [m], hodnoty RGB [0, 255] a normály X, Y, Z.
Obr. 10: Nastavení funkce Batch processing
30
4.3 Spojování skenů Další zpracování mračen bodů probíhalo v programu Cyclone. Do výsledného modelu bylo zapracováváno 91 mračen bodů, které zachycují celý hrad kromě interiéru věže. Pro vytvoření celistvého modelu hradu je nutno jednotlivá mračna spojit do jednoho celku. Mračna byla spojována pomocí identických bodů nebo na základě překrytů.
4.3.1 Import dat do softwaru Cyclone Aby bylo možno nahrát redukovaná mračna do programu Cyclone bylo potřeba přepsat koncovku jednotlivých souborů z VTX na TXT. Při spuštění programu je otevřeno okno „Cyclone Navigator“, který zobrazuje nabídku serverů. V těchto serverech jsou ukládány databáze, které obsahují soubory s naměřenými daty. Pro práci byl využit server „ZEPPELIN – PC (unshared)“. V rámci tohoto serveru byly vytvořeny dvě pracovní databáze s názvy „helfenburg“ a
„helf_palac“.
Databázi
lze
vytvořit
v
nabídce
Databases – Configure Databases – Add (Obr. 11).
Obr. 11: Zakládání nové databáze
31
pravého
tlačítka
myši:
Do první databáze „helf_palac“ bylo importováno 17 mračen, která zobrazují pouze vnitřní hrad. Do druhé databáze „helfenburg“ bylo importováno zbylých 74 mračen, která zobrazující hradby, exteriér věže, nádvoří a část exteriéru vnitřního hradu. Import byl proveden kliknutím pravým tlačítkem na příslušnou databázi a zvolením možnosti Import. Na Obr. 12 je zobrazeno dialogové okno importu. V záložce Options je třeba nastavit správné hodnoty rozsahu barev a intenzity, aby se importované mračno bodů dalo vhodně zobrazit. Dále je potřeba nastavit kolik sloupců má použitý formát dat importovaného souboru, případně kolik řádků hlavičky souboru je třeba přeskočit. Pro každý importovaný sloupec je potřeba nastavit, který druh dat z použitého formátu VTX obsahuje (souřadnice X, Y, Z; atd.).
Obr. 12: Nastavení importu dat Cyclone
32
Pro každé importované mračno se vytvoří samostatný tzv. ScanWorld (Obr. 13), který obsahuje [8]: ControlSpace (objekty vstupující do registrace mračen) ModelSpace (informace o objektech, změny v databázi ScanWorldu) ModelSpace View (pracovní prostředí, kde probíhají změny mračna a modelování) Scans (údaje o naskenovaném mračnu bodů) Images (složka obsahující případné fotografie pořízené skenerem v průběhu měření)
Obr. 13: Struktura Scanworldu
4.3.2 Modelování identických bodů Ve všech skenech bylo nutno postupně vyhledat identické body. Identické body představují naskenované kulové terče, které byly rozmisťovány během skenování. Z jednotlivého stanoviska byla vždy naskenována pouze viditelná část povrchu kulového terče. Tudíž z různých stanovisek je naskenována jiná část, jejíž body nejsou mezi stanovisky totožné. Proto bylo nutno těmito body proložit kouli a identickým bodem byl zvolen střed vložené koule. Pro vkládané koule byl zafixován jejich průměr na 200 mm, což je hodnota udávaná výrobcem. K tomuto nastavení slouží funkce Edit – Object Preferences – Fit Diameter (Obr. 14).
33
Obr. 14: Nastavení průměru koule
Modelování koulí probíhalo tak, že byl v mračnu vybrán jeden z bodů příslušející kouli, která zde má být proložena. Použila se funkce Create Object – Region Grow – Sphere. Funkce vybere automaticky další body z mračna, které kouli představují (Obr. 15). Pomocí nastavení Region Size je možno upravit velikost vybrané oblasti, možnosti nastavení jsou uvedeny na Obr. 16.
Obr. 15: Vlevo: koule viditelná v mračnu. Uprostřed: koule vybraná funkcí Region Grow. Vpravo: vymodelovaná koule
34
Obr. 16: Možnosti funkce Region Grow
Do středu takto vymodelované koule byl vložen identický bod tzv. vertex, který je potřebný k registraci. Pro vložení identického bodu byla využita funkce Create Object – Insert – Vertex. Každý vložený identický bod byl očíslován funkcí Tools – Registration – Add/Edit Registration Label. Takto vložené a označené body dokáže program při registraci sám rozpoznat jako body identické. Identické body reprezentující střed jedné koule mají v jednotlivých mračnech shodné číslo. Veškeré modelování probíhalo v prostředí ModelSpace. Proto je na závěr potřeba připravené mračno označit a vložit do ControlSpace. Vkládané identické body a jejich označení se do ControlSpace vkládají automaticky. Toto je nutné provést, jelikož program při registraci mračen pracuje pouze s daty, která jsou v ControlSpace.
4.3.3 Registrace Během registrace jsou jednotlivá mračna zachycující skenovaný objekt transformována do jednoho společného souřadnicového systému. Výsledkem registrace bude tedy ucelený model tohoto objektu. Výpočet transformačního klíče probíhá z dvojích souřadnic identických bodů nebo na základě překrytu [5]. 35
Celý proces registrace hradu byl pro značný počet zpracovávaných mračen bodů rozdělen do tří navazujících bloků.
4.3.3.1 Popis registrace Zpracování jednotlivých registrací probíhá obdobně. V dané databázi je v menu nejprve vytvořena registrace zvolením možnosti Create – Registration. Dále pracujeme v otevřeném okně vytvořené registrace. Postupně jsou prováděny tyto kroky:
Pomocí funkce ScanWorld – Add ScanWorld jsou vybrány skeny, které budou v dané registraci spojeny.
Nastavení stanoviska, do jehož souřadnicové soustavy budou ostatní mračna transformována. Nastavení se provádí funkcí ScanWorld – Set
Home
ScanWorld
(domácí
ScanWorld).
Automaticky
je za Home ScanWorld vybráno první mračno a je zvýrazněno tučně.
Kontrola mračen, zda jsou označena jako horizontovaná. Pokud je mračno takto označeno, za jeho názvem je uvedeno „Leveled“. Zpracovávaná mračna jsou již horizontovaná (skener byl vybaven dvouosým kompenzátorem). Toto správné označení je důležité pro výsledky registrace.
Vyhledání identických bodů a jejich vzájemné propojení se provede automaticky funkcí Constraint – Auto-Add Constraints.
Kontrolu, odstranění a případnou editaci vytvořených podmínek lze provést v záložce Constraint List. Je možno zde také nastavit váhu, se kterou podmínka pro identický bod bude vstupovat do výpočtu. V případě nedostatečného počtu podmínek, lze doplnit podmínky přidáním „překrytu mračen“ (podrobněji v kapitole 4.3.3.3).
Výpočet registrace se provede funkcí Registration – Register.
Kontrola spojení mračen. Pomocí funkce Registration – View Interim Result se zobrazí náhled spojení registrovaných mračen a lze provést vizuální kontrolu. Podrobněji o prováděných kontrolách v kapitole 4.3.3.5.
Ukončení registrace. Před spuštěním výpočtu je nejprve potřeba v záložce ModelSpaces označit ControlSpace jednotlivých mračen, která jsou 36
registrována. Funkcí Registration – Create ScanWorld/Freeze Registration je vytvořen ModelSpace, který obsahuje registrovaná mračna. Skeny zde nejsou spojeny do jednoho mračna, ale jsou stále reprezentovány jednotlivými mračny.
Zobrazení a uložení protokolu o provedené registraci (Registration – Show Diagnostics). V protokolu jsou vypsány dosažené odchylky na jednotlivých identických bodech, celková průměrná odchylka registrace a transformační klíč pro jednotlivá mračna (Mean Absolute Error for Enabled Constraints).
Většinu popsaných důležitých funkcí je možno spustit pomocí ikon zobrazených v pásu pod jednotlivými záložkami funkcí. Nejčastěji používané funkce jsou zvýrazněny na Obr. 17.
Obr. 17: Pracovní okno registrace
4.3.3.2 Postup registrace První blok registrace: V prvním bloku bylo vytvořeno deset registrací, které byly vždy tvořeny skupinou jednotlivých na sebe navazujících mračen. Tyto skupiny byly voleny tak, aby vybrané skeny bylo možno spojit pomocí vymodelovaných identických bodů. Registrace byly vytvářeny v databázích „helfenburg“ a „helf_palac“ dle příslušnosti registrovaných mračen.
37
Podmínky pro vzájemné propojení registrovaných mračen byly u většiny vytvářeny pouze pomocí vymodelovaných identických bodů. Váhy u všech těchto identických bodů byly nastaveny na hodnotu 1. Přehled spojených mračen v těchto registracích a jejich dosažené průměrné odchylky jsou uvedeny v Tab. 7. U dvou problémových registrací (označeny I, J) bylo potřeba pro spojení mračen rozšířit podmínky. U těchto registrací se ukázalo, že podmínky tvořené pouze dostupnými identickými body jsou zde nedostačující. První problém byl u registrace „I“, která je tvořena mračny v prostoru východní části hradeb a věže. Při provedení kontrolních řezů bylo zjištěno, že vzájemná vzdálenost registrovaných mračen v prostoru věže je větší než 6 cm. Proto byla v mračnech doplněna ploška vymodelovaná na jednom z kvádrů věže. Střed této plošky vstoupil do registrace jako další identický bod. Poté již vzdálenost v kontrolním řezu mezi registrovanými mračny klesla na zhruba 1,5 cm. Více tato chyba nebyla řešena, jelikož kompletní model věže zpracoval nezávisle Bc. Martin Toušek v rámci své diplomové práce. Registrace „J“ je tvořena skeny hradeb v prostoru severní strany hradu. Při prvotní registraci bylo zjištěno, že registrované skeny na sebe nesprávně navazují. Problém byl vyřešen ořezáním jednotlivých mračen. Ořezání problémových částí hradeb bylo provedeno od postavení skeneru zhruba 15 m dále na obě strany. Tyto ořezané vzdálenější body od skeneru byly vyhodnoceny jako chybně naskenované pod příliš ostrým úhlem dopadu. Do nové registrace vstoupila ořezaná mračna a mezi podmínky bylo přidáno pět vzájemných překrytů. Váhy pro identické body byly nastaveny na hodnotu 1 a pro překryty na hodnotu 0,7. Tyto hodnoty pro nastavení vah, byly zvoleny na základě výsledků několika zkušebních výpočtů registrace. Dále byl vyloučen z podmínek identický bod s číslem 320, který byl zhodnocen jako nespolehlivý. Tento bod představuje koule, která byla umístěna v terénu na vrstvě listí. Tato koule nejspíše v průběhu pořizování skenů změnila svoji polohu oproti původnímu postavení.
38
registrace označení název A reg1-10 B reg13-17 C reg205-210 D reg214-217 E reg219-223 F reg251-254 G reg236-238 H reg256-258 I reg229-233 J reg269-274
spojená mračna
průměrná absolutní odchylka [m]
1 - 10 13 - 17 205 - 210 214, 215, 217 219 -223 251 - 254 236 - 238 256 - 258 229 -233 269 - 274
0,0023 0,0020 0,0032 0,0024 0,0012 0,0028 0,0026 0,0013 0,0083 0,0118
Tab. 7: Přehled registrovaných mračen v prvním bloku
Druhý blok registrace: V druhém bloku bylo vytvořeno osm registrací, které byly nazvány „doregistrace“. Tyto doregistrace vznikly připojením zbylých jednotlivých mračen k již existujícím skupinám. Vazby v doregistracích jsou tvořeny kombinací překrytů a identických bodů. Pro tvorbu těchto registrací byla založena nová pracovní databáze „helf_vse“. Před tvorbou těchto registrací byl vytvořen rozvrh pro přiřazení volných skenů k již vytvořeným skupinám. V těchto registracích byly nejčastěji tvořeny vazby pomocí překrytů. Pro dosažení správnosti výpočtu bylo nejprve přistoupeno k očištění jednotlivých mračen i již vytvořených skupin od přebytečných bodů (více v kapitole 4.3.3.4). Při tvorbě registrací byl kladen důraz na vzájemné propojení mezi stávající skupinou skenů a připojovanými skeny. Váhy pro jednotlivé typy vazeb byly v doregistracích voleny dvěma způsoby. Když vazby byly tvořeny pouze překryty, váha pro překryt byla volena na hodnotu 1. V případě, že byly vazby tvořeny identickými body a překryty, identické body měly váhu 1 a překryty 0,7. Přehled spojených mračen v těchto registracích a jejich dosažené průměrné odchylky jsou uvedeny v Tab. 8.
39
doregistrace označení název K doregistrace13-17 L doregistrace205-210 M doregistrace214-217 N doregistrace219-223 O doregistrace251-254 P doregistrace236-238 Q doregistrace256-258 R doregistrace229-233
průměrná absolutní odchylka [m] B, 201, 202 0,0021 C, 204, 211 - 213, 218 0,0031 D, 216 0,0003 E, 225 - 227, 203 0,0024 F, 224, 239 - 244, 255 0,0022 G, 235, 247 - 250, 266 - 268 0,0017 H, 245, 246, 261 - 265, 259, 260 0,0028 I, 228, 234 0,0036 spojené registrace a mračna
Tab. 8: Přehled registrovaných mračen v druhém bloku
Třetí blok registrace: Ve třetím bloku vznikla poslední registrace modelu hradu, která spojila jednotlivé části do jednoho celku. Tato registrace byla provedena v databázi s názvem „helf_vse“. Důležité zde bylo vzájemně propojit jednotlivé části modelu vytvořené předchozími registracemi. Vazby mezi jednotlivými částmi byly vytvořeny kombinací překrytů a identických bodů. Váhy pro identické body a překryty byly voleny shodně, jako v případě registrací v druhém bloku. Přehled spojených registrací a dosažená průměrná odchylka je uvedena v Tab. 9. finální registrace označení název S registrace_fin
spojené registrace
průměrná absolutní odchylka [m]
A, J, K, L, M, N, O, P, Q, R
0,0107
Tab. 9: Přehled registrace třetího bloku
Veškeré protokoly z provedených registracích jsou zařazeny mezi elektronické přílohy této práce (složka č. 2: Protokoly z provedených registrací).
4.3.3.3 Tvorba podmínek na základě překrytu Při tomto způsobu spojování mračen je využita pro určení transformačního klíče metoda minimální vzdálenosti povrchů v překrytových oblastech. Oproti využití identických bodů vykazuje nižší přesnost a v některých případech selhává. Pro výpočet
40
transformačních rovnic je využíván ICP algoritmus (Iterative Closest Point). Počáteční hodnoty rotace a translace jsou zadány uživatelem [5]. Podmínky tvořené zadáním překrytu je v okně registrace možno zadat funkcí Cloud Constraint – Cloud Constraints Wizard. Spuštěním této funkce se nejprve otevře okno, ve kterém je potřeba vybrat dvojice skenů, mezi kterými bude překryt zadáván. Po odsouhlasení výběru jsou zobrazena dvě pracovní okna, ve kterých jsou mračna vybrané dvojice. V těchto zobrazených mračnech je potřeba vybrat alespoň dva páry bodů, které představují přibližné souřadnice pro začátek výpočtu. Výběr potvrdíme funkcí Preview, která nám zobrazí náhled spojení mračen pomocí zadaných bodů. Když jsme s náhledem spokojeni, tak ho uzavřeme a podmínku uložíme funkcí Constraint. Zadaná vazba se nám zobrazí v seznamu podmínek (Constraint List). Náhled pracovního okna při zadávání párů bodů je na Obr. 18. Takto vytvořenou podmínku lze případně optimalizovat v nabídce pravého tlačítka myši: Cloud Constraint – Optimize Cloud Alignment.
Obr. 18: Zadávání překrytu v programu Cyclone
4.3.3.4 Hrubé očištění skenů Během skenování byly naskenovány i jiné předměty než součásti hradu a bezprostřední terén. Proto bylo třeba z jednotlivých mračen odstranit přebytečnou
41
vegetaci, stojící osoby, měřické vybavení použité během zaměřování bodového pole apod. Tyto objekty nejsou nehybné, a na jednotlivých skenech mohou být zachyceny v jiné pozici nebo v rozdílném postavení. Přebytečné body v mračnu byly nejprve vybrány funkcí Polygonal Fence Mode a následně vymazány příkazem Fence – Delete inside.
4.3.3.5 Kontrola registrovaných mračen Při tvorbě každé registrace proběhla vždy průběžná vizuální kontrola. Kontrola byla prováděna v ModelSpace, který je vytvořen funkcí Registration – View Interim Result. Kontrola spočívala v provedení tenkých vodorovných a svislých řezů. Vždy bylo provedeno několik vybraných řezů pro pokrytí celé oblasti spojovaných mračen bodů. Řezy byly tvořeny pomocí funkce Limit Box nebo kopií výřezu mračna do nového ModelSpace (Fence – Copy Fenced to New ModelSpace). V každém řezu bylo zjišťováno, zda jsou mračna spojována správně a nevznikají mezi nimi mezery. Pomocí výběru bodů ze zkoumaných mračen lze měřit jejich vzájemnou vzdálenost (Tools – Measure – Distance – Point to Point). Na Obr. 19 je znázorněn kontrolní řez, ve kterém můžeme vidět nedoléhající mračna a jejich vzdálenost.
Obr. 19: Nedoléhající mračna v kontrolním řezu
4.4 Zhodnocení modelu na základě kontrolního měření Ověření
přesnosti
vytvořeného
3D
modelu
(registrovaných
mračen)
bylo provedeno pomocí porovnání kontrolních délek. Porovnány byly odpovídající
42
si délky získané měřením a odečtené z vytvořeného modelu. Tyto délky mají za hlavní úkol ověřit, zda vytvořený model hradu odpovídá realitě a není deformován. Jedna sada kontrolních délek je určena ze souřadnic bodů získaných kontrolním měřením (kapitola 3). Druhá sada délek je určena ze souřadnic kontrolních bodů odečtených ve vytvořeném modelu.
4.4.1 Modelové souřadnice kontrolních bodů Jak již bylo popsáno v kapitole 3, k dispozici jsou dva druhy kontrolních bodů. Jedná se o body na rozích, hranách kamenů hradu a jiných stabilních prvcích. Druhým typem kontrolních bodů jsou středy šachovnicových terčů. Kontrolní body jsou vyhodnoceny ve výsledném modelu. Středy šachovnicových terčů: Ve výsledném mračnu finální registrace byly ve středech naskenovaných terčů vytvořeny identické body s příslušnými čísly identických bodů (Obr. 20). Čísla bodů byla navýšena o konstantu 9000, aby identické body nebyly zaměněny s body, které byly takto vkládány do středů modelovaných koulí pro registraci mračen. Tyto identické body mají modelové souřadnice a jsou v samostatné vrstvě.
Obr. 20: Naskenovaný terč v mračnu
Rohy a hrany kamenů: Do mračna v místech druhého typu kontrolních bodů byly taktéž vkládány identické body (Obr. 21). Vzhledem k tomu, že stavebním materiálem hradu je pískovec a nejsou zde dostupné ostré a snadno identifikovatelné hrany, bylo tyto body v mračnu obtížnější vyhodnotit. Čísla těchto kontrolních bodů byla navýšena o konstantu 900, aby taktéž nebyly zaměněny s již vloženými body. 43
V mračnu finálního modelu bylo vyhodnoceno 66 kontrolních bodů a 19 středů šachovnicových terčů. Některé body, které bylo obtížné podle dokumentace vyhodnotit, byly vyřazeny. Vyhodnocené souřadnice kontrolních bodů byly pro další zpracování vyexportovány do textového souboru.
Obr. 21: Kontrolní bod na hradbách
4.4.2 Zhodnocení modelu Zhodnocení modelu popisuje podrobně ve své diplomové práci Bc. Alžběta Prokopová. Zde je uveden pouze postup a závěr zhodnocení. Z modelových a naměřených souřadnic kontrolních bodů byly v programu Matlab vypočteny prostorové vzdálenosti mezi kontrolními body. Pro každý soubor délek byla vytvořena matice obsahující kombinaci délek mezi všemi kontrolními body. Odečtením těchto matic byly určeny rozdíly mezi délkami odpovídající modelu a realitě. Vypočtené rozdíly vzdáleností byly porovnány s mezními odchylkami, které byly stanoveny na základě úvahy o přesnosti. Hodnoty mezních odchylek rozdílu délky pro oba typy kontrolních bodů jsou následující: Pro šachovnicové terče:
𝛿∆𝑑1 = 20 𝑚𝑚
Pro hrany kamenů:
𝛿∆𝑑2 = 40 𝑚𝑚
Pro šachovnicový terč a hranu kamenu:
𝛿∆𝑑3 = 32 𝑚𝑚
Tabulka s rozdíly všech délek je velmi rozsáhlá, proto je zařazena mezi elektronické přílohy této práce (složka č. 3: Rozdíly délek). Na základě dosažených rozdílů délek
44
a jejich porovnání s mezními odchylkami lze usoudit, že výsledný model hradu není nijak deformován. Pro zhodnocení přesnosti modelu byly vybrány nezávislé délky. Zvoleny byly nejdelší délky mezi šachovnicovými terči ve směru podélném, příčném a osy Z. Vybrané délky jsou uvedeny v následujících tabulkách (Tab. 10, Tab. 11, Tab. 12) a jsou také znázorněny na přehledce, která je mezi přílohami této práce (příloha č. 5). Dále byla vypočtena z rozdílů vybraných vzdáleností směrodatná odchylka rozdílů délek (𝜎𝑟 ): 𝜎𝑟 = √ kde:
d
rozdíl délek
n
počet délek
∑ ∆𝑑 2 𝑛
Hodnota směrodatné odchylky rozdílů vybraných délek je 15 mm [9]. d [m] Δd [m] 81,485 -0,021 100,479 -0,018 92,902 -0,012 91,421 -0,021
Číslo bodu Číslo bodu 9112 9118 9111 9125 9108 9120 9127 9122
Tab. 10: Délky v podélném směru
d [m] 60,361 47,442 24,151 23,414
Číslo bodu Číslo bodu 9112 9127 9107 9103 9116 9102 9118 9125
Δd [m] 0,015 -0,017 0,004 -0,004
Tab. 11: Délky v příčném směru
Číslo bodu Číslo bodu 9101 9001 9127 9103 9122 9003
d [m] 8,370 42,957 45,758
Δd [m] -0,007 0,002 -0,023
Tab. 12: Délky ve směru osy Z
45
|ΔZ| 6,934 8,713 17,305
5 Transformace modelu do S-JTSK a Bpv Výsledný model hradu byl transformován do souřadnicového systému S-JTSK a výškového systému Bpv pomocí identických bodů. Shodnostní 3D transformace proběhla v programu Cyclone prostřednictvím nové registrace. Dále byla k výslednému modelu připojena dvě zapomenutá mračna s čísly 11 a 12. Nejprve byl vytvořen nový ScanWorld „SJTSK_Lukas“, do kterého bylo importováno vybraných 18 identických bodů. Tyto identické body jsou zaměřené středy šachovnicových terčů, které byly rozmisťovány po hradu během skenování. Importované souřadnice bodů bylo potřeba nejprve upravit, jelikož program pracuje s kartézským souřadnicovým systémem. Souřadnice X a Y v S-JTSK bylo potřeba násobit hodnotou -1 a vzájemně zaměnit. Po úpravě máme tedy x = - Y, y = - X a z = Z. Upravené souřadnice identických bodů byly pro import uloženy do nového textového souboru. Přehled identických bodů a jejich použitých souřadnic je uveden v Tab. 13. Číslo bodu 9001 9003 9101 9102 9103 9106 9107 9108 9111 9112 9114 9116 9117 9118 9120 9122 9125 9127
Y [m] 737990,478 737990,084 737985,791 737992,972 737976,591 737983,646 737989,635 737967,281 737952,549 737982,424 737998,989 738014,235 738016,279 738020,213 738020,863 738015,617 738000,992 737936,849
S-JTSK X [m] 988443,537 988438,428 988443,514 988429,258 988437,188 988467,783 988482,621 988482,760 988492,464 988490,538 988468,007 988436,192 988432,669 988418,384 988407,093 988404,628 988405,432 988450,977
Z [m] 327,681 325,873 320,747 319,284 318,461 316,868 314,607 314,709 319,156 310,867 310,685 310,159 309,872 309,388 309,055 308,568 306,069 309,748
x [m] -737990,478 -737990,084 -737985,791 -737992,972 -737976,591 -737983,646 -737989,635 -737967,281 -737952,549 -737982,424 -737998,989 -738014,235 -738016,279 -738020,213 -738020,863 -738015,617 -738000,992 -737936,849
cyclone y [m] -988443,537 -988438,428 -988443,514 -988429,258 -988437,188 -988467,783 -988482,621 -988482,760 -988492,464 -988490,538 -988468,007 -988436,192 -988432,669 -988418,384 -988407,093 -988404,628 -988405,432 -988450,977
Tab. 13: Souřadnice použitých identických bodů
46
z [m] 327,681 325,873 320,747 319,284 318,461 316,868 314,607 314,709 319,156 310,867 310,685 310,159 309,872 309,388 309,055 308,568 306,069 309,748
V ModelSpace finální registrace mračen („registrace_fin“) byly zkontrolovány již vložené identické body v místech naskenovaných terčů. Tyto identické body mají modelové souřadnice a jsou v samostatné vrstvě. Čísla bodů jsou navýšena o konstantu 9000, aby identické body nebyly zaměněny s body, které byly takto vkládány do středů modelovaných koulí pro registraci mračen. Po těchto přípravách byla vytvořena nová registrace „SJTSK“. Do registrace byl vložen ScanWorld s identickými body („SJTSK_Lukas“), který byl označen jako domácí (Home ScanWorld), jehož souřadnicová soustava bude použita. Dále byl vložen ScanWorld s finální registrací mračen („registrace_fin“) a zbývající ScanWorldy s čísly 11, 12. Registrace proběhla již podle výše popsaného postupu v kapitole 4.3.3. Provedená registrace byla ukončena funkcí Create ScanWorld/Freeze Registration a následně byl vytvořen výsledný ModelSpace. Tento ModelSpace obsahuje všechna registrovaná mračna v S-JTSK a Bpv. V protokolu o registraci (Registration – Show Diagnostics) jsou uvedeny dosažené odchylky na jednotlivých identických bodech, průměrná odchylka celé registrace a transformační klíč. Tento protokol je zařazen mezi přílohy. Průměrná absolutní odchylka registrace dosáhla hodnoty 13 mm.
5.1 Sjednocení do jednoho mračna Vytvořený ModelSpace při registraci „SJTSK“ obsahuje jednotlivá mračna, skeny nejsou spojeny do jednoho celku. Pro další práci bylo potřeba všechna zpracovaná mračna bodů spojit do mračna jediného. Jelikož spojení mračen je nevratný proces byla vytvořena kopie výsledného ScanWorldu, která byla vložena do nově vytvořené databáze „helf_unify“. Pro spojení mračen byla použita funkce unify. V ModelSpace byla vybrána všechna mračna a spuštěna funkce Tools – Unify Clouds (Obr. 22).
47
Obr. 22: Funkce Unify Clouds
Při použití této funkce nastal během jejího průběhu problém. Vždy, když výpočet dosáhl hodnoty zhruba 50 %, proces byl náhle ukončen a zobrazilo se chybové hlášení (Obr. 23).
Obr. 23: Chybová hláška funkce unify
S obdobným problémem se setkali i další uživatelé programu, kteří uvedli, že problém nejspíše způsobuje antivirový program. Experimentálně bylo zjištěno, že tento problém opravdu způsobuje antivirový program. Problém byl odstraněn vypnutím rezidentní ochrany antivirového programu před spuštěním této funkce.
5.2 Výsledný model hradu Výsledný 3D model hradu je reprezentován jedním mračnem bodů, které zobrazuje celý hrad. Mračno obsahuje 235,5 milionů bodů a je v souřadnicovém systému
48
S-JTSK a výškovém systému Bpv. Na Obr. 24 je zobrazen náhled výsledného 3D modelu hradu v podobě mračna bodů.
Obr. 24: Výsledné mračno bodů celého hradu
49
6 Drátový model hradeb Výsledný 3D model hradu je podkladem pro vyhotovení dalších výstupů. Bude zhotoven hrubý drátový model a podrobnější síťový model hradu. Tvorba drátového modelu hradu byla rozdělena mezi studentky Bc. Alžbětu Prokopovou a Bc. Petru Dífkovou. Zpracování síťového modelu hradu je úkolem studentky Bc. Jany Poesové. Aby bylo možno s finálním modelem pracovat i v programu Geomagic, byl transformován do redukovaného souřadnicového systému S-JTSK. Tento program neumožňuje na rozdíl od programu Cyclone pracovat se souřadnicemi bodů, které jsou v řádech desetitisíců a více. Transformace modelu do redukovaného souřadnicového systému proběhla pomocí nové registrace obdobně jako v případě transformace do systému S-JTSK. Pouze byly upraveny souřadnice identických bodů, na základě kterých je transformace provedena a které jsou naimportovány do nového Scanworldu. Souřadnice identických bodů bylo potřeba upravit stejně jako v případě transformace do S-JTSK a navíc redukovat. Souřadnice Y byla redukována o hodnotu 737 000 a souřadnice X o hodnotu 988 000. V Tab. 14 jsou uvedeny čísla bodů a souřadnice použitých identických bodů pro tuto transformaci. Číslo bodu 9001 9003 9101 9102 9103 9106 9107 9108 9111
x [m]
y [m]
z [m]
-7990,478 -7990,084 -7985,791 -7992,972 -7976,591 -7983,646 -7989,635 -7967,281 -7952,549
-443,537 -438,428 -443,514 -429,258 -437,188 -467,783 -482,621 -482,760 -492,464
327,681 325,873 320,747 319,284 318,461 316,868 314,607 314,709 319,156
Číslo bodu 9112 9114 9116 9117 9118 9120 9122 9125 9127
x [m]
y [m]
z [m]
-7982,424 -7998,989 -8014,235 -8016,279 -8020,213 -8020,863 -8015,617 -8000,992 -7936,849
-490,538 -468,007 -436,192 -432,669 -418,384 -407,093 -404,628 -405,432 -450,977
310,867 310,685 310,159 309,872 309,388 309,055 308,568 306,069 309,748
Tab. 14: Redukované souřadnice identických bodů
Tato diplomová práce se dále zabývá tvorbou prostorového modelu vnějšího hradebního pásu. Bylo vyzkoušeno několik postupů. Detailně je dále popsán pouze použitý postup. 50
6.1 Tvorba CAD modelu Zhotovený drátový model hradu by měl být využit pro prezentaci hradu na internetových stránkách, kde jej mohou návštěvníci virtuálně prozkoumat. Pro účely této vizualizace je pro model dostačující nižší přesnost. Detaily menší než 10 cm nejsou zobrazeny. Výsledné mračno mělo pro vyhodnocování v programu Cyclone přílišnou hustotu bodů. V tomto programu není bohužel dostupná funkce, kterou by bylo možno hustotu bodů zředit. Pro tvorbu CAD modelu bylo jako podklad využito zregistrované mračno bodů, které vytvořila studentka Bc. Jana Poesová v programu Geomagic Studio. Tento model vznikl registrováním původních neočištěných mračen bodů zadáním transformačních klíčů (ukázka na Obr. 25) určených z jednotlivých registrací, které byly prováděny v programu Cyclone při tvorbě výsledného modelu hradu (kapitola 4.3.3). V programu Geomagic Studio byla z tohoto modelu vyříznuta pouze část mračna zobrazující vnější hradební pás. V oddělené části mračna byla hustota bodů naředěna na vzájemný rozestup 5 cm. Takto upravené mračno bodů bylo vyexportováno ve formátu VTX a naimportováno do programu Cyclone.
Obr. 25: Ukázka transformačního klíče pro mračno
V programu Cyclone proběhlo nejprve ruční vyhodnocení. Jelikož hrad není tvořen přesnými geometrickými tvary nelze využít nástroje pro automatické vyhodnocení, které je možno výhodně použít při modelování např. potrubí a rovin. Ruční vyhodnocování probíhalo po menších blocích. Jednotlivý blok obsahuje oblast se čtyřmi až osmi cimbuřími hradeb. Mezi body mračna byly vkládány linie reprezentující charakteristické hrany jednotlivých objektů hradeb (Obr.
26). Vyhodnocená liniová kresba byla
vyexportována v CAD formátu .DXF, který je možno otevřít v programu Microstation, kde byla kresba dále upravena.
51
Obr. 26: Liniová kresba – vyhodnocené hrany
V programu Microstation byly vyhodnocené linie spojeny, aby vytvořili uzavřený drátový model a uzavřené rovinné obrazce byly zaplochovány. Jednotlivé vyhodnocené bloky zde byly spojeny do jednoho celku. Modelována byla cimbuří a zdivo hradeb k rozhraní s terénem. Ruční vyhodnocení, začištění kresby a plochování bylo zdlouhavé a časově náročné. Výsledný model je uložený ve formátu .DGN mezi elektronickými přílohami této práce (složka č. 4: CAD model). Pro pouhé prohlížení je zde uložen také ve formátu .PDF. Na Obr. 27 je zobrazen náhled části tohoto modelu.
Obr. 27: CAD model
52
Závěr Cílem této diplomové práce bylo zaměření a vytvoření aktuální dokumentace hradu Helfenburk u Úštěka. Na hradě bylo nejprve vybudováno nové bodové pole, aby výsledná dokumentace mohla být v souřadnicovém systému S-JTSK a výškovém systému Bpv. Pro zaměření celého hradu byla zvolena metoda laserového skenování. Zaměření hradu bylo provedeno skenerem Trimble TX5. Pro pokrytí celého exteriéru hradu bylo pořízeno 91 skenů. V průběhu měření byly na hradu rozmísťovány kulové a šachovnicové terče. Dále bylo provedeno kontrolní měření, jehož úkolem bylo ověřit přesnost vytvářeného modelu. Pro usnadnění pozdějšího vyhledání zaměřených kontrolních bodů v modelu byla vytvořena jejich dokumentace. Dokumentace kontrolních bodů je jedním z výstupů této práce. Naměřená data byla dále zpracována ve specializovaných programech pro laserové skenování: Geomagic Studio a Cyclone. Procesem registrace bylo všech 91 pořízených mračen bodů transformováno do společného souřadnicového systému. V průběhu registrací byly prováděny kontroly, zda na sebe jednotlivá mračna bodů správně navazují. Kontroly byly realizovány vybranými tenkými řezy v podélném a příčném směru. Prvním výsledkem je 3D model hradu reprezentovaný jednotlivými mračny bodů, které byly při registraci spojeny. Přesnost modelu byla ověřena porovnáním délek z kontrolního měření. Touto kontrolou bylo prokázáno, že model není nikterak deformován. Hodnota směrodatné odchylky rozdílů vybraných délek je 15 mm. Dalším krokem byla transformace 3D modelu hradu do požadovaného souřadnicového systému S-JTSK a výškového systému Bpv. Aby v modelu byly zachovány skutečné vzdálenosti byla použita shodnostní 3D transformace. Transformace byla provedena na základě 18 identických bodů, které jsou voleny ve středech zaměřených šachovnicových terčů. Výsledný 3D model hradu je reprezentován jedním mračnem bodů a je v souřadnicovém systému S-JTSK a výškovém systému Bpv. Část mračna představující vnější pás hradeb byla z 3D modelu hradu vyříznuta. Z této části mračna byl vytvořen CAD model, který je dalším výstupem této práce. V programu Cyclone byly nejprve liniemi vyhodnoceny hrany části cimbuří a hradeb. 53
Vyhodnocená liniová kresba byla dále v programu Microstation uzavřena a zaplochována. Při zpracování dokumentace hradu jsem se podílela na budování bodového pole, zaměření hradu, zpracování naměřených dat (modelování identických koulí, registrace), doplňujícím kontrolním měření a tvorbě CAD modelu pásu vnějších hradeb. Zvolená metoda laserového skenování se pro podrobné zaměření celého hradu ukázala jako vhodná. Běžnými geodetickými metodami měření by nebylo možno vytvořit takto podrobnou dokumentaci jako je výsledný 3D model hradu v podobě mračna bodů. Nevýhodou tohoto postupu zpracování dokumentace byla velmi časově náročná registrace pořízených mračen bodů. Registrace mračen trvala několik měsíců oproti zaměření hradu, které trvalo pouze několik dní. Registrace mračen byla časově náročná z následujících důvodů. Při měření bylo dostupných pouze šest kulových terčů, které byly opakovaně rozmisťovány. Na pořízených skenech nebyl zachycen dostatečný počet identických koulí, aby následně bylo mračna bodů možné vzájemně provázat. Z tohoto důvodu bylo potřeba doplnit podmínky mezi mračny zadáním vzájemných překrytů. Tímto se proces registrace stal složitějším. Dalším problémem byla zdlouhavá domluva studentů spolupracujících na tvorbě modelu hradu. Konkrétní dělba a organizace práce byla podceněna. Tímto vznikaly dlouhé časové prodlevy mezi jednotlivými etapami práce. Tyto prodlevy vznikaly dlouhými diskuzemi o postupu práce. Dalším z aspektů průtahů jsou nezkušenosti studentů s takto rozsáhlým projektem. Hlavní výstupy práce jsou dostupné v elektronické podobě. Data jsou po domluvě s vedoucím práce uloženy na školním PC, na kterém probíhalo zpracování.
54
Použité zdroje [1] KOL. AUTORŮ, 1984. Hrady, zámky a tvrze v Čechách, na Moravě a ve Slezku: Sv. 3, Severní Čechy. Praha: Svoboda. ISBN 25-089-84. [2] Letecké snímky: Helfenburk u Úštěka [online]. [cit. 2015-03-10]. Dostupné z: http://www.flyfoto.cz/2012/06/hrad-helfenburk-u-usteka.html [3] KOTREJCH, Vladimír. 2014. Osobní poznámky o historii zaměření Hrádku. [4] VOSYKA, Lukáš. 2015. Vybudování vztažné sítě pro detailní zaměření hradu Helfenburk. Diplomová práce. Dostupné z: http://gama.fsv.cvut.cz/~cepek/proj/dp/2015/lukas-vosyka-dp-2015.pdf. [5] ŠTRONER, Martin a KOL. AUTORŮ. 3D skenovací systémy. 1. vyd. Praha: České vysoké učení technické v Praze, 2013. ISBN 978-80-01-05371-3. [6] TOUŠEK, Martin. 2013. Testování skenovacího systému Trimble TX5. Bakalářská práce. Dostupné z: http://gama.fsv.cvut.cz/~cepek/proj/bp/2013/martin-tousek-bp-2013.pdf. [7] Trimble TX5 scanner: Datasheet [online]. [cit. 2015-05-10]. Dostupné z: http://trl.trimble.com/docushare/dsweb/Get/Document-628869/022504122_Trimble_TX5_DS_1012_LR.pdf [8] Leica Cyclone software: Guide [online]. [cit. 2015-05-10]. Dostupné z: http://gmv.cast.uark.edu/scanning/software/leica-software/leicacyclone/cyclone-shortcut-guides/leica-cyclone-interface-basics-2-3/ [9] ŠTRONER, Martin a HAMPACHER, Miroslav. Zpracování a analýza měření v inženýrské geodézii. 1. vyd. Praha: České vysoké učení technické v Praze, 2011. ISBN 978-80-01-04900-6.
55
Použité zkratky Bpv
výškový systém baltský - po vyrovnání
CAD
computer aided drafting (počítačem podporované kreslení)
CCD kamera
elektronická součástka pro snímání obrazové informace
GB
gigabyte
(označení
jednotky
pro
množství
informace
v informatice) GNSS
Global Navigation Satellite Systém (globální družicový polohový systém)
ICP
Iterative Closest Point
RGB
barevný model (červená, zelená, modrá)
S-JTSK
souřadnicový systém Jednotné trigonometrické sítě katastrální
3D
trojdimenzionální
56
Seznam obrázků Obr. 1: Letecký pohled na Helfenburk [2] ............................................................ 11 Obr. 2: Vlevo: Bývalá hlavní vstupní brána s brankou pro pěší. Vpravo: Dolní brána, dnes sloužící jako vstupní. ..................................................................................... 12 Obr. 3: Výřez z výsledného tachymetrického plánu [3] ....................................... 14 Obr. 4: Ukázka použité stabilizace ........................................................................ 16 Obr. 5: Skener Trimble TX5................................................................................... 19 Obr. 6: Ukázka umístění použitých terčů ............................................................. 22 Obr. 7: Menu nastavení parametrů měření [6] .................................................... 24 Obr. 8: Ukázka pořízené fotodokumentace kontrolních bodů ............................ 26 Obr. 9: Nastavení funkce Uniform ........................................................................ 30 Obr. 10: Nastavení funkce Batch processing ........................................................ 30 Obr. 11: Zakládání nové databáze ........................................................................ 31 Obr. 12: Nastavení importu dat Cyclone .............................................................. 32 Obr. 13: Struktura Scanworldu ............................................................................. 33 Obr. 14: Nastavení průměru koule ....................................................................... 34 Obr. 15: Vlevo: koule viditelná v mračnu. Uprostřed: koule vybraná funkcí Region Grow. Vpravo: vymodelovaná koule ......................................................... 34 Obr. 16: Možnosti funkce Region Grow ............................................................... 35 Obr. 17: Pracovní okno registrace ........................................................................ 37 Obr. 18: Zadávání překrytu v programu Cyclone ................................................. 41 Obr. 19: Nedoléhající mračna v kontrolním řezu ................................................. 42 Obr. 20: Naskenovaný terč v mračnu ................................................................... 43 Obr. 21: Kontrolní bod na hradbách..................................................................... 44
57
Obr. 22: Funkce Unify Clouds ............................................................................... 48 Obr. 23: Chybová hláška funkce unify .................................................................. 48 Obr. 24: Výsledné mračno bodů celého hradu .................................................... 49 Obr. 25: Ukázka transformačního klíče pro mračno ............................................ 51 Obr. 26: Liniová kresba – vyhodnocené hrany ..................................................... 52 Obr. 27: CAD model .............................................................................................. 52
58
Seznam tabulek Tab. 1: Rozpis provedených prací ......................................................................... 15 Tab. 2: Rozpis provedeného skenování ................................................................ 20 Tab. 3: Přehled pořízených skenů ........................................................................ 23 Tab. 4: Použité parametry skenování ................................................................... 24 Tab. 5: Ukázka z přehledové tabulky .................................................................... 27 Tab. 6: Konfigurace školního PC ........................................................................... 28 Tab. 7: Přehled registrovaných mračen v prvním bloku ...................................... 39 Tab. 8: Přehled registrovaných mračen v druhém bloku ..................................... 40 Tab. 9: Přehled registrace třetího bloku............................................................... 40 Tab. 10: Délky v podélném směru ........................................................................ 45 Tab. 11: Délky v příčném směru ........................................................................... 45 Tab. 12: Délky ve směru osy Z .............................................................................. 45 Tab. 14: Souřadnice použitých identických bodů ................................................. 46 Tab. 15: Redukované souřadnice identických bodů............................................. 50
59
Seznam příloh 1) Přehledka s rozmístěním stanovisek a terčů 2) Přehledka s umístěním kontrolních bodů 3) Seznam souřadnic použitých bodů bodového pole a jejich použitá čísla 4) Seznam souřadnic kontrolních bodů 5) Přehledka kontrolní délky Elektronické přílohy na přiloženém CD: 1) Kontrolní měření: a) Protokol o výpočtu souřadnic kontrolních bodů b) Tabulka s popisem kontrolních bodů c) Fotodokumentace kontrolních bodů 2) Protokoly z provedených registrací mračen 3) Rozdíly kontrolních délek 4) CAD model pásu hradeb a) ve formátu .DGN b) Ve formátu .PDF 5) Text diplomové práce
Elektronické přílohy uložené na školním PC: 1) Naměřená data 2) Protokoly z provedených registrací mračen 3) Výsledné mračno bodů zobrazující celý hrad (v samostatné vrstvě jsou uloženy vyhodnocené kontrolní body) 4) Výsledný 3D model v S-JTSK, Bpv 5) Výsledný CAD model pásu hradeb ve formátu .DGN
60
Příloha č. 3: Seznam souřadnic a použitá čísla bodů bodového pole Číslo bodu (měření) 4001 4002 4003 4004 4005 4006 4007 4008 4009 4010
Y [m]
X [m]
Z [m]
737932,028 737964,757 737951,479 737946,324 737984,759 737994,010 737980,917 738001,072 737978,133 738017,712
988 469,142 988488,338 988462,533 988448,668 988428,006 988411,937 988441,580 988456,949 988481,945 988417,055
313,130 312,910 317,610 311,980 318,590 310,530 317,700 312,950 313,100 310,030
Číslo bodu (místopisy) 506 507 701 501 703 505 702 503 502 504
Příloha č. 4: Seznam souřadnic kontrolních bodů Číslo bodu 901 902 903 904 905 906 907 908 909 910 911 912 913 914 916 917 918 919 920 921 922 923
Y [m] 737947,915 737948,266 737947,192 737948,530 737942,011 737941,198 737943,074 737937,861 737937,425 737936,910 737937,006 737939,126 737938,029 737959,789 737948,548 737946,554 737936,372 737983,453 737981,906 737981,693 737979,325 737975,811
souřadnice X [m] 988466,529 988467,839 988465,357 988466,224 988455,686 988454,970 988457,408 988452,185 988452,110 988451,893 988451,905 988452,647 988452,253 988480,101 988470,777 988469,134 988473,172 988417,410 988418,601 988418,784 988420,661 988423,210
Z [m] 321,187 322,564 322,550 318,464 315,477 316,830 316,086 313,813 314,055 312,975 311,752 312,198 311,674 315,914 313,860 314,064 313,754 315,311 315,504 315,488 317,372 315,650
Číslo bodu 924 925 926 927 929 930 933 934 935 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967
Y [m] 737975,671 737973,600 737983,853 737984,959 737983,561 737980,985 737976,725 737986,590 737989,751 737992,863 737992,343 737991,547 737965,308 737977,136 737978,386 738006,391 738005,724 738005,460 738000,106 737999,792 738001,772 738002,517 738004,156 737994,356 737992,392 737987,268 737987,778 737973,318 737984,032 738018,617 738017,679 738018,795 738018,556 737942,299 737944,505 737944,386 737945,567 737948,402 737948,406 737955,340
souřadnice X [m] 988423,311 988424,808 988416,768 988416,247 988418,650 988420,605 988423,834 988416,294 988435,323 988429,308 988432,372 988432,958 988442,863 988437,755 988439,466 988418,768 988418,089 988417,910 988462,169 988462,673 988458,830 988457,392 988454,177 988442,110 988447,671 988462,939 988464,754 988487,635 988489,109 988416,627 988410,574 988412,424 988419,012 988453,506 988456,379 988456,331 988457,840 988463,716 988463,714 988459,433
Z [m] 317,278 317,347 317,193 316,279 313,224 312,800 313,521 313,443 318,462 318,544 325,748 324,014 327,982 319,963 319,506 321,826 322,497 323,829 314,632 315,399 314,270 314,547 314,553 331,110 320,469 320,556 316,383 316,847 316,709 311,229 311,679 312,975 312,880 315,479 315,532 319,370 319,379 318,704 319,549 318,478
Číslo bodu 968 970 971 972 9001 9003 9005 9013 9015 9017 9101 9102 9103 9104 9106 9107 9108 9109 9110 9111 9112 9114 9116 9117 9118 9120 9122 9125 9126 9127
Y [m] 737960,797 737956,497 737956,444 737955,331 737990,478 737990,084 737988,446 737993,230 737996,527 737995,043 737985,791 737992,972 737976,591 737964,159 737983,646 737989,635 737967,281 737963,757 737963,825 737952,549 737982,424 737998,989 738014,235 738016,279 738020,213 738020,863 738015,617 738000,992 737940,053 737936,849
souřadnice X [m] 988455,635 988449,450 988469,525 988470,561 988443,537 988438,428 988447,163 988433,147 988436,034 988431,690 988443,514 988429,258 988437,188 988441,948 988467,783 988482,621 988482,760 988477,855 988473,216 988492,464 988490,538 988468,007 988436,192 988432,669 988418,384 988407,093 988404,628 988405,432 988450,195 988450,977
Z [m] 318,503 317,301 328,991 329,212 327,681 325,873 327,002 324,551 326,667 332,152 320,747 319,284 318,461 333,723 316,868 314,607 314,709 328,582 316,033 319,156 310,867 310,685 310,159 309,872 309,388 309,055 308,568 306,069 311,276 309,748