Vysoké učení technické v Brně Brno University of Technology
Fakulta strojního inženýrství Ústav konstruování / Odbor metodiky konstruování Faculty of Mechanical Engineering Institute of Machine and Industrial Design / Department of Design Theory
Metodika vyhodnocování tvarových odchylek výkovků pomocí analýzy obrazu
[Pojednání ke státní doktorské zkoušce] [Discourse on the Dissertation]
Autor práce: Ing. Aneta Zatočilová Author
Brno 2012
Vysoké učení technické v Brně Brno University of Technology
Fakulta strojního inženýrství Ústav konstruování / Odbor metodiky konstruování Faculty of Mechanical Engineering Institute of Machine and Industrial Design / Department of Design Theory
Metodika vyhodnocování tvarových odchylek výkovků pomocí analýzy obrazu
[Pojednání ke státní doktorské zkoušce] [Discourse on the Dissertation]
Autor práce: Ing. Aneta Zatočilová Author
Vedoucí práce: doc. Ing. Jan Brandejs, CSc. Supervisor
Brno 2012
TATO PRÁCE JE SPOLUFINANCOVÁNA EVROPSKÝM SOCIÁLNÍM FONDEM A STÁTNÍM ROZPOČTEM ČESKÉ REPUBLIKY.
OBSAH
1 ÚVOD ....................................................................................................................... 5 2 VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE.......................................................................................... 9 3 SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ ................................................ 10 3.1 Metody měření a vyhodnocování tvaru objektů vysoké teploty....................... 10 3.1.1 Aktivní systémy založené na laserovém skenování povrchu .................... 11 3.1.2 Aktivní optické systémy s xenonovým světlem ........................................ 14 3.1.3 Pasivní systémy založené fotogrammetrii a analýze obrazu ..................... 17 3.2 Průmyslové aplikace blízké fotogrammetrie .................................................... 21 3.2.1 Rozdělení fotogrammetrických systémů ................................................... 21 3.2.2 Měření polohy více cílových bodů současně ............................................. 22 3.2.3 Sledování pozice a orientace dotykové sondy ........................................... 23 3.2.4 Měření ploch .............................................................................................. 24 3.3 Dílčí kroky fotogrammetrického měření .......................................................... 25 3.3.1 Rozpoznání klíčových primitiv ................................................................. 25 3.3.2 Lokalizace polohy elips a výpočet jejich parametrů ................................. 26 3.3.3Vnitřní kalibrace fotoaparátu ...................................................................... 26 3.3.4 Vnější kalibrace fotoaparátu ...................................................................... 28 4 ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE ................................... 30 5 VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ................................................................................................................. 32 6 SOUČASNÝ STAV ŘEŠENÍ............................................................................... 33 6.1Předběžné měření v reálných podmínkách provozu ...................................... 33 6.2 Metodika ....................................................................................................... 35 6.3 Návrh a realizace kalibrace........................................................................... 38 6.3 Návrh a realizace detekce hran ..................................................................... 42 7 ZÁVĚR .................................................................................................................. 44 8 SEZNAM LITERATURY ................................................................................... 45 9 SEZNAM OBRÁZKŮ A GRAFŮ ....................................................................... 51 10 SEZNAM TABULEK ........................................................................................ 53
strana
3
strana
4
ÚVOD
1
1 ÚVOD Základním prvkem fotogrammetrických systémů je obrazový snímač. Vývoj CCD1 čipu, společně se zdokonalováním počítačů, elektro-optických komponent, laseru a dalších optických zdrojů umožnil úspěšné zavedení optických metod měření a vyhodnocování v mnoha vědních oborech. V posledních patnácti letech došlo ke značnému pokroku ve vývoji technologií měření tvaru prostorových objektů a to zejména díky sjednocení poznatků z oblastí počítačového vidění a počítačové grafiky. Důsledkem vývoje CCD čipu bylo masivní rozšíření digitální fotografie v druhé polovině 90. let. Souběžně s vývojem digitální fotografie byly vyvíjeny metody jejího zpracování pomocí počítače. Měření tvaru prostorových objektů zpravidla vyžaduje implementaci matematických metod, které zajistí zpracování dat z jednotlivých snímků v jednotném souřadném systému [1].
Obr.1-1:Vztah mezi dosahovanou přesností měření a velikostí měřeného objektu u různých bezdotykových měřicích metod [2] 2.
V současné době existuje množství metod měření trojrozměrných objektů (obr. 1-1). Metody se liší kromě principu funkce, zejména dosahovanou přesností a pracovní vzdáleností, z čehož vyplývá také oblast jejich aplikace. Optické metody měření jsou často rozdělovány na metody pasivní, aktivní a kombinované (obr.1-2).
CCD (charge-coupled device) - snímač obrazové informace. Interferometry – interferometrie; industrial metrology – průmyslová metrologie; industrial photogrammetry – průmyslová fotogrammetrie; teodolite – teodolit; laser scanning – laserové skenování; architectural and engineering photogrammetry – fotogrammetrie aplikovaná na architektonické památky a na speciální technické oblasti v inženýrské geodesii); tachymetric surveying – tachymetrický průzkum (rychlé měření délek a vodorovných a svislých úhlů); DGPS (differential global positioning system) - diferenciální GPS; aerial photogrammetry – letecká fotogrammetrie; remote sensing - dálkový průzkum Země. 1 2
strana
5
ÚVOD
Obr. 1-2: Rozdělení optických metod 3D měření, podle [3, 4].
Pasivní metody nevyžadují interakci mezi měřeným objektem nebo scénou. Aktivní optické systémy pracují na principu interakce mezi objektem či jeho okolím a nějakým snímačem. Interakce se snímaným objektem může nastat[3]: • fyzickým kontaktem mezi snímačem a objektem (kontaktní 3D skenery), • emisí světla (strukturovaného nebo bílého) a následným snímáním světelné stopy vzniklé na objektu optickým systémem, • emisí elektromagnetické vlny[3].
Obr. 1-3: Dosahované rozlišení a pracovní vzdálenost tří metod vyhodnocování třetího rozměru 3[5]
Time-of-flight (TOF) - měření doby letu vysílaného impulzu; triangulation – triangulace; interferometry – interferometre. 3
strana
6
ÚVOD
Bezdotykové optické metody (tj. aktivní i pasivní) se dále dělí podle principu získávání třetího rozměru (tj. vzdálenosti měřeného bodu od obrazové roviny snímače). Třetí rozměr je vyhodnocován jednou z následujících metod: • triangulace, • měření doby letu vysílaného impulzu (TOF – „time-of-flight“) • a interferometrie [5]. Metody využívající el. magnetické záření (tj. pasivní i aktivní) můžeme rozdělit podle vlnové délky, která je snímačem vyhodnocována [1.1], jak ukazuje obr. 1-4.
4
Obr. 1-4: Rozdělení bezkontaktních optických metod měření dle vyhodnocované vlnové délky 4 [2].
Je tedy zřejmé, že fotogrammetrie patří mezi pasivní metody vyhodnocování elektromagnetického záření metodou triangulace. Nabídka trhu dokazuje, že mezi pasivními systémy využívanými v průmyslu jednoznačně převládají právě systémy fotogrammetrické. Americká fotogrammetrická společnost (ASPRS) definovala v roce 1980 pojem fotogrammetrie: „Fotogrammetrie je umění, věda a technika, která se zabývá získáváním informací o fyzických objektech a prostředí na základě procesu zaznamenávání, měření a interpretace fotografických snímků, vzorů elektromagnetického záření a dalších jevů“ [6]. Systémy se dělí na do dvou skupin podle vzdálenosti fotoaparátu od měřeného subjektu: • fotogrammetrie letecká (aerial) - vzdálenost > 300 m, • fotogrammetrie pozemní (close range nebo terrestrial) - vzdálenost < 300 m) [6]. Aktivní optické metody měření se začaly v průmyslu využívat koncem 20. století. Technologie využívající tyto metody (kontaktní i bezkontaktní) získaly označení 3D skenery. Podle doslovné citace je 3D skener „zařízení, které analyzuje reálné objekty nebo prostředí za účelem sběru dat o jejich tvaru, případně i vzhledu (např. barvě). Data mohou být využita k vytvoření digitálních trojrozměrných modelů“ [3].
Laser tracking – stanovení polohy cílových bodů metodu interferometrie paprsku koherentního světla.
4
strana
7
ÚVOD
3D skenování je založeno na vyhodnocování kontrastního světelného vzoru promítaného na povrch objektu [3]. Podle použitého typu světleného zdroje se dělí bezdotykové 3D skenery na dvě základní kategorie: • •
laserové 3D skenery (laser scanner), optické 3D skenery (optical 3D scanner, structured-light 3D scanner, fringe projection scanner)
V průmyslu našly uplatnění především bezkontaktní laserové a optické 3D skenery a kontaktní souřadné měřící stroje CMM 5 (viz kap. 3.2.3). V současné době umožňují 3D skenery získat kompletní informace o povrchu trojrozměrného objektu ve formě mračna bodů s přesností na desetiny až setiny milimetrů. Rozsah měřených objektů je od několika centimetrů až metrů dle typu zařízení. Omezení těchto metod jsou především dána pracovním prostředím a typem měřeného objektu. Přehled optických metod používaných pro měření prostorových objektů publikoval v roce 2000 Chen at al. [7]. Technologie trojrozměrné digitalizace nejsou zcela univerzální. Ačkoli je v současné době na trhu široký výběr bezdotykových optických systémů, stále existují aplikace, které vyžadují vývoj specifického systému. Takovou aplikací je také měření tvarových odchylek výkovků. Fotogrammetrické on-line systémy umožňují automatizované měření rozměrných objektů s odchylkou v řádech mm. Nicméně schopnost měřit pomocí fotogrammetrie také objekty vysoké teploty není příliš dobře popsána. Této oblasti je však v posledních letech věnována pozornost, což dokazuje potenciál metody.
CMM (Coordinate measuring machine) - kontaktní souřadný měřící stroj. Poloha bodů na součásti se měří formou souřadnicových rozměrů v osách X, Y, Z vůči definovanému souřadnému systému. 5
strana
8
VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE
2 VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE V úvodu pojednání byly představeny metody bezdotykového měření tvaru objektů. Kapitola Shrnutí současného stavu poznaní je zaměřena na metody měření tvaru výkovků o vysoké teplotě. Dále se pojednání zabývá metodami a aplikacemi průmyslové fotogrammetrie, neboť základní principy fotogrammetrie budou využity při návrhu systému. Výsledkem práce by měl být jednoduchý avšak robustní systém, který umožní vyhodnocování přímosti osy rozměrných výkovků on-line (tzn. za provozu výrobní linky). Motivace k vývoji systému vyplývá z reálných problémů, které řeší např. společnost Žďas a.s. Řešený problém: Vývoj systému pro měření přímosti osy válcových výkovků za tepla i za studena. Uvažované reálné podmínky pro měření: Při vývoji systému je nutné uvažovat reálné podmínky, ve kterých by měl pracovat. To jsou především omezené světelné podmínky, reálná velikost měřených dílů (válcové výkovky o délce max. 20 m a průměru max. 60 cm), nutnost on-line měření a „real-time“ vyhodnocování a další. Modelový systém pro experimentální testování a vývoj Za účelem vývoje a optimalizace systému bude sestavena jednoduchá měřicí stanice. Měřeným vzorkem bude válcová tyč o délce 60 cm a průměru max. 6 cm. Měření bude prováděno pomocí dvou digitálních fotoaparátů s vysokým rozlišením. Testování bude probíhat za přesně vymezených a kontrolovaných podmínek s vyloučením dalších intervenujících proměnných. Při provádění laboratorních experimentů však budou uvažovány také podmínky reálné, k nimž celý experiment směřuje. Přesto se jedná pouze o model reálné situace, proto jsou při vyvozování obecných závěrů uvažována možná modelová zkreslení. Návrh metodiky vývoje systému: Hlavní část práce bude zaměřena na návrh softwaru a jeho optimalizaci pro dané pracovní podmínky. Optimalizace softwaru bude probíhat v laboratorních podmínkách experimentálním měřením pomocí navržené aparatury. Testovaný vzorek bude nejdříve měřen za běžné pokojové teploty. Výsledky budou porovnány s měřením pomocí certifikovaného optického 3D skeneru Atos III Triple Scan 6. Po odladění softwaru pro měření za studena, bude modifikován pro měření vzorku vysoké teploty. Vyplývající předběžné cíle: • návrh měřicí stanice, • vývoj měřicího softwaru, • porovnání odchylky měření s měřením pomocí průmyslového 3D skeneru Atos, • testování a optimalizace systému při měření žhavých výkovků. Atos III Triple Scan – nejvýkonnější průmyslový optický 3D skener od firmy GOM. Dosahuje hustoty naměřených bodů až 0,01 mm. Viz. http://www.gom.com/metrology-systems/systemoverview/atos-triple-scan.html 6
strana
9
2
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3 SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ 3.1 Metody měření a vyhodnocování tvaru objektů vysoké teploty Rozměrné výkovky jsou výchozími polotovary, které slouží k výrobě komponent zvláště pro lodní, jaderný nebo petrochemický průmysl. Polotovar se kove při teplotě 850 až 1300 °C v otevřené formě pomocí hydraulických lisů s tlakem až 10 000 tun. Jedná se o cenově velmi nákladnou, často malosériovou až kusovou produkci. S rostoucími nároky na efektivitu výroby, rostou také nároky na tvarovou a rozměrovou přesnost výkovků. Z toho vyplývají také nároky na přesné, rychlé a bezpečné měření tvaru a rozměrů. Měření musí být součástí výrobního cyklu, jednak z důvodu eliminace opakovaného ohřevu, ale také z důvodu bezprostřední opravy zjištěných nedostatků. Dosud se běžně řeší měření rozměrných výkovků vysoké teploty pomocí jednoduchých přípravků ručně. Obsluha se tak musí pohybovat v blízkosti žhavého výkovku, což může být příčinou úrazů. Navíc jsou získány pouze neúplná a nepřesná data. Kromě ručních dotykových měřidel bylo také vyvinuto několik automatizovaných dotykových měřících systémů. Např. Nye at al. [8] v roce 2001 publikoval studii, která popisuje specifický dotykový měřicí systém, který byl zabudován do manipulátoru. Systém však nebyl flexibilní a neumožňoval získat komplexní rozměry výkovku. Zajímavý dotykový měřicí systém prezentoval již v roce 1986 Siemer [9]. Jeho systém byl také závislý na otáčivém pohybu výkovku při kování, vykazoval dobré výsledky v laboratorních podmínkách, nicméně byl špatně implementovatelný do průmyslové výroby. Dotykové měření výkovků neposkytuje dostatečnou flexibilitu, přesnost ani rychlost měření. Na druhou stranu, běžné bezdotykové optické měřicí systémy, jsou pro tento účel stěží aplikovatelné, zejména z důvodu vysoké teploty a velkých rozměrů výkovku. Z toho vyplývá nutnost vývoje specifického bezdotykového systému. Pokud je autorce známo, dosud bylo publikováno pouze málo studií zabývajících se vlivem teploty na přesnost, efektivitu a funkčnost optických měřicích systémů nebo optických snímačů. Yamauchi [10] publikoval článek zaměřený na studium chyb měření tvaru objektu vysoké teploty optickou metodou. Porovnává velikost chyb stanovených pomocí matematického modelu s chybami zjištěnými experimentálním měřením. Základní informace o vlivu teploty na laserové měřiče vzdálenosti a profilometry uvádí firma Omron [11], která se specializuje na výrobu snímačů pro širokou škálu aplikací. Zmiňuje tři základní problémy, které mohou způsobit chybu měření. Z nich vyplývají omezení, která jsou obdobná pro optické snímače vzdálenosti, laserové profilometry a optické či laserové 3D skenery: • Optický snímač, může vykazovat chybu způsobenou detekcí IR 7 paprsků, které emituje žhavý objekt. V případě skenerů s bílým světlem navíc dochází nedostatečnému kontrastu promítaného vzoru na povrchu žhavého objektu. Z toho důvodu je naprostá většina 3D skenerů určena k měření za stálých světlených i teplotních podmínek, které jsou omezeny minimem i maximem [11].
7
IR (infra red) – elektromagn. záření v infračervené části spektra o vlnové délce 760 nm až 1 mm.
strana
10
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Vysoká teplota při měření může způsobit chyby měření vlivem dilatace komponent. Tento problém je však snadné eliminovat požitím skeneru s vyšší pracovní vzdáleností [11]. • Dalším obtížně řešitelným problémem je chyba měření vlivem rozdílného lomu světla v oblasti kolem snímače a tepelně ovlivněném okolí žhavého objektu [11]. Ze zmíněných problémů a nároků na reálný provoz vyplývají kritéria, které by měl splňovat specifický bezdotykový měřící systém: • Dostatečná vzdálenost optického snímače (eliminace vlivu IR emisí výkovku). • Rychlé měření bez nutnosti zásahu do výrobního cyklu výkovku. • Online systém – automatické „real-time“ vyhodnocování (viz také kap. 3.2). Dosud bylo publikováno několik studií, které se zabývají vývojem a implementací specifických bezdotykových optických systémů určených pro měření průmyslových polotovarů vysoké teploty. Z těchto studií je patrné, že pro tyto aplikace byly vyvíjeny především tři typy bezdotykových měřicích systémů: • aktivní systémy založené na laserovém skenování povrchu, • aktivní optické systémy s xenonovým světlem, • pasivní systémy založené na fotogrammetrii, detekci entit a analýze obrazu. 3.1.1 Aktivní systémy založené na laserovém skenování povrchu Systémy založené vyhodnocování metodou TOF: • Zhisong Tian, Feng Ga, Zhenlin Jin, Xianchao Zhao 2009 - Dimension measurement of hot large forgings with a novel time-of-flight system[12]. Systém popsaný v tomto článku je založený na pulzním LADARU 8 . LADAR je laserový dálkoměr, který pracuje na principu TOF, což znamená, že je vyhodnocována doba, za kterou urazí laserový pulzní paprsek od emitoru k cílovému bodu a zpět. Popsaný systém se skládá z vyhodnocovací stanice TOF, skenovacího zařízení se dvěma stupni volnosti a sféricko-paralelního mechanismu se dvěma motory (SPM), který muže být také považován za dva teodolity. Systém je určen k vyhodnocování průměru a délky rozměrných výkovků při maximální teplotě 1600 °C ze vzdálenosti 20 m. Dva teodolity kontrolují naklápění skenovacího zařízení v ose X a Y. Tímto způsobem je získána 3D geometrie části výkovku (viz obr. 2-1).
Obr. 2-1:3D model části povrchu výkovku, získaný měřením pomocí laserového TOF systému [12].
LADAR (LaserDetection And Ranging) - laserový dálkoměr pracující na principu vyhodnocování doby letu impulzu. 8
strana
11
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
V článku je popsána konstrukce a funkce systému, ale také princip vyhodnocování dat. Systém byl testován nejdříve v laboratorních podmínkách a následně v reálných podmínkách výrobního závodu. Zjištěná přesnost systému se pohybuje v řádech milimetrů. Článek však neuvádí nedostatky systému, nebo oblasti potencionálního vylepšení. Článek navazuje na studii publikovanou v roce 1993: Profiling of hot surfaces by pulsed time-of-flight laser range finder techniques. Kari Määtta, Juha Kostamovaara, and Risto Myllylä[13]. • LaCam Forge od firmy Ferrotron a Menteq Technologies Inc.[14] Tento systém je rovněž určený přímo pro měření žhavých výkovků při otevřeném kování. Je založený na stejném principu jako systém, který publikoval Tian at al. využívá skenování pomocí pulzního polovodičového laseru a vyhodnocování metodou TOF. Před kováním je možné kontrolovat rozložení teploty na povrchu výkovku pomocí teplotní mapy. Při kování je možné měřit průběžně délku výkovku bez nutnosti přerušení výrobního procesu a softwarově optimalizovat proces kování. Po kování je možné zobrazit finální rozměry, přímost osy a válcovitost. Standardní proces kování s využitím měření pomocí LaCam systému je zobrazen na obr. 2-2. Obr. 2-3 zobrazuje výstupní protokol měření.
Obr. 2-2: Proces optimalizace tvaru výkovku pomocí laserového systému LaCam od firmy Ferrotron [14].
strana
12
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-3:Protokol z měření pomocí systémuLaCam. Z výsledků je patrné, že maximální průhyb překračuje povolenou odchylku 10 mm [14].
Systémy založené na vyhodnocování metodou triangulace: • TopScan od firmy Tecnogamma [15]. Systém určený přímo pro měření výkovků vysoké teploty je založen na vyhodnocování laserového paprsku metodou triangulace. Je určen k měření výkovků o teplotě 350 až 1150 °C a maximálních rozměrech 7 m. Přesnost dosahuje ± 5 mm, ve sférických souřadnicích ± 3 mm. Výstupem je obrys výkovku se základními rozměry.
Obr. 2-4:Výstupní data získaná měřením žhavého výkovku laserovým systémem TopScan, který nabízí firma Tecnogamma [15].
strana
13
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Laserový systém vyvinutý ve spolupráci firem IMEGO, Mefos a SINTEF[16] Tento systém byl vyvíjen od roku 2005 díky spolupráci tří švédských ocelárenských společností. Skládá se z laserového projektoru a spektrálně - selektivní videokamery. Modrý laser promítá na povrch výkovku kontrastní pruhy, jejichž poloha je vyhodnocována metodou triangulace. Při měření manipulátor otáčí s výkovkem a kamera vytvoří několik překrývajících se snímků s vysokým rozlišením. Snímky slouží k sestavení 3D modelu (viz obr. 2-5). Systém umožňuje měření průměru v libovolné vzdálenosti, měření délky, křivosti osy a odchylek od požadovaného tvaru.
Obr. 2-5:Výsledky měření pomocí systému od firemIMEGO, Mefos a SINTEF. Vlevo: výsledný 3D model výkovku s odchylkami od požadovaného tvaru, vpravo řez výkovkem [16].
3.1.2 Aktivní optické systémy s xenonovým světlem • Zhenyuan Jia, Bangguo Wang, Wei Liu, Yuwen Sun. 2009 - An improved image acquiring method for machine vision measurement of hot formed parts [17]. Tento článek poprvé popisuje využití spektrálně selektivní metody pro účely měření žhavých výkovků. Celý výkovek je osvětlen pomocí xenonové lampy. Systém se skládá ze dvou monochromatických CCD snímačů s filtrem s dolní propustí, xenonové lampy a obslužného počítače (viz obr. 2-6).
Obr. 2-6: Prvky optického systému s xenonovým světlem [17].
strana
14
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-7: Princip spektrálně-selektivní metody. Levá svislá osa udává rozsah citlivosti CCD snímače. Pravá osa udává sílu radiace dokonale černého tělesa. Vodorovná osa udává vlnovou délku. Šedá oblast znázorňuje spektrum, které bude snímačem detekováno [17].
Princip spočívá v selekci spektra, ve kterém nejsou CCD ovlivněny radiací výkovku a současně jsou citlivé na xenonové světlo (viz obr. 2-7). Na povrch výkovku je promítáno xenonové světlo, které se odráží k CCD snímačům, přičemž IR záření emitované výkovkem je odfiltrováno. Díky tomu mohou být na povrchu výkovku extrahovány body, které slouží ke složení výsledného 3D modelu ze stereosnímků (obr. 2-8).
Obr. 2-8: Výsledek extrakce bodů na povrchu výkovku [17].
Tento článek poprvé uvedl metodu extrakce bodů na povrchu žhavého výkovku. Přesnost měření rozměrů výkovku pomocí tohoto systému však nebyla popsána. Studie z následujících let uvádí, že extrakce tak velkého množství bodů vede k nízké rychlosti systému, z důvodu výpočetní náročnosti. Na tento výzkum navazuje Liu [18], který prezentoval vývoj „fringe projection“ 3D skeneru s xenonovým světlem založený na spektrálně selektivní metodě (viz další odstavec).
strana
15
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Wei Liu, Xinghua Jia, Zhengyuan Jia, Shuangjun Liu, Bangguo Wang, Jian Du. 2011 - Fast dimensional measurement method and experiment of the forgings under high temperature[18]. Tento článek vychází ze studie Jia at al. [17]. Na rozdíl od metod publikovaných v [17] využívá tento systém vyhodnocování promítaného světelného vzoru. Systém se skládá ze dvou CCD kamer s IR filtry, projektoru a obslužného počítače, jak ukazuje následující obrázek.
Obr 2-9: Vlevo: schéma optického stereo systému, vpravo jeho reálná podoba [18].
K detekci promítaných pruhů byl použit modifikovaný Gaussův model a spektrálně selektivní metoda popsaná již v roce 2009 [17]. Intenzita světelných pruhů v kolmém směru nemá na povrchu žhavého objektu typické Gaussovo rozložení. Oproti běžnému průběhu křivky, má křivka za podmínek vysoké teploty zploštělý tvar vrcholu, jak ukazuje obr. 2-10 vpravo. Z toho důvodu byl použit modifikovaný Gaussův model, který je blíže popsaný v článku.
Obr.k 2-10: Modifikovaný Gausův model aplikovaný na světelné pruhy promítané na žhavý výkovek. Vlevo: světelné pruhy na výkovku, modifikovaný gausův model v místě pruhu [18].
Pomocí první derivace intenzity byla získána funkce, která je zobrazena na obr. 2-11 vpravo. Následně byla druhá derivace této funkce položena rovno nule. Výsledkem jsou kořeny funkce („ridge points“ na obr. 2-11). Poloha středových pixelů promítaných pruhů („bottom point“ na obr. 2-11) byla na závěr vypočítána Newtonovou iterativní metodou se subpixelovou přesností pomocí dvou kořenů funkce.
strana
16
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-11: Aplikace první derivace na obraz výkovku. Vlevo: detail světelných pruhů po konvoluci, vpravo: výsledná funkce představující profil světelného pruhu [18].
Dále bylo nutné určit polohu hraničních bodů na výkovku. V článku jsou popsány dva algoritmy, které slouží k výpočtu polohy koncových bodů na povrchu výkovku. První algoritmus je založen na skutečnosti, že intenzita šedi je v pozadí výkovku značně odlišná, druhý ze skutečnosti, že na rozhraní výkovku a pozadí mají promítané pruhy významnou změnu křivosti. Aby mohly být extrahovány rozměry výkovku, musejí být pruhy promítány v horizontálním i vertikálním směru (viz obr. 2-12). Systém byl testován v laboratorních i reálných podmínkách. Přesnost rozměrů u čtyř měřených vzorků byla nižší než 1%.
Obr. 2-12: Hraniční body detekované na ocelové desce. Vlevo: hraniční body v horizontálním směru, vpravo: hraniční body ve vertikálním směru [18].
3.1.3 Pasivní systémy založené fotogrammetrii a analýze obrazu • S. Chaudhuria, A. Mukherjeeb, P. K. Duttab, A. Patrab. 2003 - Automatic measurement of frontal area and volume of fluid ingots using image-based instrumentation [19]. Automatizovaný fotogrammetrický systém byl v provozu aplikován např. k predikci místa řezu ingotu dle požadované hmotnosti dílčích částí. Ingoty mají 12 hran v průřezu a kónický tvar. Vytvořený software automaticky snímkuje ingot, předzpracuje snímky a pomocí jejich analýzy počítá objem, hmotnost a délku. Článek uvádí možné řešení výpočtu obsahu čelní plochy: • Výpočtem množství pixelů uvnitř hraniční oblasti.
strana
17
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Integrací 0 až 2π radiálního vektoru, který vychází z těžiště hraničních bodů. • Detekcí rohových bodů nebo vrcholů objektu a následnou derivací ze známé geometrie objektu. Poslední dvě metody jsou vhodné spíše pro objekty jednoduchého geometrického tvaru, z toho důvodu byla zvolena první metoda. Článek obsahuje popis použitého algoritmu, jehož hlavními částmi je prahování a afinitní transformace. Průměrná odchylka hmotnosti dosahuje ± 0.5%. Systém byl zaveden v průmyslovém závodě v Indii.
• S. B. Dworkin, T. J. Nye. 2004 - Image processing for machine vision measurement of hot formed parts [20]. Studie popisuje metodu měření základních rozměrů výkovku vysoké teploty pomocí adaptivního prahování a detekce hran. Přitom je testována přesnost měření u různých typů vstupních snímků: • monochromatické snímky ve viditelném spektru, • barevné snímky ve viditelném spektru, • monochromatické snímky s IR filtrem 9. Z článku vyplývá, že nejlepší vlastnosti pro prahování a detekci hran mají monochromatické snímky s IR filtrem. Článek uvádí, že běžné CCD snímače bez filtru nejsou schopné vytvořit vhodné snímky pro detekci hran prahováním a přesné měření rozměrů výkovku. Na následujícím obrázku je monochromatický snímek výkovku o teplotě 900 až 1300 ̊C, který byl pořízen s IR filtrem. Na obr. 2-14 jsou zobrazeny výsledky prahování.
Obr. 2-13:Monochromatický snímek žhavého výkovku při použití IR filtru [20].
Obr. 2-14: Výsledné prahování snímků. Vlevo monochromatický bez IR filtru, vpravo s IR filtrem [20].
9
IR filtr - propouští pouze IR paprsky o vlnové délce vyšší než specifikace daného filtru.
strana
18
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Ch. Hu, B. Liu, X. Song. 2008 - A novel edge detection approach used for online dimensional measurement of heavy forging[21] Článek představuje novou metodu detekce hran, která překonává problémy ostatních detektorů při aplikaci na výkovky vysoké teploty. Jednotlivé kroky zpracování obrazu byly navrženy tak, aby byly co nejlépe splněny podmínky Cannyho algoritmu (viz kap. 3.3.1). U snímků výkovků o vysoké teplotě byly navíc aplikovány kroky pro předzpracování a finálního zpracování obrazu. Barevný 24 bitový obraz byl převeden na 8mi bitový. Intenzita pixelů byla rozšířena na celou oblast 0 až 255 a světlost obrazu byla zvýšena. Následně byl obraz upraven adaptivním prahováním, které poprvé publikoval Daxing Z. at al. V roce 2008 [22]. Výsledky předzpracování obrazu jsou na obr. 2-15. Některé segmenty pixelů v blízkosti objektu byly nesprávně označeny jako hrany, nicméně hranice okují byly eliminovány, zatímco většina hraničních pixelů bylo zachováno, beze změny polohy. Proces předzpracování vstupního obrazu není standardní krok detekce hran, ale slouží ke zvýšení její účinnosti při použití na žhavé objekty.
Obr. 2-15:Snímky žhavého výkovku ve stupních šedi. Vlevo původní snímek, vpravo předzpracovaný snímek [21].
Na vstupní předzpracovaný obraz byla gradientní metoda detekce hran. Za účelem zlepšení výsledků detekce hran, byl nejdříve aplikován filtr s dolní propustí pro odstranění aditivního šumu. V první fázi tedy byla provedena konvoluce obrazu ve dvou směrech pomocí jádra s Gaussovým rozložením. Ve druhé fázi byl aplikován filtr s horní propustí – diferenciální operátor slouží k nalezení hran v obraze. V poslední fázi byla provedena validace hran, za účelem zachování skutečných hrana eliminace redundantních informací.
Obr. 2-16:Výsledky experimentu s různými metodami. Vlevo: detekce pomocí Cannyho hranového detektoru, vpravo: detekce pomocí prezentované metody bez předzpracování snímku [21].
strana
19
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Provedený experiment dokazuje, že předzpracování a finální úprava obrazu jsou nepostradatelnými kroky při detekci hran výkovků vysoké teploty. Výsledky experimentu a porovnání s výsledky detekce hran pomocí Cannyho detektoru jsou na obr. 2-16 a 2-17.
Obr. 2-17:Výsledky experimentu s různými metodami. Vlevo: detekce pomocí prezentované metody bez finálního zpracování snímku, vpravo: detekce pomocí prezentované metody [21].
strana
20
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.2 Průmyslové aplikace blízké fotogrammetrie 3.2.1 Rozdělení fotogrammetrických systémů Pozemní fotogrammetrické systémy používané v průmyslu se zpravidla dělí na dvě kategorie – off-line a on-line systémy [23].
Obr.2-18: Proces on-line a off-line fotogrammetrického měření v průmyslu[8].
Nejvýraznějším rozdílem je schopnost měření a vyhodnocení „real-time“ nebo téměř „real-time“, což je výsadou on-line systémů. Typický on-line systém se skládá ze dvou nebo více kalibrovaných fotoaparátů s danou polohou a orientací. Tyto systémy jsou určeny k automatizovanému měření v konkrétním prostoru. Bývají například součástí výrobní linky, kde slouží ke kontrole kvality. Přesnost on-line systémů bývá obvykle nižší než přesnost off-line systémů, neboť jsou omezeny počtem snímků [23, 24]. Off-line systémy jsou určeny pro univerzální měření mnoha různých průmyslových aplikací. Jejich univerzalita spočívá zejména v možnosti použití fotoaparátu a objektivu dle požadované přesnosti a velikosti objektu. Fáze měření a vyhodnocování je oddělená. Off-line systémy dosahují vyšší přesnosti, podle [23] okolo 0,05 mm. Na rozdíl do on-line systémů mohou být použity také analogové fotoaparáty a klasické fotografie. Již studie z roku 1988 [24] uvádí rozdělení průmyslových fotogrammetrických systémů podle použitých technologií.
strana
21
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-19: Rozdělení průmyslových fotogrammetrických systémů [24].
Nejběžnějšími aplikacemi fotogrammetrie jsou: • měření polohy většího množství cílových bodů, • sledování pozice a orientace dotykové sondy, • měření tvaru plochy [23, 24]. Firmy AICON [27] a Nikon Metrology [28] nabízí kompletní sortiment fotogrammetrických systémů, které pokrývají všechny tři skupiny průmyslových aplikací. 3.2.2 Měření polohy více cílových bodů současně Přesná lokalizace a digitalizace diskrétních bodů je patrně nejběžnější aplikací fotogrammetrie v průmyslu. Obvykle se jedná o systémy off-line měření využívající jeden fotoaparát s vysokým rozlišením, speciální značky a algoritmus „bundle adjustment“ (viz kap. 3.2.4). Na povrch měřeného objektu se aplikují dva typy značek – kódované a nekódované. Oba typy značek mívají obvykle kruhový tvar s bílým středem a černým okrajovým pásem. Kódované značky slouží přímo k digitalizaci diskrétních bodů na povrchu objektu. Nekódované značky mají navíc po obvodě umístěný unikátní znak, který umožňuje automatické rozpoznání každého bodu softwarem. Tyto značky slouží k nalezení geometrických vztahů mezi jednotlivými obrazovými rovinami. Díky tomu je možné určit polohu detekovaných nekódovaných bodů ze všech snímků vůči jednomu souřadnému systému. V okolí objektu musí být umístěna také jedna až dvě měřítkové tyče, které slouží k transformaci obrazu do správného měřítka. Přenos dat do počítače bývá řešen buď paměťovou kartou, nebo bezdrátově. Naměřená data se zpracovávají pomocí specializovaného softwaru ihned po měření nebo s libovolným odstupem času [30]. Přesnost měření je možné ovlivnit hustotou rozmístění značek a počtem pořízených snímků. Digitalizace každé značky deseti až dvanácti snímky zajistí přesnost RMS 1-sigma 10 až 0,025 mm [23]. RMS 1-sigma (root mean square nebo 1 standard deviation) – střední kvadratická chyba s úrovní spolehlivosti 68%. 10
strana
22
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr.2-20: Vlevo: fotogrammetrické měření pomocí tzv. kódovaných a nekódovaných značek, vpravo: příklad nekódovaných fotogrammetrických značek [29].
3.2.3 Sledování pozice a orientace dotykové sondy Měření pomocí dotykové sondy se provádí ručním nebo strojním polohováním. Systémy se strojním polohováním se nazývají 3D souřadnicové měřicí stroje, podle anglického „coordinate measurement machine“ (CMM). Tyto systémy se dále dělí na analogové a diskrétní. Diskrétní CMM měří jednotlivé body na povrchu objektu postupně přemisťováním sondy, přičemž je měřena poloha sondy mechanicky. Analogové CMM měří tvar objektu v celé ploše - sonda plynule kopíruje celý měřený povrch. Pohyb sondy je v tomto případě, stejně jako v případě ručního polohování sondou, snímán optickým systémem. Ten zaměřuje orientaci a polohu hrotu sondy a kontrolních bodů, které jsou umístěny na jejím těle. Fotogrammetrický systém může být složen pouze z jednoho optického snímače, ale i ze dvou nebo více snímačů. Obvykle platí, že čím více snímačů sleduje pohyb sondy, tím přesnějších výsledků měřicí systém dosahuje.
Obr. 2-21: Dotykové sondy: vlevo ruční dotyková sonda od firmy Metronor [26], vpravo nahoře ruční dotyková sonda od firmy AICON [27], vpravo dole CMM od firmy Reinshaw [25].
strana
23
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Mezi největší firmy nabízející CMM patří Reinshaw [25], Metronor (dnes spojen s Nikon Metrology) [26, 28] a AICON [27]. Systémy s ručním polohováním nabízí firma Metronor [26], AICON [27] a GOM [29]. 3.2.4 Měření ploch Fotogrammetrické metody měření plochy vyžadují použití minimálně dvou snímků, které jsou složeny do jednoho souřadného systému pomocí detekovaných bodů na povrchu objektu. Měřená plocha musí být opatřena vhodnou texturou, která umožní detekci a skládání vlícovacích bodů. Textura na měřeném objektu může být přirozená nebo uměle vytvořená například přilepením kontrastních pásek nebo projekcí. Hlavní výhodou fotogrammetrického měření ploch je možnost měření dynamické scény. Příkladem je měření, při němž se pohybují snímače, nebo dochází k pohybu či deformaci měřené plochy, jak ukazuje např. obr. 2-23. Častými aplikacemi jsou např. měření deformací karoserie při „crash testech“ nebo při měření funkce airbagů při nárazu. Dynamické změny tvaru plochy není možné měřit žádnou jinou optickou technologií včetně 3D skenování [23]. Fotogrammetrické systémy určené k měření dynamické scény nabízí např. firma GOM [29] nebo AICON [27].
Obr. 2-22: Vlevo: projekce vzoru a fotogrammetrické měření systémem PRO-SPOT [31], vpravo: fotogrammetrický systém určený k měření deformací od firmy GOM [29].
Obr. 2-23:Měření dynamické scény pomocí fotogrammetrického systému od firmy AICON[27].
strana
24
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.3 Dílčí kroky fotogrammetrického měření Klíčové informace v digitální podobě jsou výsledkem řešení komplexního algoritmu. Přesnost výsledku a rychlost výpočtu závisí na mnoha parametrech. Nedílnou součástí návrhu algoritmu je jeho testování a ladění. Neustále se zvyšující se nároky na přesnost optických systémů, jsou důvodem stále aktivního výzkumu těchto dílčích algoritmů. Cílem výzkumu bývá především zpřesnění výsledné rekonstrukce, ale také snížení výpočetní náročnosti, nebo přizpůsobení algoritmů pro specifické podmínky měření. Kompletní rekonstrukce 3D informací založená na fotogrammetrii obvykle obsahuje tyto dílčí procesy: • • • •
rozpoznání klíčových primitiv (značek), lokalizace polohy klíčových primitiv a výpočet jejich parametrů, kalibrace vnitřních parametrů fotoaparátu, kalibrace vnějších parametrů fotoaparátu - nalezení přesné orientace všech snímků vůči jednomu souřadnému systému.
3.3.1 Rozpoznání klíčových primitiv Metody zpracování digitálního obrazu byly vyvíjeny již v 60. letech 20. století, zejména pro účely satelitního a lékařského snímkování, ale také za účelem zlepšování vizuálních vlastností digitální fotografie. Do 70. lech bránila vývoji těchto metod špatná dostupnost počítačů i jejich výkonnost. Značná pozornost byla již od počátku věnována rozpoznávání objektů v digitálním obraze. První fází rozpoznávání objektů v obraze je detekce hran nebo segmentace obrazu. V oblasti průmyslové fotogrammetrie se využívá výhradně detekce hran [32]. Přehled metod detekce hran v digitálním obraze publikoval v roce 1975Davis [33], v roce 1982 Peli a Malah [34] a v roce 1986Zioua Tabbone [35]. Problematice detekce hran v digitálním obraze je stále věnována pozornost. Studie se v posledních letech zabývají především optimalizací algoritmů pro konkrétní úlohy detekce [36]. Mechanismy detekce je možné rozdělit podle technik, které používají: založené na první derivaci, založené na druhé derivaci, nepracující s derivací, nebo detektory pracující ve frekvenční oblasti. Velká skupina metod aproximuje derivaci pomocí konvoluce s vhodným jádrem. Mezi nejznámější operátory patří: Sobel, Prewitt, Roberts a Laplacian. V současné době je stále považováno za optimální řešení detekce hran Cannyho algoritmus [37], který nejlépe splňuje tři základní kritéria: • Minimální chybovost - poměr mezi množstvím detekovaných hran a množstvím šumu musí být co nejmenší. • Lokalizace – poloha detekované hrany musí být co nejblíže její skutečné poloze v obraze. • Jednoznačná odezva – operátor musí detekovat pouze jednu odezvu na každou hranu. Optimalizacím Cannyho algoritmu se v posledních letech věnoval např. Ding a Goshtasby [38] nebo R. Medina-Carnicer at al. [39]. Detekce hran pomocí Cannyho detektoru dosahuje pouze pixelové přesnosti. V posledních letech jsou vyvíjeny metody, jejichž cílem je sub-pixelová detekce, např. [40].
strana
25
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.3.2 Lokalizace polohy elips a výpočet jejich parametrů Fotogrammetrické systémy využívají detekci hran ke zjištění přesné polohy středů kruhových značek. Je zřejmé, že kruhové značky se v obrazové rovině perspektivní projekcí transformují do tvaru elipsy. Klíčovou oblastí fotogrammetrie je přesná lokalizace těchto geometrických primitiv na detekovaných hranách. Metody pro lokalizaci a výpočet parametrů elipsy je možné rozdělit do dvou hlavních skupin – metody shlukovací („voting/clustering“) a optimalizační. Studii zabývající srovnáním metod publikoval v roce 1994 Shortis at al. [41]. K metodám spadajícím do první skupiny patří algoritmy založené na Haughové transformaci [42-45], RANSAC [46, 47] a fuzzy shlukování [48, 49]. Tyto metody jsou odolné proti vlivu odlehlých pixelů, a umožňují detekci několika primitiv současně. Nevýhodou těchto metod je především výpočetní složitost a nižší rychlost. Druhá skupina metod je založena na optimalizaci cílové funkce. Počítá se optimální poloha elipsy vůči detekovaným hraničním bodům metodu nejmenších čtverců [50-57]. Hlavní výhodou těchto metod je jejich přesnost a rychlost. Nicméně citlivost k odlehlým pixelům je vyšší než v případě shlukovacích metod. V posledních letech je věnována pozornost výhradně této skupině metod. Jiandong at al. publikoval v roce 2009 studii popisující základní kroky algoritmu fotogrammetrické extrakce diskrétních bodů. Studie se hlouběji zabývá algoritmem přesné lokalizace fotogrammetrických značek [58]. 3.3.3Vnitřní kalibrace fotoaparátu Vnitřní kalibrace slouží k nalezení vztahů mezi body 3D scény a jejich průměty do 2D souřadnic snímku. Vnitřní kalibraci fotoaparátu je možné rozdělit na lineární a nelineární. Lineární kalibrace se provádí za účelem určení parametrů zobrazovacího systému, což jsou [59]: • konstanta kamery – ohnisková vzdálenost („principal distance“), • střed zobrazení („principal point“), • formát snímku [59]. Pět lineárních parametrů kamery je vyjádřeno v horní trojúhelníkové matici 3x3: 𝛼𝛼𝑥𝑥 𝐴𝐴 = � 0 0
𝛾𝛾 𝛼𝛼𝑦𝑦 0
𝑢𝑢0 𝑣𝑣0 � 1
Parametry αx= f . mx a αy = f . my reprezentují ohniskovou vzdálenost f v pixelech, kde mx a my jsou konstanty měřítka mezi pixely a vzdáleností. γ představuje faktor úhlové vychýlení osy x a y („skew factor“). A u0, v0 představují vzdálenost středu zobrazení od středu snímku [59]. Nelineárními vnitřními parametry fotoaparátu jsou parametry zkreslení objektivu: • radiální zkreslení způsobené nedokonalou křivostí čoček (posun bodu na snímku o radiální vzdálenost), obsahuje složku zkreslení způsobeného excentricitou čoček („decentring distorsion“),
strana
26
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• tečné (tangenciální) zkreslení způsobené převážně asymetrií radiálního zkreslení u soustavy čoček (často se zanedbává), • a „thin prism“ zkreslení, které je způsobené naklopením čoček vůči obrazovému snímači [59]. Do dnes bylo vyvinuto a úspěšně implementováno mnoho metod, které teoreticky řeší všechny problémy vnitřní kalibrace. I když je lineární a nelineární kalibrace obvykle řešena současně, pozornost je v posledních dvaceti letech věnována výhradně studiu nelineárního zkreslení čoček, zejména za účelem zpřesnění a zjednodušení kalibrace. Důvodem je to, že najít přesný model nelineárního zkreslení čoček je mnohem komplikovanějším problémem, než hledání lineárních parametrů zobrazovací soustavy. Conrady v roce 1919 poprvé popsal zkreslení způsobené excentrickou polohou čoček. V roce 1966 byl představen model s názvem „plumb line method“, který dělí zkreslení na 3 základní typy (Brown-Conrady model [60, 61]). Tento model umožňoval určit konstantu fotoaparátu a střed zobrazení, nicméně eliminace zkreslení objektivu byla ošetřena pouze částečně. I přesto našel od té doby uplatnění v mnoha aplikacích. Do dnes bylo prezentováno několik modifikovaných metod, nicméně tyto modifikace se zaměřovaly zejména na matematickou interpretaci. Kalibrace čoček byla dosud převážně řešena jako součet zkreslení radiálního, tečného a „thin prism“ – takové metody jsou označovány jako konvenční. Studie [62] z roku 2008 prezentuje model, který řeší jednotlivé složky zkreslení a jejich fyzické příčiny odděleně a uvažuje jejich vzájemné působení. Tento model dosahuje stejného korekčního efektu jako konvenční modely, obsahuje však méně koeficientů a tak zjednodušuje kalibraci. Fotoaparát je považován za vnitřně kalibrovaný, je-li známá jeho ohnisková vzdálenost, střed obrazové roviny a zkreslení objektivu. Podle toho, zda kalibrace využívá kalibrační element, rozlišujeme: • metrická kalibrace – kalibrace založená na digitalizaci kalibračního pole. • nemetrická kalibrace nebo také samokalibrace – kalibrace bez použití kalibračního pole, je řešena přímo pomocí entit ve scéně. Metrická i nemetrická kalibrace je založena na srovnávání polohy daného typu prvků ve zkresleném obraze s jejich skutečnou polohou. Kontrastní prvky jsou ve scéně rozmístěny tak, že tvoří pole prvků. V případě metrické kalibrace je nejjednodušším kalibračním polem čtvercová mřížka (kalibrace pomocí bodů v průsečících), jeho obdobou jsou pravidelně rozmístěné kruhové značky. Používají se také kalibrační pole s trojúhelníkovými nebo čtvercovými značkami. Pole se snímkuje kamerou z různých úhlů ze vzdálenosti, která bude použita při měření. Během kalibrace nesmí být měněna ohnisková vzdálenost a zaostření. Zaostření musí být nastaveno na vzdálenost, ve které bude kamera následně pracovat [80]. Sun a Cooperstock prezentovali v roce 2006 studii zabývající se empirickým vyhodnocením faktorů, které mají vliv na přesnost kalibrace [63]. 3D skenery často využívají před samotným měřením kalibraci metrickou, která zajistí vyšší přesnost systému, než kalibrace nemetrická. Nicméně do průmyslových fotogrammetrických systémů byly úspěšně implementovány metody samoklaibrace pomocí značek rozmístěných na scéně [64 - 67].
strana
27
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.3.4 Vnější kalibrace fotoaparátu Kalibrace vnějších parametrů fotoaparátu se provádí se za účelem: • určení transformačních parametrů (tři rotace, tři posunutí, obecně až16 parametrů) - popisují polohu fotoaparátu v prostoru a vyjadřují vzájemný vztah světového souřadného systému a souřadného systému fotoaparátu, • určení relativní orientace fotoaparátů, • určení báze (vzdálenost středů zobrazení dvou fotoaparátů)[59]. Kalibrace vnějších parametrů je stejně jako kalibrace vnitřních parametrů založena na perspektivním promítání, které je zjednodušeno na model dírkového fotoaparátu. 2D pozice bodu v digitálním obraze je často vyjadřována vektorem [u v 1]T, zatímco pozice bodu ve 3D světových souřadnicích rozšířeným vektorem [xw yw zw 1]T. Díky tomu je možné zapsat průmět světových souřadnic do souřadnic digitálního obrazu následovně [59]: 𝑥𝑥𝑤𝑤 𝑢𝑢 𝑧𝑧𝑐𝑐 �𝑣𝑣 � = 𝐴𝐴[𝑅𝑅 𝑇𝑇] �𝑦𝑦𝑤𝑤 � 1 1 A je matice lineárních vnitřních parametrů fotoaparátu, a R, T jsou transformační matice vnějších parametrů kalibrace (rotační matice a translační vektor). Pro určení těchto transformačních parametrů se využívají dva základní modely – lineální a nelineární model [59].
Lineární model Vzájemný vztah mezi fotografiemi je určen epipolární geometrií. Nejdříve se analyzuje vzájemná poloha konkrétních bodů v prostoru a jejich průmět do obrazové roviny. Tento vztah popisuje matice fotoaparátu P. Ta závisí jak na vnitřních, tak i na vnějších parametrech. Je více způsobů, jak získat průmět P a P´,v závislosti na tom, jaká data jsou k dispozici. Nejpoužívanější je však metoda využívající tzv. esenciální matici E nebo metoda přímé lineární transformace („direct linear transformation“DLT), kterou publikoval již v roce 1971 Abdel-Aziz & Karara [68]. Esenciální matice je speciální druh fundamentální matice, který vyžaduje znalost interních kalibračních dat. Esenciální i fundamentální matice byly poprvé popsány Longuet-Higginsem [69]. Jednoznačné řešení pro určení této matice zajišťuje průmět osmi bodů – podle toho název normalizovaný osmibodový algoritmus. Existuje však mnoho dalších algoritmů, z nichž většina v některé fázi využívá osmi-bodový algoritmus [70, 71, 72]. V posledních letech byly prezentovány také metody, které řeší lineární kalibraci pouze pomocí sedmi nebo šesti bodů. Pokud je určena esenciální matice je možné vypočítat transformační parametry matice R a T metodou singulárního rozkladu. Studii, která se zabývá komplexním problémem určení orientace kamer pomocí algoritmů fundamentální a esenciální matice prezentoval v roce 1996 Luongaa Faugeras [73]. Výhodou lineárního modelu je možnost kalibrace s proměnnou ohniskovou vzdáleností. Kalibrace je poměrně rychlá a jednoduchá, nicméně vyžaduje znalost více bodů oproti kalibraci nelineární.
strana
28
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Nelineární model V oblasti blízké fotogrammetrie je často využívána nelineární kalibrace pomocí algoritmu, který je známý jako „Bundle Adjustment“. Název vychází z představy paprsků světla („bundle“), které se odrážejí od 3D struktury a konvergují v obrazovém středu každého fotoaparátu. Klasické „bundle adjustment“ algoritmy jsou iteračními metodami, které se řeší jako nelineální problém nejmenších čtverců. K jednoznačnému řešení tohoto modelu stačí znalost minimálně pěti bodů [74, 75]. Cílem celého algoritmu je nalézt globálně konzistentní skupinu parametrů, které minimalizují chybu registrace mezi všemi páry snímků. Minimalizace se často řeší metodou Lavenberg-Marquardt [76] nebo Sampsonovou aproximací [74].
Obr.2-24: Orientace snímků pomocí algoritmu „Bundle adjustment“ [79].
Algoritmus „Bundle adjustment“ byl představen v 70. letech 20. století – Brown, 1974 [77], Granshaw, 1980 [78]. Triggs publikoval v roce 2000 studii, která se zabývá přehledem algoritmů „bundle adjustment“ a popisuje také nové přístupy založené na distribučních modelech [79]. Výhodou tohoto algoritmu je nutnost znalosti pouze pěti bodů. „Bundle adjustment“ bývá navíc často kombinován se samokalibrací interních parametrů, a je tedy vhodný pro rychlé fotogrammetrické měření v průmyslu. Srovnání metod kalibrace prezentoval v roce 2006 Remondino a Fraser [80].
strana
29
ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE
4 ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE V současné době na trhu množství sofistikovaných automatizovaných systémů bezdotykového měření tvaru a rozměrů objektů. Přesto je výzkum v této oblasti stále velmi aktivní. Z rešerše bylo zjištěno, že se výzkum zaměřuje především na zvyšování přesnosti, spolehlivosti a rychlosti již zavedených systémů a jejich úpravou pro specifické pracovní podmínky. Podstata výzkumu stojí na vývoji nových matematických modelů a komplexních algoritmů, jejich testování, optimalizaci a implementaci. Na základě rešerše bylo také zjištěno, že v poslední dekádě je věnována pozornost vývoji specifických bezdotykových systémů určených pro účely měření objektů o vysoké teplotě. Zejména v těžkém ocelárenském průmyslu jsou stále rezervy v efektivitě výroby a to především z důvodu absence pokročilých technologií měření tvaru a rozměrů za tepla. Důvody nutnosti a potřebnosti řešení dané problematiky v daném čase jsou zmíněny v následujících bodech: • S rostoucími nároky na efektivitu výroby, rostou také nároky na tvarovou a rozměrovou přesnost výkovků. Z toho vyplývají také nároky na přesné, rychlé a bezpečné měření tvaru a rozměrů. Výzkumný záměr koresponduje se současnými problémy efektivity výroby v ocelárenském průmyslu (např. společnost Žďas a.s., která se zabývá výrobou rozměrných výkovků). • Dosud se běžně řeší měření rozměrných výkovků vysoké teploty pomocí etalonů či velkých posuvných měřidel ručně. Obsluha se tak musí pohybovat v blízkosti žhavého výkovku, což může být příčinou úrazů. Navíc jsou získány pouze neúplná a nepřesná data. Dotykové měření výkovků neposkytuje dostatečnou flexibilitu, přesnost ani rychlost měření. • Běžné bezdotykové optické měřicí systémy, jsou pro účel měření výkovků stěží aplikovatelné, zejména z důvodu vysoké teploty a velkých rozměrů výkovků. Z toho vyplývá nutnost vývoje specifického bezdotykového systému. • Možnosti měření žhavých objektů pomocí fotogrammetrie a analýzy obrazu jsou dosud málo prozkoumány a komerční systémy, pokud je autorce známo, zatím neexistují, na rozdíl od systémů laserových. • Výhodami fotogrammetrického systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu výkovku a v neposlední řadě také cena. Nároky na specifické systémy měření objektů vysoké teploty jsou především: • dostatečná vzdálenost optického snímače (eliminace vlivu IR emisí výkovku), • rychlé měření bez nutnosti zásahu do výrobního cyklu výkovku, • online systém – automatické „real-time“ vyhodnocování (viz také kap. 4.2). Specifické systémy pro zmíněné aplikace jsou obvykle založeny: • na laserovém skenování, • na optickém skenování a promítání xenonového světla, nebo • na principu fotogrammetrie a analýzy obrazu.
strana
30
ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE
Laserové systémy dosahují vysoké přesnosti a komplexnosti měřených dat, jejich hlavní nevýhodou velmi vysoká cena. Možnosti měření žhavých objektů pomocí fotogrammetrie a analýzy obrazu jsou dosud málo prozkoumány a komerční systémy, pokud je autorce známo, zatím neexistují, na rozdíl od systémů aktivních. Pro řešení problému měření válcovitosti výkovku není třeba získat komplexní informace o jeho tvaru, což je jedním z důvodů vývoje systému založeném na principu fotogrammetrie a analýzy obrazu. Výhodami takového systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu.
strana
31
VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ
5 VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ Cílem doktorské disertační práce je navrhnout systém, který bude automaticky vyhodnocovat odchylky osy válcových výkovků od ideální přímé osy. Přitom musí být splněny nároky na dosahovanou přesnost – maximální odchylka ±5mm při měření výkovku o délce do 10 m. Cílem práce je také experimentálně stanovit rozdíl dosahované přesnosti při měření vzorků pokojové a vysoké teploty. Pro účely experimentálního měření a ladění algoritmu bude zkonstruován jednoduchý měřící stůl. Na základě globálních cílů byly stanoveny dílčí cíle řešení: 1) Navrhnout experimentální měřící aparaturu. 2) Navrhnout a realizovat software pro automatizované měření odchylek výkovků pokojové teploty. 3) Porovnat odchylky měření systému s měřením pomocí certifikovaného optickým 3D skenerem Atos Triple Scan. 4) Modifikovat software pro měření odchylek výkovků vysoké teploty. 5) Ověřit přesnost systému při měření vzorku pokojové a vysoké teploty. 6) Popsat hlavní vlivy na dosahovanou přesnost. Na základě rešerše bude provedena změna názvu dizertační práce. Nový název: Vývoj metodiky vyhodnocování tvarových odchylek výkovků pomocí fotogrammetrie a analýzy obrazu. V průběhu řešení se předpokládá spolupráce s dalšími institucemi. Protože téma disertační práce je z části složitý matematický problém, autorka by chtěla využít spolupráci ve formě konzultací s ústavem matematiky naší fakulty. Spolupráce při řešení algoritmů a samotného programu se předpokládá s Ing. Radkem Poliščukem, Ph.D.
strana
32
SOUČASNÝ STAV ŘEŠENÍ 6
6 SOUČASNÝ STAV ŘEŠENÍ 6.1 Předběžné měření v reálných podmínkách provozu Před návrhem metodiky automatizovaného měření bylo provedeno testovací měření v reálných podmínkách provozu. Ve společnosti Ždas a.s. byly vytvořeny snímky žhavých i chladných výkovků o délce 3000 m a průměru 60 cm. Fotoaparáty byly umístěny na připraveném stativu se dvěma rameny (viz obr. 6-1).
Obr. 6-1: Navržený stativ pro testovací měření v reálných podmínkách.
Ve snímcích byly manuálně detekovány hrany výkovků (software Adobe Photoshop). Následně byly hrany extrahovány a přeneseny na bílý podklad. Binární obraz byl následně vektorizován. Vzniklé 2D vektorové křivky byly načteny do 3D modeláře. Díky známé orientaci kamer byly křivky z obou snímků složeny do jednoho souřadného systému. Ze čtyř hran výkovku byla vypočítána 3D křivka jeho osy. Ta byla srovnána s ideální osou vedoucí koncovými body 3D křivky. Na závěr bylo manuálně nalezeno místo maximálního průhybu, určena jeho velikost, poloha a orientace.
Obr. 6-2: Měření výkovku, detekce hran v rovině A-A.
Obr. 6-3: Měření výkovku, detekce hran v rovině B-B.
strana
33
SOUČASNÝ STAV ŘEŠENÍ
Obr. 6-3: Graf měření osy v rovině A-A. Maximální odchylka D= – 9 mm, ve vzdálenosti z=1551 mm.
Obr. 6-4: Graf měření osy v rovině B-B. Maximální odchylka D= – 5,7 mm, ve vzdálenosti z=1726 mm.
Obrazová rovina A-A
Obrazová rovina B-B
Obr. 6-5: Čelní pohled na měřené osy výkovku v obou rovinách. Fialová křivka je vypočtená 3D osa s maximálním průhybem 10,5 mm ve vzdálenosti z = 1711 mm pod úhlem α = - 76°.
Při testovacím měření byly ověřeny podmínky reálného provozu. Problémem při snímkování byl zejména nestabilní a nerovný povrch, na kterém byl stativ postaven. Další problémy vznikaly při rychlém manuálním vyhodnocení – zejména nedostatečné informace o geometrii scény (nebyla provedena kalibrace), z toho vyplývající problémy při skládání snímků a vyrovnání perspektivy obrazů. Předběžné měření v reálných podmínkách potvrdilo, že fotogrammetrické měření a zpracování obrazu je možným řešením vyhodnocování tvarových odchylek výkovků za tepla i za studena. K dosažení přesných a správných výsledků je nutné dodržet několik základních předpokladů: 1) kontrastní pozadí výkovku ( v místě prizmat mohou metody detekce selhávat, proto je vhodné použít prizmata s malou šířkou), 2) zajistit dostatečně stabilní a tuhý stativ, pro uchycení kamer a nastavení jejich orientace, 3) provedení kalibrace systému přímo na místě.
strana
34
SOUČASNÝ STAV ŘEŠENÍ
Výsledky měření však nelze považovat za směrodatné, protože nebyly porovnány přesnějším měřením. Vývoj metodiky měření musí byt provedeno za přesně vymezených a kontrolovaných podmínek s vyloučením intervenujících proměnných. Proto bude vývoj měřicího systému probíhat v laboratorních podmínkách, kde bude systém testován a výsledky srovnávány s výsledky měření pomocí dalších bezdotykových metod. 6.2 Metodika Na základě rešerše byla navržena předběžná metodika řešení problému. Metodika se bude dále rozvíjet a modifikovat dle výsledků laboratorních experimentů. Hlavní část práce bude zaměřena na návrh algoritmu a jeho optimalizaci pro dané pracovní podmínky. Testovaný vzorek bude nejdříve měřen za běžné pokojové teploty. Po odladění algoritmu pro měření za studena, bude modifikován pro měření vzorku vysoké teploty.
Obr. 6-6: Diagram navržené metodiky řešení.
strana
35
SOUČASNÝ STAV ŘEŠENÍ
6.2 Návrh a realizace měřicí stanice K testování a optimalizaci metodiky v laboratorních podmínkách byla vytvořena měřicí stanice. Stanice se skládá ze tří základních částí – pracovního stolu, ramen pro uchycení fotoaparátů a prizmat pro vymezení polohy měřených vzorků. Měření je prováděno pomocí dvou fotoaparátů s vysokým rozlišením. Poloha odnímatelných fotoaparátů bude pevně dána konstrukcí ramen na robustním stole, v jehož ose budou umístěna prizmata s měřeným vzorkem. Pracovní stůl musí mít v celé ploše pod měřeným vzorkem kontrastní barvu, tzn. bílou. Vyrobená měřicí stanice je zobrazena na obr. 6-7.
y
x
z
Obr. 6-7: Vyrobená měřicí stanice s naznačeným objektovým souřadným systémem.
Konfigurace měřicího stolu a parametry měřených objektů: • • • • • •
úhel mezi obrazovými rovinami fotoaparátu cca 90°, poloha měřeného objektu ve středu obrazové roviny, zvolené fotoaparáty: Canon EOS 500D, použité ohnisko objektivu 18 mm, měřené vzorky – válcové tyče o délce do 700 mm a průměru 50 až 100 mm, požadovaná velikost zorného pole v rovině měřeného objektu 700 mm.
Parametry fotoaparátu Canon EOS 500D: • • • • •
rozlišení: 15,1 Mpix - 4752×3168 pixel2 velikost čipu (h×v): 22,3x14,9 mm2 velikost pixelu: 4,69×4,70 μm2 crop faktor: 1,6 Objektiv s proměnnou ohniskovou vzdáleností: 18 až 55 mm
strana
36
SOUČASNÝ STAV ŘEŠENÍ
Při návrhu měřicí stanice byly uvažovány nároky na její maximální velikost. Ta se řídí velikostí měřeného vzorku, vzdáleností fotoaparátů s daným zorným úhlem při zvolené ohniskové vzdálenosti. Z důvodu úspory prostoru byla volena nejkratší dostupná ohnisková vzdálenost daného objektivu f18 mm. Z velikosti měřeného vzorku byla určena požadovaná velikost zorného pole v horizontálním směru zxo_pož ≥ 720 mm a vypočtena minimální potřebná vzdálenost obrazové roviny Dmin: 𝐷𝐷𝑚𝑚𝑚𝑚𝑚𝑚 = kde:
𝑧𝑧𝑧𝑧𝑜𝑜_𝑝𝑝𝑝𝑝 ž × 𝐹𝐹 720 × 18 − 𝐹𝐹 = − 18 = 563 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 22,3
Dmin
(mm)
- minimální potřebná obrazové roviny od měřeného vzorku
f
(mm)
- zvolené ohnisko
sxo
(mm)
- velikost čipu v horizontálním směru
zXo_pož
(mm)
- požadovaná velikost zorného pole v horizontálním směru
Měřicí stanice byla vyrobena s rameny o délce 400x400 mm, vzdálenost obrazové roviny od měřeného objektu D = 570 mm. Zorné pole je při F18 v této vzdálenosti 486x728 mm2. Teoreticky dosažitelná přesnost pro skutečnou velikost zorného pole a rozlišení fotoaparátů bude porovnána se skutečně dosaženou přesností po kalibraci systému: 𝑃𝑃𝑣𝑣_𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 = 𝑃𝑃ℎ_𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 =
kde:
𝑝𝑝𝑝𝑝𝑜𝑜 × 𝑧𝑧𝑧𝑧𝑜𝑜 0,00470 × 486 = = 0,1533 𝑚𝑚𝑚𝑚 ≅ 0,15 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 14,9
𝑝𝑝𝑝𝑝𝑜𝑜 × 𝑧𝑧𝑧𝑧𝑜𝑜 0,00469 × 728 = = 0,1531 𝑚𝑚𝑚𝑚 ≅ 0,15 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 22,3
pyo_Dexp (mm)
- velikost pixelu ve vertikálním směru ve vzdálenosti Dexp
pxo_Dexp (mm)
- velikost pixelu v horizontálním směru ve vzdálenosti Dexp
pyo
(μm)
- velikost pixelu ve vertikálním směru
pXo
(μm)
- velikost pixelu v horizontálním směru
syo
(mm)
- velikost čipu ve vertikálním směru
sxo
(mm)
- velikost čipu v horizontálním směru
zyo
(mm)
-velikost zorného pole ve vertikálním směru
zxo
(mm)
- velikost zorného pole v horizontálním směru
Pro reálné podmínky musí být voleny průmyslové kamery, které jsou teplotně odolnější, robustnější konstrukce, monochromatické, bez zabudovaného IR filtru. Průmyslové kamery s vysokým rozlišením používají čipy Kodak KAI. V současné
strana
37
SOUČASNÝ STAV ŘEŠENÍ
době je čipem s nejvyšším rozlišením v této řadě čip KAI-29050 11 s 28,8 Mpix. Například průmyslové monochromatické kamery JAI AM-1600 GE 12, nebo Allied Vision Technology Pike F-1600 13 s čipem Kodak KAI-16000 mají rozlišení 16 Mpix (4872x3248 pixelů, velikost pixelu 7,4×7,4 µm). Mohou pracovat při teplotách -5 °C až 45 °C.
Obr. 6-8: Průmyslová monochromatická kamera AM-1600 GE s rozlišením 4872×3248 pixelů od společnosti JAI.
Obr. 6-9: Průmyslová monochromatická kamera Pike F-1600 s rozlišením 4872×3248 pixelů od společnosti Allied Vision Technology.
6.3 Návrh a realizace kalibrace Byla zvolena metrická kalibrace (viz kap. 3.3.3). Postup při kalibraci je možné rozdělit do dvou základních bodů: • snímkování, • zpracování. Snímkování Před samotným snímkováním je třeba zvolit vhodné kalibrační pole podle vzdálenosti snímkování a typu objektivu. Při měření výkovků v reálných podmínkách by musela být provedena kalibrace na vzdálenost 10 m. Je tedy nutné zvolit takové kalibrační pole, které je možné použít i pro velké rozměry zorného pole. Proto byl zvolen jako základní kalibrační prvek lanko s pravidelně rozloženými ryskami s roztečí 50 mm v horizontálním směru (osa x). Ve scéně byly také umístěny měrky ve zbývajících dvou směrech (v ose y a z). Orientace fotoaparátů byla zvolena tak, aby střed obrazové roviny ležel přibližně na značce ve středu scény. Scéna byla snímkována vždy jednou kamerou pouze z jednoho úhlu. Pole v ose x sloužilo k určení radiálního zkreslení obrazu (viz. kap. 3.3.3), měrky v ose y a z, sloužily k určení orientace obrazové roviny. CCD čipy Kodak KAI, viz http://www.truesenseimaging.com/markets/machine-vision monochromatická kamera JAI AM-1600 GE, viz http://eshop.integrys.com/imaging-and-video/category-4/model/16mp-en.html 13 Allied Vision Technology Pike F-1600, viz http://www.alliedvisiontec.com/us/products/cameras/firewire/pike.html 11 12
strana
38
SOUČASNÝ STAV ŘEŠENÍ
Obr. 6-10: Snímek pro kalibraci s kalibračními měrkami a obrazovým souřadným systémem. Zobrazená mřížka kopíruje deformaci obrazu způsobenou radiálním zkreslením.
Zpracování Nejdříve byla určena deformace obrazu radiální zkreslením. Ostatní komponenty nelineárního zkreslení (tangential, thin prism) nebyly řešeny z důvodu jejich zanedbatelného vlivu na přesnost. Střed obrazové roviny byl přiřazen středové značce na lanku. Dále byl určen počet pixelů mezi jednotlivými značkami v horizontálním směru a srovnán s reálnou roztečí značek. Odchylky způsobené radiálním zkreslením jsou patrné z tabulky 1. - počet pixelů mezi jednotlivými značkami Δpx klesá s rostoucí vzdáleností od optické osy. Na obr. 6-5 je zobrazena mřížka odpovídající radiálnímu zkreslení obrazu. Z tabulky je také patrné, že v blízkosti optické osy odpovídá velikost pixelu teoretické vypočtené hodnotě. Tab. 1: Vyhodnocení odchylek způsobených radiálním zkreslením: mm – skutečná poloha značek v (mm), Δpx počet pixelů mezi jednotlivými značkami, Δmm - skutečná rozteč značek v (mm), Δmm/Δpx – velikost pixelu mezi jednotlivými značkami v (mm).
Kalibrace kamera A Jednotky / značka Č.1 Č.2 Č.3 Č.4 Č.5 Č.6 Č.7 Č.8
Kalibrace kamera B
px
mm
Δpx
Δmm Δmm/Δpx px
-70 254 566 879 1183 1480 1773 2062
-10,9382 39,0618 89,0618 139,0618 189,0618 239,0618 289,0618 339,0618
324 312 313 304 297 293 289
50 50 50 50 50 50 50
0,1543 0,1603 0,1597 0,1645 0,1684 0,1706 0,1730
-4 322 641 959 1275 1581 1887 2190
mm
Δpx Δmm
Δmm/Δpx
-0,6250 49,3750 99,3750 149,3750 199,3750 249,3750 299,3750 349,3750
0 326 319 318 316 306 306 303
0 0,1534 0,1567 0,1572 0,1582 0,1634 0,1634 0,1650
0 50 50 50 50 50 50 50
strana
39
SOUČASNÝ STAV ŘEŠENÍ
Dalším krokem byla vnitřní kalibrace lineárních parametrů zobrazovacího systému, tj. výpočet polohy a orientace obrazové roviny. V obraze byla odečtena velikost některých prvků scény s kalibračním lankem - viz obr. 6-6 a tab. 2.
Obr. 6-11: Schéma geometrie systému v rovině zx (zOxO). Červeně je značen obrazový souřadný systém, modře objektový. Tab. 2: Velikost některých prvků scény odečtená ze snímků.
kamera A
kamera B
55,3 87,0 86,0
53,5 82,0 86,0
L (mm) M (mm) H (mm)
Průmět středu stolu v optických rovinách Průmět středové značky kalibračního lanka na stole Výška lanka nad stolem v klidu
Ze známých rozměrů ve scéně byl dopočítán úhel α natočení obrazové osy zo od objektové osy z a další pomocné úhly ve scéně: 𝐻𝐻 180 𝛼𝛼 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 𝑀𝑀 𝜋𝜋
86 180 ≐ 44,67° 𝛼𝛼𝐴𝐴 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 87 𝜋𝜋
86 180 ≐ 46,36° 𝛼𝛼𝐵𝐵 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 82 𝜋𝜋
𝛾𝛾 = 90 − 𝛼𝛼 − 𝛽𝛽
𝛾𝛾𝐴𝐴 = 90 − 44,67 − 39,79 = 5,54°
𝛾𝛾𝐵𝐵 = 90 − 39,79 − 38,27 = 5,37°
strana
40
(𝐿𝐿 × 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠) 180 𝛽𝛽 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � �× (𝐻𝐻 − 𝐿𝐿 ∗ 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐) 𝜋𝜋
(55,3 × 𝑠𝑠𝑠𝑠𝑠𝑠 44,67) 180 𝛽𝛽𝐴𝐴 = � �× ≐ 39,79° 86,0 − 55,3 × 𝑐𝑐𝑐𝑐𝑐𝑐 44,67 𝜋𝜋
(53,5 × 𝑠𝑠𝑠𝑠𝑠𝑠 46,36) 180 𝛽𝛽𝐵𝐵 = � �× ≐ 38,27° 86,0 − 53,5 × 𝑐𝑐𝑐𝑐𝑐𝑐 46,36 𝜋𝜋
SOUČASNÝ STAV ŘEŠENÍ Tab. 3: Výsledky výpočtů základních úhlů geometrie systému v rovině zx.
Jednotky/ úhel α β γ
Kamera A (rad) (deg) 0,78 44,67 0,69 39,79 0,10 5,54
Kamera B (rad) (deg) 0,81 46,36 0,67 38,27 0,09 5,37
popis sklon osy zo od osy z objektového souřadného systému sklon osy zo od osy y objektového souřadného systému úhel mezi optickou osou zo a průmětem středu stolu
Vzdálenost mezi průnikem optických os a obrazovými rovinami, byla vypočítána podle vzorce (viz obr. 6-11): 𝜋𝜋 − 𝛾𝛾 2 𝜋𝜋 𝐷𝐷𝐴𝐴 = 55,3 × 𝑡𝑡𝑡𝑡𝑡𝑡 − 5,54 ≐ 570,38 𝑚𝑚𝑚𝑚 2 𝐷𝐷 = 𝐿𝐿 × 𝑡𝑡𝑡𝑡𝑡𝑡
𝐷𝐷𝐵𝐵 = 53,5 × 𝑡𝑡𝑡𝑡𝑡𝑡
𝜋𝜋 − 5,37 ≐ 569,56 𝑚𝑚𝑚𝑚 2
Kontrolní měření Na závěr bylo provedeno měření k předběžnému stanovení dosahované přesnosti. Kalibrační lanko bylo zatíženo, přičemž průhyb změřený ve směru objektové osy y byl přibližně 7,5 mm. K nastavení průhybu lanka byla využita měrka s jemnou stupnicí (viz obr. 6-12). Nicméně umístění lanka na dané rysce bylo ověřeno pouhým okem. Z toho vyplývá, že z tohoto předběžného měření nelze vyvozovat obecné závěry o přesnosti, slouží tedy pouze orientačně.
Obr. 6-12: Snímek scény se zatíženým lankem.
Ze snímku se zatíženým lankem byla odečtena poloha jeho průmětu do optické osy yo (viz tab. 4). Pomocí následujících vztahů byla vypočtena poloha lanka v objektové ose y a odchylka ω vypočtené hodnoty od skutečné:
strana
41
SOUČASNÝ STAV ŘEŠENÍ
𝑌𝑌𝑜𝑜 cos 𝛼𝛼 5,3 = ≐ 7,5 𝑚𝑚𝑚𝑚 cos 44,67 5,1 = ≐ 7,4 𝑚𝑚𝑚𝑚 cos 46,36
𝜔𝜔 = 𝑌𝑌 − 𝑌𝑌𝑣𝑣
𝑌𝑌𝑣𝑣 = 𝑌𝑌𝑣𝑣𝑣𝑣
𝑌𝑌𝑣𝑣𝑣𝑣
𝜔𝜔𝐴𝐴 = −7,5 − (−7,45) = 0,05 ≐ 0 𝑚𝑚𝑚𝑚
𝜔𝜔𝐵𝐵 = −7,5 − (−7,39) = 0,11 ≐ 0,1 𝑚𝑚𝑚𝑚
Tab. 4: Výsledky předběžné kontroly přesnosti po kalibraci systému.
Jednotky/ označení Y Yo Yv ω
Kamera A (mm)
Kamera B (mm)
-7,5 -5,3 -7,45 0
-7,5 -5,1 -7,39 0,1
popis skutečný průhyb lanka průmět průhybu lanka do optické osy y vypočtený průhyb lanka v objektové ose y odchylka měření zaokrouhlená na desetiny
6.3 Návrh a realizace detekce hran Pro automatizovanou detekci hran výkovků byla s ohledem na charakter analyzovaného obrazu navržena gradientní metoda, validující jasovou hodnotu každého bodu analyzovaného řezu příčným mediánovým filtrem. Motivací pro toto řešení byla snaha o snížení chybovosti, dané vlivem impulzního šumu v obraze. Ten je zde představován zejména přítomností různých kontrolních značek a povrchových nedokonalostí na zkoumaném objektu (rýhy, okuje) a v jeho okolí. Nastavení vhodné šířky plovoucího filtračního okna a prahovacích parametrů bude podrobně popsáno v disertační práci. V současné době probíhá testování metody na snímcích zdeformované tyče. V případě selhání automatického vyhledávání je možné chybně identifikované řezy v uživatelském rozhraní programu buď vyřadit, nebo opravit manuálně.
Obr. 6-13: Snímek deformované válcové tyče po automatické detekci hran.
strana
42
SOUČASNÝ STAV ŘEŠENÍ
Výsledkem detekce hran je soubor hraničních bodů, z nichž je vypočítána 2D osa válcové tyče. Díky známé orientaci obou obrazových rovin bude proveden přepočet dvou 2D křivek na 3D křivku osy. Následně bude provedena manuální kontrola přesnosti výsledné 3D křivky porovnáním s výsledky měření pomocí systému Atos III triple scan. Výsledek rekonstrukce tvaru válcové tyče, která byla měřena pomocí zmíněného průmyslového skeneru je zobrazen na následujícím obrázku. Skutečná osa výkovku byla porovnána s ideální osou vedenou koncovými body skutečné osy.
Obr. 6-14:Výsledek měření a softwarové rekonstrukce deformované tyče pomocí systému Atos s určeným místem maximálního průhybu.
strana
43
ZÁVĚR
7 ZÁVĚR Pojednání ke státní doktorské zkoušce shrnuje současný stav poznání v oblasti bezdotykových optických metod měření se zaměřením na metody fotogrammetrie, jejich aplikace a základní kroky fotogrammetrického měření. Zaměřuje se blíže na možnosti měření tvaru rozměrných válcových výkovků za vysoké teploty. Z literatury vyplívá, že metody měření žhavých objektů jsou založeny na laserovém skenování, optickém skenování a promítání xenonového světla, nebo na principu fotogrammetrie a analýzy obrazu. Přitom v současné době se běžně provádí měření pomocí etalonů či velkých posuvných měřidel ručně. Problémem je vysoká cena a dostupnost specifických laserových systémů, které jsou v současné době jedinou komerčně vyráběnou technologií pro měření žhavých výkovků. Tyto systémy zajistí měření komplexního tvaru výkovku. Pro účely měření základních geometrických odchylek tvaru válcových výkovků (přímost osy, válcovitost) však není nutné znát komplexní informace o geometrii celého objektu. Měření založené na principech fotogrammetrie a analýzy obrazu je adekvátním řešením problému, což vyplývá také z rešerše. Výhodami fotogrammetrického systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu výkovku a v neposlední řadě také cena. Vývoj a testování metodiky probíhá v laboratorních podmínkách z důvodů zajištění vymezených a kontrolovaných podmínek měření s vyloučením intervenujících proměnných. V současné době je řešení práce ve fázi testování algoritmu detekce hran u chladných válcových tyčí. V další fázi proběhne vyhodnocení přesnosti navrženého systému srovnáním výsledků s měřením pomocí průmyslového skeneru Atos III Triple Scan. Navržená metodika se bude dále rozvíjet a modifikovat dle výsledků laboratorních experimentů. Předpokládá se měření tyčí různých rozměrů s různými velikostmi deformací, dále bude testována opakovatelnost měření. Po odladění softwaru pro měření výkovků za studena bude probíhat jeho modifikace pro měření za tepla. Výsledky měření žhavých výkovků budou porovnány s výsledky měření pomocí laserového systému.
strana
44
SEZNAM LITERATURY 8
8 SEZNAM LITERATURY ÚVOD [1] GIROD, B., GREINER, G.ANDNIEMANN, H. Principles of 3D Image Analysis and Synthesis.USA: Kluwer Academic Publishers, 2000. p. 466, ISBN 0-7923-7850-4. [2] LUHMANN, Thomas, et al. Close Range Photogrammetry: Principles, techniques and applications. Dunbeath, Caithness KW6 6EG, Scotland, UK: Whittles Publishing, 2011. Fundamental methods, p. 510. ISBN 978-184995057-2. [3] BERNARD, A. Reverse engineering for rapid product development: a state of the art. In K. G. HARDING.Three-Dimensional Imaging, Optical Metrology, and Inspection V. 1999, vol. 3835, p. 50-63. [4] REMONDINO, F.; EL-HAKIM, S. Image-Based 3D Modeling: A review. The photogrammetric Record. September 2006, vol. 115, no. 21, p. 269–291. [5] VOJÁČEK, Antonín. PMD senzor & 3D měření vzdálenosti - 1. část "Princip" [online]. 2006-07-03 [cit. 2011-06-23].
. [6] MATTHEWS, Neffra; NOBLE, Tom. Aerial and Close-Range Photogrammetric Technology: Providing Resource Documentation, Interpretation and Preservation. In [online]. Denver Federal Center: BLM, Technical Note 428, 2008 [cit. 2011-06-23]. < http://www.blm.gov/nstc/library/pdf/TN428.pdf>. [7] CHEN, F., BROWN, G.M. AND SONG, M. Overview of three-dimensional shape measurement using optical methods. Optical Engineering, 2000, vol. 39, no. 1, p. 10-22. METODY MĚŘENÍ ROZMĚRNÝCH VÝKOVKŮ [8] NYE, T.J., ELBADAN, A.M., BONE, G.M. Real-time process characterization of open die forging for adaptive control. ASME Journal of Engineering Materials andTechnology, 2001, vol. 123, no. 4, p. 511–516. [9] SIEMER, E., NEISCHWITZ, P., KOPP, R. Quality optimized process control in open die forging. Stahl und Eisen, 1986. Vol. 106, no. 8, p. 383–388. [10] YAMAUCHI, M. Errors in optical shape measurement caused by a hightemperature atmosphere.Opt. Eng., Sep 02, 2009, vol. 48, no. 9. [11]OMRON. OMRON.[online]. c2011 [cit. 2011-07-25]. [12] TIAN, Z. at al. Dimension measurement of hot large forgings with a novel timeof-flight system. The International Journal of Advanced Manufacturing Technology, 2009, vol. 44, no. 1, p. 125-132. [13] MÄTTA, K., KOSTAMOVAARA, J. AND MYLLYL, R. Profiling of hot surfaces by pulsed time-of-flight laser range finder techniques. Applied optics,
strana
45
SEZNAM LITERATURY
1993, vol. 32, no. 27, p. 5334-5347. [14] MITEQ. FERROTRON A MINTEQ DIVISION [online]. c2011 [cit. 2011-07-25]. LaCam Forge.. [15] MERMEC GROUP. MERMEC GROUP [online]. c2009 [cit. 2011-07-25]. Open Die Forgings.. [16] SINTEF. SINTEF [online]. November 19 2008 [cit. 2011-07-25]. Measuring 3D shape of steel shafts during forging. . [17] JIA, Z. at al. An improved image acquiring method for machine vision measurement of hot formed parts. Journal of Materials Processing Technology, 2010, vol. 210, no. 2, p. 267-271. [18] LIU, W. at al. Fast dimensional measurement method and experiment of the forgings under high temperature. Journal of Materials Processing Technology, 2011. vol. 211, Issue 2, p. 237–244. [19] CHAUDHURI, S. at al. Automatic measurement of frontal area and volume of fluted ingots using image-based instrumentation. Measurement, 2003, vol. 33, no. 4, p. 325-332. [20] DWORKIN, S.B. AND NYE, T.J. Image processing for machine vision measurement of hot formed parts. Journal of materials processing technology, 2006, vol. 174, no. 1-3, p. 1-6. [21] Ch. Hu, B. Liu, X. Song. A novel edge detection approach used for online dimensional measurement of heavy forging.2008 International Conference on Optical Instruments and Technology: Optoelectronic Measurement Technology and Applications, Proc. of SPIE,vol. 7160. FOTOGRAMMETRIE – APLIKACEV PRŮMYSLU [23] LUHMANN, T. Close range photogrammetry for industrial applications. ISPRS Journal of Photogrammetry and Remote Sensing, 2010, vol. 65, no. 6, p. 558569. [24] FRASER, C.S. State of the art in industrial photogrammetry. IAPRS, Kyoto, 1988, vol. 27, p. 166-181. [25] Reinshow. Reinshow [online]. 2001-2011 [cit. 2011-06-28]. Dostupné z WWW: . [26] Metronor. Metronor [online]. 2007 [cit. 2011-06-28]. Dostupné z WWW: < http://www.metronor.com/ >. [27] AICON. AICON [online]. 2011 [cit. 2011-06-28]. Dostupné z WWW:
strana
46
SEZNAM LITERATURY
[28] Nikon Metrology. Nikon Metrology NV [online]. 2011 [cit. 2011-06-30]. Dostupné z WWW: . [29] GOM. GOM: Optical Measurement Techniques [online]. 2011 [cit. 2011-06-28]. Dostupné z WWW: < http://www.gom.com/>. [30] JIANDONG, Z., LIYAN, Z. AND XIAOYU, D. Accurate 3D Target Positioning in Close Range Photogrammetry with Implicit Image Correction. Chinese Journal of Aeronautics, 2009, vol. 22, no. 6, p. 649-657. [31] Géodésie Maintenance Services. Géodésie Maintenance Services [online]. 1998 [cit. 2011-06-30]. Dostupné z WWW: . DETEKCE HRAN [32] ROSENFELD, A. Picture processing by computer. ACM Computing Surveys (CSUR). 1969, vol. 1,no. 3, p. 147-176. [33] DAVIS, L. S. A survey of edge detection techniques. Computer Graphics and Image Processing, 1975, vol. 4, no. 3, p. 248-270. [34] PELI, T. AND MALAH, D. A study of edge detection algorithms* 1. Computer Graphics and Image Processing, 1982, vol. 20, no. 1, p. 1-21. [35] ZIOU, D. AND TABBONE, S. Edge detection techniques-an overview. International Journal of Pattern Recognition and Image Analysis, 1998. [36] PAPARI, G. AND PETKOV, N. Edge and line oriented contour detection: State of the art. Image andVision Computing, 2011, vol. 29, no. 2-3, p. 79-103. [37] CANNY J. A computational approach to edge detection. IEEE Transactions on Pattern Analysis and Machine intelligence. 1986; 8(6): 679-698. [38] MEDINA-CARNICER, R. at al. A novel method to look for the hysteresis thresholds for the Canny edge detector. Pattern Recognition, Jun 2011, vol. 44, no. 6, p. 1201-1211. [39] DING, L. AND GOSHTASBY, A. On the Canny edge detector. Pattern Recognition, 2001, vol. 34, no. 3, p. 721-725. [40] DA, F.P. AND ZHANG, H. Sub-pixel edge detection based on an improved moment. Image and Vision Computing, Dec 2010, vol. 28, no. 12, p. 1645-1658. LOKALIZACE POLOHY ELIPS A VÝPOČET JEJÍCH PARAMETRŮ [41] SHORTIS, M.R., CLARKE, T.A. AND SHORT, T. A comparison of some techniques for the subpixel location of discrete target images. SPIE Videometrics, 1994, vol. 2350, no. III, p. 239-249. [42] YUEN, H.K., ILLINGWORTH, J. AND KITTLER, J. Detecting partially occluded ellipses using the Hough transform. Image and Vision Computing, 1989, vol. 7, no. 1, p. 31-37. [43] YIP, R.K.K., TAM, P.K.S. AND LEUNG, D.N.K. Modification of hough transform for circles and ellipses detection using a 2-dimensional array. Pattern Recognition, 1992, vol. 25, no. 9, p. 1007-1022.
strana
47
SEZNAM LITERATURY
[44] WU, W.-Y. AND WANG, M.-J.J. Elliptical object detection by using its geometric properties. Pattern Recognition, 1993, vol. 26, no. 10, p. 1499-1509. [45] NAIR, P.S. AND SAUNDERS, A.T. Hough transform based ellipse detection algorithm. Pattern Recognition Letters, 1996, vol. 17, no. 7, p. 777-784. [46] ROSIN, P.L. Ellipse fitting by accumulating five-point fits. Pattern Recognition Letters, 1993, vol. 14, no. 8, p. 661-669 [47] WERMAN, M. AND GEYZEL, Z. Fitting a second degree curve in the presence of error. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 1995, vol. 17, no. 2, p. 207-211. [48] DAVE, R.N. AND BHASWAN, K. Adaptive fuzzy shells clustering and detection of ellipses. Neural Networks, IEEE Transactions on, 1992, vol. 3, no. 5, p. 643-662. [49] GATH, I. AND HOORY, D. Fuzzy clustering of elliptic ring-shaped clusters. Pattern Recognition Letters, 1995, vol. 16, no. 7, p. 727-741. [50] GANDER, W., GOLUB, G.H. AND STREBEL, R. Least-squares fitting of circles and ellipses. BIT Numerical Mathematics,1994, vol. 34, no. 4, p. 558578. [51] FITZGIBBON, A.W., PILU, M. AND FISHER, R.B. Direct least squares fitting of ellipses. Pattern Recognition, Proceedings of the 13th International Conference on, 1996, vol. 1, p. 253-257. [52] HAL R, R. AND FLUSSER, J. Numerically stable direct least squares fitting of ellipses. The Sixth International Conference in Central Europe on Computer Graphics and Visualization, 1998, p. 59-108. [53] ANDREW, F. Direct Least Square Fitting of Ellipses. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1999, vol. 21, p. 476-480. [54] MAINI, E.S. Enhanced direct least square fitting of ellipses. International Journal of Pattern Recognition & Artificial Intelligence, 2006, vol. 20, no. 6, p. 939-953. [55] RAY, A. AND SRIVASTAVA, D.C. Non-linear least squares ellipse fitting using the genetic algorithm with applications to strain analysis. Journal of Structural Geology, 2008, vol. 30, no. 12, p. 1593-1602. [56] CHAUDHURI, D. A simple least squares method for fitting of ellipses and circles depends on border points of a two-tone image and their 3-D extensions. Pattern Recognition Letters, 2010, vol. 31, no. 9, p. 818-829. [57] KANATANI, K. AND RANGARAJAN, P. Hyper least squares fitting of circles and ellipses. Computational Statistics & Data Analysis, 2011, vol. 55, no. 6, p. 2197-2208. [58] JIANDONG, Z., LIYAN, Z. AND XIAOYU, D. Accurate 3D Target Positioning in Close Range Photogrammetry with Implicit Image Correction. Chinese Journal of Aeronautics, 2009, vol. 22, no. 6, p. 649-657.
strana
48
SEZNAM LITERATURY
VNITŘNÍ KALIBRACE [59] HARTLEY, R. AND ZISSERMAN, A. Multiple view geometry in Computer Vision. Edtion ed.: Cambridge university press, 2000. ISBN 0-521-54051-8. [60] D. C. BROWN, Decentering distortion of lenses, Photogrammetric Engineering and Remote Sensing, 1966, no. 24, p. 555–566. [61] D. C. BROWN, Close-range camera calibration, Photogrammetric Engineering and Remote Sensing, 1971, no. 42, p. 855–866. [62] WANG, J. at al. A new calibration model of camera lens distortion. Pattern Recognition, 2008, vol. 41, no. 2, p. 607-615. [63] SUN, W. AND COOPERSTOCK, J. R. An empirical evaluation of factors influencing camera calibration accuracy using three publicly available techniques. Machine Vision and Applications, 2006, vol. 17, no. 1, p. 51-67. [64] Z. ZHANG, On the epipolar geometry between two images with lens distortion, Proceedings of International Conference on Computer Vision and Pattern Recognition, 1996, vol. 1, no. I, p. 407-411. [65] STEIN, G. P. Lens distortion calibration using point correspondences. Proceedings of the International Conference on Computer Vision and Pattern Recognition, 1997, p. 602-608. [66] FITZGIBBON, A. W. Simultaneous linear estimation of multiple view geometry and lens distortion. IEEE International Conference on Computer Vision and Pattern Recognition, 2001, vol. 1, p. I-125. [67] RICOLFE-VIALA, C. AND SÁNCHEZ-SALMERÓN, A.-J. Robust metric calibration of non-linear camera lens distortion. Pattern Recognition, 2010, vol. 43, no. 4, p. 1688-1699. VNĚJŠÍ KALIBRACE [68] ABDEL-AZIZ, Y. I. and KARARA, H. M., 1971: Direct linear transformation into object space coordinates in close-range photogrammetry. Proc. Symposium on Close-Range Photogrammetry, p. 1-18 [69] LONGUET-HIGGINS, H. C. A Computer Algorithm for Reconstructing a Scene from Two Projections, Nature, Septemebr 1981, vol. 293, p. 133-135. [70] HARTLEY, R.I. In defense of the eight-point algorithm. Ieee Transactions on Pattern Analysis and Machine Intelligence, Jun 1997, vol. 19, no. 6, p. 580-593. [71] CHOJNACKI, W. at al. Revisiting Hartley's normalized eight-point algorithm. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2003, p. 1172-1177. [72] NISTÉR, D. An efficient solution to the five-point relative pose problem. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2004, vol. 26, no. 6, p. 756-770. [73] LUONG, Q.T. AND FAUGERAS, O.D. The fundamental matrix: Theory, algorithms, and stability analysis. International Journal of Computer Vision, 1996, vol. 17, no. 1, p. 43-75.
strana
49
SEZNAM LITERATURY
[74] STEWÉNIUS, H. at al. A minimal solution for relative pose with unknown focal length. Computer Vision and Pattern Recognition, 2005. CVPR 2005. IEEE Computer Society Conference on, 2005, vol. 2, p. 789-794. [75] STEWENIUS, H., ENGELS, C. AND NISTER, D. Recent developments on direct relative orientation. Isprs Journal of Photogrammetry and Remote Sensing, Jun 2006, vol. 60, no. 4, p. 284-294. [76] LOURAKIS, M.I.A. AND ARGYROS, A.A. Is Levenberg-Marquardt the most efficient optimization algorithm for implementing bundle adjustment? 2005. [77] D. C. BROWN. The bundle adjustment-progress and prospects. Int. Archives Photogrammetry, 1976. vol. 21, no. 3, p. 3–33 [78] S. GRANSHAW. Bundle adjustment methods in engineering photogrammetry. Photogrammetric Record, 1980, vol. 10, no. 56, p. 181–207. [79] TRIGGS, B. at al. Bundle adjustment--a modern synthesis. Vision algorithms: theory and practice, 2000, p. 153-177. [80] REMONDINO, F. AND FRASER, C. Digital camera calibration methods: considerations and comparisons. IPRS International Archives of Photogrammetry, Remote Sensing and Spatial Information Sciences, 2006, vol. 36, no. 5, p. 266-272.
strana
50
SEZNAM OBRÁZKŮ A GRAFŮ 9
9 SEZNAM OBRÁZKŮ A GRAFŮ Obr. 2-1
Vztah mezi dosahovanou přesností měření a velikostí měřeného objektu u různých bezdotykových měřicích metod [2]
Obr. 1-2
Rozdělení optických metod 3D měření, podle [3, 4].
Obr. 1-3
Dosahované rozlišení a pracovní vzdálenost tří metod vyhodnocování třetího rozměru [5].
Obr. 1-4
Rozdělení bezkontaktních optických metod měření dle vyhodnocované vlnové délky [2].
Obr. 2-1
3D model části povrchu výkovku, získaný měřením pomocí laserového TOF systému [12].
Obr. 2-2
Proces optimalizace tvaru výkovku pomocí laserového systému LaCam od firmy Ferrotron [14].
Obr. 2-3
Protokol z měření pomocí systému LaCam. Z výsledků je patrné, že maximální průhyb překračuje povolenou odchylku 10 mm[14].
Obr. 2-4
Výstupní data získaná měřením žhavého výkovku laserovým systémem TopScan, který nabízí firma Tecnogamma [15].
Obr. 2-5
Výsledky měření pomocí systému od firemIMEGO, Mefos a SINTEF[16].
Obr. 2-6
Prvky optického systému s xenonovým světlem [17].
Obr. 2-7
Princip spektrálně-selektivní metody[17].
Obr. 2-8
Výsledek extrakce bodů na povrchu výkovku [17].
Obr. 2-9
Schéma optického stereo systému s xenonovou lampou a jeho reálná podoba [18].
Obr. 2-10 Modifikovaný Gausův model aplikovaný na světelné pruhy promítané na žhavý výkovek[18]. Obr. 2-11 Aplikace první derivace na obraz výkovku[18]. Obr. 2-12 Hraniční body detekované na ocelové desce[18]. Obr. 2-13 Monochromatický snímek žhavého výkovku při použití IR filtru [20]. Obr. 2-14 Výsledné prahování snímků, vlevo monochromatický bez IR filtru, vpravo s IR filtrem [20]. Obr. 2-15 Snímky žhavého výkovku ve stupních šedi[21]. Obr. 2-16 Výsledky experimentu s různými metodami[21]. Obr. 2-17 Výsledky experimentu s různými metodami[21]. Obr. 2-18 Proces on-line a off-line fotogrammetrického měření v průmyslu [8]. Obr. 2-19 Rozdělení průmyslových fotogrammetrických systémů [24]. Obr. 2-20 Fotogrammetrické měření pomocí tzv. kódovaných a nekódovaných značek [29].
strana
51
SEZNAM OBRÁZKŮ A GRAFŮ
Obr. 2-21 Dotykové sondy: ruční dotyková sonda od firmy Metronor [26], ruční dotyková sonda od firmy AICON [27], CMM od firmy Reinshaw [25]. Obr. 2-22 Projekce vzoru a fotogrammetrické měření systémem PRO-SPOT [31], fotogrammetrický systém určený k měření deformací od firmy GOM [29]. Obr. 2-23 Měření dynamické scény pomocí fotogrammetrického systému od firmy AICON [27]. Obr. 2-24 Orientace snímků pomocí algoritmu „Bundle adjustment“ [79]. Obr. 6-1
Navržený stativ pro testovací měření v reálných podmínkách.
Obr. 6-2
Měření výkovku, detekce hran v rovině A-A.
Obr. 6-3
Měření výkovku, detekce hran v rovině B-B.
Obr. 6-3
Graf měření osy v rovině A-A.
Obr. 6-4
Graf měření osy v rovině B-B.
Obr. 6-5
Čelní pohled na měřené osy výkovku v obou rovinách.
Obr. 6-6
Diagram navržené metodiky řešení.
Obr. 6-7
Vyrobená měřicí stanice s naznačeným objektovým souřadným systémem.
Obr. 6-8
Průmyslová monochromatická kameraAM-1600 GE s rozlišením 4872×3248 pixelů od společnosti JAI.
Obr. 6-9
Průmyslová monochromatickákamera Pike F-1600 s rozlišením 4872×3248 pixelů od společnosti Allied Vision Technology.
Obr. 6-10 Snímek měřené scény pro kalibraci s kalibračními měrkami a obrazovým souřadným systémem. Obr. 6-11 Schéma geometrie systému v rovině zx (zoxo). Obr. 6-12 Snímek scény se zatíženým lankem. Obr. 6-13 Snímek deformované válcové tyče, určený k testování navržené metody detekce hran. Obr. 6-14 Výsledek měření a softwarové rekonstrukcedeformované tyče pomocí systému Atos s určeným místem maximálního průhybu.
strana
52
SEZNAM TABULEK
10
10 SEZNAM TABULEK Tab. 1 Vyhodnocení odchylek způsobených radiálním zkreslením
39
Tab. 2 Velikost některých prvků scény odečtená ze snímků
40
Tab. 3 Výsledky výpočtů základních úhlů geometrie systému v rovině zx
41
Tab. 4 Výsledky předběžné kontroly přesnosti po kalibraci systému
42
strana
53
OBSAH
1 ÚVOD ....................................................................................................................... 5 2 VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE.......................................................................................... 9 3 SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ ................................................ 10 3.1 Metody měření a vyhodnocování tvaru objektů vysoké teploty....................... 10 3.1.1 Aktivní systémy založené na laserovém skenování povrchu .................... 11 3.1.2 Aktivní optické systémy s xenonovým světlem ........................................ 14 3.1.3 Pasivní systémy založené fotogrammetrii a analýze obrazu ..................... 17 3.2 Průmyslové aplikace blízké fotogrammetrie .................................................... 21 3.2.1 Rozdělení fotogrammetrických systémů ................................................... 21 3.2.2 Měření polohy více cílových bodů současně ............................................. 22 3.2.3 Sledování pozice a orientace dotykové sondy ........................................... 23 3.2.4 Měření ploch .............................................................................................. 24 3.3 Dílčí kroky fotogrammetrického měření .......................................................... 25 3.3.1 Rozpoznání klíčových primitiv ................................................................. 25 3.3.2 Lokalizace polohy elips a výpočet jejich parametrů ................................. 26 3.3.3Vnitřní kalibrace fotoaparátu ...................................................................... 26 3.3.4 Vnější kalibrace fotoaparátu ...................................................................... 28 4 ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE ................................... 30 5 VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ................................................................................................................. 32 6 SOUČASNÝ STAV ŘEŠENÍ............................................................................... 33 6.1Předběžné měření v reálných podmínkách provozu ...................................... 33 6.2 Metodika ....................................................................................................... 35 6.3 Návrh a realizace kalibrace........................................................................... 38 6.3 Návrh a realizace detekce hran ..................................................................... 42 7 ZÁVĚR .................................................................................................................. 44 8 SEZNAM LITERATURY ................................................................................... 45 9 SEZNAM OBRÁZKŮ A GRAFŮ ....................................................................... 51 10 SEZNAM TABULEK ........................................................................................ 53
strana
3
strana
4
ÚVOD
1
1 ÚVOD Základním prvkem fotogrammetrických systémů je obrazový snímač. Vývoj CCD1 čipu, společně se zdokonalováním počítačů, elektro-optických komponent, laseru a dalších optických zdrojů umožnil úspěšné zavedení optických metod měření a vyhodnocování v mnoha vědních oborech. V posledních patnácti letech došlo ke značnému pokroku ve vývoji technologií měření tvaru prostorových objektů a to zejména díky sjednocení poznatků z oblastí počítačového vidění a počítačové grafiky. Důsledkem vývoje CCD čipu bylo masivní rozšíření digitální fotografie v druhé polovině 90. let. Souběžně s vývojem digitální fotografie byly vyvíjeny metody jejího zpracování pomocí počítače. Měření tvaru prostorových objektů zpravidla vyžaduje implementaci matematických metod, které zajistí zpracování dat z jednotlivých snímků v jednotném souřadném systému [1].
Obr.1-1:Vztah mezi dosahovanou přesností měření a velikostí měřeného objektu u různých bezdotykových měřicích metod [2] 2.
V současné době existuje množství metod měření trojrozměrných objektů (obr. 1-1). Metody se liší kromě principu funkce, zejména dosahovanou přesností a pracovní vzdáleností, z čehož vyplývá také oblast jejich aplikace. Optické metody měření jsou často rozdělovány na metody pasivní, aktivní a kombinované (obr.1-2).
CCD (charge-coupled device) - snímač obrazové informace. Interferometry – interferometrie; industrial metrology – průmyslová metrologie; industrial photogrammetry – průmyslová fotogrammetrie; teodolite – teodolit; laser scanning – laserové skenování; architectural and engineering photogrammetry – fotogrammetrie aplikovaná na architektonické památky a na speciální technické oblasti v inženýrské geodesii); tachymetric surveying – tachymetrický průzkum (rychlé měření délek a vodorovných a svislých úhlů); DGPS (differential global positioning system) - diferenciální GPS; aerial photogrammetry – letecká fotogrammetrie; remote sensing - dálkový průzkum Země. 1 2
strana
5
ÚVOD
Obr. 1-2: Rozdělení optických metod 3D měření, podle [3, 4].
Pasivní metody nevyžadují interakci mezi měřeným objektem nebo scénou. Aktivní optické systémy pracují na principu interakce mezi objektem či jeho okolím a nějakým snímačem. Interakce se snímaným objektem může nastat[3]: • fyzickým kontaktem mezi snímačem a objektem (kontaktní 3D skenery), • emisí světla (strukturovaného nebo bílého) a následným snímáním světelné stopy vzniklé na objektu optickým systémem, • emisí elektromagnetické vlny[3].
Obr. 1-3: Dosahované rozlišení a pracovní vzdálenost tří metod vyhodnocování třetího rozměru 3[5]
Time-of-flight (TOF) - měření doby letu vysílaného impulzu; triangulation – triangulace; interferometry – interferometre. 3
strana
6
ÚVOD
Bezdotykové optické metody (tj. aktivní i pasivní) se dále dělí podle principu získávání třetího rozměru (tj. vzdálenosti měřeného bodu od obrazové roviny snímače). Třetí rozměr je vyhodnocován jednou z následujících metod: • triangulace, • měření doby letu vysílaného impulzu (TOF – „time-of-flight“) • a interferometrie [5]. Metody využívající el. magnetické záření (tj. pasivní i aktivní) můžeme rozdělit podle vlnové délky, která je snímačem vyhodnocována [1.1], jak ukazuje obr. 1-4.
4
Obr. 1-4: Rozdělení bezkontaktních optických metod měření dle vyhodnocované vlnové délky 4 [2].
Je tedy zřejmé, že fotogrammetrie patří mezi pasivní metody vyhodnocování elektromagnetického záření metodou triangulace. Nabídka trhu dokazuje, že mezi pasivními systémy využívanými v průmyslu jednoznačně převládají právě systémy fotogrammetrické. Americká fotogrammetrická společnost (ASPRS) definovala v roce 1980 pojem fotogrammetrie: „Fotogrammetrie je umění, věda a technika, která se zabývá získáváním informací o fyzických objektech a prostředí na základě procesu zaznamenávání, měření a interpretace fotografických snímků, vzorů elektromagnetického záření a dalších jevů“ [6]. Systémy se dělí na do dvou skupin podle vzdálenosti fotoaparátu od měřeného subjektu: • fotogrammetrie letecká (aerial) - vzdálenost > 300 m, • fotogrammetrie pozemní (close range nebo terrestrial) - vzdálenost < 300 m) [6]. Aktivní optické metody měření se začaly v průmyslu využívat koncem 20. století. Technologie využívající tyto metody (kontaktní i bezkontaktní) získaly označení 3D skenery. Podle doslovné citace je 3D skener „zařízení, které analyzuje reálné objekty nebo prostředí za účelem sběru dat o jejich tvaru, případně i vzhledu (např. barvě). Data mohou být využita k vytvoření digitálních trojrozměrných modelů“ [3].
Laser tracking – stanovení polohy cílových bodů metodu interferometrie paprsku koherentního světla.
4
strana
7
ÚVOD
3D skenování je založeno na vyhodnocování kontrastního světelného vzoru promítaného na povrch objektu [3]. Podle použitého typu světleného zdroje se dělí bezdotykové 3D skenery na dvě základní kategorie: • •
laserové 3D skenery (laser scanner), optické 3D skenery (optical 3D scanner, structured-light 3D scanner, fringe projection scanner)
V průmyslu našly uplatnění především bezkontaktní laserové a optické 3D skenery a kontaktní souřadné měřící stroje CMM 5 (viz kap. 3.2.3). V současné době umožňují 3D skenery získat kompletní informace o povrchu trojrozměrného objektu ve formě mračna bodů s přesností na desetiny až setiny milimetrů. Rozsah měřených objektů je od několika centimetrů až metrů dle typu zařízení. Omezení těchto metod jsou především dána pracovním prostředím a typem měřeného objektu. Přehled optických metod používaných pro měření prostorových objektů publikoval v roce 2000 Chen at al. [7]. Technologie trojrozměrné digitalizace nejsou zcela univerzální. Ačkoli je v současné době na trhu široký výběr bezdotykových optických systémů, stále existují aplikace, které vyžadují vývoj specifického systému. Takovou aplikací je také měření tvarových odchylek výkovků. Fotogrammetrické on-line systémy umožňují automatizované měření rozměrných objektů s odchylkou v řádech mm. Nicméně schopnost měřit pomocí fotogrammetrie také objekty vysoké teploty není příliš dobře popsána. Této oblasti je však v posledních letech věnována pozornost, což dokazuje potenciál metody.
CMM (Coordinate measuring machine) - kontaktní souřadný měřící stroj. Poloha bodů na součásti se měří formou souřadnicových rozměrů v osách X, Y, Z vůči definovanému souřadnému systému. 5
strana
8
VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE
2 VYMEZENÍ ŘEŠENÉ PROBLEMATIKY A PŘEDBĚŽNÉ CÍLE DISERTAČNÍ PRÁCE V úvodu pojednání byly představeny metody bezdotykového měření tvaru objektů. Kapitola Shrnutí současného stavu poznaní je zaměřena na metody měření tvaru výkovků o vysoké teplotě. Dále se pojednání zabývá metodami a aplikacemi průmyslové fotogrammetrie, neboť základní principy fotogrammetrie budou využity při návrhu systému. Výsledkem práce by měl být jednoduchý avšak robustní systém, který umožní vyhodnocování přímosti osy rozměrných výkovků on-line (tzn. za provozu výrobní linky). Motivace k vývoji systému vyplývá z reálných problémů, které řeší např. společnost Žďas a.s. Řešený problém: Vývoj systému pro měření přímosti osy válcových výkovků za tepla i za studena. Uvažované reálné podmínky pro měření: Při vývoji systému je nutné uvažovat reálné podmínky, ve kterých by měl pracovat. To jsou především omezené světelné podmínky, reálná velikost měřených dílů (válcové výkovky o délce max. 20 m a průměru max. 60 cm), nutnost on-line měření a „real-time“ vyhodnocování a další. Modelový systém pro experimentální testování a vývoj Za účelem vývoje a optimalizace systému bude sestavena jednoduchá měřicí stanice. Měřeným vzorkem bude válcová tyč o délce 60 cm a průměru max. 6 cm. Měření bude prováděno pomocí dvou digitálních fotoaparátů s vysokým rozlišením. Testování bude probíhat za přesně vymezených a kontrolovaných podmínek s vyloučením dalších intervenujících proměnných. Při provádění laboratorních experimentů však budou uvažovány také podmínky reálné, k nimž celý experiment směřuje. Přesto se jedná pouze o model reálné situace, proto jsou při vyvozování obecných závěrů uvažována možná modelová zkreslení. Návrh metodiky vývoje systému: Hlavní část práce bude zaměřena na návrh softwaru a jeho optimalizaci pro dané pracovní podmínky. Optimalizace softwaru bude probíhat v laboratorních podmínkách experimentálním měřením pomocí navržené aparatury. Testovaný vzorek bude nejdříve měřen za běžné pokojové teploty. Výsledky budou porovnány s měřením pomocí certifikovaného optického 3D skeneru Atos III Triple Scan 6. Po odladění softwaru pro měření za studena, bude modifikován pro měření vzorku vysoké teploty. Vyplývající předběžné cíle: • návrh měřicí stanice, • vývoj měřicího softwaru, • porovnání odchylky měření s měřením pomocí průmyslového 3D skeneru Atos, • testování a optimalizace systému při měření žhavých výkovků. Atos III Triple Scan – nejvýkonnější průmyslový optický 3D skener od firmy GOM. Dosahuje hustoty naměřených bodů až 0,01 mm. Viz. http://www.gom.com/metrology-systems/systemoverview/atos-triple-scan.html 6
strana
9
2
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3 SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ 3.1 Metody měření a vyhodnocování tvaru objektů vysoké teploty Rozměrné výkovky jsou výchozími polotovary, které slouží k výrobě komponent zvláště pro lodní, jaderný nebo petrochemický průmysl. Polotovar se kove při teplotě 850 až 1300 °C v otevřené formě pomocí hydraulických lisů s tlakem až 10 000 tun. Jedná se o cenově velmi nákladnou, často malosériovou až kusovou produkci. S rostoucími nároky na efektivitu výroby, rostou také nároky na tvarovou a rozměrovou přesnost výkovků. Z toho vyplývají také nároky na přesné, rychlé a bezpečné měření tvaru a rozměrů. Měření musí být součástí výrobního cyklu, jednak z důvodu eliminace opakovaného ohřevu, ale také z důvodu bezprostřední opravy zjištěných nedostatků. Dosud se běžně řeší měření rozměrných výkovků vysoké teploty pomocí jednoduchých přípravků ručně. Obsluha se tak musí pohybovat v blízkosti žhavého výkovku, což může být příčinou úrazů. Navíc jsou získány pouze neúplná a nepřesná data. Kromě ručních dotykových měřidel bylo také vyvinuto několik automatizovaných dotykových měřících systémů. Např. Nye at al. [8] v roce 2001 publikoval studii, která popisuje specifický dotykový měřicí systém, který byl zabudován do manipulátoru. Systém však nebyl flexibilní a neumožňoval získat komplexní rozměry výkovku. Zajímavý dotykový měřicí systém prezentoval již v roce 1986 Siemer [9]. Jeho systém byl také závislý na otáčivém pohybu výkovku při kování, vykazoval dobré výsledky v laboratorních podmínkách, nicméně byl špatně implementovatelný do průmyslové výroby. Dotykové měření výkovků neposkytuje dostatečnou flexibilitu, přesnost ani rychlost měření. Na druhou stranu, běžné bezdotykové optické měřicí systémy, jsou pro tento účel stěží aplikovatelné, zejména z důvodu vysoké teploty a velkých rozměrů výkovku. Z toho vyplývá nutnost vývoje specifického bezdotykového systému. Pokud je autorce známo, dosud bylo publikováno pouze málo studií zabývajících se vlivem teploty na přesnost, efektivitu a funkčnost optických měřicích systémů nebo optických snímačů. Yamauchi [10] publikoval článek zaměřený na studium chyb měření tvaru objektu vysoké teploty optickou metodou. Porovnává velikost chyb stanovených pomocí matematického modelu s chybami zjištěnými experimentálním měřením. Základní informace o vlivu teploty na laserové měřiče vzdálenosti a profilometry uvádí firma Omron [11], která se specializuje na výrobu snímačů pro širokou škálu aplikací. Zmiňuje tři základní problémy, které mohou způsobit chybu měření. Z nich vyplývají omezení, která jsou obdobná pro optické snímače vzdálenosti, laserové profilometry a optické či laserové 3D skenery: • Optický snímač, může vykazovat chybu způsobenou detekcí IR 7 paprsků, které emituje žhavý objekt. V případě skenerů s bílým světlem navíc dochází nedostatečnému kontrastu promítaného vzoru na povrchu žhavého objektu. Z toho důvodu je naprostá většina 3D skenerů určena k měření za stálých světlených i teplotních podmínek, které jsou omezeny minimem i maximem [11].
7
IR (infra red) – elektromagn. záření v infračervené části spektra o vlnové délce 760 nm až 1 mm.
strana
10
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Vysoká teplota při měření může způsobit chyby měření vlivem dilatace komponent. Tento problém je však snadné eliminovat požitím skeneru s vyšší pracovní vzdáleností [11]. • Dalším obtížně řešitelným problémem je chyba měření vlivem rozdílného lomu světla v oblasti kolem snímače a tepelně ovlivněném okolí žhavého objektu [11]. Ze zmíněných problémů a nároků na reálný provoz vyplývají kritéria, které by měl splňovat specifický bezdotykový měřící systém: • Dostatečná vzdálenost optického snímače (eliminace vlivu IR emisí výkovku). • Rychlé měření bez nutnosti zásahu do výrobního cyklu výkovku. • Online systém – automatické „real-time“ vyhodnocování (viz také kap. 3.2). Dosud bylo publikováno několik studií, které se zabývají vývojem a implementací specifických bezdotykových optických systémů určených pro měření průmyslových polotovarů vysoké teploty. Z těchto studií je patrné, že pro tyto aplikace byly vyvíjeny především tři typy bezdotykových měřicích systémů: • aktivní systémy založené na laserovém skenování povrchu, • aktivní optické systémy s xenonovým světlem, • pasivní systémy založené na fotogrammetrii, detekci entit a analýze obrazu. 3.1.1 Aktivní systémy založené na laserovém skenování povrchu Systémy založené vyhodnocování metodou TOF: • Zhisong Tian, Feng Ga, Zhenlin Jin, Xianchao Zhao 2009 - Dimension measurement of hot large forgings with a novel time-of-flight system[12]. Systém popsaný v tomto článku je založený na pulzním LADARU 8 . LADAR je laserový dálkoměr, který pracuje na principu TOF, což znamená, že je vyhodnocována doba, za kterou urazí laserový pulzní paprsek od emitoru k cílovému bodu a zpět. Popsaný systém se skládá z vyhodnocovací stanice TOF, skenovacího zařízení se dvěma stupni volnosti a sféricko-paralelního mechanismu se dvěma motory (SPM), který muže být také považován za dva teodolity. Systém je určen k vyhodnocování průměru a délky rozměrných výkovků při maximální teplotě 1600 °C ze vzdálenosti 20 m. Dva teodolity kontrolují naklápění skenovacího zařízení v ose X a Y. Tímto způsobem je získána 3D geometrie části výkovku (viz obr. 2-1).
Obr. 2-1:3D model části povrchu výkovku, získaný měřením pomocí laserového TOF systému [12].
LADAR (LaserDetection And Ranging) - laserový dálkoměr pracující na principu vyhodnocování doby letu impulzu. 8
strana
11
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
V článku je popsána konstrukce a funkce systému, ale také princip vyhodnocování dat. Systém byl testován nejdříve v laboratorních podmínkách a následně v reálných podmínkách výrobního závodu. Zjištěná přesnost systému se pohybuje v řádech milimetrů. Článek však neuvádí nedostatky systému, nebo oblasti potencionálního vylepšení. Článek navazuje na studii publikovanou v roce 1993: Profiling of hot surfaces by pulsed time-of-flight laser range finder techniques. Kari Määtta, Juha Kostamovaara, and Risto Myllylä[13]. • LaCam Forge od firmy Ferrotron a Menteq Technologies Inc.[14] Tento systém je rovněž určený přímo pro měření žhavých výkovků při otevřeném kování. Je založený na stejném principu jako systém, který publikoval Tian at al. využívá skenování pomocí pulzního polovodičového laseru a vyhodnocování metodou TOF. Před kováním je možné kontrolovat rozložení teploty na povrchu výkovku pomocí teplotní mapy. Při kování je možné měřit průběžně délku výkovku bez nutnosti přerušení výrobního procesu a softwarově optimalizovat proces kování. Po kování je možné zobrazit finální rozměry, přímost osy a válcovitost. Standardní proces kování s využitím měření pomocí LaCam systému je zobrazen na obr. 2-2. Obr. 2-3 zobrazuje výstupní protokol měření.
Obr. 2-2: Proces optimalizace tvaru výkovku pomocí laserového systému LaCam od firmy Ferrotron [14].
strana
12
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-3:Protokol z měření pomocí systémuLaCam. Z výsledků je patrné, že maximální průhyb překračuje povolenou odchylku 10 mm [14].
Systémy založené na vyhodnocování metodou triangulace: • TopScan od firmy Tecnogamma [15]. Systém určený přímo pro měření výkovků vysoké teploty je založen na vyhodnocování laserového paprsku metodou triangulace. Je určen k měření výkovků o teplotě 350 až 1150 °C a maximálních rozměrech 7 m. Přesnost dosahuje ± 5 mm, ve sférických souřadnicích ± 3 mm. Výstupem je obrys výkovku se základními rozměry.
Obr. 2-4:Výstupní data získaná měřením žhavého výkovku laserovým systémem TopScan, který nabízí firma Tecnogamma [15].
strana
13
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Laserový systém vyvinutý ve spolupráci firem IMEGO, Mefos a SINTEF[16] Tento systém byl vyvíjen od roku 2005 díky spolupráci tří švédských ocelárenských společností. Skládá se z laserového projektoru a spektrálně - selektivní videokamery. Modrý laser promítá na povrch výkovku kontrastní pruhy, jejichž poloha je vyhodnocována metodou triangulace. Při měření manipulátor otáčí s výkovkem a kamera vytvoří několik překrývajících se snímků s vysokým rozlišením. Snímky slouží k sestavení 3D modelu (viz obr. 2-5). Systém umožňuje měření průměru v libovolné vzdálenosti, měření délky, křivosti osy a odchylek od požadovaného tvaru.
Obr. 2-5:Výsledky měření pomocí systému od firemIMEGO, Mefos a SINTEF. Vlevo: výsledný 3D model výkovku s odchylkami od požadovaného tvaru, vpravo řez výkovkem [16].
3.1.2 Aktivní optické systémy s xenonovým světlem • Zhenyuan Jia, Bangguo Wang, Wei Liu, Yuwen Sun. 2009 - An improved image acquiring method for machine vision measurement of hot formed parts [17]. Tento článek poprvé popisuje využití spektrálně selektivní metody pro účely měření žhavých výkovků. Celý výkovek je osvětlen pomocí xenonové lampy. Systém se skládá ze dvou monochromatických CCD snímačů s filtrem s dolní propustí, xenonové lampy a obslužného počítače (viz obr. 2-6).
Obr. 2-6: Prvky optického systému s xenonovým světlem [17].
strana
14
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-7: Princip spektrálně-selektivní metody. Levá svislá osa udává rozsah citlivosti CCD snímače. Pravá osa udává sílu radiace dokonale černého tělesa. Vodorovná osa udává vlnovou délku. Šedá oblast znázorňuje spektrum, které bude snímačem detekováno [17].
Princip spočívá v selekci spektra, ve kterém nejsou CCD ovlivněny radiací výkovku a současně jsou citlivé na xenonové světlo (viz obr. 2-7). Na povrch výkovku je promítáno xenonové světlo, které se odráží k CCD snímačům, přičemž IR záření emitované výkovkem je odfiltrováno. Díky tomu mohou být na povrchu výkovku extrahovány body, které slouží ke složení výsledného 3D modelu ze stereosnímků (obr. 2-8).
Obr. 2-8: Výsledek extrakce bodů na povrchu výkovku [17].
Tento článek poprvé uvedl metodu extrakce bodů na povrchu žhavého výkovku. Přesnost měření rozměrů výkovku pomocí tohoto systému však nebyla popsána. Studie z následujících let uvádí, že extrakce tak velkého množství bodů vede k nízké rychlosti systému, z důvodu výpočetní náročnosti. Na tento výzkum navazuje Liu [18], který prezentoval vývoj „fringe projection“ 3D skeneru s xenonovým světlem založený na spektrálně selektivní metodě (viz další odstavec).
strana
15
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Wei Liu, Xinghua Jia, Zhengyuan Jia, Shuangjun Liu, Bangguo Wang, Jian Du. 2011 - Fast dimensional measurement method and experiment of the forgings under high temperature[18]. Tento článek vychází ze studie Jia at al. [17]. Na rozdíl od metod publikovaných v [17] využívá tento systém vyhodnocování promítaného světelného vzoru. Systém se skládá ze dvou CCD kamer s IR filtry, projektoru a obslužného počítače, jak ukazuje následující obrázek.
Obr 2-9: Vlevo: schéma optického stereo systému, vpravo jeho reálná podoba [18].
K detekci promítaných pruhů byl použit modifikovaný Gaussův model a spektrálně selektivní metoda popsaná již v roce 2009 [17]. Intenzita světelných pruhů v kolmém směru nemá na povrchu žhavého objektu typické Gaussovo rozložení. Oproti běžnému průběhu křivky, má křivka za podmínek vysoké teploty zploštělý tvar vrcholu, jak ukazuje obr. 2-10 vpravo. Z toho důvodu byl použit modifikovaný Gaussův model, který je blíže popsaný v článku.
Obr.k 2-10: Modifikovaný Gausův model aplikovaný na světelné pruhy promítané na žhavý výkovek. Vlevo: světelné pruhy na výkovku, modifikovaný gausův model v místě pruhu [18].
Pomocí první derivace intenzity byla získána funkce, která je zobrazena na obr. 2-11 vpravo. Následně byla druhá derivace této funkce položena rovno nule. Výsledkem jsou kořeny funkce („ridge points“ na obr. 2-11). Poloha středových pixelů promítaných pruhů („bottom point“ na obr. 2-11) byla na závěr vypočítána Newtonovou iterativní metodou se subpixelovou přesností pomocí dvou kořenů funkce.
strana
16
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-11: Aplikace první derivace na obraz výkovku. Vlevo: detail světelných pruhů po konvoluci, vpravo: výsledná funkce představující profil světelného pruhu [18].
Dále bylo nutné určit polohu hraničních bodů na výkovku. V článku jsou popsány dva algoritmy, které slouží k výpočtu polohy koncových bodů na povrchu výkovku. První algoritmus je založen na skutečnosti, že intenzita šedi je v pozadí výkovku značně odlišná, druhý ze skutečnosti, že na rozhraní výkovku a pozadí mají promítané pruhy významnou změnu křivosti. Aby mohly být extrahovány rozměry výkovku, musejí být pruhy promítány v horizontálním i vertikálním směru (viz obr. 2-12). Systém byl testován v laboratorních i reálných podmínkách. Přesnost rozměrů u čtyř měřených vzorků byla nižší než 1%.
Obr. 2-12: Hraniční body detekované na ocelové desce. Vlevo: hraniční body v horizontálním směru, vpravo: hraniční body ve vertikálním směru [18].
3.1.3 Pasivní systémy založené fotogrammetrii a analýze obrazu • S. Chaudhuria, A. Mukherjeeb, P. K. Duttab, A. Patrab. 2003 - Automatic measurement of frontal area and volume of fluid ingots using image-based instrumentation [19]. Automatizovaný fotogrammetrický systém byl v provozu aplikován např. k predikci místa řezu ingotu dle požadované hmotnosti dílčích částí. Ingoty mají 12 hran v průřezu a kónický tvar. Vytvořený software automaticky snímkuje ingot, předzpracuje snímky a pomocí jejich analýzy počítá objem, hmotnost a délku. Článek uvádí možné řešení výpočtu obsahu čelní plochy: • Výpočtem množství pixelů uvnitř hraniční oblasti.
strana
17
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Integrací 0 až 2π radiálního vektoru, který vychází z těžiště hraničních bodů. • Detekcí rohových bodů nebo vrcholů objektu a následnou derivací ze známé geometrie objektu. Poslední dvě metody jsou vhodné spíše pro objekty jednoduchého geometrického tvaru, z toho důvodu byla zvolena první metoda. Článek obsahuje popis použitého algoritmu, jehož hlavními částmi je prahování a afinitní transformace. Průměrná odchylka hmotnosti dosahuje ± 0.5%. Systém byl zaveden v průmyslovém závodě v Indii.
• S. B. Dworkin, T. J. Nye. 2004 - Image processing for machine vision measurement of hot formed parts [20]. Studie popisuje metodu měření základních rozměrů výkovku vysoké teploty pomocí adaptivního prahování a detekce hran. Přitom je testována přesnost měření u různých typů vstupních snímků: • monochromatické snímky ve viditelném spektru, • barevné snímky ve viditelném spektru, • monochromatické snímky s IR filtrem 9. Z článku vyplývá, že nejlepší vlastnosti pro prahování a detekci hran mají monochromatické snímky s IR filtrem. Článek uvádí, že běžné CCD snímače bez filtru nejsou schopné vytvořit vhodné snímky pro detekci hran prahováním a přesné měření rozměrů výkovku. Na následujícím obrázku je monochromatický snímek výkovku o teplotě 900 až 1300 ̊C, který byl pořízen s IR filtrem. Na obr. 2-14 jsou zobrazeny výsledky prahování.
Obr. 2-13:Monochromatický snímek žhavého výkovku při použití IR filtru [20].
Obr. 2-14: Výsledné prahování snímků. Vlevo monochromatický bez IR filtru, vpravo s IR filtrem [20].
9
IR filtr - propouští pouze IR paprsky o vlnové délce vyšší než specifikace daného filtru.
strana
18
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• Ch. Hu, B. Liu, X. Song. 2008 - A novel edge detection approach used for online dimensional measurement of heavy forging[21] Článek představuje novou metodu detekce hran, která překonává problémy ostatních detektorů při aplikaci na výkovky vysoké teploty. Jednotlivé kroky zpracování obrazu byly navrženy tak, aby byly co nejlépe splněny podmínky Cannyho algoritmu (viz kap. 3.3.1). U snímků výkovků o vysoké teplotě byly navíc aplikovány kroky pro předzpracování a finálního zpracování obrazu. Barevný 24 bitový obraz byl převeden na 8mi bitový. Intenzita pixelů byla rozšířena na celou oblast 0 až 255 a světlost obrazu byla zvýšena. Následně byl obraz upraven adaptivním prahováním, které poprvé publikoval Daxing Z. at al. V roce 2008 [22]. Výsledky předzpracování obrazu jsou na obr. 2-15. Některé segmenty pixelů v blízkosti objektu byly nesprávně označeny jako hrany, nicméně hranice okují byly eliminovány, zatímco většina hraničních pixelů bylo zachováno, beze změny polohy. Proces předzpracování vstupního obrazu není standardní krok detekce hran, ale slouží ke zvýšení její účinnosti při použití na žhavé objekty.
Obr. 2-15:Snímky žhavého výkovku ve stupních šedi. Vlevo původní snímek, vpravo předzpracovaný snímek [21].
Na vstupní předzpracovaný obraz byla gradientní metoda detekce hran. Za účelem zlepšení výsledků detekce hran, byl nejdříve aplikován filtr s dolní propustí pro odstranění aditivního šumu. V první fázi tedy byla provedena konvoluce obrazu ve dvou směrech pomocí jádra s Gaussovým rozložením. Ve druhé fázi byl aplikován filtr s horní propustí – diferenciální operátor slouží k nalezení hran v obraze. V poslední fázi byla provedena validace hran, za účelem zachování skutečných hrana eliminace redundantních informací.
Obr. 2-16:Výsledky experimentu s různými metodami. Vlevo: detekce pomocí Cannyho hranového detektoru, vpravo: detekce pomocí prezentované metody bez předzpracování snímku [21].
strana
19
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Provedený experiment dokazuje, že předzpracování a finální úprava obrazu jsou nepostradatelnými kroky při detekci hran výkovků vysoké teploty. Výsledky experimentu a porovnání s výsledky detekce hran pomocí Cannyho detektoru jsou na obr. 2-16 a 2-17.
Obr. 2-17:Výsledky experimentu s různými metodami. Vlevo: detekce pomocí prezentované metody bez finálního zpracování snímku, vpravo: detekce pomocí prezentované metody [21].
strana
20
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ 3.2
3.2 Průmyslové aplikace blízké fotogrammetrie 3.2.1 Rozdělení fotogrammetrických systémů Pozemní fotogrammetrické systémy používané v průmyslu se zpravidla dělí na dvě kategorie – off-line a on-line systémy [23].
Obr.2-18: Proces on-line a off-line fotogrammetrického měření v průmyslu[8].
Nejvýraznějším rozdílem je schopnost měření a vyhodnocení „real-time“ nebo téměř „real-time“, což je výsadou on-line systémů. Typický on-line systém se skládá ze dvou nebo více kalibrovaných fotoaparátů s danou polohou a orientací. Tyto systémy jsou určeny k automatizovanému měření v konkrétním prostoru. Bývají například součástí výrobní linky, kde slouží ke kontrole kvality. Přesnost on-line systémů bývá obvykle nižší než přesnost off-line systémů, neboť jsou omezeny počtem snímků [23, 24]. Off-line systémy jsou určeny pro univerzální měření mnoha různých průmyslových aplikací. Jejich univerzalita spočívá zejména v možnosti použití fotoaparátu a objektivu dle požadované přesnosti a velikosti objektu. Fáze měření a vyhodnocování je oddělená. Off-line systémy dosahují vyšší přesnosti, podle [23] okolo 0,05 mm. Na rozdíl do on-line systémů mohou být použity také analogové fotoaparáty a klasické fotografie. Již studie z roku 1988 [24] uvádí rozdělení průmyslových fotogrammetrických systémů podle použitých technologií.
strana
21
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr. 2-19: Rozdělení průmyslových fotogrammetrických systémů [24].
Nejběžnějšími aplikacemi fotogrammetrie jsou: • měření polohy většího množství cílových bodů, • sledování pozice a orientace dotykové sondy, • měření tvaru plochy [23, 24]. Firmy AICON [27] a Nikon Metrology [28] nabízí kompletní sortiment fotogrammetrických systémů, které pokrývají všechny tři skupiny průmyslových aplikací. 3.2.2 Měření polohy více cílových bodů současně Přesná lokalizace a digitalizace diskrétních bodů je patrně nejběžnější aplikací fotogrammetrie v průmyslu. Obvykle se jedná o systémy off-line měření využívající jeden fotoaparát s vysokým rozlišením, speciální značky a algoritmus „bundle adjustment“ (viz kap. 3.2.4). Na povrch měřeného objektu se aplikují dva typy značek – kódované a nekódované. Oba typy značek mívají obvykle kruhový tvar s bílým středem a černým okrajovým pásem. Kódované značky slouží přímo k digitalizaci diskrétních bodů na povrchu objektu. Nekódované značky mají navíc po obvodě umístěný unikátní znak, který umožňuje automatické rozpoznání každého bodu softwarem. Tyto značky slouží k nalezení geometrických vztahů mezi jednotlivými obrazovými rovinami. Díky tomu je možné určit polohu detekovaných nekódovaných bodů ze všech snímků vůči jednomu souřadnému systému. V okolí objektu musí být umístěna také jedna až dvě měřítkové tyče, které slouží k transformaci obrazu do správného měřítka. Přenos dat do počítače bývá řešen buď paměťovou kartou, nebo bezdrátově. Naměřená data se zpracovávají pomocí specializovaného softwaru ihned po měření nebo s libovolným odstupem času [30]. Přesnost měření je možné ovlivnit hustotou rozmístění značek a počtem pořízených snímků. Digitalizace každé značky deseti až dvanácti snímky zajistí přesnost RMS 1-sigma 10 až 0,025 mm [23]. RMS 1-sigma (root mean square nebo 1 standard deviation) – střední kvadratická chyba s úrovní spolehlivosti 68%. 10
strana
22
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Obr.2-20: Vlevo: fotogrammetrické měření pomocí tzv. kódovaných a nekódovaných značek, vpravo: příklad nekódovaných fotogrammetrických značek [29].
3.2.3 Sledování pozice a orientace dotykové sondy Měření pomocí dotykové sondy se provádí ručním nebo strojním polohováním. Systémy se strojním polohováním se nazývají 3D souřadnicové měřicí stroje, podle anglického „coordinate measurement machine“ (CMM). Tyto systémy se dále dělí na analogové a diskrétní. Diskrétní CMM měří jednotlivé body na povrchu objektu postupně přemisťováním sondy, přičemž je měřena poloha sondy mechanicky. Analogové CMM měří tvar objektu v celé ploše - sonda plynule kopíruje celý měřený povrch. Pohyb sondy je v tomto případě, stejně jako v případě ručního polohování sondou, snímán optickým systémem. Ten zaměřuje orientaci a polohu hrotu sondy a kontrolních bodů, které jsou umístěny na jejím těle. Fotogrammetrický systém může být složen pouze z jednoho optického snímače, ale i ze dvou nebo více snímačů. Obvykle platí, že čím více snímačů sleduje pohyb sondy, tím přesnějších výsledků měřicí systém dosahuje.
Obr. 2-21: Dotykové sondy: vlevo ruční dotyková sonda od firmy Metronor [26], vpravo nahoře ruční dotyková sonda od firmy AICON [27], vpravo dole CMM od firmy Reinshaw [25].
strana
23
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Mezi největší firmy nabízející CMM patří Reinshaw [25], Metronor (dnes spojen s Nikon Metrology) [26, 28] a AICON [27]. Systémy s ručním polohováním nabízí firma Metronor [26], AICON [27] a GOM [29]. 3.2.4 Měření ploch Fotogrammetrické metody měření plochy vyžadují použití minimálně dvou snímků, které jsou složeny do jednoho souřadného systému pomocí detekovaných bodů na povrchu objektu. Měřená plocha musí být opatřena vhodnou texturou, která umožní detekci a skládání vlícovacích bodů. Textura na měřeném objektu může být přirozená nebo uměle vytvořená například přilepením kontrastních pásek nebo projekcí. Hlavní výhodou fotogrammetrického měření ploch je možnost měření dynamické scény. Příkladem je měření, při němž se pohybují snímače, nebo dochází k pohybu či deformaci měřené plochy, jak ukazuje např. obr. 2-23. Častými aplikacemi jsou např. měření deformací karoserie při „crash testech“ nebo při měření funkce airbagů při nárazu. Dynamické změny tvaru plochy není možné měřit žádnou jinou optickou technologií včetně 3D skenování [23]. Fotogrammetrické systémy určené k měření dynamické scény nabízí např. firma GOM [29] nebo AICON [27].
Obr. 2-22: Vlevo: projekce vzoru a fotogrammetrické měření systémem PRO-SPOT [31], vpravo: fotogrammetrický systém určený k měření deformací od firmy GOM [29].
Obr. 2-23:Měření dynamické scény pomocí fotogrammetrického systému od firmy AICON[27].
strana
24
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ 3.3
3.3 Dílčí kroky fotogrammetrického měření Klíčové informace v digitální podobě jsou výsledkem řešení komplexního algoritmu. Přesnost výsledku a rychlost výpočtu závisí na mnoha parametrech. Nedílnou součástí návrhu algoritmu je jeho testování a ladění. Neustále se zvyšující se nároky na přesnost optických systémů, jsou důvodem stále aktivního výzkumu těchto dílčích algoritmů. Cílem výzkumu bývá především zpřesnění výsledné rekonstrukce, ale také snížení výpočetní náročnosti, nebo přizpůsobení algoritmů pro specifické podmínky měření. Kompletní rekonstrukce 3D informací založená na fotogrammetrii obvykle obsahuje tyto dílčí procesy: • • • •
rozpoznání klíčových primitiv (značek), lokalizace polohy klíčových primitiv a výpočet jejich parametrů, kalibrace vnitřních parametrů fotoaparátu, kalibrace vnějších parametrů fotoaparátu - nalezení přesné orientace všech snímků vůči jednomu souřadnému systému.
3.3.1 Rozpoznání klíčových primitiv Metody zpracování digitálního obrazu byly vyvíjeny již v 60. letech 20. století, zejména pro účely satelitního a lékařského snímkování, ale také za účelem zlepšování vizuálních vlastností digitální fotografie. Do 70. lech bránila vývoji těchto metod špatná dostupnost počítačů i jejich výkonnost. Značná pozornost byla již od počátku věnována rozpoznávání objektů v digitálním obraze. První fází rozpoznávání objektů v obraze je detekce hran nebo segmentace obrazu. V oblasti průmyslové fotogrammetrie se využívá výhradně detekce hran [32]. Přehled metod detekce hran v digitálním obraze publikoval v roce 1975Davis [33], v roce 1982 Peli a Malah [34] a v roce 1986Zioua Tabbone [35]. Problematice detekce hran v digitálním obraze je stále věnována pozornost. Studie se v posledních letech zabývají především optimalizací algoritmů pro konkrétní úlohy detekce [36]. Mechanismy detekce je možné rozdělit podle technik, které používají: založené na první derivaci, založené na druhé derivaci, nepracující s derivací, nebo detektory pracující ve frekvenční oblasti. Velká skupina metod aproximuje derivaci pomocí konvoluce s vhodným jádrem. Mezi nejznámější operátory patří: Sobel, Prewitt, Roberts a Laplacian. V současné době je stále považováno za optimální řešení detekce hran Cannyho algoritmus [37], který nejlépe splňuje tři základní kritéria: • Minimální chybovost - poměr mezi množstvím detekovaných hran a množstvím šumu musí být co nejmenší. • Lokalizace – poloha detekované hrany musí být co nejblíže její skutečné poloze v obraze. • Jednoznačná odezva – operátor musí detekovat pouze jednu odezvu na každou hranu. Optimalizacím Cannyho algoritmu se v posledních letech věnoval např. Ding a Goshtasby [38] nebo R. Medina-Carnicer at al. [39]. Detekce hran pomocí Cannyho detektoru dosahuje pouze pixelové přesnosti. V posledních letech jsou vyvíjeny metody, jejichž cílem je sub-pixelová detekce, např. [40].
strana
25
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.3.2 Lokalizace polohy elips a výpočet jejich parametrů Fotogrammetrické systémy využívají detekci hran ke zjištění přesné polohy středů kruhových značek. Je zřejmé, že kruhové značky se v obrazové rovině perspektivní projekcí transformují do tvaru elipsy. Klíčovou oblastí fotogrammetrie je přesná lokalizace těchto geometrických primitiv na detekovaných hranách. Metody pro lokalizaci a výpočet parametrů elipsy je možné rozdělit do dvou hlavních skupin – metody shlukovací („voting/clustering“) a optimalizační. Studii zabývající srovnáním metod publikoval v roce 1994 Shortis at al. [41]. K metodám spadajícím do první skupiny patří algoritmy založené na Haughové transformaci [42-45], RANSAC [46, 47] a fuzzy shlukování [48, 49]. Tyto metody jsou odolné proti vlivu odlehlých pixelů, a umožňují detekci několika primitiv současně. Nevýhodou těchto metod je především výpočetní složitost a nižší rychlost. Druhá skupina metod je založena na optimalizaci cílové funkce. Počítá se optimální poloha elipsy vůči detekovaným hraničním bodům metodu nejmenších čtverců [50-57]. Hlavní výhodou těchto metod je jejich přesnost a rychlost. Nicméně citlivost k odlehlým pixelům je vyšší než v případě shlukovacích metod. V posledních letech je věnována pozornost výhradně této skupině metod. Jiandong at al. publikoval v roce 2009 studii popisující základní kroky algoritmu fotogrammetrické extrakce diskrétních bodů. Studie se hlouběji zabývá algoritmem přesné lokalizace fotogrammetrických značek [58]. 3.3.3Vnitřní kalibrace fotoaparátu Vnitřní kalibrace slouží k nalezení vztahů mezi body 3D scény a jejich průměty do 2D souřadnic snímku. Vnitřní kalibraci fotoaparátu je možné rozdělit na lineární a nelineární. Lineární kalibrace se provádí za účelem určení parametrů zobrazovacího systému, což jsou [59]: • konstanta kamery – ohnisková vzdálenost („principal distance“), • střed zobrazení („principal point“), • formát snímku [59]. Pět lineárních parametrů kamery je vyjádřeno v horní trojúhelníkové matici 3x3: 𝛼𝛼𝑥𝑥 𝐴𝐴 = � 0 0
𝛾𝛾 𝛼𝛼𝑦𝑦 0
𝑢𝑢0 𝑣𝑣0 � 1
Parametry αx= f . mx a αy = f . my reprezentují ohniskovou vzdálenost f v pixelech, kde mx a my jsou konstanty měřítka mezi pixely a vzdáleností. γ představuje faktor úhlové vychýlení osy x a y („skew factor“). A u0, v0 představují vzdálenost středu zobrazení od středu snímku [59]. Nelineárními vnitřními parametry fotoaparátu jsou parametry zkreslení objektivu: • radiální zkreslení způsobené nedokonalou křivostí čoček (posun bodu na snímku o radiální vzdálenost), obsahuje složku zkreslení způsobeného excentricitou čoček („decentring distorsion“),
strana
26
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
• tečné (tangenciální) zkreslení způsobené převážně asymetrií radiálního zkreslení u soustavy čoček (často se zanedbává), • a „thin prism“ zkreslení, které je způsobené naklopením čoček vůči obrazovému snímači [59]. Do dnes bylo vyvinuto a úspěšně implementováno mnoho metod, které teoreticky řeší všechny problémy vnitřní kalibrace. I když je lineární a nelineární kalibrace obvykle řešena současně, pozornost je v posledních dvaceti letech věnována výhradně studiu nelineárního zkreslení čoček, zejména za účelem zpřesnění a zjednodušení kalibrace. Důvodem je to, že najít přesný model nelineárního zkreslení čoček je mnohem komplikovanějším problémem, než hledání lineárních parametrů zobrazovací soustavy. Conrady v roce 1919 poprvé popsal zkreslení způsobené excentrickou polohou čoček. V roce 1966 byl představen model s názvem „plumb line method“, který dělí zkreslení na 3 základní typy (Brown-Conrady model [60, 61]). Tento model umožňoval určit konstantu fotoaparátu a střed zobrazení, nicméně eliminace zkreslení objektivu byla ošetřena pouze částečně. I přesto našel od té doby uplatnění v mnoha aplikacích. Do dnes bylo prezentováno několik modifikovaných metod, nicméně tyto modifikace se zaměřovaly zejména na matematickou interpretaci. Kalibrace čoček byla dosud převážně řešena jako součet zkreslení radiálního, tečného a „thin prism“ – takové metody jsou označovány jako konvenční. Studie [62] z roku 2008 prezentuje model, který řeší jednotlivé složky zkreslení a jejich fyzické příčiny odděleně a uvažuje jejich vzájemné působení. Tento model dosahuje stejného korekčního efektu jako konvenční modely, obsahuje však méně koeficientů a tak zjednodušuje kalibraci. Fotoaparát je považován za vnitřně kalibrovaný, je-li známá jeho ohnisková vzdálenost, střed obrazové roviny a zkreslení objektivu. Podle toho, zda kalibrace využívá kalibrační element, rozlišujeme: • metrická kalibrace – kalibrace založená na digitalizaci kalibračního pole. • nemetrická kalibrace nebo také samokalibrace – kalibrace bez použití kalibračního pole, je řešena přímo pomocí entit ve scéně. Metrická i nemetrická kalibrace je založena na srovnávání polohy daného typu prvků ve zkresleném obraze s jejich skutečnou polohou. Kontrastní prvky jsou ve scéně rozmístěny tak, že tvoří pole prvků. V případě metrické kalibrace je nejjednodušším kalibračním polem čtvercová mřížka (kalibrace pomocí bodů v průsečících), jeho obdobou jsou pravidelně rozmístěné kruhové značky. Používají se také kalibrační pole s trojúhelníkovými nebo čtvercovými značkami. Pole se snímkuje kamerou z různých úhlů ze vzdálenosti, která bude použita při měření. Během kalibrace nesmí být měněna ohnisková vzdálenost a zaostření. Zaostření musí být nastaveno na vzdálenost, ve které bude kamera následně pracovat [80]. Sun a Cooperstock prezentovali v roce 2006 studii zabývající se empirickým vyhodnocením faktorů, které mají vliv na přesnost kalibrace [63]. 3D skenery často využívají před samotným měřením kalibraci metrickou, která zajistí vyšší přesnost systému, než kalibrace nemetrická. Nicméně do průmyslových fotogrammetrických systémů byly úspěšně implementovány metody samoklaibrace pomocí značek rozmístěných na scéně [64 - 67].
strana
27
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
3.3.4 Vnější kalibrace fotoaparátu Kalibrace vnějších parametrů fotoaparátu se provádí se za účelem: • určení transformačních parametrů (tři rotace, tři posunutí, obecně až16 parametrů) - popisují polohu fotoaparátu v prostoru a vyjadřují vzájemný vztah světového souřadného systému a souřadného systému fotoaparátu, • určení relativní orientace fotoaparátů, • určení báze (vzdálenost středů zobrazení dvou fotoaparátů)[59]. Kalibrace vnějších parametrů je stejně jako kalibrace vnitřních parametrů založena na perspektivním promítání, které je zjednodušeno na model dírkového fotoaparátu. 2D pozice bodu v digitálním obraze je často vyjadřována vektorem [u v 1]T, zatímco pozice bodu ve 3D světových souřadnicích rozšířeným vektorem [xw yw zw 1]T. Díky tomu je možné zapsat průmět světových souřadnic do souřadnic digitálního obrazu následovně [59]: 𝑥𝑥𝑤𝑤 𝑢𝑢 𝑧𝑧𝑐𝑐 �𝑣𝑣 � = 𝐴𝐴[𝑅𝑅 𝑇𝑇] �𝑦𝑦𝑤𝑤 � 1 1 A je matice lineárních vnitřních parametrů fotoaparátu, a R, T jsou transformační matice vnějších parametrů kalibrace (rotační matice a translační vektor). Pro určení těchto transformačních parametrů se využívají dva základní modely – lineální a nelineární model [59].
Lineární model Vzájemný vztah mezi fotografiemi je určen epipolární geometrií. Nejdříve se analyzuje vzájemná poloha konkrétních bodů v prostoru a jejich průmět do obrazové roviny. Tento vztah popisuje matice fotoaparátu P. Ta závisí jak na vnitřních, tak i na vnějších parametrech. Je více způsobů, jak získat průmět P a P´,v závislosti na tom, jaká data jsou k dispozici. Nejpoužívanější je však metoda využívající tzv. esenciální matici E nebo metoda přímé lineární transformace („direct linear transformation“DLT), kterou publikoval již v roce 1971 Abdel-Aziz & Karara [68]. Esenciální matice je speciální druh fundamentální matice, který vyžaduje znalost interních kalibračních dat. Esenciální i fundamentální matice byly poprvé popsány Longuet-Higginsem [69]. Jednoznačné řešení pro určení této matice zajišťuje průmět osmi bodů – podle toho název normalizovaný osmibodový algoritmus. Existuje však mnoho dalších algoritmů, z nichž většina v některé fázi využívá osmi-bodový algoritmus [70, 71, 72]. V posledních letech byly prezentovány také metody, které řeší lineární kalibraci pouze pomocí sedmi nebo šesti bodů. Pokud je určena esenciální matice je možné vypočítat transformační parametry matice R a T metodou singulárního rozkladu. Studii, která se zabývá komplexním problémem určení orientace kamer pomocí algoritmů fundamentální a esenciální matice prezentoval v roce 1996 Luongaa Faugeras [73]. Výhodou lineárního modelu je možnost kalibrace s proměnnou ohniskovou vzdáleností. Kalibrace je poměrně rychlá a jednoduchá, nicméně vyžaduje znalost více bodů oproti kalibraci nelineární.
strana
28
SHRNUTÍ SOUČASNÉHO STAVU POZNÁNÍ
Nelineární model V oblasti blízké fotogrammetrie je často využívána nelineární kalibrace pomocí algoritmu, který je známý jako „Bundle Adjustment“. Název vychází z představy paprsků světla („bundle“), které se odrážejí od 3D struktury a konvergují v obrazovém středu každého fotoaparátu. Klasické „bundle adjustment“ algoritmy jsou iteračními metodami, které se řeší jako nelineální problém nejmenších čtverců. K jednoznačnému řešení tohoto modelu stačí znalost minimálně pěti bodů [74, 75]. Cílem celého algoritmu je nalézt globálně konzistentní skupinu parametrů, které minimalizují chybu registrace mezi všemi páry snímků. Minimalizace se často řeší metodou Lavenberg-Marquardt [76] nebo Sampsonovou aproximací [74].
Obr.2-24: Orientace snímků pomocí algoritmu „Bundle adjustment“ [79].
Algoritmus „Bundle adjustment“ byl představen v 70. letech 20. století – Brown, 1974 [77], Granshaw, 1980 [78]. Triggs publikoval v roce 2000 studii, která se zabývá přehledem algoritmů „bundle adjustment“ a popisuje také nové přístupy založené na distribučních modelech [79]. Výhodou tohoto algoritmu je nutnost znalosti pouze pěti bodů. „Bundle adjustment“ bývá navíc často kombinován se samokalibrací interních parametrů, a je tedy vhodný pro rychlé fotogrammetrické měření v průmyslu. Srovnání metod kalibrace prezentoval v roce 2006 Remondino a Fraser [80].
strana
29
ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE
4 ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE V současné době na trhu množství sofistikovaných automatizovaných systémů bezdotykového měření tvaru a rozměrů objektů. Přesto je výzkum v této oblasti stále velmi aktivní. Z rešerše bylo zjištěno, že se výzkum zaměřuje především na zvyšování přesnosti, spolehlivosti a rychlosti již zavedených systémů a jejich úpravou pro specifické pracovní podmínky. Podstata výzkumu stojí na vývoji nových matematických modelů a komplexních algoritmů, jejich testování, optimalizaci a implementaci. Na základě rešerše bylo také zjištěno, že v poslední dekádě je věnována pozornost vývoji specifických bezdotykových systémů určených pro účely měření objektů o vysoké teplotě. Zejména v těžkém ocelárenském průmyslu jsou stále rezervy v efektivitě výroby a to především z důvodu absence pokročilých technologií měření tvaru a rozměrů za tepla. Důvody nutnosti a potřebnosti řešení dané problematiky v daném čase jsou zmíněny v následujících bodech: • S rostoucími nároky na efektivitu výroby, rostou také nároky na tvarovou a rozměrovou přesnost výkovků. Z toho vyplývají také nároky na přesné, rychlé a bezpečné měření tvaru a rozměrů. Výzkumný záměr koresponduje se současnými problémy efektivity výroby v ocelárenském průmyslu (např. společnost Žďas a.s., která se zabývá výrobou rozměrných výkovků). • Dosud se běžně řeší měření rozměrných výkovků vysoké teploty pomocí etalonů či velkých posuvných měřidel ručně. Obsluha se tak musí pohybovat v blízkosti žhavého výkovku, což může být příčinou úrazů. Navíc jsou získány pouze neúplná a nepřesná data. Dotykové měření výkovků neposkytuje dostatečnou flexibilitu, přesnost ani rychlost měření. • Běžné bezdotykové optické měřicí systémy, jsou pro účel měření výkovků stěží aplikovatelné, zejména z důvodu vysoké teploty a velkých rozměrů výkovků. Z toho vyplývá nutnost vývoje specifického bezdotykového systému. • Možnosti měření žhavých objektů pomocí fotogrammetrie a analýzy obrazu jsou dosud málo prozkoumány a komerční systémy, pokud je autorce známo, zatím neexistují, na rozdíl od systémů laserových. • Výhodami fotogrammetrického systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu výkovku a v neposlední řadě také cena. Nároky na specifické systémy měření objektů vysoké teploty jsou především: • dostatečná vzdálenost optického snímače (eliminace vlivu IR emisí výkovku), • rychlé měření bez nutnosti zásahu do výrobního cyklu výkovku, • online systém – automatické „real-time“ vyhodnocování (viz také kap. 4.2). Specifické systémy pro zmíněné aplikace jsou obvykle založeny: • na laserovém skenování, • na optickém skenování a promítání xenonového světla, nebo • na principu fotogrammetrie a analýzy obrazu.
strana
30
ANALÝZA A INTERPRETACE A ZHODNOCENÍ POZNATKŮ ZÍSKANÝCH NA ZÁKLADĚ KRITICKÉ REŠERŠE
Laserové systémy dosahují vysoké přesnosti a komplexnosti měřených dat, jejich hlavní nevýhodou velmi vysoká cena. Možnosti měření žhavých objektů pomocí fotogrammetrie a analýzy obrazu jsou dosud málo prozkoumány a komerční systémy, pokud je autorce známo, zatím neexistují, na rozdíl od systémů aktivních. Pro řešení problému měření válcovitosti výkovku není třeba získat komplexní informace o jeho tvaru, což je jedním z důvodů vývoje systému založeném na principu fotogrammetrie a analýzy obrazu. Výhodami takového systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu.
strana
31
VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ
5 VYMEZENÍ CÍLE DISERTAČNÍ PRÁCE A NÁVRH ZPŮSOBU JEJÍHO ŘEŠENÍ Cílem doktorské disertační práce je navrhnout systém, který bude automaticky vyhodnocovat odchylky osy válcových výkovků od ideální přímé osy. Přitom musí být splněny nároky na dosahovanou přesnost – maximální odchylka ±5mm při měření výkovku o délce do 10 m. Cílem práce je také experimentálně stanovit rozdíl dosahované přesnosti při měření vzorků pokojové a vysoké teploty. Pro účely experimentálního měření a ladění algoritmu bude zkonstruován jednoduchý měřící stůl. Na základě globálních cílů byly stanoveny dílčí cíle řešení: 1) Navrhnout experimentální měřící aparaturu. 2) Navrhnout a realizovat software pro automatizované měření odchylek výkovků pokojové teploty. 3) Porovnat odchylky měření systému s měřením pomocí certifikovaného optickým 3D skenerem Atos Triple Scan. 4) Modifikovat software pro měření odchylek výkovků vysoké teploty. 5) Ověřit přesnost systému při měření vzorku pokojové a vysoké teploty. 6) Popsat hlavní vlivy na dosahovanou přesnost. Na základě rešerše bude provedena změna názvu dizertační práce. Nový název: Vývoj metodiky vyhodnocování tvarových odchylek výkovků pomocí fotogrammetrie a analýzy obrazu. V průběhu řešení se předpokládá spolupráce s dalšími institucemi. Protože téma disertační práce je z části složitý matematický problém, autorka by chtěla využít spolupráci ve formě konzultací s ústavem matematiky naší fakulty. Spolupráce při řešení algoritmů a samotného programu se předpokládá s Ing. Radkem Poliščukem, Ph.D.
strana
32
SOUČASNÝ STAV ŘEŠENÍ 6
6 SOUČASNÝ STAV ŘEŠENÍ 6.1 Předběžné měření v reálných podmínkách provozu Před návrhem metodiky automatizovaného měření bylo provedeno testovací měření v reálných podmínkách provozu. Ve společnosti Ždas a.s. byly vytvořeny snímky žhavých i chladných výkovků o délce 3000 m a průměru 60 cm. Fotoaparáty byly umístěny na připraveném stativu se dvěma rameny (viz obr. 6-1).
Obr. 6-1: Navržený stativ pro testovací měření v reálných podmínkách.
Ve snímcích byly manuálně detekovány hrany výkovků (software Adobe Photoshop). Následně byly hrany extrahovány a přeneseny na bílý podklad. Binární obraz byl následně vektorizován. Vzniklé 2D vektorové křivky byly načteny do 3D modeláře. Díky známé orientaci kamer byly křivky z obou snímků složeny do jednoho souřadného systému. Ze čtyř hran výkovku byla vypočítána 3D křivka jeho osy. Ta byla srovnána s ideální osou vedoucí koncovými body 3D křivky. Na závěr bylo manuálně nalezeno místo maximálního průhybu, určena jeho velikost, poloha a orientace.
Obr. 6-2: Měření výkovku, detekce hran v rovině A-A.
Obr. 6-3: Měření výkovku, detekce hran v rovině B-B.
strana
33
SOUČASNÝ STAV ŘEŠENÍ
Obr. 6-3: Graf měření osy v rovině A-A. Maximální odchylka D= – 9 mm, ve vzdálenosti z=1551 mm.
Obr. 6-4: Graf měření osy v rovině B-B. Maximální odchylka D= – 5,7 mm, ve vzdálenosti z=1726 mm.
Obrazová rovina A-A
Obrazová rovina B-B
Obr. 6-5: Čelní pohled na měřené osy výkovku v obou rovinách. Fialová křivka je vypočtená 3D osa s maximálním průhybem 10,5 mm ve vzdálenosti z = 1711 mm pod úhlem α = - 76°.
Při testovacím měření byly ověřeny podmínky reálného provozu. Problémem při snímkování byl zejména nestabilní a nerovný povrch, na kterém byl stativ postaven. Další problémy vznikaly při rychlém manuálním vyhodnocení – zejména nedostatečné informace o geometrii scény (nebyla provedena kalibrace), z toho vyplývající problémy při skládání snímků a vyrovnání perspektivy obrazů. Předběžné měření v reálných podmínkách potvrdilo, že fotogrammetrické měření a zpracování obrazu je možným řešením vyhodnocování tvarových odchylek výkovků za tepla i za studena. K dosažení přesných a správných výsledků je nutné dodržet několik základních předpokladů: 1) kontrastní pozadí výkovku ( v místě prizmat mohou metody detekce selhávat, proto je vhodné použít prizmata s malou šířkou), 2) zajistit dostatečně stabilní a tuhý stativ, pro uchycení kamer a nastavení jejich orientace, 3) provedení kalibrace systému přímo na místě.
strana
34
SOUČASNÝ STAV ŘEŠENÍ
Výsledky měření však nelze považovat za směrodatné, protože nebyly porovnány přesnějším měřením. Vývoj metodiky měření musí byt provedeno za přesně vymezených a kontrolovaných podmínek s vyloučením intervenujících proměnných. Proto bude vývoj měřicího systému probíhat v laboratorních podmínkách, kde bude systém testován a výsledky srovnávány s výsledky měření pomocí dalších bezdotykových metod. 6.2 Metodika Na základě rešerše byla navržena předběžná metodika řešení problému. Metodika se bude dále rozvíjet a modifikovat dle výsledků laboratorních experimentů. Hlavní část práce bude zaměřena na návrh algoritmu a jeho optimalizaci pro dané pracovní podmínky. Testovaný vzorek bude nejdříve měřen za běžné pokojové teploty. Po odladění algoritmu pro měření za studena, bude modifikován pro měření vzorku vysoké teploty.
Obr. 6-6: Diagram navržené metodiky řešení.
strana
35
SOUČASNÝ STAV ŘEŠENÍ
6.2 Návrh a realizace měřicí stanice K testování a optimalizaci metodiky v laboratorních podmínkách byla vytvořena měřicí stanice. Stanice se skládá ze tří základních částí – pracovního stolu, ramen pro uchycení fotoaparátů a prizmat pro vymezení polohy měřených vzorků. Měření je prováděno pomocí dvou fotoaparátů s vysokým rozlišením. Poloha odnímatelných fotoaparátů bude pevně dána konstrukcí ramen na robustním stole, v jehož ose budou umístěna prizmata s měřeným vzorkem. Pracovní stůl musí mít v celé ploše pod měřeným vzorkem kontrastní barvu, tzn. bílou. Vyrobená měřicí stanice je zobrazena na obr. 6-7.
y
x
z
Obr. 6-7: Vyrobená měřicí stanice s naznačeným objektovým souřadným systémem.
Konfigurace měřicího stolu a parametry měřených objektů: • • • • • •
úhel mezi obrazovými rovinami fotoaparátu cca 90°, poloha měřeného objektu ve středu obrazové roviny, zvolené fotoaparáty: Canon EOS 500D, použité ohnisko objektivu 18 mm, měřené vzorky – válcové tyče o délce do 700 mm a průměru 50 až 100 mm, požadovaná velikost zorného pole v rovině měřeného objektu 700 mm.
Parametry fotoaparátu Canon EOS 500D: • • • • •
rozlišení: 15,1 Mpix - 4752×3168 pixel2 velikost čipu (h×v): 22,3x14,9 mm2 velikost pixelu: 4,69×4,70 μm2 crop faktor: 1,6 Objektiv s proměnnou ohniskovou vzdáleností: 18 až 55 mm
strana
36
SOUČASNÝ STAV ŘEŠENÍ
Při návrhu měřicí stanice byly uvažovány nároky na její maximální velikost. Ta se řídí velikostí měřeného vzorku, vzdáleností fotoaparátů s daným zorným úhlem při zvolené ohniskové vzdálenosti. Z důvodu úspory prostoru byla volena nejkratší dostupná ohnisková vzdálenost daného objektivu f18 mm. Z velikosti měřeného vzorku byla určena požadovaná velikost zorného pole v horizontálním směru zxo_pož ≥ 720 mm a vypočtena minimální potřebná vzdálenost obrazové roviny Dmin: 𝐷𝐷𝑚𝑚𝑚𝑚𝑚𝑚 = kde:
𝑧𝑧𝑧𝑧𝑜𝑜_𝑝𝑝𝑝𝑝 ž × 𝐹𝐹 720 × 18 − 𝐹𝐹 = − 18 = 563 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 22,3
Dmin
(mm)
- minimální potřebná obrazové roviny od měřeného vzorku
f
(mm)
- zvolené ohnisko
sxo
(mm)
- velikost čipu v horizontálním směru
zXo_pož
(mm)
- požadovaná velikost zorného pole v horizontálním směru
Měřicí stanice byla vyrobena s rameny o délce 400x400 mm, vzdálenost obrazové roviny od měřeného objektu D = 570 mm. Zorné pole je při F18 v této vzdálenosti 486x728 mm2. Teoreticky dosažitelná přesnost pro skutečnou velikost zorného pole a rozlišení fotoaparátů bude porovnána se skutečně dosaženou přesností po kalibraci systému: 𝑃𝑃𝑣𝑣_𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 = 𝑃𝑃ℎ_𝐷𝐷𝐷𝐷𝐷𝐷𝐷𝐷 =
kde:
𝑝𝑝𝑝𝑝𝑜𝑜 × 𝑧𝑧𝑧𝑧𝑜𝑜 0,00470 × 486 = = 0,1533 𝑚𝑚𝑚𝑚 ≅ 0,15 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 14,9
𝑝𝑝𝑝𝑝𝑜𝑜 × 𝑧𝑧𝑧𝑧𝑜𝑜 0,00469 × 728 = = 0,1531 𝑚𝑚𝑚𝑚 ≅ 0,15 𝑚𝑚𝑚𝑚 𝑠𝑠𝑠𝑠𝑜𝑜 22,3
pyo_Dexp (mm)
- velikost pixelu ve vertikálním směru ve vzdálenosti Dexp
pxo_Dexp (mm)
- velikost pixelu v horizontálním směru ve vzdálenosti Dexp
pyo
(μm)
- velikost pixelu ve vertikálním směru
pXo
(μm)
- velikost pixelu v horizontálním směru
syo
(mm)
- velikost čipu ve vertikálním směru
sxo
(mm)
- velikost čipu v horizontálním směru
zyo
(mm)
-velikost zorného pole ve vertikálním směru
zxo
(mm)
- velikost zorného pole v horizontálním směru
Pro reálné podmínky musí být voleny průmyslové kamery, které jsou teplotně odolnější, robustnější konstrukce, monochromatické, bez zabudovaného IR filtru. Průmyslové kamery s vysokým rozlišením používají čipy Kodak KAI. V současné
strana
37
SOUČASNÝ STAV ŘEŠENÍ
době je čipem s nejvyšším rozlišením v této řadě čip KAI-29050 11 s 28,8 Mpix. Například průmyslové monochromatické kamery JAI AM-1600 GE 12, nebo Allied Vision Technology Pike F-1600 13 s čipem Kodak KAI-16000 mají rozlišení 16 Mpix (4872x3248 pixelů, velikost pixelu 7,4×7,4 µm). Mohou pracovat při teplotách -5 °C až 45 °C.
Obr. 6-8: Průmyslová monochromatická kamera AM-1600 GE s rozlišením 4872×3248 pixelů od společnosti JAI.
Obr. 6-9: Průmyslová monochromatická kamera Pike F-1600 s rozlišením 4872×3248 pixelů od společnosti Allied Vision Technology.
6.3 Návrh a realizace kalibrace Byla zvolena metrická kalibrace (viz kap. 3.3.3). Postup při kalibraci je možné rozdělit do dvou základních bodů: • snímkování, • zpracování. Snímkování Před samotným snímkováním je třeba zvolit vhodné kalibrační pole podle vzdálenosti snímkování a typu objektivu. Při měření výkovků v reálných podmínkách by musela být provedena kalibrace na vzdálenost 10 m. Je tedy nutné zvolit takové kalibrační pole, které je možné použít i pro velké rozměry zorného pole. Proto byl zvolen jako základní kalibrační prvek lanko s pravidelně rozloženými ryskami s roztečí 50 mm v horizontálním směru (osa x). Ve scéně byly také umístěny měrky ve zbývajících dvou směrech (v ose y a z). Orientace fotoaparátů byla zvolena tak, aby střed obrazové roviny ležel přibližně na značce ve středu scény. Scéna byla snímkována vždy jednou kamerou pouze z jednoho úhlu. Pole v ose x sloužilo k určení radiálního zkreslení obrazu (viz. kap. 3.3.3), měrky v ose y a z, sloužily k určení orientace obrazové roviny. CCD čipy Kodak KAI, viz http://www.truesenseimaging.com/markets/machine-vision monochromatická kamera JAI AM-1600 GE, viz http://eshop.integrys.com/imaging-and-video/category-4/model/16mp-en.html 13 Allied Vision Technology Pike F-1600, viz http://www.alliedvisiontec.com/us/products/cameras/firewire/pike.html 11 12
strana
38
SOUČASNÝ STAV ŘEŠENÍ
Obr. 6-10: Snímek pro kalibraci s kalibračními měrkami a obrazovým souřadným systémem. Zobrazená mřížka kopíruje deformaci obrazu způsobenou radiálním zkreslením.
Zpracování Nejdříve byla určena deformace obrazu radiální zkreslením. Ostatní komponenty nelineárního zkreslení (tangential, thin prism) nebyly řešeny z důvodu jejich zanedbatelného vlivu na přesnost. Střed obrazové roviny byl přiřazen středové značce na lanku. Dále byl určen počet pixelů mezi jednotlivými značkami v horizontálním směru a srovnán s reálnou roztečí značek. Odchylky způsobené radiálním zkreslením jsou patrné z tabulky 1. - počet pixelů mezi jednotlivými značkami Δpx klesá s rostoucí vzdáleností od optické osy. Na obr. 6-5 je zobrazena mřížka odpovídající radiálnímu zkreslení obrazu. Z tabulky je také patrné, že v blízkosti optické osy odpovídá velikost pixelu teoretické vypočtené hodnotě. Tab. 1: Vyhodnocení odchylek způsobených radiálním zkreslením: mm – skutečná poloha značek v (mm), Δpx počet pixelů mezi jednotlivými značkami, Δmm - skutečná rozteč značek v (mm), Δmm/Δpx – velikost pixelu mezi jednotlivými značkami v (mm).
Kalibrace kamera A Jednotky / značka Č.1 Č.2 Č.3 Č.4 Č.5 Č.6 Č.7 Č.8
Kalibrace kamera B
px
mm
Δpx
Δmm Δmm/Δpx px
-70 254 566 879 1183 1480 1773 2062
-10,9382 39,0618 89,0618 139,0618 189,0618 239,0618 289,0618 339,0618
324 312 313 304 297 293 289
50 50 50 50 50 50 50
0,1543 0,1603 0,1597 0,1645 0,1684 0,1706 0,1730
-4 322 641 959 1275 1581 1887 2190
mm
Δpx Δmm
Δmm/Δpx
-0,6250 49,3750 99,3750 149,3750 199,3750 249,3750 299,3750 349,3750
0 326 319 318 316 306 306 303
0 0,1534 0,1567 0,1572 0,1582 0,1634 0,1634 0,1650
0 50 50 50 50 50 50 50
strana
39
SOUČASNÝ STAV ŘEŠENÍ
Dalším krokem byla vnitřní kalibrace lineárních parametrů zobrazovacího systému, tj. výpočet polohy a orientace obrazové roviny. V obraze byla odečtena velikost některých prvků scény s kalibračním lankem - viz obr. 6-6 a tab. 2.
Obr. 6-11: Schéma geometrie systému v rovině zx (zOxO). Červeně je značen obrazový souřadný systém, modře objektový. Tab. 2: Velikost některých prvků scény odečtená ze snímků.
kamera A
kamera B
55,3 87,0 86,0
53,5 82,0 86,0
L (mm) M (mm) H (mm)
Průmět středu stolu v optických rovinách Průmět středové značky kalibračního lanka na stole Výška lanka nad stolem v klidu
Ze známých rozměrů ve scéně byl dopočítán úhel α natočení obrazové osy zo od objektové osy z a další pomocné úhly ve scéně: 𝐻𝐻 180 𝛼𝛼 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 𝑀𝑀 𝜋𝜋
86 180 ≐ 44,67° 𝛼𝛼𝐴𝐴 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 87 𝜋𝜋
86 180 ≐ 46,36° 𝛼𝛼𝐵𝐵 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � � × 82 𝜋𝜋
𝛾𝛾 = 90 − 𝛼𝛼 − 𝛽𝛽
𝛾𝛾𝐴𝐴 = 90 − 44,67 − 39,79 = 5,54°
𝛾𝛾𝐵𝐵 = 90 − 39,79 − 38,27 = 5,37°
strana
40
(𝐿𝐿 × 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠) 180 𝛽𝛽 = 𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎 � �× (𝐻𝐻 − 𝐿𝐿 ∗ 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐) 𝜋𝜋
(55,3 × 𝑠𝑠𝑠𝑠𝑠𝑠 44,67) 180 𝛽𝛽𝐴𝐴 = � �× ≐ 39,79° 86,0 − 55,3 × 𝑐𝑐𝑐𝑐𝑐𝑐 44,67 𝜋𝜋
(53,5 × 𝑠𝑠𝑠𝑠𝑠𝑠 46,36) 180 𝛽𝛽𝐵𝐵 = � �× ≐ 38,27° 86,0 − 53,5 × 𝑐𝑐𝑐𝑐𝑐𝑐 46,36 𝜋𝜋
SOUČASNÝ STAV ŘEŠENÍ Tab. 3: Výsledky výpočtů základních úhlů geometrie systému v rovině zx.
Jednotky/ úhel α β γ
Kamera A (rad) (deg) 0,78 44,67 0,69 39,79 0,10 5,54
Kamera B (rad) (deg) 0,81 46,36 0,67 38,27 0,09 5,37
popis sklon osy zo od osy z objektového souřadného systému sklon osy zo od osy y objektového souřadného systému úhel mezi optickou osou zo a průmětem středu stolu
Vzdálenost mezi průnikem optických os a obrazovými rovinami, byla vypočítána podle vzorce (viz obr. 6-11): 𝜋𝜋 − 𝛾𝛾 2 𝜋𝜋 𝐷𝐷𝐴𝐴 = 55,3 × 𝑡𝑡𝑡𝑡𝑡𝑡 − 5,54 ≐ 570,38 𝑚𝑚𝑚𝑚 2 𝐷𝐷 = 𝐿𝐿 × 𝑡𝑡𝑡𝑡𝑡𝑡
𝐷𝐷𝐵𝐵 = 53,5 × 𝑡𝑡𝑡𝑡𝑡𝑡
𝜋𝜋 − 5,37 ≐ 569,56 𝑚𝑚𝑚𝑚 2
Kontrolní měření Na závěr bylo provedeno měření k předběžnému stanovení dosahované přesnosti. Kalibrační lanko bylo zatíženo, přičemž průhyb změřený ve směru objektové osy y byl přibližně 7,5 mm. K nastavení průhybu lanka byla využita měrka s jemnou stupnicí (viz obr. 6-12). Nicméně umístění lanka na dané rysce bylo ověřeno pouhým okem. Z toho vyplývá, že z tohoto předběžného měření nelze vyvozovat obecné závěry o přesnosti, slouží tedy pouze orientačně.
Obr. 6-12: Snímek scény se zatíženým lankem.
Ze snímku se zatíženým lankem byla odečtena poloha jeho průmětu do optické osy yo (viz tab. 4). Pomocí následujících vztahů byla vypočtena poloha lanka v objektové ose y a odchylka ω vypočtené hodnoty od skutečné:
strana
41
SOUČASNÝ STAV ŘEŠENÍ
𝑌𝑌𝑜𝑜 cos 𝛼𝛼 5,3 = ≐ 7,5 𝑚𝑚𝑚𝑚 cos 44,67 5,1 = ≐ 7,4 𝑚𝑚𝑚𝑚 cos 46,36
𝜔𝜔 = 𝑌𝑌 − 𝑌𝑌𝑣𝑣
𝑌𝑌𝑣𝑣 = 𝑌𝑌𝑣𝑣𝑣𝑣
𝑌𝑌𝑣𝑣𝑣𝑣
𝜔𝜔𝐴𝐴 = −7,5 − (−7,45) = 0,05 ≐ 0 𝑚𝑚𝑚𝑚
𝜔𝜔𝐵𝐵 = −7,5 − (−7,39) = 0,11 ≐ 0,1 𝑚𝑚𝑚𝑚
Tab. 4: Výsledky předběžné kontroly přesnosti po kalibraci systému.
Jednotky/ označení Y Yo Yv ω
Kamera A (mm)
Kamera B (mm)
-7,5 -5,3 -7,45 0
-7,5 -5,1 -7,39 0,1
popis skutečný průhyb lanka průmět průhybu lanka do optické osy y vypočtený průhyb lanka v objektové ose y odchylka měření zaokrouhlená na desetiny
6.3 Návrh a realizace detekce hran Pro automatizovanou detekci hran výkovků byla s ohledem na charakter analyzovaného obrazu navržena gradientní metoda, validující jasovou hodnotu každého bodu analyzovaného řezu příčným mediánovým filtrem. Motivací pro toto řešení byla snaha o snížení chybovosti, dané vlivem impulzního šumu v obraze. Ten je zde představován zejména přítomností různých kontrolních značek a povrchových nedokonalostí na zkoumaném objektu (rýhy, okuje) a v jeho okolí. Nastavení vhodné šířky plovoucího filtračního okna a prahovacích parametrů bude podrobně popsáno v disertační práci. V současné době probíhá testování metody na snímcích zdeformované tyče. V případě selhání automatického vyhledávání je možné chybně identifikované řezy v uživatelském rozhraní programu buď vyřadit, nebo opravit manuálně.
Obr. 6-13: Snímek deformované válcové tyče po automatické detekci hran.
strana
42
SOUČASNÝ STAV ŘEŠENÍ
Výsledkem detekce hran je soubor hraničních bodů, z nichž je vypočítána 2D osa válcové tyče. Díky známé orientaci obou obrazových rovin bude proveden přepočet dvou 2D křivek na 3D křivku osy. Následně bude provedena manuální kontrola přesnosti výsledné 3D křivky porovnáním s výsledky měření pomocí systému Atos III triple scan. Výsledek rekonstrukce tvaru válcové tyče, která byla měřena pomocí zmíněného průmyslového skeneru je zobrazen na následujícím obrázku. Skutečná osa výkovku byla porovnána s ideální osou vedenou koncovými body skutečné osy.
Obr. 6-14:Výsledek měření a softwarové rekonstrukce deformované tyče pomocí systému Atos s určeným místem maximálního průhybu.
strana
43
ZÁVĚR
7 ZÁVĚR Pojednání ke státní doktorské zkoušce shrnuje současný stav poznání v oblasti bezdotykových optických metod měření se zaměřením na metody fotogrammetrie, jejich aplikace a základní kroky fotogrammetrického měření. Zaměřuje se blíže na možnosti měření tvaru rozměrných válcových výkovků za vysoké teploty. Z literatury vyplívá, že metody měření žhavých objektů jsou založeny na laserovém skenování, optickém skenování a promítání xenonového světla, nebo na principu fotogrammetrie a analýzy obrazu. Přitom v současné době se běžně provádí měření pomocí etalonů či velkých posuvných měřidel ručně. Problémem je vysoká cena a dostupnost specifických laserových systémů, které jsou v současné době jedinou komerčně vyráběnou technologií pro měření žhavých výkovků. Tyto systémy zajistí měření komplexního tvaru výkovku. Pro účely měření základních geometrických odchylek tvaru válcových výkovků (přímost osy, válcovitost) však není nutné znát komplexní informace o geometrii celého objektu. Měření založené na principech fotogrammetrie a analýzy obrazu je adekvátním řešením problému, což vyplývá také z rešerše. Výhodami fotogrammetrického systému by měla být snadná instalace a implementace bez nutnosti zásahu do výrobního cyklu výkovku a v neposlední řadě také cena. Vývoj a testování metodiky probíhá v laboratorních podmínkách z důvodů zajištění vymezených a kontrolovaných podmínek měření s vyloučením intervenujících proměnných. V současné době je řešení práce ve fázi testování algoritmu detekce hran u chladných válcových tyčí. V další fázi proběhne vyhodnocení přesnosti navrženého systému srovnáním výsledků s měřením pomocí průmyslového skeneru Atos III Triple Scan. Navržená metodika se bude dále rozvíjet a modifikovat dle výsledků laboratorních experimentů. Předpokládá se měření tyčí různých rozměrů s různými velikostmi deformací, dále bude testována opakovatelnost měření. Po odladění softwaru pro měření výkovků za studena bude probíhat jeho modifikace pro měření za tepla. Výsledky měření žhavých výkovků budou porovnány s výsledky měření pomocí laserového systému.
strana
44
SEZNAM LITERATURY 8
8 SEZNAM LITERATURY ÚVOD [1] GIROD, B., GREINER, G.ANDNIEMANN, H. Principles of 3D Image Analysis and Synthesis.USA: Kluwer Academic Publishers, 2000. p. 466, ISBN 0-7923-7850-4. [2] LUHMANN, Thomas, et al. Close Range Photogrammetry: Principles, techniques and applications. Dunbeath, Caithness KW6 6EG, Scotland, UK: Whittles Publishing, 2011. Fundamental methods, p. 510. ISBN 978-184995057-2. [3] BERNARD, A. Reverse engineering for rapid product development: a state of the art. In K. G. HARDING.Three-Dimensional Imaging, Optical Metrology, and Inspection V. 1999, vol. 3835, p. 50-63. [4] REMONDINO, F.; EL-HAKIM, S. Image-Based 3D Modeling: A review. The photogrammetric Record. September 2006, vol. 115, no. 21, p. 269–291. [5] VOJÁČEK, Antonín. PMD senzor & 3D měření vzdálenosti - 1. část "Princip" [online]. 2006-07-03 [cit. 2011-06-23]. . [6] MATTHEWS, Neffra; NOBLE, Tom. Aerial and Close-Range Photogrammetric Technology: Providing Resource Documentation, Interpretation and Preservation. In [online]. Denver Federal Center: BLM, Technical Note 428, 2008 [cit. 2011-06-23]. < http://www.blm.gov/nstc/library/pdf/TN428.pdf>. [7] CHEN, F., BROWN, G.M. AND SONG, M. Overview of three-dimensional shape measurement using optical methods. Optical Engineering, 2000, vol. 39, no. 1, p. 10-22. METODY MĚŘENÍ ROZMĚRNÝCH VÝKOVKŮ [8] NYE, T.J., ELBADAN, A.M., BONE, G.M. Real-time process characterization of open die forging for adaptive control. ASME Journal of Engineering Materials andTechnology, 2001, vol. 123, no. 4, p. 511–516. [9] SIEMER, E., NEISCHWITZ, P., KOPP, R. Quality optimized process control in open die forging. Stahl und Eisen, 1986. Vol. 106, no. 8, p. 383–388. [10] YAMAUCHI, M. Errors in optical shape measurement caused by a hightemperature atmosphere.Opt. Eng., Sep 02, 2009, vol. 48, no. 9. [11]OMRON. OMRON.[online]. c2011 [cit. 2011-07-25]. [12] TIAN, Z. at al. Dimension measurement of hot large forgings with a novel timeof-flight system. The International Journal of Advanced Manufacturing Technology, 2009, vol. 44, no. 1, p. 125-132. [13] MÄTTA, K., KOSTAMOVAARA, J. AND MYLLYL, R. Profiling of hot surfaces by pulsed time-of-flight laser range finder techniques. Applied optics,
strana
45
SEZNAM LITERATURY
1993, vol. 32, no. 27, p. 5334-5347. [14] MITEQ. FERROTRON A MINTEQ DIVISION [online]. c2011 [cit. 2011-07-25]. LaCam Forge.. [15] MERMEC GROUP. MERMEC GROUP [online]. c2009 [cit. 2011-07-25]. Open Die Forgings.. [16] SINTEF. SINTEF [online]. November 19 2008 [cit. 2011-07-25]. Measuring 3D shape of steel shafts during forging. . [17] JIA, Z. at al. An improved image acquiring method for machine vision measurement of hot formed parts. Journal of Materials Processing Technology, 2010, vol. 210, no. 2, p. 267-271. [18] LIU, W. at al. Fast dimensional measurement method and experiment of the forgings under high temperature. Journal of Materials Processing Technology, 2011. vol. 211, Issue 2, p. 237–244. [19] CHAUDHURI, S. at al. Automatic measurement of frontal area and volume of fluted ingots using image-based instrumentation. Measurement, 2003, vol. 33, no. 4, p. 325-332. [20] DWORKIN, S.B. AND NYE, T.J. Image processing for machine vision measurement of hot formed parts. Journal of materials processing technology, 2006, vol. 174, no. 1-3, p. 1-6. [21] Ch. Hu, B. Liu, X. Song. A novel edge detection approach used for online dimensional measurement of heavy forging.2008 International Conference on Optical Instruments and Technology: Optoelectronic Measurement Technology and Applications, Proc. of SPIE,vol. 7160. FOTOGRAMMETRIE – APLIKACEV PRŮMYSLU [23] LUHMANN, T. Close range photogrammetry for industrial applications. ISPRS Journal of Photogrammetry and Remote Sensing, 2010, vol. 65, no. 6, p. 558569. [24] FRASER, C.S. State of the art in industrial photogrammetry. IAPRS, Kyoto, 1988, vol. 27, p. 166-181. [25] Reinshow. Reinshow [online]. 2001-2011 [cit. 2011-06-28]. Dostupné z WWW: . [26] Metronor. Metronor [online]. 2007 [cit. 2011-06-28]. Dostupné z WWW: < http://www.metronor.com/ >. [27] AICON. AICON [online]. 2011 [cit. 2011-06-28]. Dostupné z WWW:
strana
46
SEZNAM LITERATURY
[28] Nikon Metrology. Nikon Metrology NV [online]. 2011 [cit. 2011-06-30]. Dostupné z WWW: . [29] GOM. GOM: Optical Measurement Techniques [online]. 2011 [cit. 2011-06-28]. Dostupné z WWW: < http://www.gom.com/>. [30] JIANDONG, Z., LIYAN, Z. AND XIAOYU, D. Accurate 3D Target Positioning in Close Range Photogrammetry with Implicit Image Correction. Chinese Journal of Aeronautics, 2009, vol. 22, no. 6, p. 649-657. [31] Géodésie Maintenance Services. Géodésie Maintenance Services [online]. 1998 [cit. 2011-06-30]. Dostupné z WWW: . DETEKCE HRAN [32] ROSENFELD, A. Picture processing by computer. ACM Computing Surveys (CSUR). 1969, vol. 1,no. 3, p. 147-176. [33] DAVIS, L. S. A survey of edge detection techniques. Computer Graphics and Image Processing, 1975, vol. 4, no. 3, p. 248-270. [34] PELI, T. AND MALAH, D. A study of edge detection algorithms* 1. Computer Graphics and Image Processing, 1982, vol. 20, no. 1, p. 1-21. [35] ZIOU, D. AND TABBONE, S. Edge detection techniques-an overview. International Journal of Pattern Recognition and Image Analysis, 1998. [36] PAPARI, G. AND PETKOV, N. Edge and line oriented contour detection: State of the art. Image andVision Computing, 2011, vol. 29, no. 2-3, p. 79-103. [37] CANNY J. A computational approach to edge detection. IEEE Transactions on Pattern Analysis and Machine intelligence. 1986; 8(6): 679-698. [38] MEDINA-CARNICER, R. at al. A novel method to look for the hysteresis thresholds for the Canny edge detector. Pattern Recognition, Jun 2011, vol. 44, no. 6, p. 1201-1211. [39] DING, L. AND GOSHTASBY, A. On the Canny edge detector. Pattern Recognition, 2001, vol. 34, no. 3, p. 721-725. [40] DA, F.P. AND ZHANG, H. Sub-pixel edge detection based on an improved moment. Image and Vision Computing, Dec 2010, vol. 28, no. 12, p. 1645-1658. LOKALIZACE POLOHY ELIPS A VÝPOČET JEJÍCH PARAMETRŮ [41] SHORTIS, M.R., CLARKE, T.A. AND SHORT, T. A comparison of some techniques for the subpixel location of discrete target images. SPIE Videometrics, 1994, vol. 2350, no. III, p. 239-249. [42] YUEN, H.K., ILLINGWORTH, J. AND KITTLER, J. Detecting partially occluded ellipses using the Hough transform. Image and Vision Computing, 1989, vol. 7, no. 1, p. 31-37. [43] YIP, R.K.K., TAM, P.K.S. AND LEUNG, D.N.K. Modification of hough transform for circles and ellipses detection using a 2-dimensional array. Pattern Recognition, 1992, vol. 25, no. 9, p. 1007-1022.
strana
47
SEZNAM LITERATURY
[44] WU, W.-Y. AND WANG, M.-J.J. Elliptical object detection by using its geometric properties. Pattern Recognition, 1993, vol. 26, no. 10, p. 1499-1509. [45] NAIR, P.S. AND SAUNDERS, A.T. Hough transform based ellipse detection algorithm. Pattern Recognition Letters, 1996, vol. 17, no. 7, p. 777-784. [46] ROSIN, P.L. Ellipse fitting by accumulating five-point fits. Pattern Recognition Letters, 1993, vol. 14, no. 8, p. 661-669 [47] WERMAN, M. AND GEYZEL, Z. Fitting a second degree curve in the presence of error. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 1995, vol. 17, no. 2, p. 207-211. [48] DAVE, R.N. AND BHASWAN, K. Adaptive fuzzy shells clustering and detection of ellipses. Neural Networks, IEEE Transactions on, 1992, vol. 3, no. 5, p. 643-662. [49] GATH, I. AND HOORY, D. Fuzzy clustering of elliptic ring-shaped clusters. Pattern Recognition Letters, 1995, vol. 16, no. 7, p. 727-741. [50] GANDER, W., GOLUB, G.H. AND STREBEL, R. Least-squares fitting of circles and ellipses. BIT Numerical Mathematics,1994, vol. 34, no. 4, p. 558578. [51] FITZGIBBON, A.W., PILU, M. AND FISHER, R.B. Direct least squares fitting of ellipses. Pattern Recognition, Proceedings of the 13th International Conference on, 1996, vol. 1, p. 253-257. [52] HAL R, R. AND FLUSSER, J. Numerically stable direct least squares fitting of ellipses. The Sixth International Conference in Central Europe on Computer Graphics and Visualization, 1998, p. 59-108. [53] ANDREW, F. Direct Least Square Fitting of Ellipses. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1999, vol. 21, p. 476-480. [54] MAINI, E.S. Enhanced direct least square fitting of ellipses. International Journal of Pattern Recognition & Artificial Intelligence, 2006, vol. 20, no. 6, p. 939-953. [55] RAY, A. AND SRIVASTAVA, D.C. Non-linear least squares ellipse fitting using the genetic algorithm with applications to strain analysis. Journal of Structural Geology, 2008, vol. 30, no. 12, p. 1593-1602. [56] CHAUDHURI, D. A simple least squares method for fitting of ellipses and circles depends on border points of a two-tone image and their 3-D extensions. Pattern Recognition Letters, 2010, vol. 31, no. 9, p. 818-829. [57] KANATANI, K. AND RANGARAJAN, P. Hyper least squares fitting of circles and ellipses. Computational Statistics & Data Analysis, 2011, vol. 55, no. 6, p. 2197-2208. [58] JIANDONG, Z., LIYAN, Z. AND XIAOYU, D. Accurate 3D Target Positioning in Close Range Photogrammetry with Implicit Image Correction. Chinese Journal of Aeronautics, 2009, vol. 22, no. 6, p. 649-657.
strana
48
SEZNAM LITERATURY
VNITŘNÍ KALIBRACE [59] HARTLEY, R. AND ZISSERMAN, A. Multiple view geometry in Computer Vision. Edtion ed.: Cambridge university press, 2000. ISBN 0-521-54051-8. [60] D. C. BROWN, Decentering distortion of lenses, Photogrammetric Engineering and Remote Sensing, 1966, no. 24, p. 555–566. [61] D. C. BROWN, Close-range camera calibration, Photogrammetric Engineering and Remote Sensing, 1971, no. 42, p. 855–866. [62] WANG, J. at al. A new calibration model of camera lens distortion. Pattern Recognition, 2008, vol. 41, no. 2, p. 607-615. [63] SUN, W. AND COOPERSTOCK, J. R. An empirical evaluation of factors influencing camera calibration accuracy using three publicly available techniques. Machine Vision and Applications, 2006, vol. 17, no. 1, p. 51-67. [64] Z. ZHANG, On the epipolar geometry between two images with lens distortion, Proceedings of International Conference on Computer Vision and Pattern Recognition, 1996, vol. 1, no. I, p. 407-411. [65] STEIN, G. P. Lens distortion calibration using point correspondences. Proceedings of the International Conference on Computer Vision and Pattern Recognition, 1997, p. 602-608. [66] FITZGIBBON, A. W. Simultaneous linear estimation of multiple view geometry and lens distortion. IEEE International Conference on Computer Vision and Pattern Recognition, 2001, vol. 1, p. I-125. [67] RICOLFE-VIALA, C. AND SÁNCHEZ-SALMERÓN, A.-J. Robust metric calibration of non-linear camera lens distortion. Pattern Recognition, 2010, vol. 43, no. 4, p. 1688-1699. VNĚJŠÍ KALIBRACE [68] ABDEL-AZIZ, Y. I. and KARARA, H. M., 1971: Direct linear transformation into object space coordinates in close-range photogrammetry. Proc. Symposium on Close-Range Photogrammetry, p. 1-18 [69] LONGUET-HIGGINS, H. C. A Computer Algorithm for Reconstructing a Scene from Two Projections, Nature, Septemebr 1981, vol. 293, p. 133-135. [70] HARTLEY, R.I. In defense of the eight-point algorithm. Ieee Transactions on Pattern Analysis and Machine Intelligence, Jun 1997, vol. 19, no. 6, p. 580-593. [71] CHOJNACKI, W. at al. Revisiting Hartley's normalized eight-point algorithm. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2003, p. 1172-1177. [72] NISTÉR, D. An efficient solution to the five-point relative pose problem. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2004, vol. 26, no. 6, p. 756-770. [73] LUONG, Q.T. AND FAUGERAS, O.D. The fundamental matrix: Theory, algorithms, and stability analysis. International Journal of Computer Vision, 1996, vol. 17, no. 1, p. 43-75.
strana
49
SEZNAM LITERATURY
[74] STEWÉNIUS, H. at al. A minimal solution for relative pose with unknown focal length. Computer Vision and Pattern Recognition, 2005. CVPR 2005. IEEE Computer Society Conference on, 2005, vol. 2, p. 789-794. [75] STEWENIUS, H., ENGELS, C. AND NISTER, D. Recent developments on direct relative orientation. Isprs Journal of Photogrammetry and Remote Sensing, Jun 2006, vol. 60, no. 4, p. 284-294. [76] LOURAKIS, M.I.A. AND ARGYROS, A.A. Is Levenberg-Marquardt the most efficient optimization algorithm for implementing bundle adjustment? 2005. [77] D. C. BROWN. The bundle adjustment-progress and prospects. Int. Archives Photogrammetry, 1976. vol. 21, no. 3, p. 3–33 [78] S. GRANSHAW. Bundle adjustment methods in engineering photogrammetry. Photogrammetric Record, 1980, vol. 10, no. 56, p. 181–207. [79] TRIGGS, B. at al. Bundle adjustment--a modern synthesis. Vision algorithms: theory and practice, 2000, p. 153-177. [80] REMONDINO, F. AND FRASER, C. Digital camera calibration methods: considerations and comparisons. IPRS International Archives of Photogrammetry, Remote Sensing and Spatial Information Sciences, 2006, vol. 36, no. 5, p. 266-272.
strana
50
SEZNAM OBRÁZKŮ A GRAFŮ 9
9 SEZNAM OBRÁZKŮ A GRAFŮ Obr. 2-1
Vztah mezi dosahovanou přesností měření a velikostí měřeného objektu u různých bezdotykových měřicích metod [2]
Obr. 1-2
Rozdělení optických metod 3D měření, podle [3, 4].
Obr. 1-3
Dosahované rozlišení a pracovní vzdálenost tří metod vyhodnocování třetího rozměru [5].
Obr. 1-4
Rozdělení bezkontaktních optických metod měření dle vyhodnocované vlnové délky [2].
Obr. 2-1
3D model části povrchu výkovku, získaný měřením pomocí laserového TOF systému [12].
Obr. 2-2
Proces optimalizace tvaru výkovku pomocí laserového systému LaCam od firmy Ferrotron [14].
Obr. 2-3
Protokol z měření pomocí systému LaCam. Z výsledků je patrné, že maximální průhyb překračuje povolenou odchylku 10 mm[14].
Obr. 2-4
Výstupní data získaná měřením žhavého výkovku laserovým systémem TopScan, který nabízí firma Tecnogamma [15].
Obr. 2-5
Výsledky měření pomocí systému od firemIMEGO, Mefos a SINTEF[16].
Obr. 2-6
Prvky optického systému s xenonovým světlem [17].
Obr. 2-7
Princip spektrálně-selektivní metody[17].
Obr. 2-8
Výsledek extrakce bodů na povrchu výkovku [17].
Obr. 2-9
Schéma optického stereo systému s xenonovou lampou a jeho reálná podoba [18].
Obr. 2-10 Modifikovaný Gausův model aplikovaný na světelné pruhy promítané na žhavý výkovek[18]. Obr. 2-11 Aplikace první derivace na obraz výkovku[18]. Obr. 2-12 Hraniční body detekované na ocelové desce[18]. Obr. 2-13 Monochromatický snímek žhavého výkovku při použití IR filtru [20]. Obr. 2-14 Výsledné prahování snímků, vlevo monochromatický bez IR filtru, vpravo s IR filtrem [20]. Obr. 2-15 Snímky žhavého výkovku ve stupních šedi[21]. Obr. 2-16 Výsledky experimentu s různými metodami[21]. Obr. 2-17 Výsledky experimentu s různými metodami[21]. Obr. 2-18 Proces on-line a off-line fotogrammetrického měření v průmyslu [8]. Obr. 2-19 Rozdělení průmyslových fotogrammetrických systémů [24]. Obr. 2-20 Fotogrammetrické měření pomocí tzv. kódovaných a nekódovaných značek [29].
strana
51
SEZNAM OBRÁZKŮ A GRAFŮ
Obr. 2-21 Dotykové sondy: ruční dotyková sonda od firmy Metronor [26], ruční dotyková sonda od firmy AICON [27], CMM od firmy Reinshaw [25]. Obr. 2-22 Projekce vzoru a fotogrammetrické měření systémem PRO-SPOT [31], fotogrammetrický systém určený k měření deformací od firmy GOM [29]. Obr. 2-23 Měření dynamické scény pomocí fotogrammetrického systému od firmy AICON [27]. Obr. 2-24 Orientace snímků pomocí algoritmu „Bundle adjustment“ [79]. Obr. 6-1
Navržený stativ pro testovací měření v reálných podmínkách.
Obr. 6-2
Měření výkovku, detekce hran v rovině A-A.
Obr. 6-3
Měření výkovku, detekce hran v rovině B-B.
Obr. 6-3
Graf měření osy v rovině A-A.
Obr. 6-4
Graf měření osy v rovině B-B.
Obr. 6-5
Čelní pohled na měřené osy výkovku v obou rovinách.
Obr. 6-6
Diagram navržené metodiky řešení.
Obr. 6-7
Vyrobená měřicí stanice s naznačeným objektovým souřadným systémem.
Obr. 6-8
Průmyslová monochromatická kameraAM-1600 GE s rozlišením 4872×3248 pixelů od společnosti JAI.
Obr. 6-9
Průmyslová monochromatickákamera Pike F-1600 s rozlišením 4872×3248 pixelů od společnosti Allied Vision Technology.
Obr. 6-10 Snímek měřené scény pro kalibraci s kalibračními měrkami a obrazovým souřadným systémem. Obr. 6-11 Schéma geometrie systému v rovině zx (zoxo). Obr. 6-12 Snímek scény se zatíženým lankem. Obr. 6-13 Snímek deformované válcové tyče, určený k testování navržené metody detekce hran. Obr. 6-14 Výsledek měření a softwarové rekonstrukcedeformované tyče pomocí systému Atos s určeným místem maximálního průhybu.
strana
52
SEZNAM TABULEK
10
10 SEZNAM TABULEK Tab. 1 Vyhodnocení odchylek způsobených radiálním zkreslením
39
Tab. 2 Velikost některých prvků scény odečtená ze snímků
40
Tab. 3 Výsledky výpočtů základních úhlů geometrie systému v rovině zx
41
Tab. 4 Výsledky předběžné kontroly přesnosti po kalibraci systému
42
strana
53